Tick-borne encephalitis virus sequenced directly from questing and blood-feeding ticks reveals quasispecies variance

Eilat virus isolated from Culex univittatus mosquitoes from the Namibian Zambezi Region influences in vitro superinfection with alpha- and flaviviruses in a virus-species-dependent manner

The genus Alphavirus harbors arboviruses of great concern, such as the Chikungunya virus and the equine encephalitis viruses. Transmission of pathogenic alphaviruses by mosquitoes could be influenced by insect-specific alphaviruses such as Eilat virus (EILV). However, insect-specific alphaviruses are rarely found in wild mosquitoes and only a few have been described in the literature. Here, we report the detection of EILV in a Culex univittatus mosquito from the north-eastern Namibian Zambezi region. Full genome analysis of MP458-NA-2018 showed 94.5% nucleotide identity to an EILV isolate from Israel. MP458-NA-2018 grouped with EILV in phylogenetic analysis and was placed within the clade of insect-specific alphaviruses. The virus was isolated in mosquito cells and shown to be restricted to insects as hosts by the inability to infect different vertebrate cell lines and a complete block of virus replication at 34°C. We further showed that infection of cells with EILV MP458-NA-2018 reduced production of infectious particles of Sindbis virus by 2000-fold over the entire course of infection, whereas reduction rates of Bagaza and Middleburg virus were approximately 3-10-fold and dependent on time after infection. While production of infectious particles of cells superinfected with the Chikungunya virus were approximately 30-fold reduced and more pronounced at a lower multiplicity of infection of 0.01, EILV seemed to enhance production of West Nile virus infectious particles by >5-fold when superinfected at a multiplicitiy of infection of 0.1. In summary, EILV from the Namibian Zambezi region influences in vitro replication of endemic flavi- and alphaviruses.

Differential susceptibility of geographically distinct Ixodes ricinus populations to tick-borne encephalitis virus and louping ill virus

Tick-borne encephalitis virus (TBEV) is an emerging pathogen in the Netherlands. Multiple divergent viral strains are circulating and the focal distribution of TBEV remains poorly understood. This may, however, be explained by differences in the susceptibility of tick populations for specific viruses and viral strains, and by viral strains having higher infection success in their local tick population. We investigated this hypothesis by exposing Dutch Ixodes ricinus ticks to two different TBEV strains: TBEV-NL from the Netherlands and TBEV-Neudoerfl from Austria. In addition, we exposed ticks to louping Ill virus (LIV), which is endemic to large parts of the United Kingdom and Ireland, but has not been reported in the Netherlands. Ticks were collected from two locations in the Netherlands: one location without evidence of TBEV circulation and one location endemic for the TBEV-NL strain. Ticks were infected in a biosafety level 3 laboratory using an artificial membrane feeding system. Ticks collected from the region without evidence of TBEV circulation had lower infection rates for TBEV-NL as compared to TBEV-Neudoerfl. Vice versa, ticks collected from the TBEV-NL endemic region had higher infection rates for TBEV-NL compared to TBEV-Neudoerfl. In addition, LIV infection rates were much lower in Dutch ticks compared to TBEV, which may explain why LIV is not present in the Netherlands. Our findings show that ticks from two distinct geographical populations differ in their susceptibility to TBEV strains, which could be the result of differences in the genetic background of the tick populations.

Rescue and in vitro characterization of a divergent TBEV-Eu strain from the Netherlands

