High-Throughput Sequencing Reveals Three Rhabdoviruses Persisting in the IRE/CTVM19 Cell Line.

Cell cultures derived from ticks have become a commonly used tool for the isolation and study of tick-borne pathogens and tick biology. The IRE/CTVM19 cell line, originating from embryos of Ixodes ricinus, is one such line. Previously, reovirus-like particles, as well as sequences with similarity to rhabdoviruses and iflaviruses, were detected in the IRE/CTVM19 cell line, suggesting the presence of multiple persisting viruses. Subsequently, the full genome of an IRE/CTVM19-associated rhabdovirus was recovered from a cell culture during the isolation of the Alongshan virus. In the current work, we used high-throughput sequencing to describe a virome of the IRE/CTVM19 cell line. In addition to the previously detected IRE/CTVM19-associated rhabdovirus, two rhabdoviruses were detected: Chimay rhabdovirus and Norway mononegavirus 1. In the follow-up experiments, we were able to detect both positive and negative RNA strands of the IRE/CTVM19-associated rhabdovirus and Norway mononegavirus 1 in the IRE/CTVM19 cells, suggesting their active replication in the cell line. Passaging attempts in cell lines of mammalian origin failed for all three discovered rhabdoviruses.

Viromes of Tabanids from Russia.

Advances in sequencing technologies and bioinformatics have greatly enhanced our knowledge of virus biodiversity. Currently, the viromes of hematophagous invertebrates, such as mosquitoes and ixodid ticks, are being actively studied. Tabanidae (Diptera) are a widespread family, with members mostly known for their persistent hematophagous behavior. They transmit viral, bacterial, and other pathogens, both biologically and mechanically. However, tabanid viromes remain severely understudied. In this study, we used high-throughput sequencing to describe the viromes of several species in the Hybomitra, Tabanus, Chrysops, and Haematopota genera, which were collected in two distant parts of Russia: the Primorye Territory and Ryazan Region. We assembled fourteen full coding genomes of novel viruses, four partial coding genomes, as well as several fragmented viral sequences, which presumably belong to another twelve new viruses. All the discovered viruses were tested for their ability to replicate in mammalian porcine embryo kidney (PEK), tick HAE/CTVM8, and mosquito C6/36 cell lines. In total, 16 viruses were detected in at least one cell culture after three passages (for PEK and C6/36) or 3 weeks of persistence in HAE/CTVM8. However, in the majority of cases, qPCR showed a decline in virus load over time.

Quantitative Polymerase Chain Reaction System for Alongshan Virus Detection.

The recently discovered Jingmenvirus group includes viruses with a segmented genome, RNA of a positive polarity, and several proteins with distant homology to the proteins of the members of the genus Orthoflavivirus. Some Jingmenvirus group members, namely the Alongshan virus (ALSV) and Jingmen tick virus, are reported to be tick-borne human pathogens that can cause a wide variety of symptoms. The ALSV is widely distributed in Eurasia, yet no reliable assay that can detect it exists. We describe a qPCR system for ALSV detection. Our data showed that this system can detect as little as 104 copies of the ALSV in a sample. The system showed no amplification of the common tick-borne viruses circulating in Eurasia, i.e., the Yanggou tick virus-which is another Jingmenvirus group member-or some known members of the genus Orthoflavivirus. The qPCR system was tested and had no nonspecific signal for the Ixodes ricinus, I. persulcatus, Dermacentor reticulatus, D. marginatus, Haemaphysalis concinna, and H. japonica ticks. The qPCR system had no nonspecific signal for human and sheep serum as well. Overall, the qPCR system described here can be used for reliable and quantitative ALSV detection.

Conserved Sequences in the 5' and 3' Untranslated Regions of Jingmenvirus Group Representatives.

The Jingmenvirus group (JVG), with members such as Jingmen tick virus (JMTV), Alongshan virus (ALSV), Yanggou tick virus (YGTV), and Takachi virus (TAKV), is drawing attention due to evidence of it causing disease in humans and its unique genome architecture. In the current work, complete untranslated regions (UTRs) of four strains of ALSV and eight strains of YGTV were obtained. An analysis of these sequences, as well as JVG sequences from GenBank, uncovered several regions within viral UTRs that were highly conserved for all the segments and viruses. Bioinformatics predictions suggested that the UTRs of all the segments of YGTV, ALSV, and JMTV could form similar RNA structures. The most notable feature of these structures was a stable stem-loop with one (5' UTR) or two (3' UTR) AAGU tetraloops on the end of a hairpin.

