Functional RNA Structures in the 3'UTR of Tick-Borne, Insect-Specific and No-Known-Vector Flaviviruses

Tick-borne flavivirus exoribonuclease-resistant RNAs contain a double loop structure

Viruses from the Flaviviridae family contain human relevant pathogens that generate subgenomic noncoding RNAs during infection using structured exoribonuclease resistant RNAs (xrRNAs). These xrRNAs block progression of host cell's 5' to 3' exoribonucleases. The structures of several xrRNAs from mosquito-borne and insect-specific flaviviruses reveal a conserved fold in which a ring-like motif encircles the 5' end of the xrRNA. However, the xrRNAs found in tick-borne and no known vector flaviviruses have distinct characteristics, and their 3-D fold was unsolved. Here, we verify the presence of xrRNAs in the encephalitis-causing tick-borne Powassan Virus. We characterize their secondary structure and obtain a mid-resolution map of one of these xrRNAs using cryo-EM, revealing a unique double-loop ring element. Integrating these results with covariation analysis, biochemical data, and existing high-resolution structural information yields a model in which the core of the fold matches the previously solved xrRNA fold, but the expanded double loop ring is remodeled upon encountering the exoribonuclease. These results are representative of a broad class of xrRNAs and reveal a conserved strategy of structure-based exoribonuclease resistance achieved through a unique topology across a viral family of importance to global health.

Inter-species dynamics of non-coding RNAs: Impact on host immunomodulation and pathogen survival

Non-coding RNAs (ncRNAs) are composed of nucleotides that do not encode proteins but instead serve as guides. It interacts with amino acids at precise genomic sites, influencing chromatin structure and gene expression. These ncRNAs contribute to numerous inter-species dynamics, including those within the vector-host-pathogen triad. Vector-associated ncRNAs are released into hosts to combat the host immune system and sustain arthropod viability. Conversely, hosts may utilize specific ncRNAs as part of their defences to counteract pathogen-carrying vectors. Moreover, pathogens transmitted through vectors' saliva into hosts carry ncRNAs that enhances their virulence. While recent investigations have primarily focused on vector-associated ncRNAs in animal hosts, only a few have explored the functions of pathogen-associated ncRNAs and their role in initiating infections. Our review delves into the historical prospects of ncRNAs, mechanisms by which pathogen-derived ncRNAs influence host-pathogen interactions, regulate gene expression, and evade host defences. Ultimately, it underscores the importance ncRNAs mediated regulatory network in vector-host-pathogen dynamics, offering new strategies to combat vector-borne diseases and enhance public health outcomes.

Identification of Three Novel O-Linked Glycans in the Envelope Protein of Tick-Borne Encephalitis Virus

The tick-borne encephalitis virus is a pathogen endemic to northern Europe and Asia, transmitted through bites from infected ticks. It is a member of the Flaviviridae family and possesses a positive-sense, single-stranded RNA genome encoding a polypeptide that is processed into seven non-structural and three structural proteins, including the envelope (E) protein. The glycosylation of the E protein, involving a single N-linked glycan at position N154, plays a critical role in viral infectivity and pathogenesis. Here, we dissected the entire glycosylation profile of the E protein using liquid chromatography-tandem mass spectrometry and identified three novel O-linked glycans, which were found at relatively low frequency. One of the O-linked glycans was positioned close to the highly conserved N-linked glycan site, and structural analysis suggested that it may be relevant for the function of the E 150-loop. The N154 site was found to be glycosylated with a high frequency, containing oligomannose or complex-type structures, some of which were fucosylated. An unusually high portion of oligomannose N-linked glycan structures exhibited compositions that are normally observed on proteins when they are translocated from the endoplasmic reticulum to the trans-Golgi network, suggesting disruption of the glycan processing pathway in the infected cells from which the E protein was obtained.

Pan-flavivirus analysis reveals sfRNA-independent, 3' UTR-biased siRNA production from an insect-specific flavivirus

RNA interference (RNAi) plays an essential role in mosquito antiviral immunity, but it is not known whether viral small interfering RNA (siRNA) profiles differ between mosquito-borne and mosquito-specific viruses. A pan-Orthoflavivirus analysis in Aedes albopictus cells revealed that viral siRNAs were evenly distributed across the viral genome of most representatives of the Flavivirus genus. In contrast, siRNA production was biased toward the 3' untranslated region (UTR) of the genomes of classical insect-specific flaviviruses (cISF), which was most pronounced for Kamiti River virus (KRV), a virus with a unique, 1.2 kb long 3' UTR. KRV-derived siRNAs were produced in high quantities and almost exclusively mapped to the 3' UTR. We mapped the 5' end of KRV subgenomic flavivirus RNAs (sfRNAs), products of the 5'-3' exoribonuclease XRN1/Pacman stalling on secondary RNA structures in the 3' UTR of the viral genome. We found that KRV produces high copy numbers of a long, 1,017 nt sfRNA1 and a short, 421 nt sfRNA2, corresponding to two predicted XRN1-resistant elements. Expression of both sfRNA1 and sfRNA2 was reduced in Pacman-deficient Aedes albopictus cells; however, this did not correlate with a shift in viral siRNA profiles. We suggest that cISFs, particularly KRV, developed a unique mechanism to produce high amounts of siRNAs as a decoy for the antiviral RNAi response in an sfRNA-independent manner.IMPORTANCEThe Flavivirus genus contains diverse mosquito viruses ranging from insect-specific viruses circulating exclusively in mosquito populations to mosquito-borne viruses that cause disease in humans and animals. Studying the mechanisms of virus replication and antiviral immunity in mosquitoes is important to understand arbovirus transmission and may inform the development of disease control strategies. In insects, RNA interference (RNAi) provides broad antiviral activity and constitutes a major immune response against viruses. Comparing diverse members of the Flavivirus genus, we found that all flaviviruses are targeted by RNAi. However, the insect-specific Kamiti River virus was unique in that small interfering RNAs are highly skewed toward its uniquely long 3' untranslated region. These results suggest that mosquito-specific viruses have evolved unique mechanisms for genome replication and immune evasion.