Tick-borne encephalitis virus (TBEV) may cause tick-borne encephalitis (TBE), a potential life-threatening infection of the central nervous system in humans. Phylogenetically, TBEVs can be subdivided into three main subtypes, which differ in endemic region and pathogenic potential. In 2016, TBEV was first detected in the Netherlands. One of two detected strains, referred to as Salland, belonged to the TBEV-Eu subtype, yet diverged ≥ 2% on amino acid level from other members of this subtype. Here, we report the successful rescue of this strain using infectious subgenomic amplicons and its subsequent in vitro characterization by comparison to two well-characterized TBEV-Eu strains; Neudoerfl and Hypr. In the human alveolar epithelial cell line A549, growth kinetics of Salland were comparable to the high pathogenicity TBEV-Eu strain Hypr, and both strains grew considerably faster than the mildly pathogenic strain Neudoerfl. In the human neuroblastoma cell line SK-N-SH, Salland replicated faster and to higher infectious titers than both reference strains. All three TBEV strains infected primary human monocyte-derived dendritic cells to a similar extent and interacted with the type I interferon system in a similar manner. The current study serves as the first in vitro characterization of the novel, divergent TBEV-Eu strain Salland.

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.

Complete Genome Sequencing of Tick-Borne Encephalitis Virus Directly from Clinical Samples: Comparison of Shotgun Metagenomic and Targeted Amplicon-Based Sequencing

The clinical presentation of tick-borne encephalitis virus (TBEV) infection varies from asymptomatic to severe meningoencephalitis or meningoencephalomyelitis. The TBEV subtype has been suggested as one of the most important risk factors for disease severity, but TBEV genetic characterization is difficult. Infection is usually diagnosed in the post-viremic phase, and so relevant clinical samples of TBEV are extremely rare and, when present, are associated with low viral loads. To date, only two complete TBEV genomes sequenced directly from patient clinical samples are publicly available. The aim of this study was to develop novel protocols for the direct sequencing of the TBEV genome, enabling studies of viral genetic determinants that influence disease severity. We developed a novel oligonucleotide primer scheme for amplification of the complete TBEV genome. The primer set was tested on 21 clinical samples with various viral loads and collected over a 15-year period using the two most common sequencing platforms. The amplicon-based strategy was compared to direct shotgun sequencing. Using the novel primer set, we successfully obtained nearly complete TBEV genomes (>90% of genome) from all clinical samples, including those with extremely low viral loads. Comparison of consensus sequences of the TBEV genome generated using the novel amplicon-based strategy and shotgun sequencing showed no difference. We conclude that the novel primer set is a powerful tool for future studies on genetic determinants of TBEV that influence disease severity and will lead to a better understanding of TBE pathogenesis.

Analysing the genomes of two tick-borne encephalitis viruses isolated in Hungary in 1952 and 2019

The genomes of two Tick-borne encephalitis virus (TBEV) strains were fully sequenced and compared to those of known Hungarian strains. One was a laboratory strain (KEM-1) isolated in 1952, which had gone through hundreds of passages both on Vero cell cultures and in laboratory mice, while the other was a recent isolate (2019) from questing female ticks. The laboratory strain formed a monophyletic group with the already published 4 Hungarian strains on the evolutionary tree, located relatively close to Finnish (Kumlinge) and Russian (Absettarov) strains. This KEM-1 strain was phylogenetically distantly related both to the geographically close reference strain Neudörfl and the chronologically close Czech isolates from 1953. The 2019 isolate, KEM-195 was related to TBEV isolates from Southern Slovakia and Styria, and had the longest (328 nucleotides) deletion in its 3'-non-coding region among published sequences of strains of European subtype. Our results show that decades of laboratory passage have not altered the viral genome too much and that at least two distinct branches of TBEV strains circulate in Hungary.

Targeting the NS2B-NS3 protease of tick-borne encephalitis virus with pan-flaviviral protease inhibitors