Distribution and Characterisation of Tick-Borne Flavi-, Flavi-like, and Phenuiviruses in the Chelyabinsk Region of Russia.

In this work, we presented data from a two-year study of flavi-, flavi-like, and phenuiviruses circulation in the population of ixodid ticks in the Chelyabinsk region. We isolated three tick-borne encephalitis virus (TBEV) strains from I. persulcatus, which was not detected in the ticks of the genus Dermacentor. The virus prevalence ranged from 0.66% to 2.28%. The Yanggou tick virus (YGTV) is widespread in steppe and forest-steppe zones and is mainly associated with ticks of the genus Dermacentor. We isolated 26 strains from D. reticulatus, D. marginatus, and I. persulcatus ticks in the HAE/CTVM8 tick cell line. The virus prevalence ranged from 1.58% to 4.18% in D. reticulatus, ranged from 0.78% to 3.93% in D. marginatus, and was 0.66% in I. persulcatus. There was combined focus of TBEV and YGTV in the territory of the Chelyabinsk region. The Alongshan virus (ALSV) was found to be associated with I. persulcatus ticks and is spread in forest zone. We detected 12 amplicons and isolated 7 strains of ALSV in tick cells. The virus prevalence ranged from 1.13% to 6.00%. The phlebovirus Gomselga and unclassified phenuivirus Stavropol were associated with I. persulcatus and D. reticulatus ticks, respectively. Virus prevalence of the unclassified phenuivirus Stavropol in the Chelyabinsk region is lower than that in neighbouring regions.

Integrated Jingmenvirus Polymerase Gene in Ixodes ricinus Genome.

Members of the jingmenviruses group have been found in arthropods and mammals on all continents except Australia and Antarctica. Two viruses of this group were isolated from patients with fever after a tick bite. Using a nested RT-PCR assay targeting a jingmenvirus polymerase gene fragment, we screened ticks collected in seven regions of Russia and found that the abundant jingmenvirus-positive were of Ixodes ricinus species, with the prevalence ranging from 19.8% to 34.3%. In all cases, DNase/RNase treatment suggested that the detected molecule was DNA and subsequent next generation sequencing (NGS) proved that the viral polymerase gene was integrated in the I. ricinus genome. The copy number of the integrated polymerase gene was quantified by qPCR relative to the ITS2 gene and estimated as 1.32 copies per cell. At least three different genetic variants of the integrated polymerase gene were found in the territory of Russia. Phylogenetic analysis of the integrated jingmenvirus polymerase gene showed the highest similarity with the sequence of the correspondent gene obtained in Serbia from I. ricinus.

Differentiation of Laboratory-Obtained Ixodes ricinus × Ixodes persulcatus Hybrid Ticks: Selection of Suitable Genes.

Ixodes ricinus and Ixodes persulcatus ticks are the main vectors of tick-borne encephalitis virus and some bacterial pathogens. The regions where these tick species live overlap, forming large sympatric areas. It has previously been shown that these tick species have no morphological barrier, and interspecies crossing is possible with the appearance of sterile hybrids. It has also been shown that hybrid larvae and nymphs can be differentiated using discriminant functions based on a set of morphological features. However, such an approach is laborious and rather ineffective with adult ticks, making a molecular approach necessary. In the current work, we tested the ability of different systems to differentiate laboratory-obtained hybrid ticks. Our data suggest that commonly used primer sets that target rRNA are unsuitable for hybrid tick determination, likely due to the rRNA region being linked with the X chromosome in I. ricinus and I. persulcatus ticks. We tested several primer sets targeting different non rRNA genes to assess their ability to determine hybrids. The best primer set, Toll_R, targeting the putative Toll gene, showed little to no bias when used for DNA amplification from hybrid ticks. Thus, Toll gene can be further used for hybrid detection.

Virome of Three Termite Species from Southern Vietnam.