Xinyang flavivirus, from Haemaphysalis flava ticks in Henan Province, China, defines a basal, likely tick-only Orthoflavivirus clade

Tick-borne orthoflaviviruses (TBFs) are classified into three conventional groups based on genetics and ecology: mammalian, seabird and probable-TBF group. Recently, a fourth basal group has been identified in Rhipicephalus ticks from Africa: Mpulungu flavivirus (MPFV) in Zambia and Ngoye virus (NGOV) in Senegal. Despite attempts, isolating these viruses in vertebrate and invertebrate cell lines or intracerebral injection of newborn mice with virus-containing homogenates has remained unsuccessful. In this study, we report the discovery of Xinyang flavivirus (XiFV) in Haemaphysalis flava ticks from Xìnyáng, Henan Province, China. Phylogenetic analysis shows that XiFV was most closely related to MPFV and NGOV, marking the first identification of this tick orthoflavivirus group in Asia. We developed a reverse transcriptase quantitative PCR assay to screen wild-collected ticks and egg clutches, with absolute infection rates of 20.75 % in adult females and 15.19 % in egg clutches, suggesting that XiFV could be potentially spread through transovarial transmission. To examine potential host range, dinucleotide composition analyses revealed that XiFV, MPFV and NGOV share a closer composition to classical insect-specific orthoflaviviruses than to vertebrate-infecting TBFs, suggesting that XiFV could be a tick-only orthoflavivirus. Additionally, both XiFV and MPFV lack a furin cleavage site in the prM protein, unlike other TBFs, suggesting these viruses might exist towards a biased immature particle state. To examine this, chimeric Binjari virus with XIFV-prME (bXiFV) was generated, purified and analysed by SDS-PAGE and negative-stain transmission electron microscopy, suggesting prototypical orthoflavivirus size (~50 nm) and bias towards uncleaved prM. In silico structural analyses of the 3'-untranslated regions show that XiFV forms up to five pseudo-knot-containing stem-loops and a prototypical orthoflavivirus dumbbell element, suggesting the potential for multiple exoribonuclease-resistant RNA structures.

Viral intra-host evolutionary dynamics revealed via serial passage of Japanese encephalitis virus in vitro

Analyses of viral inter- and intra-host mutations could better guide the prevention and control of infectious diseases. For a long time, studies on viral evolution have focused on viral inter-host variations. Next-generation sequencing has accelerated the investigations of viral intra-host diversity. However, the theoretical basis and dynamic characteristics of viral intra-host mutations remain unknown. Here, using serial passages of the SA14-14-2 vaccine strain of Japanese encephalitis virus (JEV) as the in vitro model, the distribution characteristics of 1,788 detected intra-host single-nucleotide variations (iSNVs) and their mutated frequencies from 477 deep-sequenced samples were analyzed. Our results revealed that in adaptive (baby hamster kidney (BHK)) cells, JEV is under a nearly neutral selection pressure, and both non-synonymous and synonymous mutations represent an S-shaped growth trend over time. A higher positive selection pressure was observed in the nonadaptive (C6/36) cells, and logarithmic growth in non-synonymous iSNVs and linear growth in synonymous iSNVs were observed over time. Moreover, the mutation rates of the NS4B protein and the untranslated region (UTR) of the JEV are significantly different between BHK and C6/36 cells, suggesting that viral selection pressure is regulated by different cellular environments. In addition, no significant difference was detected in the distribution of mutated frequencies of iSNVs between BHK and C6/36 cells.

Strukturierte RNAs in Viren

Evolutionarily conserved RNAs in untranslated regions are key regulators of the viral life cycle. Exoribonuclease-resistant RNAs (xrRNAs) are particularly interesting examples of structurally conserved elements because they actively dysregulate the messenger RNA (mRNA) degradation machinery of host cells, thereby mediating viral pathogenicity. We review the principles of RNA structure conservation in viruses and discuss potential applications of xrRNAs in synthetic biology and future mRNA vaccines.

Interactions of host miRNAs in the flavivirus 3´UTR genome: From bioinformatics predictions to practical approaches

The genus Flavivirus of the Flaviviridae family includes important viruses, such as Dengue, Zika, West Nile, Japanese encephalitis, Murray Valley encephalitis, tick-borne encephalitis, Yellow fever, Saint Louis encephalitis, and Usutu viruses. They are transmitted by mosquitoes or ticks, and they can infect humans, causing fever, encephalitis, or haemorrhagic fever. The treatment resources for these diseases and the number of vaccines available are limited. It has been discovered that eukaryotic cells synthesize small RNA molecules that can bind specifically to sequences present in messenger RNAs to inhibit the translation process, thus regulating gene expression. These small RNAs have been named microRNAs, and they have an important impact on viral infections. In this review, we compiled the available information on miRNAs that can interact with the 3’ untranslated region (3’UTR) of the flavivirus genome, a conserved region that is important for viral replication and translation.

Comparative analysis of genome-encoded viral sequences reveals the evolutionary history of flavivirids (family Flaviviridae)

Flavivirids (family Flaviviridae) are a group of positive-strand ribonucleic acid (RNA) viruses that pose serious risks to human and animal health on a global scale. Here, we use flavivirid-derived deoxyribonucleic acid (DNA) sequences, identified in animal genomes, to reconstruct the long-term evolutionary history of family Flaviviridae. We demonstrate that flavivirids are >100 million years old and show that this timing can be combined with dates inferred from co-phyletic analysis to produce a cohesive overview of their evolution, distribution, and diversity wherein the main flavivirid subgroups originate in early animals and broadly co-diverge with major animal phyla. In addition, we reveal evidence that the 'classical flaviviruses' of vertebrates, most of which are transmitted via blood-feeding arthropod vectors, originally evolved in haematophagous arachnids and later acquired the capacity to be transmitted by insects. Our findings imply that the biological properties of flavivirids have been acquired gradually over the course of animal evolution. Thus, broad-scale comparative analysis will likely reveal fundamental insights into their biology. We therefore published our results via an open, extensible, database (Flavivirid-GLUE), which we constructed to facilitate the wider utilisation of genomic data and evolution-related domain knowledge in flavivirid research.