Tick-borne encephalitis (TBE) is a severe neurological disorder caused by tick-borne encephalitis virus (TBEV), a member of the Flavivirus genus. Currently, two vaccines are available in Europe against TBEV. However, TBE cases have been rising in Sweden for the past twenty years, and thousands of cases are reported in Europe, emphasizing the need for antiviral treatments against this virus. The NS2B-NS3 protease is essential for flaviviral life cycle and has been studied as a target for the design of inhibitors against several well-known flaviviruses, but not TBEV. In the present study, Compound 86, a known tripeptidic inhibitor of dengue (DENV), West Nile (WNV) and Zika (ZIKV) proteases, was predicted to be active against TBEV protease using a combination of in silico techniques. Further, Compound 86 was found to inhibit recombinant TBEV protease with an IC₅₀ = 0.92 μM in the in vitro enzymatic assay. Additionally, two more peptidic analogues were synthetized and they displayed inhibitory activities against both TBEV and ZIKV proteases. In particular, Compound 104 inhibited ZIKV protease with an IC₅₀ = 0.25 μM. These compounds represent the first reported inhibitors of TBEV protease to date and provides valuable information for the further development of TBEV as well as pan-flavivirus protease inhibitors.

High-throughput sequencing of two European strains of tick-borne encephalitis virus (TBEV), Hochosterwitz and 1993/783

Tick-borne encephalitis virus (TBEV) is a medically important arbovirus, widespread in Europe and Asia. The virus is primarily transmitted to humans and animals by bites from ticks and, in rare cases, by consumption of unpasteurized dairy products. The aim of this study was to sequence and characterize two TBEV strains with amplicon sequencing by designing overlapping primers. The amplicon sequencing, via Illumina MiSeq, covering nearly the entire TBEV genome, was successful: We retrieved and characterized the complete polyprotein sequence of two TBEV strains, Hochosterwitz and 1993/783 from Austria and Sweden, respectively. In this study the previous phylogenetic analysis of both strains was confirmed to be of the European subtypes of TBEV (TBEV-Eu) by whole genome sequencing. The Hochosterwitz strain clustered with the two strains KrM 93 and KrM 213 from South Korea, and the 1993/783 strain clustered together with the NL/UH strain from the Netherlands. Our study confirms the suitability and rapidness of the high-throughput sequencing method used to produce complete TBEV genomes from TBEV samples of high viral load giving high-molecular-weight cDNA with large overlapping amplicons.

Tick-Borne Encephalitis Virus Adaptation in Different Host Environments and Existence of Quasispecies

A highly virulent strain (Hypr) of tick-borne encephalitis virus (TBEV) was serially subcultured in the mammalian porcine kidney stable (PS) and Ixodes ricinus tick (IRE/CTVM19) cell lines, producing three viral variants. These variants exhibited distinct plaque sizes and virulence in a mouse model. Comparing the full-genome sequences of all variants, several nucleotide changes were identified in different genomic regions. Furthermore, different sequential variants were revealed to co-exist within one sample as quasispecies. Interestingly, the above-mentioned nucleotide changes found within the whole genome sequences of the new variants were present alongside the nucleotide sequence of the parental strain, which was represented as a minority quasispecies. These observations further imply that TBEV exists as a heterogeneous population that contains virus variants pre-adapted to reproduction in different environments, probably enabling virus survival in ticks and mammals.

Tick-borne encephalitis virus inhibits rRNA synthesis and host protein production in human cells of neural origin

Tick-borne encephalitis virus (TBEV), a member of the genus Flavivirus (Flaviviridae), is a causative agent of a severe neuroinfection. Recently, several flaviviruses have been shown to interact with host protein synthesis. In order to determine whether TBEV interacts with this host process in its natural target cells, we analysed de novo protein synthesis in a human cell line derived from cerebellar medulloblastoma (DAOY HTB-186). We observed a significant decrease in the rate of host protein synthesis, including the housekeeping genes HPRT1 and GAPDH and the known interferon-stimulated gene viperin. In addition, TBEV infection resulted in a specific decrease of RNA polymerase I (POLR1) transcripts, 18S and 28S rRNAs and their precursor, 45-47S pre-rRNA, but had no effect on the POLR3 transcribed 5S rRNA levels. To our knowledge, this is the first report of flavivirus-induced decrease of specifically POLR1 rRNA transcripts accompanied by host translational shut-off.