Modern metagenomic approaches enable the effective discovery of novel viruses in previously unexplored organisms. Termites are significant ecosystem converters and influencers. As with the majority of tropical forest insects, termites are studied insufficiently, and termite virome remains especially understudied. Here, we studied the virome of lichenophagous and mycophagous termites (Hospitalitermes bicolor, Macrotermes carbonarius and Odontotermes wallonensis) collected in the Cat Tien National Park (Vietnam). We assembled four full genomes of novel viruses related to Solemoviridae, Lispiviridae, Polycipiviridae and Kolmioviridae. We also found several contigs with relation to Chuviridae and Deltaflexiviridae that did not correspond to complete virus genomes. All the novel viruses clustered phylogenetically with previously identified viruses of the termites. Deltaflexi-like contigs were identified in the fungi-cultivating M. carbonarius and showed homology with viruses recently discovered in the edible basidiomycete mushrooms.

Possible Arbovirus Found in Virome of Melophagus ovinus.

Members of the Lipopteninae subfamily are blood-sucking ectoparasites of mammals. The sheep ked (Melophagus ovinus) is a widely distributed ectoparasite of sheep. It can be found in most sheep-rearing areas and can cause skin irritation, restlessness, anemia, weight loss and skin injuries. Various bacteria and some viruses have been detected in M. ovinus; however, the virome of this ked has never been studied using modern approaches. Here, we study the virome of M. ovinus collected in the Republic of Tuva, Russia. In our research, we were able to assemble full genomes for five novel viruses, related to the Rhabdoviridae (Sigmavirus), Iflaviridae, Reoviridae and Solemoviridae families. Four viruses were found in all five of the studied pools, while one virus was found in two pools. Phylogenetically, all of the novel viruses clustered together with various recently described arthropod viruses. All the discovered viruses were tested on their ability to replicate in the mammalian porcine embryo kidney (PEK) cell line. Aksy-Durug Melophagus sigmavirus RNA was detected in the PEK cell line cultural supernate after the first, second and third passages. Such data imply that this virus might be able to replicate in mammalian cells, and thus, can be considered as a possible arbovirus.

Geographical and Tick-Dependent Distribution of Flavi-Like Alongshan and Yanggou Tick Viruses in Russia.

The genus Flavivirus includes related, unclassified segmented flavi-like viruses, two segments of which have homology with flavivirus RNA-dependent RNA polymerase NS5 and RNA helicase-protease NS3. This group includes such viruses as Jingmen tick virus, Alongshan virus, Yanggou tick virus and others. We detected the Yanggou tick virus in Dermacentor nuttalli and Dermacentor marginatus ticks in two neighbouring regions of Russia. The virus prevalence ranged from 0.5% to 8.0%. We detected RNA of the Alongshan virus in 44 individuals or pools of various tick species in eight regions of Russia. The virus prevalence ranged from 0.6% to 7.8%. We demonstrated the successful replication of the Yanggou tick virus and Alongshan virus in IRE/CTVM19 and HAE/CTVM8 tick cell lines without a cytopathic effect. According to the phylogenetic analysis, we divided the Alongshan virus into two groups: an Ixodes persulcatus group and an Ixodes ricinus group. In addition, the I. persulcatus group can be divided into European and Asian subgroups. We found amino acid signatures specific to the I. ricinus and I. persulcatus groups and also distinguished between the European and Asian subgroups of the I. persulcatus group.

Isolation and Characterisation of Alongshan Virus in Russia.

In recent decades, many new flavi-like viruses have been discovered predominantly in different invertebrates and, as was recently shown, some of them may cause disease in humans. The Jingmenvirus (JMV) group holds a special place among flaviviruses and flavi-like viruses because they have a segmented ssRNA(+) genome. We detected Alongshan virus (ALSV), which is a representative of the JMV group, in ten pools of adult Ixodes persulcatus ticks collected in two geographically-separated Russian regions. Three of the ten strains were isolated in the tick cell line IRE/CTVM19. One of the strains persisted in the IRE/CTVM19 cells without cytopathic effect for three years. Most ALSV virions purified from tick cells were spherical with a diameter of approximately 40.5 nm. In addition, we found smaller particles of approximately 13.1 nm in diameter. We obtained full genome sequences of all four segments of two of the isolated ALSV strains, and partial sequences of one segment from the third strain. Phylogenetic analysis on genome segment 2 of the JMV group clustered our novel strains with other ALSV strains. We found evidence for the existence of a novel upstream open reading frame in the glycoprotein-coding segment of ALSV and other members of the JMV group.

Evaluation of the population heterogeneity of TBEV laboratory variants using high-throughput sequencing.