Evolutionary traits of Tick-borne encephalitis virus: Pervasive non-coding RNA structure conservation and molecular epidemiology

Tick-borne encephalitis virus (TBEV) is the aetiological agent of tick-borne encephalitis, an infectious disease of the central nervous system that is often associated with severe sequelae in humans. While TBEV is typically classified into three subtypes, recent evidence suggests a more varied range of TBEV subtypes and lineages that differ substantially in the architecture of their 3ʹ untranslated region (3ʹUTR). Building on comparative genomic approaches and thermodynamic modelling, we characterize the TBEV UTR structureome diversity and propose a unified picture of pervasive non-coding RNA structure conservation. Moreover, we provide an updated phylogeny of TBEV, building on more than 220 publicly available complete genomes, and investigate the molecular epidemiology and phylodynamics with Nextstrain, a web-based visualization framework for real-time pathogen evolution.

Recent advances in RNA structurome

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

In Vitro and In Vivo Characterization of a New Strain of Mosquito Flavivirus Derived from Culicoides

Mosquito-specific flaviviruses comprise a group of insect-specific viruses with a single positive RNA, which can affect the duplication of mosquito-borne viruses and the life growth of mosquitoes, and which have the potential to be developed as a vaccine platform for mosquito-borne viruses. In this study, a strain of mosquito flavivirus (MFV) YN15-283-02 was detected in Culicoides collected from Yunnan, China. The isolation of the purified MFV YN15-283-02 from cell culture failed, and the virus was then rescued by an infectious clone. To study the biological features of MFV YN15-283-02 in vitro and in vivo, electron microscopy, phylogenetic tree, and viral growth kinetic analyses were performed in both cell lines and mosquitoes. The rescued MFV (rMFV) YN15-283-02 duplicated and reached a peak in C6/36 cells at 6 d.p.i. with approximately 2 × 10⁶ RNA copies/μL (RNA to cell ratio of 0.1), but without displaying a cytopathic effect. In addition, the infection rate for the rMFV in Ae.aegypti show a low level in both larvae (≤15%) and adult mosquitoes (≤12%).

Rational design of West Nile virus vaccine through large replacement of 3' UTR with internal poly(A)

The genus Flavivirus comprises numerous emerging and re-emerging arboviruses causing human illness. Vaccines are the best approach to prevent flavivirus diseases. But pathogen diversities are always one of the major hindrances for timely development of new vaccines when confronting unpredicted flavivirus outbreaks. We used West Nile virus (WNV) as a model to develop a new live-attenuated vaccine (LAV), WNV-poly(A), by replacing 5' portion (corresponding to SL and DB domains in WNV) of 3'-UTR with internal poly(A) tract. WNV-poly(A) not only propagated efficiently in Vero cells, but also was highly attenuated in mouse model. A single-dose vaccination elicited robust and long-lasting immune responses, conferring full protection against WNV challenge. Such "poly(A)" vaccine strategy may be promising for wide application in the development of flavivirus LAVs because of its general target regions in flaviviruses.

A Mathematical Analysis of HDV Genotypes: From Molecules to Cells

Hepatitis D virus (HDV) is classified according to eight genotypes. The various genotypes are included in the HDVdb database, where each HDV sequence is specified by its genotype. In this contribution, a mathematical analysis is performed on RNA sequences in HDVdb. The RNA folding predicted structures of the Genbank HDV genome sequences in HDVdb are classified according to their coarse-grain tree-graph representation. The analysis allows discarding in a simple and efficient way the vast majority of the sequences that exhibit a rod-like structure, which is important for the virus replication, to attempt to discover other biological functions by structure consideration. After the filtering, there remain only a small number of sequences that can be checked for their additional stem-loops besides the main one that is known to be responsible for virus replication. It is found that a few sequences contain an additional stem-loop that is responsible for RNA editing or other possible functions. These few sequences are grouped into two main classes, one that is well-known experimentally belonging to genotype 3 for patients from South America associated with RNA editing, and the other that is not known at present belonging to genotype 7 for patients from Cameroon. The possibility that another function besides virus replication reminiscent of the editing mechanism in HDV genotype 3 exists in HDV genotype 7 has not been explored before and is predicted by eigenvalue analysis. Finally, when comparing native and shuffled sequences, it is shown that HDV sequences belonging to all genotypes are accentuated in their mutational robustness and thermodynamic stability as compared to other viruses that were subjected to such an analysis.

Rapid structure-function insights via hairpin-centric analysis of big RNA structure probing datasets

The functions of RNA are often tied to its structure, hence analyzing structure is of significant interest when studying cellular processes. Recently, large-scale structure probing (SP) studies have enabled assessment of global structure-function relationships via standard data summarizations or local folding. Here, we approach structure quantification from a hairpin-centric perspective where putative hairpins are identified in SP datasets and used as a means to capture local structural effects. This has the advantage of rapid processing of big (e.g. transcriptome-wide) data as RNA folding is circumvented, yet it captures more information than simple data summarizations. We reformulate a statistical learning algorithm we previously developed to significantly improve precision of hairpin detection, then introduce a novel nucleotide-wise measure, termed the hairpin-derived structure level (HDSL), which captures local structuredness by accounting for the presence of likely hairpin elements. Applying HDSL to data from recent studies recapitulates, strengthens and expands on their findings which were obtained by more comprehensive folding algorithms, yet our analyses are orders of magnitude faster. These results demonstrate that hairpin detection is a promising avenue for global and rapid structure-function analysis, furthering our understanding of RNA biology and the principal features which drive biological insights from SP data.