Tick-borne encephalitis virus (TBEV) is a causative agent of a severe human neuroinfection that threatens Europe and Asia. Little is known about the interaction of this neurotropic virus with neural cells, even though this may be important to better understand why or how TBEV can cause high pathogenicity in humans, especially following neural cell infection. Here, we showed that TBEV induced host translational shut-off in cells of neural origin. In addition, TBEV interfered also with the expression of host ribosomal RNAs. Interestingly, the transcriptional shut-off was documented for rRNA species transcribed by RNA polymerase I (18S rRNA, 28S rRNA and their precursor 45-47S pre-rRNA), but not for RNA polymerase III rRNA transcripts (5S rRNA). Artificial inhibition of host translation using cycloheximide resulted in the decrease of all rRNA species. Based on these data, TBEV seems to specifically target transcription of RNA polymerase I. These new findings further increase our understanding of TBEV interactions with a key target cell type.

Comparison of Schmallenberg virus sequences isolated from mammal host and arthropod vector

Schmallenberg virus (SBV) is the member of Peribunyaviridae family, which comprises pathogens of importance for human and veterinary medicine. The virus is transmitted only between animals and mainly by biting midges of the genus Culicoides. This study was performed in order to determine SBV genetic diversity and elucidate the host–vector adaptation. All three viral segments were analysed for sequence variability and phylogenetic relations. The Polish SBV strains obtained from acute infections of cattle, congenital cases in sheep, and from Culicoides midges were sequenced using Sanger and next-generation sequencing (NGS) methods. The obtained sequences were genetically similar (99.2–100% identity) to the first-detected strain BH80/11—4 from German cattle. The sampling year and origin of Polish sequences had no effect on molecular diversity of SBV. Considering all analysed Polish as well as European sequences, ovine-derived sequences were the most variable, while the midge ones were more conserved and encompassed unique substitutions located mainly in nonstructural protein S. SBV sequences isolated from Culicoides are the first submitted to GenBank and reported.

Viral Determinants of Virulence in Tick-Borne Flaviviruses

Tick-borne flaviviruses have a global distribution and cause significant human disease, including encephalitis and hemorrhagic fever, and often result in neurologic sequelae. There are two distinct properties that determine the neuropathogenesis of a virus. The ability to invade the central nervous system (CNS) is referred to as the neuroinvasiveness of the agent, while the ability to infect and damage cells within the CNS is referred to as its neurovirulence. Examination of laboratory variants, cDNA clones, natural isolates with varying pathogenicity, and virally encoded immune evasion strategies have contributed extensively to our understanding of these properties. Here we will review the major viral determinants of virulence that contribute to pathogenesis and influence both neuroinvasiveness and neurovirulence properties of tick-borne flaviviruses, focusing particularly on the envelope protein (E), nonstructural protein 5 (NS5), and the 3′ untranslated region (UTR).

Hungarian tick-borne encephalitis viruses isolated from a 0.5-ha focus are closely related to Finnish strains

Four tick-borne encephalitis virus strains were isolated from a small 0.5-ha focus over a six-year-long period (2011-2016) in Hungary. Two strains with identical genomes were isolated from Ixodes ricinus and Haemaphysalis concinna two months apart, which shows that the virus had not evolved separately in these tick species. Whole-genome sequencing of the virus revealed that the isolates differed from each other in 4 amino acids and 9 nucleotides. The calculated substitution rates indicated that the speed of genome evolution differs from habitat to habitat, and continuously changes even within the same focus. The amino acid changes affected the capsid, envelope, NS2a and NS5 genes, and one mutation each occurred in the 5' and 3' NCR as well as the premembrane, NS2a and NS5 genes. Phylogenetic analyses based on complete coding ORF sequences showed that the isolates belong to the European subtype of the virus and are closely related to the Finnish Kumlinge strains, the Bavarian isolate Leila and two isolates of Russian origin, but more distantly related to viruses from the neighbouring Central European countries. These isolates obviously have a common origin and are probably connected by migrating birds. These are the first published complete Hungarian TBEV sequences.