We studied minor variants within two tick-borne encephalitis virus (TBEV) populations with a common ancestor: the mouse brain-adapted variant EK-328c and the tick-adapted variant M. High-throughput sequencing with custom amplicons from RT-PCR viral RNA was performed on Illumina MiSeq 2*250 paired-end v2 chemistry. Using the LowFreq program (default settings) and Sanger-sequenced consensus as a reference, variants with an abundance of 1 % and above within the studied populations were identified. Using the obtained data in the context of our previous studies, we concluded that TBEV variants, which are different from the major population phenotype and can become a major part of the viral population under favourable environmental conditions, can exist at abundances of less than 1 % in the long-term. The comparison of our data with the literature allowed us to conclude that the laboratory variant EK-328c and variant M have similar SNV counts to TBEV variants from natural populations and some fast-evolving RNA viruses.

Properties of the tick-borne encephalitis virus population during persistent infection of ixodid ticks and tick cell lines.

Tick-borne encephalitis virus (TBEV) is the causative agent of tick-borne encephalitis (TBE), a vector-borne zoonotic neuroinfection. For successful circulation in natural foci the virus has to survive in the vector for a long period of time. Information about the effect of long-term infection of ticks on properties of the viral population is of great importance. In recent years, changes in the eco-epidemiology of TBEV due to changes in distribution of ixodid ticks have been observed. These changes in TBEV-endemic areas could result in a shift of the main tick vector species, which in turn may lead to changes in properties of the virus. In the present study we evaluated the selective pressure on the TBEV population during persistent infection of various species of ticks and tick cell lines. TBEV effectively replicated and formed persistent infection in ticks and tick cell lines of the vector species (Ixodes spp.), potential vectors (Dermacentor spp.) and non-vector ticks (Hyalomma spp.). During TBEV persistence in Ixodes and Dermacentor ticks, properties of the viral population remained virtually unchanged. In contrast, persistent TBEV infection of tick cell lines from both vector and non-vector ticks favoured selection of viral variants with low neuroinvasiveness for laboratory mice and substitutions in the E protein that increased local positive charge of the virion. Thus, selective pressure on viral population may differ in ticks and tick cell lines during persistent infection. Nevertheless, virus variants with properties of the original strain adapted to mouse CNS were not eliminated from the viral population during long-term persistence of TBEV in ticks and tick cell lines.

Lethal Experimental Tick-Borne Encephalitis Infection: Influence of Two Strains with Similar Virulence on the Immune Response.

Tick-borne encephalitis virus (TBEV) is a tick-transmitted arbovirus that causes serious diseases in humans in Europe and Northern Asia. About 6000-10,000 cases are registered annually, and one-third of them lead to sequela with different degrees of severity. Two TBEV strains (Absettarov and EK-328) similar in virulence rate in laboratory mice were used to study pathogenesis and immune response upon lethal infection in mice. The strains differed in the dynamics of appearance of virus, IFNs and other cytokines in blood of mice, and ability to induce a cytokine storm in the terminal stages of disease and a non-sterile immunity. Moreover, the TBEV strains differed in characteristics of their interactions with DCs: level of reproduction in these cells, virus dose triggering IFN-α production, and impact on DCs' maturation. Infection of DCs with Absettarov strain led to IFN-α induction only at high multiplicity of infection (MOI), and an increased amount of the mature DCs with high adhesion activity and low-level of MHCII positive cells. While reproduction of the EK-328 strain in DCs was less efficient, a low dose of the virus induced IFN-α production and stimulated maturation of DCs with relatively low adhesive capacity, but with the high percentage of cells expressing MHCII molecules. Thus, the studied strains differed significantly in the impact on DCs' maturation and antigen presentation to CD4+ lymphocytes. Injection of low (103 PFU) and high (106 PFU) doses of both TBEV strains caused a lethal infection in mice. At the same time, the dose of the virus in the inoculum, regardless of the strain properties, affected the following virulence characteristics: the time of virus appearance in brain (day 4-5 vs. day 1 p.i.), time of IFN-α appearance in blood (10 h vs. 5 h p.i.), concentration of IFN-α in blood, and induction of IFN-α during infection of DCs. Therefore, virulent TBEV strains during lethal infection can interact differently with the host immune system, and the infectious dose has an impact on both: virus spread in the infected organism and immune response activation.