Supramolecular Cylinders Target Bulge Structures in the 5' UTR of the RNA Genome of SARS-CoV-2 and Inhibit Viral Replication

The untranslated regions (UTRs) of viral genomes contain a variety of conserved yet dynamic structures crucial for viral replication, providing drug targets for the development of broad spectrum anti-virals. We combine in vitro RNA analysis with molecular dynamics simulations to build the first 3D models of the structure and dynamics of key regions of the 5' UTR of the SARS-CoV-2 genome. Furthermore, we determine the binding of metallo-supramolecular helicates (cylinders) to this RNA structure. These nano-size agents are uniquely able to thread through RNA junctions and we identify their binding to a 3-base bulge and the central cross 4-way junction located in stem loop 5. Finally, we show these RNA-binding cylinders suppress SARS-CoV-2 replication, highlighting their potential as novel anti-viral agents.

Supramolecular Cylinders Target Bulge Structures in the 5' UTR of the RNA Genome of SARS-CoV-2 and Inhibit Viral Replication*

The untranslated regions (UTRs) of viral genomes contain a variety of conserved yet dynamic structures crucial for viral replication, providing drug targets for the development of broad spectrum anti-virals. We combine in vitro RNA analysis with molecular dynamics simulations to build the first 3D models of the structure and dynamics of key regions of the 5' UTR of the SARS-CoV-2 genome. Furthermore, we determine the binding of metallo-supramolecular helicates (cylinders) to this RNA structure. These nano-size agents are uniquely able to thread through RNA junctions and we identify their binding to a 3-base bulge and the central cross 4-way junction located in stem loop 5. Finally, we show these RNA-binding cylinders suppress SARS-CoV-2 replication, highlighting their potential as novel anti-viral agents.

Shared properties and singularities of exoribonuclease-resistant RNAs in viruses

What viral RNA genomes lack in size, they make up for in intricacy. Elaborate RNA structures embedded in viral genomes can hijack essential cellular mechanisms aiding virus propagation. Exoribonuclease-resistant RNAs (xrRNAs) are an emerging class of viral elements, which resist degradation by host cellular exoribonucleases to produce viral RNAs with diverse roles during infection. Detailed three-dimensional structural studies of xrRNAs from flaviviruses and a subset of plant viruses led to a mechanistic model in which xrRNAs block enzymatic digestion using a ring-like structure that encircles the 5' end of the resistant structure. In this mini-review, we describe the state of our understanding of the phylogenetic distribution of xrRNAs, their structures, and their conformational dynamics. Because xrRNAs have now been found in several major superfamilies of RNA viruses, they may represent a more widely used strategy than currently appreciated. Could xrRNAs represent a 'molecular clock' that would help us understand virus evolution and pathogenicity? The more we study xrRNAs in viruses, the closer we get to finding xrRNAs within cellular RNAs.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions

Motivation