Identification and analysis of host proteins that interact with the 3'-untranslated region of tick-borne encephalitis virus genomic RNA

Tick-borne encephalitis virus (TBEV) causes severe neurological disease, but the pathogenetic mechanism is unclear. The conformational structure of the 3'-untranslated region (UTR) of TBEV is associated with its virulence. We tried to identify host proteins interacting with the 3'-UTR of TBEV. Cellular proteins of HEK293T cells were co-precipitated with biotinylated RNAs of the 3'-UTR of low- and high-virulence TBEV strains and subjected to mass spectrometry analysis. Fifteen host proteins were found to bind to the 3'-UTR of TBEV, four of which-cold shock domain containing-E1 (CSDE1), spermatid perinuclear RNA binding protein (STRBP), fragile X mental retardation protein (FMRP), and interleukin enhancer binding factor 3 (ILF3)-bound specifically to that of the low-virulence strain. An RNA immunoprecipitation and pull-down assay confirmed the interactions of the complete 3'-UTRs of TBEV genomic RNA with CSDE1, FMRP, and ILF3. Partial deletion of the stem loop (SL) 3 to SL 5 structure of the variable region of the 3'-UTR did not affect interactions with the host proteins, but the interactions were markedly suppressed by deletion of the complete SL 3, 4, and 5 structures, as in the high-virulence TBEV strain. Further analysis of the roles of host proteins in the neurologic pathogenicity of TBEV is warranted.

The 5' and 3' Untranslated Regions of the Flaviviral Genome

Flaviviruses are enveloped arthropod-borne viruses with a single-stranded, positive-sense RNA genome that can cause serious illness in humans and animals. The 11 kb 5′ capped RNA genome consists of a single open reading frame (ORF), and is flanked by 5′ and 3′ untranslated regions (UTR). The ORF is a polyprotein that is processed into three structural and seven non-structural proteins. The UTRs have been shown to be important for viral replication and immune modulation. Both of these regions consist of elements that are essential for genome cyclization, resulting in initiation of RNA synthesis. Genome mutation studies have been employed to investigate each component of the essential elements to show the necessity of each component and its role in viral RNA replication and growth. Furthermore, the highly structured 3′UTR is responsible for the generation of subgenomic flavivirus RNA (sfRNA) that helps the virus evade host immune response, thereby affecting viral pathogenesis. In addition, changes within the 3′UTR have been shown to affect transmissibility between vector and host, which can influence the development of vaccines.

Deep sequencing analysis of tick-borne encephalitis virus from questing ticks at natural foci reveals similarities between quasispecies pools of the virus

Every year, tick-borne encephalitis virus (TBEV) causes severe central nervous system infection in 10 000 to 15 000 people in Europe and Asia. TBEV is maintained in the environment by an enzootic cycle that requires a tick vector and a vertebrate host, and the adaptation of TBEV to vertebrate and invertebrate environments is essential for TBEV persistence in nature. This adaptation is facilitated by the error-prone nature of the virus’s RNA-dependent RNA polymerase, which generates genetically distinct virus variants called quasispecies. TBEV shows a focal geographical distribution pattern where each focus represents a TBEV hotspot. Here, we sequenced and characterized two TBEV genomes, JP-296 and JP-554, from questing Ixodes ricinus ticks at a TBEV focus in central Sweden. Phylogenetic analysis showed geographical clustering among the newly sequenced strains and three previously sequenced Scandinavian strains, Toro-2003, Saringe-2009 and Mandal-2009, which originated from the same ancestor. Among these five Scandinavian TBEV strains, only Mandal-2009 showed a large deletion within the 3′ non-coding region (NCR), similar to the highly virulent TBEV strain Hypr. Deep sequencing of JP-296, JP-554 and Mandal-2009 revealed significantly high quasispecies diversity for JP-296 and JP-554, with intact 3′NCRs, compared to the low diversity in Mandal-2009, with a truncated 3′NCR. Single-nucleotide polymorphism analysis showed that 40 % of the single-nucleotide polymorphisms were common between quasispecies populations of JP-296 and JP-554, indicating a putative mechanism for how TBEV persists and is maintained within its natural foci.

The host cell response to tick-borne encephalitis virus

Tick-borne encephalitis virus is the most prevalent autochthonous arbovirus in Europe and an important travel-associated virus. Complications of the infection could lead to lethal encephalitis in susceptible individuals. However, despite its clinical relevance and expanding geographical distribution, most of our knowledge on its pathogenesis is inferred from studies on other flaviviruses. Molecular details of the host cell response to infection are scarce leading to a poor understanding of the antiviral pathways and viral countermeasures that are critical to determine the outcome of the infection. In this work the relevant literature is reviewed and the key elements of tick-borne encephalitis virus infection of human cells are identified, which requires further investigation.

The role of the poly(A) tract in the replication and virulence of tick-borne encephalitis virus

The tick-borne encephalitis virus (TBEV) is a flavivirus transmitted to humans, usually via tick bites. The virus causes tick-borne encephalitis (TBE) in humans, and symptoms range from mild flu-like symptoms to severe and long-lasting sequelae, including permanent brain damage. It has been suggested that within the population of viruses transmitted to the mammalian host, quasispecies with neurotropic properties might become dominant in the host resulting in neurological symptoms. We previously demonstrated the existence of TBEV variants with variable poly(A) tracts within a single blood-fed tick. To characterize the role of the poly(A) tract in TBEV replication and virulence, we generated infectious clones of Torö-2003 with the wild-type (A)₃C(A)₆ sequence (Torö-6A) or with a modified (A)₃C(A)₃₈ sequence (Torö-38A). Torö-38A replicated poorly compared to Torö-6A in cell culture, but Torö-38A was more virulent than Torö-6A in a mouse model of TBE. Next-generation sequencing of TBEV genomes after passaging in cell culture and/or mouse brain revealed mutations in specific genomic regions and the presence of quasispecies that might contribute to the observed differences in virulence. These data suggest a role for quasispecies development within the poly(A) tract as a virulence determinant for TBEV in mice.

A novel virus from Macrosiphum euphorbiae with similarities to members of the family Flaviviridae

A virus with a large genome was identified in the transcriptome of the potato aphid (Macrosiphum euphorbiae) and was named Macrosiphum euphorbiae virus 1 (MeV-1). The MeV-1 genome is 22 780 nt in size, including 3' and 5' non-coding regions, with a single large ORF encoding a putative polyprotein of 7333 aa. The C-terminal region of the predicted MeV-1 polyprotein contained sequences with similarities to helicase, methyltransferase and RNA-dependent RNA polymerase (RdRp) motifs, while the N-terminal region lacked any motifs including structural proteins. Phylogenetic analysis of the helicase placed MeV-1 close to pestiviruses, while the RdRp region placed it close to pestiviruses and flaviviruses, suggesting MeV-1 has a positive-polarity ssRNA genome and is a member of the family Flaviviridae. Since the MeV-1 genome is predicted to contain a methyltransferase, a gene present typically in flaviviruses but not pestiviruses, MeV-1 is likely a member of the genus Flavivirus. MeV-1 was present in nymphal and adult stages of the aphid, aphid saliva and plant tissues fed upon by aphids. However, the virus was unable to multiply and spread in tomato plants. In addition, dsRNA, the replication intermediate of RNA viruses, was isolated from virus-infected M. euphorbiae and not from tomato plants infested with the aphid. Furthermore, nymphs laid without exposure to infected plants harboured the virus, indicating that MeV-1 is an aphid-infecting virus likely transmitted transovarially. The virus was present in M. euphorbiae populations from Europe but not from North America and was absent in all other aphid species tested.