The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3′ untranslated regions (3′UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies.

Results

We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3′UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3′UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses.

Availability and implementation

The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Novel Virus of Flaviviridae Associated with Sexual Precocity in Macrobrachium rosenbergii

Since 2010, sexual precocity, a typical sign of the iron prawn syndrome (IPS), resulting in the reduced size of farmed giant freshwater prawns Macrobrachium rosenbergii, has caused substantial production losses. However, the cause of IPS was not clear. We ran tests for eight major shrimp pathogens, but none were detected from IPS-affected prawns. We performed the histopathological examination of tissues and identified an eosinophilic inclusion in the perinuclear cytoplasm of cells in various tissues associated with nervous and endocrinal functions in the compound eyes. A subsequent bioassay with viral extracts of IPS-affected samples reproduced the gross signs of IPS. Metatranscriptomic sequencing identified a novel virus of Flaviviridae in all IPS-affected M. rosenbergii prawns, which was not found in samples without IPS. This virus contains a positive-sense, single-stranded RNA genome of 12,630 nucleotides (nt). Phylogenetic analysis of the conserved RdRp and NS3 domains showed that it may belong to a new genus between Jingmenvirus and Flavivirus. Under transmission electron microscopy (TEM), putative virus particles showed as spherical with a diameter of 40 to 60 nm. In situ hybridization found hybridization signals consistent with the histopathology in the compound eyes from IPS-affected M. rosenbergii. We provisionally name this virus infectious precocity virus (IPV) and propose the binominal Latin name Crustaflavivirus infeprecoquis gen. nov., sp. nov. We developed a nested reverse transcription-PCR diagnostic assay and confirmed that all IPS-affected prawns tested IPV positive but normal prawns tested negative. Collectively, our study revealed a novel virus of Flaviviridae associated with sexual precocity in M. rosenbergii. IMPORTANCE The iron prawn syndrome (IPS), also described as sexual precocity, results in the reduced size of farmed prawns at harvest and significant economic losses. IPS has been frequently reported in populations of farmed Macrobrachium rosenbergii since 2010, but the cause was heretofore unknown. Here, we reported a novel virus identified from prawns with IPS using infection experiments, metatranscriptomic sequencing, and transmission electron microscopy and provisionally named it infectious precocity virus (IPV). Phylogenetic analysis showed that IPV represents a new genus, proposed as Crustaflavivirus gen. nov., in the family Flaviviridae. This study provides novel insight that a viral infection may cause pathological change and sexual maturation and subsequently affect crustacean growth. Therefore, we call for quarantine inspection of IPV in transboundary trade of live M. rosenbergii and enhanced surveillance of IPV in aquaculture in the region and globally.

Mechanisms Underlying Host Range Variation in Flavivirus: From Empirical Knowledge to Predictive Models

Preventing and controlling epidemics caused by vector-borne viruses are particularly challenging due to their diverse pool of hosts and highly adaptive nature. Many vector-borne viruses belong to the Flavivirus genus, whose members vary greatly in host range and specificity. Members of the Flavivirus genus can be categorized to four main groups: insect-specific viruses that are maintained solely in arthropod populations, mosquito-borne viruses and tick-borne viruses that are transmitted to vertebrate hosts by mosquitoes or ticks via blood feeding, and those with no-known vector. The mosquito-borne group encompasses the yellow fever, dengue, and West Nile viruses, all of which are globally spread and cause severe morbidity in humans. The Flavivirus genus is genetically diverse, and its members are subject to different host-specific and vector-specific selective constraints, which do not always align. Thus, understanding the underlying genetic differences that led to the diversity in host range within this genus is an important aspect in deciphering the mechanisms that drive host compatibility and can aid in the constant arms-race against viral threats. Here, we review the phylogenetic relationships between members of the genus, their infection bottlenecks, and phenotypic and genomic differences. We further discuss methods that utilize these differences for prediction of host shifts in flaviviruses and can contribute to viral surveillance efforts.

Information Encoded by the Flavivirus Genomes beyond the Nucleotide Sequence

The genus Flavivirus comprises numerous, small, single positive-stranded RNA viruses, many of which are important human pathogens. To store all the information required for their successful propagation, flaviviruses use discrete structural genomic RNA elements to code for functional information by the establishment of dynamic networks of long-range RNA-RNA interactions that promote specific folding. These structural elements behave as true cis-acting, non-coding RNAs (ncRNAs) and have essential regulatory roles in the viral cycle. These include the control of the formation of subgenomic RNAs, known as sfRNAs, via the prevention of the complete degradation of the RNA genome. These sfRNAs are important in ensuring viral fitness. This work summarizes our current knowledge of the functions performed by the genome conformations and the role of RNA-RNA interactions in these functions. It also reviews the role of RNA structure in the production of sfRNAs across the genus Flavivirus, and their existence in related viruses.

An African tick flavivirus forming an independent clade exhibits unique exoribonuclease-resistant RNA structures in the genomic 3'-untranslated region

Tick-borne flaviviruses (TBFVs) infect mammalian hosts through tick bites and can cause various serious illnesses, such as encephalitis and hemorrhagic fevers, both in humans and animals. Despite their importance to public health, there is limited epidemiological information on TBFV infection in Africa. Herein, we report that a novel flavivirus, Mpulungu flavivirus (MPFV), was discovered in a Rhipicephalus muhsamae tick in Zambia. MPFV was found to be genetically related to Ngoye virus detected in ticks in Senegal, and these viruses formed a unique lineage in the genus Flavivirus. Analyses of dinucleotide contents of flaviviruses indicated that MPFV was similar to those of other TBFVs with a typical vertebrate genome signature, suggesting that MPFV may infect vertebrate hosts. Bioinformatic analyses of the secondary structures in the 3′-untranslated regions (UTRs) revealed that MPFV exhibited unique exoribonuclease-resistant RNA (xrRNA) structures. Utilizing biochemical approaches, we clarified that two xrRNA structures of MPFV in the 3′-UTR could prevent exoribonuclease activity. In summary, our findings provide new information regarding the geographical distribution of TBFV and xrRNA structures in the 3′-UTR of flaviviruses.

Guapiaçu virus, a new insect-specific flavivirus isolated from two species of Aedes mosquitoes from Brazil

Classical insect-flaviviruses (cISFVs) and dual host-related insect-specific flavivirus (dISFV) are within the major group of insect-specific flavivirus. Remarkably dISFV are evolutionarily related to some of the pathogenic flavivirus, such as Zika and dengue viruses. The Evolutionary relatedness of dISFV to flavivirus allowed us to investigate the evolutionary principle of host adaptation. Additionally, dISFV can be used for the development of flavivirus vaccines and to explore underlying principles of mammalian pathogenicity. Here we describe the genetic characterization of a novel putative dISFV, termed Guapiaçu virus (GUAPV). Distinct strains of GUAPV were isolated from pools of Aedes terrens and Aedes scapularis mosquitoes. Additionally, we also detected viral GUAPV RNA in a plasma sample of an individual febrile from the Amazon region (North of Brazil). Although GUAPV did not replicate in tested mammalian cells, 3'UTR secondary structures duplication and codon usage index were similar to pathogenic flavivirus.

Dynamic Molecular Epidemiology Reveals Lineage-Associated Single-Nucleotide Variants That Alter RNA Structure in Chikungunya Virus

Chikungunya virus (CHIKV) is an emerging Alphavirus which causes millions of human infections every year. Outbreaks have been reported in Africa and Asia since the early 1950s, from three CHIKV lineages: West African, East Central South African, and Asian Urban. As new outbreaks occurred in the Americas, individual strains from the known lineages have evolved, creating new monophyletic groups that generated novel geographic-based lineages. Building on a recently updated phylogeny of CHIKV, we report here the availability of an interactive CHIKV phylodynamics dataset, which is based on more than 900 publicly available CHIKV genomes. We provide an interactive view of CHIKV molecular epidemiology built on Nextstrain, a web-based visualization framework for real-time tracking of pathogen evolution. CHIKV molecular epidemiology reveals single nucleotide variants that change the stability and fold of locally stable RNA structures. We propose alternative RNA structure formation in different CHIKV lineages by predicting more than a dozen RNA elements that are subject to perturbation of the structure ensemble upon variation of a single nucleotide.

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.

Different tertiary interactions create the same important 3D features in a distinct flavivirus xrRNA

During infection by a flavivirus (FV), cells accumulate noncoding subgenomic flavivirus RNAs (sfRNAs) that interfere with several antiviral pathways. These sfRNAs are formed by structured RNA elements in the 3' untranslated region (UTR) of the viral genomic RNA, which block the progression of host cell exoribonucleases that have targeted the viral RNA. Previous work on these exoribonuclease-resistant RNAs (xrRNAs) from mosquito-borne FVs revealed a specific three-dimensional fold with a unique topology in which a ring-like structure protectively encircles the 5' end of the xrRNA. Conserved nucleotides make specific tertiary interactions that support this fold. Examination of more divergent FVs reveals differences in their 3' UTR sequences, raising the question of whether they contain xrRNAs and if so, how they fold. To answer this, we demonstrated the presence of an authentic xrRNA in the 3' UTR of the Tamana bat virus (TABV) and solved its structure by X-ray crystallography. The structure reveals conserved features from previously characterized xrRNAs, but in the TABV version these features are created through a novel set of tertiary interactions not previously seen in xrRNAs. This includes two important A-C interactions, four distinct backbone kinks, several ordered Mg²⁺ ions, and a C⁺-G-C base triple. The discovery that the same overall architecture can be achieved by very different sequences and interactions in distantly related flaviviruses provides insight into the diversity of this type of RNA and will inform searches for undiscovered xrRNAs in viruses and beyond.

Superinfection Exclusion in Mosquitoes and Its Potential as an Arbovirus Control Strategy

The continuing emergence of arbovirus disease outbreaks around the world, despite the use of vector control strategies, warrants the development of new strategies to reduce arbovirus transmission. Superinfection exclusion, a phenomenon whereby a primary virus infection prevents the replication of a second closely related virus, has potential to control arbovirus disease emergence and outbreaks. This phenomenon has been observed for many years in plants, insects and mammalian cells. In this review, we discuss the significance of identifying novel vector control strategies, summarize studies exploring arbovirus superinfection exclusion and consider the potential for this phenomenon to be the basis for novel arbovirus control strategies.

A New Subclass of Exoribonuclease-Resistant RNA Found in Multiple Genera of Flaviviridae

Viruses have developed innovative strategies to exploit the cellular machinery and overcome the antiviral defenses of the host, often using specifically structured RNA elements. Examples are found in the Flavivirus genus (in the family Flaviviridae), where during flaviviral infection, pathogenic subgenomic flaviviral RNAs (sfRNAs) accumulate in the cell. These sfRNAs are formed when a host cell 5' to 3' exoribonuclease degrades the viral genomic RNA but is blocked by an exoribonuclease-resistant RNA structure (xrRNA) located in the viral genome's 3' untranslated region (UTR). Although known to exist in several Flaviviridae genera, the full distribution and diversity of xrRNAs in this family were unknown. Using the recently solved high-resolution structure of an xrRNA from the divergent flavivirus Tamana bat virus (TABV) as a reference, we used bioinformatic searches to identify xrRNAs in the remaining three genera of Flaviviridae: Pegivirus, Pestivirus, and Hepacivirus We biochemically and structurally characterized several examples, determining that they are genuine xrRNAs with a conserved fold. These new xrRNAs look superficially similar to the previously described xrRNAs but possess structural differences making them distinct from previous classes of xrRNAs. Overall, we have identified the presence of xrRNA in all four genera of Flaviviridae, but not in all species. Our findings thus require adjustments of previous xrRNA classification schemes and expand the previously known distribution of xrRNA in Flaviviridae.IMPORTANCE The members of the Flaviviridae comprise one of the largest families of positive-sense single-stranded RNA (+ssRNA) and are divided into the Flavivirus, Pestivirus, Pegivirus, and Hepacivirus genera. The genus Flavivirus contains many medically relevant viruses such as Zika virus, dengue virus, and Powassan virus. In these, a part of the RNA of the virus twists up into a distinct three-dimensional shape called an exoribonuclease-resistant RNA (xrRNA) that blocks the ability of the cell to "chew up" the viral RNA. Hence, part of the RNA of the virus remains intact, and this protected part is important for viral infection. These xrRNAs were known to occur in flaviviruses, but whether they existed in the other members of the family was not known. In this study, we identified a new subclass of xrRNA found not only in flaviviruses but also in the remaining three genera. The fact that these structured viral RNAs exist throughout the Flaviviridae family suggests they are important parts of the infection strategy of diverse pathogens, which could lead to new avenues of research.

A New Subclass of Exoribonuclease-Resistant RNA Found in Multiple Genera of Flaviviridae

The members of the Flaviviridae comprise one of the largest families of positive-sense single-stranded RNA (+ssRNA) and are divided into the Flavivirus, Pestivirus, Pegivirus, and Hepacivirus genera. The genus Flavivirus contains many medically relevant viruses such as Zika virus, dengue virus, and Powassan virus. In these, a part of the RNA of the virus twists up into a distinct three-dimensional shape called an exoribonuclease-resistant RNA (xrRNA) that blocks the ability of the cell to “chew up” the viral RNA. Hence, part of the RNA of the virus remains intact, and this protected part is important for viral infection. These xrRNAs were known to occur in flaviviruses, but whether they existed in the other members of the family was not known. In this study, we identified a new subclass of xrRNA found not only in flaviviruses but also in the remaining three genera. The fact that these structured viral RNAs exist throughout the Flaviviridae family suggests they are important parts of the infection strategy of diverse pathogens, which could lead to new avenues of research.

Viruses have developed innovative strategies to exploit the cellular machinery and overcome the antiviral defenses of the host, often using specifically structured RNA elements. Examples are found in the Flavivirus genus (in the family Flaviviridae), where during flaviviral infection, pathogenic subgenomic flaviviral RNAs (sfRNAs) accumulate in the cell. These sfRNAs are formed when a host cell 5′ to 3′ exoribonuclease degrades the viral genomic RNA but is blocked by an exoribonuclease-resistant RNA structure (xrRNA) located in the viral genome’s 3′ untranslated region (UTR). Although known to exist in several Flaviviridae genera, the full distribution and diversity of xrRNAs in this family were unknown. Using the recently solved high-resolution structure of an xrRNA from the divergent flavivirus Tamana bat virus (TABV) as a reference, we used bioinformatic searches to identify xrRNAs in the remaining three genera of Flaviviridae: Pegivirus, Pestivirus, and Hepacivirus. We biochemically and structurally characterized several examples, determining that they are genuine xrRNAs with a conserved fold. These new xrRNAs look superficially similar to the previously described xrRNAs but possess structural differences making them distinct from previous classes of xrRNAs. Overall, we have identified the presence of xrRNA in all four genera of Flaviviridae, but not in all species. Our findings thus require adjustments of previous xrRNA classification schemes and expand the previously known distribution of xrRNA in Flaviviridae.

Kyasanur Forest Disease and Alkhurma Hemorrhagic Fever Virus-Two Neglected Zoonotic Pathogens

Kyasanur Forest disease virus (KFDV) and Alkhurma hemorrhagic fever virus (AHFV) are tick-borne flaviviruses that cause life-threatening hemorrhagic fever in humans with case fatality rates of 3-5% for KFDV and 1-20% for AHFV, respectively. Both viruses are biosafety level 4 pathogens due to the severity of disease they cause and the lack of effective countermeasures. KFDV was discovered in India and is restricted to parts of the Indian subcontinent, whereas AHFV has been found in Saudi Arabia and Egypt. In recent years, both viruses have spread beyond their original endemic zones and the potential of AHFV to spread through ticks on migratory birds is a public health concern. While there is a vaccine with limited efficacy for KFDV used in India, there is no vaccine for AHFV nor are there any therapeutic concepts to combat infections with these viruses. In this review, we summarize the current knowledge about pathogenesis, vector distribution, virus spread, and infection control. We aim to bring attention to the potential public health threats posed by KFDV and AHFV and highlight the urgent need for the development of effective countermeasures.

Discoveries of Exoribonuclease-Resistant Structures of Insect-Specific Flaviviruses Isolated in Zambia

To monitor the arthropod-borne virus transmission in mosquitoes, we have attempted both to detect and isolate viruses from 3304 wild-caught female mosquitoes in the Livingstone (Southern Province) and Mongu (Western Province) regions in Zambia in 2017. A pan-flavivirus RT-PCR assay was performed to identify flavivirus genomes in total RNA extracted from mosquito lysates, followed by virus isolation and full genome sequence analysis using next-generation sequencing and rapid amplification of cDNA ends. We isolated a newly identified Barkedji virus (BJV Zambia) (10,899 nt) and a novel flavivirus, tentatively termed Barkedji-like virus (BJLV) (10,885 nt) from Culex spp. mosquitoes which shared 96% and 75% nucleotide identity with BJV which has been isolated in Israel, respectively. These viruses could replicate in C6/36 cells but not in mammalian and avian cell lines. In parallel, a comparative genomics screening was conducted to study evolutionary traits of the 5'- and 3'-untranslated regions (UTRs) of isolated viruses. Bioinformatic analyses of the secondary structures in the UTRs of both viruses revealed that the 5'-UTRs exhibit canonical stem-loop structures, while the 3'-UTRs contain structural homologs to exoribonuclease-resistant RNAs (xrRNAs), SL-III, dumbbell, and terminal stem-loop (3'SL) structures. The function of predicted xrRNA structures to stop RNA degradation by Xrn1 exoribonuclease was further proved by the in vitro Xrn1 resistance assay.

Understanding the Mechanisms Underlying Host Restriction of Insect-Specific Viruses

Arthropod-borne viruses contribute significantly to global mortality and morbidity in humans and animals. These viruses are mainly transmitted between susceptible vertebrate hosts by hematophagous arthropod vectors, especially mosquitoes. Recently, there has been substantial attention for a novel group of viruses, referred to as insect-specific viruses (ISVs) which are exclusively maintained in mosquito populations. Recent discoveries of novel insect-specific viruses over the past years generated a great interest not only in their potential use as vaccine and diagnostic platforms but also as novel biological control agents due to their ability to modulate arbovirus transmission. While arboviruses infect both vertebrate and invertebrate hosts, the replication of insect-specific viruses is restricted in vertebrates at multiple stages of virus replication. The vertebrate restriction factors include the genetic elements of ISVs (structural and non-structural genes and the untranslated terminal regions), vertebrate host factors (agonists and antagonists), and the temperature-dependent microenvironment. A better understanding of these bottlenecks is thus warranted. In this review, we explore these factors and the complex interplay between ISVs and their hosts contributing to this host restriction phenomenon.

Structures and Functions of the 3' Untranslated Regions of Positive-Sense Single-Stranded RNA Viruses Infecting Humans and Animals

The 3′ untranslated region (3′ UTR) of positive-sense single-stranded RNA [ssRNA(+)] viruses is highly structured. Multiple elements in the region interact with other nucleotides and proteins of viral and cellular origin to regulate various aspects of the virus life cycle such as replication, translation, and the host-cell response. This review attempts to summarize the primary and higher order structures identified in the 3′UTR of ssRNA(+) viruses and their functional roles.

Different Degrees of 5'-to-3' DAR Interactions Modulate Zika Virus Genome Cyclization and Host-Specific Replication

Mosquito-borne flaviviruses, which include many important human pathogens, such as West Nile virus (WNV), dengue virus (DENV), and Zika virus (ZIKV), have caused numerous emerging epidemics in recent years. Details of the viral genome functions necessary for effective viral replication in mosquito and vertebrate hosts remain obscure. Here, using ZIKV as a model, we found that the conserved "downstream of AUG region" (DAR), which is known to be an essential element for genome cyclization, is involved in viral replication in a host-specific manner. Mutational analysis of the DAR element showed that a single-nucleotide mismatch between the 5' DAR and the 3' DAR had little effect on ZIKV replication in mammalian cells but dramatically impaired viral propagation in mosquito cells. The revertant viruses passaged in mosquito cells generated compensatory mutations restoring the base pairing of the DAR, further confirming the importance of the complementarity of the DAR in mosquito cells. We demonstrate that a single-nucleotide mutation in the DAR is sufficient to destroy long-range RNA interaction of the ZIKV genome and affects de novo RNA synthesis at 28°C instead of 37°C, resulting in the different replication efficiencies of the mutant viruses in mosquito and mammalian cells. Our results reveal a novel function of the circular form of the flavivirus genome in host-specific viral replication, providing new ideas to further explore the functions of the viral genome during host adaptation.IMPORTANCE Flaviviruses naturally cycle between the mosquito vector and vertebrate hosts. The disparate hosts provide selective pressures that drive virus genome evolution to maintain efficient replication during host alteration. Host adaptation may occur at different stages of the viral life cycle, since host-specific viral protein processing and virion conformations have been reported in the individual hosts. However, the viral determinants and the underlying mechanisms associated with host-specific functions remain obscure. In this study, using Zika virus, we found that the DAR-mediated genome cyclization regulates viral replication differently and is under different selection pressures in mammalian and mosquito cells. A more constrained complementarity of the DAR is required in mosquito cells than in mammalian cells. Since the DAR element is stably maintained among mosquito-borne flaviviruses, our findings could provide new information for understanding the role of flavivirus genome cyclization in viral adaptation and RNA evolution in the two hosts.

Nanoscale Structure Determination of Murray Valley Encephalitis and Powassan Virus Non-Coding RNAs

Viral infections are responsible for numerous deaths worldwide. Flaviviruses, which contain RNA as their genetic material, are one of the most pathogenic families of viruses. There is an increasing amount of evidence suggesting that their 5' and 3' non-coding terminal regions are critical for their survival. Information on their structural features is essential to gain detailed insights into their functions and interactions with host proteins. In this study, the 5' and 3' terminal regions of Murray Valley encephalitis virus and Powassan virus were examined using biophysical and computational modeling methods. First, we used size exclusion chromatography and analytical ultracentrifuge methods to investigate the purity of in-vitro transcribed RNAs. Next, we employed small-angle X-ray scattering techniques to study solution conformation and low-resolution structures of these RNAs, which suggest that the 3' terminal regions are highly extended as compared to the 5' terminal regions for both viruses. Using computational modeling tools, we reconstructed 3-dimensional structures of each RNA fragment and compared them with derived small-angle X-ray scattering low-resolution structures. This approach allowed us to reinforce that the 5' terminal regions adopt more dynamic structures compared to the mainly double-stranded structures of the 3' terminal regions.

Updated Phylogeny of Chikungunya Virus Suggests Lineage-Specific RNA Architecture

Chikungunya virus (CHIKV), a mosquito-borne alphavirus of the family Togaviridae, has recently emerged in the Americas from lineages from two continents: Asia and Africa. Historically, CHIKV circulated as at least four lineages worldwide with both enzootic and epidemic transmission cycles. To understand the recent patterns of emergence and the current status of the CHIKV spread, updated analyses of the viral genetic data and metadata are needed. Here, we performed phylogenetic and comparative genomics screens of CHIKV genomes, taking advantage of the public availability of many recently sequenced isolates. Based on these new data and analyses, we derive a revised phylogeny from nucleotide sequences in coding regions. Using this phylogeny, we uncover the presence of several distinct lineages in Africa that were previously considered a single one. In parallel, we performed thermodynamic modeling of CHIKV untranslated regions (UTRs), which revealed evolutionarily conserved structured and unstructured RNA elements in the 3'UTR. We provide evidence for duplication events in recently emerged American isolates of the Asian CHIKV lineage and propose the existence of a flexible 3'UTR architecture among different CHIKV lineages.

Effects of Arbovirus Multi-Host Life Cycles on Dinucleotide and Codon Usage Patterns

Arthropod-borne viruses (arboviruses) of vertebrates including dengue, zika, chikungunya, Rift Valley fever, and blue tongue viruses cause extensive morbidity and mortality in humans, agricultural animals, and wildlife across the globe. As obligate intercellular pathogens, arboviruses must be well adapted to the cellular and molecular environment of both their arthropod (invertebrate) and vertebrate hosts, which are vastly different due to hundreds of millions of years of separate evolution. Here we discuss the comparative pressures on arbovirus RNA genomes as a result of a dual host life cycle, focusing on pressures that do not alter amino acids. We summarize what is currently known about arboviral genetic composition, such as dinucleotide and codon usage, and how cyclical infection of vertebrate and invertebrate hosts results in different genetic profiles compared with single-host viruses. To serve as a comparison, we compile what is known about arthropod tRNA, dinucleotide, and codon usages and compare this with vertebrates. Additionally, we discuss the potential roles of genetic robustness in arboviral evolution and how it may vary from other viruses. Overall, both arthropod and vertebrate hosts influence the resulting genetic composition of arboviruses, but a great deal remains to be investigated.

Musashi binding elements in Zika and related Flavivirus 3'UTRs: A comparative study in silico

Zika virus (ZIKV) belongs to a class of neurotropic viruses that have the ability to cause congenital infection, which can result in microcephaly or fetal demise. Recently, the RNA-binding protein Musashi-1 (Msi1), which mediates the maintenance and self-renewal of stem cells and acts as a translational regulator, has been associated with promoting ZIKV replication, neurotropism, and pathology. Msi1 predominantly binds to single-stranded motifs in the 3' untranslated region (UTR) of RNA that contain a UAG trinucleotide in their core. We systematically analyzed the properties of Musashi binding elements (MBEs) in the 3'UTR of flaviviruses with a thermodynamic model for RNA folding. Our results indicate that MBEs in ZIKV 3'UTRs occur predominantly in unpaired, single-stranded structural context, thus corroborating experimental observations by a biophysical model of RNA structure formation. Statistical analysis and comparison with related viruses show that ZIKV MBEs are maximally accessible among mosquito-borne flaviviruses. Our study addresses the broader question of whether other emerging arboviruses can cause similar neurotropic effects through the same mechanism in the developing fetus by establishing a link between the biophysical properties of viral RNA and teratogenicity. Moreover, our thermodynamic model can explain recent experimental findings and predict the Msi1-related neurotropic potential of other viruses.