The deduced evolution history of Omsk hemorrhagic fever virus

Inferring Evolutionary Timescale of Omsk Hemorrhagic Fever Virus

Until 2020, there were only three original complete genome (CG) nucleotide sequences of Omsk hemorrhagic fever virus (OHFV) in GenBank. For this reason, the evolutionary rate and divergence time assessments reported in the literature were based on the E gene sequences, but notably without temporal signal evaluation, such that their reliability is unclear. As of July 2022, 47 OHFV CG sequences have been published, which enables testing of temporal signal in the data and inferring unbiased and reliable substitution rate and divergence time values. Regression analysis in the TempEst software demonstrated a stronger clocklike behavior in OHFV samples for the complete open reading frame (ORF) data set (R² = 0.42) than for the E gene data set (R² = 0.11). Bayesian evaluation of temporal signal indicated very strong evidence, with a log Bayes factor of more than 5, in favor of temporal signal in all data sets. Our results based on the complete ORF sequences showed a more precise OHFV substitution rate (95% highest posterior density (HPD) interval, 9.1 × 10^-5-1.8 × 10^-4 substitutions per site per year) and tree root height (416-896 years ago) compared with previous assessments. The rate obtained is significantly higher than tick-borne encephalitis virus by at least 3.8-fold. The phylogenetic analysis and past population dynamics reconstruction revealed the declining trend of OHFV genetic diversity, but there was phylogenomic evidence that implicit virus subpopulations evolved locally and underwent an exponential growth phase.

Dates and Rates of Tick-Borne Encephalitis Virus-The Slowest Changing Tick-Borne Flavivirus

We evaluated the temporal signal and substitution rate of tick-borne encephalitis virus (TBEV) using 276 complete open reading frame (ORF) sequences with known collection dates. According to a permutation test, the TBEV Siberian subtype (TBEV-S) data set has no temporal structure and cannot be applied for substitution rate estimation without other TBEV subtypes. The substitution rate obtained suggests that the common clade of TBEV (TBEV-common), including all TBEV subtypes and louping-ill virus (LIV), is characterized by the lowest rate (1.87 × 10^-5 substitutions per site per year (s/s/y) or 1 nucleotide substitution per ORF per 4.9 years; 95% highest posterior density (HPD) interval, 1.3-2.4 × 10^-5 s/s/y) among all tick-borne flaviviruses previously assessed. Within TBEV-common, the TBEV European subtype (TBEV-E) has the lowest substitution rate (1.3 × 10^-5 s/s/y or 1 nucleotide substitution per ORF per 7.5 years; 95% HPD, 1.0-1.8 × 10^-5 s/s/y) as compared with TBEV Far-Eastern subtype (3.0 × 10^-5 s/s/y or 1 nucleotide substitution per ORF per 3.2 years; 95% HPD, 1.6-4.5 × 10^-5 s/s/y). TBEV-common representing the species tick-borne encephalitis virus diverged 9623 years ago (95% HPD interval, 6373-13,208 years). The TBEV Baikalian subtype is the youngest one (489 years; 95% HPD, 291-697 years) which differs significantly by age from TBEV-E (848 years; 95% HPD, 596-1112 years), LIV (2424 years; 95% HPD, 1572-3400 years), TBEV-FE (1936 years, 95% HPD, 1344-2598 years), and the joint clade of TBEV-S (2505 years, 95% HPD, 1700-3421 years) comprising Vasilchenko, Zausaev, and Baltic lineages.

OMSK HEMORRHAGIC FEVER VIRUS IS A TICK-BORNE ENCEPHALITIS VIRUS ADAPTED TO MUSKRAT THROUGH HOST-JUMPING

Omsk hemorrhagic fever was first described in the early 1940s and is a natural focal infection, spread exclusively in four regions of Western Siberia and associated with muskrat (Ondatra zibethicus). The etiological agent of this disease is the Omsk hemorrhagic fever virus (OHFV) which is closely related to the tick-borne encephalitis virus (TBEV), and its range entirely lies within the TBEV area. OHFV belongs to the mammalian tick-borne flaviviruses and the ecological group of arboviruses. The problem concerning the origin of OHFV remains unresolved to date. This work analyzed all nucleotide sequences of the OHFV genome obtained in the present study and available in GenBank, including the E gene fragment and the amino acid sequences of the surface glycoprotein encoded by it. The conclusions, based on the clusteron approach, suggest that OHFV originated directly from the TBEV of the Far Eastern subtype due to the host-jump phenomenon, that is, through a rapid change from an arthropod host, Ixodes persulcatus, to a rodent, O. zibethicus. The muskrat was introduced to Western Siberia in the second half of the 1930s. The peculiarities of the biology and ecology of the muskrat in the new habitat became the reason for the TBEV cross-species transmission. Calculations show that host-jumping occurred between 1931 and 1947 and accompanied a cascade of adaptive amino acid substitutions in protein E. As a result, the virus changed its transmission to contact, alimentary, and airborne routes. Based on the data obtained, OHFV would be more correctly attributed to zoonotic viruses transmitted by rodents and, accordingly, to the ecological group of roboviruses. This article is protected by copyright. All rights reserved.

Alternation between taxonomically divergent hosts is not the major determinant of flavivirus evolution

Flaviviruses display diverse epidemiological and ecological features. Tick-borne and mosquito-borne flaviviruses (TBFV and MBFV, respectively) are important human pathogens that alternate replication in invertebrate vectors and vertebrate hosts. The Flavivirus genus also includes insect-specific viruses (ISFVs) and viruses with unknown invertebrate hosts. It is generally accepted that viruses that alternate between taxonomically different hosts evolve slowly and that the evolution of MBFVs and TBFVs is dominated by strong constraints, with limited episodes of positive selection. We exploited the availability of flavivirus genomes to test these hypotheses and to compare their rates and patterns of evolution. We estimated the substitution rates of CFAV and CxFV (two ISFVs) and, by taking into account the time-frame of measurement, compared them with those of other flaviviruses. Results indicated that CFAV and CxFV display relatively different substitution rates. However, these data, together with estimates for single-host members of the Flaviviridae family, indicated that MBFVs do not display relatively slower evolution. Conversely, TBFVs displayed some of lowest substitution rates among flaviviruses. Analysis of selective patterns over longer evolutionary time-frames confirmed that MBFVs evolve under strong purifying selection. Interestingly, TBFVs and ISFVs did not show extremely different levels of constraint, although TBFVs alternate among hosts, whereas ISFVs do not. Additional results showed that episodic positive selection drove the evolution of MBFVs, despite their high constraint. Positive selection was also detected on two branches of the TBFVs phylogeny that define the seabird clade. Thus, positive selection was much more common during the evolution of arthropod-borne flaviviruses than previously thought. Overall, our data indicate that flavivirus evolutionary patterns are complex and most likely determined by multiple factors, not limited to the alternation between taxonomically divergent hosts. The frequency of both positive and purifying selection, especially in MBFVs, suggests that a minority of sites in the viral polyprotein experience weak constraint and can evolve to generate new viral phenotypes and possibly promote adaptation to new hosts.

Molecular variability and genetic structure of Omsk hemorrhagic fever virus, based on analysis of the complete genome sequences

Omsk hemorrhagic fever virus (OHFV) is the etiological agent of Omsk hemorrhagic fever, a disease described in the 1940s in Western Siberia. However, until now, it has been represented in GenBank by just four complete genome sequences, which do not reflect the real genetic diversity of the virus in nature. In this study, we analyzed the molecular variability and genetic structure of OHFV based on 20 complete genome sequences, fifteen of which were obtained for the first time. All these sequences belong to virus strains isolated at different times from three regions of Western Siberia. The results suggest that the genetic diversity of OHFV is significantly wider than previously thought and is represented by at least three subtypes, rather than two. This broadens our understanding of the evolutionary history of OHFV. Also, it is argued that the OHFV reference strain Bogoluvovska (NC_005062) is actually a Kubrin strain and that either cross-contamination or a laboratory error was the cause of this.

Complete mitogenome of the ixodid tick Dermacentor reticulatus (Acari: Ixodida)

Here, we present the complete mitochondrial DNA sequence of Dermacentor reticulatus. The mitogenome is 14,806 bp and contains 13 protein-coding, 2 rRNA, and 22 tRNA genes, along with 2 control regions. Dermacentor reticulatus mitogenome has the common mitochondrial gene order of Metastriata ticks. It is phylogenetically close to the mitogenomes of Dermacentor ticks, of which D. everestanus mitogenome is the closest with 85.7% similarity. These data provide insights into the phylogenetic relations among Dermacentor ticks.

Exploiting the Legacy of the Arbovirus Hunters

In recent years, it has become evident that a generational gap has developed in the community of arbovirus research. This apparent gap is due to the dis-investment of training for the next generation of arbovirologists, which threatens to derail the rich history of virus discovery, field epidemiology, and understanding of the richness of diversity that surrounds us. On the other hand, new technologies have resulted in an explosion of virus discovery that is constantly redefining the virosphere and the evolutionary relationships between viruses. This paradox presents new challenges that may have immediate and disastrous consequences for public health when yet to be discovered arboviruses emerge. In this review we endeavor to bridge this gap by providing a historical context for the work being conducted today and provide continuity between the generations. To this end, we will provide a narrative of the thrill of scientific discovery and excitement and the challenges lying ahead.

Tick-Borne Viruses

Ticks are important vectors for the transmission of pathogens including viruses. The viruses carried by ticks also known as tick-borne viruses (TBVs), contain a large group of viruses with diverse genetic properties and are concluded in two orders, nine families, and at least 12 genera. Some members of the TBVs are notorious agents causing severe diseases with high mortality rates in humans and livestock, while some others may pose risks to public health that are still unclear to us. Herein, we review the current knowledge of TBVs with emphases on the history of virus isolation and identification, tick vectors, and potential pathogenicity to humans and animals, including assigned species as well as the recently discovered and unassigned species. All these will promote our understanding of the diversity of TBVs, and will facilitate the further investigation of TBVs in association with both ticks and vertebrate hosts.

Dermacentor reticulatus: a vector on the rise

Dermacentor reticulatus is a hard tick species with extraordinary biological features. It has a high reproduction rate, a rapid developmental cycle, and is also able to overcome years of unfavourable conditions. Dermacentor reticulatus can survive under water for several months and is cold-hardy even compared to other tick species. It has a wide host range: over 60 different wild and domesticated hosts are known for the three active developmental stages. Its high adaptiveness gives an edge to this tick species as shown by new data on the emergence and establishment of D. reticulatus populations throughout Europe. The tick has been the research focus of a growing number of scientists, physicians and veterinarians. Within the Web of Science database, more than a fifth of the over 700 items published on this species between 1897 and 2015 appeared in the last three years (2013–2015). Here we attempt to synthesize current knowledge on the systematics, ecology, geographical distribution and recent spread of the species and to highlight the great spectrum of possible veterinary and public health threats it poses. Canine babesiosis caused by Babesia canis is a severe leading canine vector-borne disease in many endemic areas. Although less frequently than Ixodes ricinus, D. reticulatus adults bite humans and transmit several Rickettsia spp., Omsk haemorrhagic fever virus or Tick-borne encephalitis virus. We have not solely collected and reviewed the latest and fundamental scientific papers available in primary databases but also widened our scope to books, theses, conference papers and specialists colleagues’ experience where needed. Besides the dominant literature available in English, we also tried to access scientific literature in German, Russian and eastern European languages as well. We hope to inspire future research projects that are necessary to understand the basic life-cycle and ecology of this vector in order to understand and prevent disease threats. We conclude that although great strides have been made in our knowledge of the eco-epidemiology of this species, several gaps still need to be filled with basic research, targeting possible reservoir and vector roles and the key factors resulting in the observed geographical spread of D. reticulatus.

Comparative Metagenomic Profiling of Symbiotic Bacterial Communities Associated with Ixodes persulcatus, Ixodes pavlovskyi and Dermacentor reticulatus Ticks

Ixodes persulcatus, Ixodes pavlovskyi, and Dermacentor reticulatus ticks inhabiting Western Siberia are responsible for the transmission of a number of etiological agents that cause human and animal tick-borne diseases. Because these ticks are abundant in the suburbs of large cities, agricultural areas, and popular tourist sites and frequently attack people and livestock, data regarding the microbiomes of these organisms are required. Using metagenomic 16S profiling, we evaluate bacterial communities associated with I. persulcatus, I. pavlovskyi, and D. reticulatus ticks collected from the Novosibirsk region of Russia. A total of 1214 ticks were used for this study. DNA extracted from the ticks was pooled according to tick species and sex. Sequencing of the V3-V5 domains of 16S rRNA genes was performed using the Illumina Miseq platform. The following bacterial genera were prevalent in the examined communities: Acinetobacter (all three tick species), Rickettsia (I. persulcatus and D. reticulatus) and Francisella (D. reticulatus). B. burgdorferi sensu lato and B. miyamotoi sequences were detected in I. persulcatus and I. pavlovskyi but not in D. reticulatus ticks. The pooled samples of all tick species studied contained bacteria from the Anaplasmataceae family, although their occurrence was low. DNA from A. phagocytophilum and Candidatus Neoehrlichia mikurensis was first observed in I. pavlovskyi ticks. Significant inter-species differences in the number of bacterial taxa as well as intra-species diversity related to tick sex were observed. The bacterial communities associated with the I. pavlovskyi ticks displayed a higher biodiversity compared with those of the I. persulcatus and D. reticulatus ticks. Bacterial community structure was also diverse across the studied tick species, as shown by permutational analysis of variance using the Bray-Curtis dissimilarity metric (p = 0.002). Between-sex variation was confirmed by PERMANOVA testing in I. persulcatus (p = 0.042) and I. pavlovskyi (p = 0.042) ticks. Our study indicated that 16S metagenomic profiling could be used for rapid assessment of the occurrence of medically important bacteria in tick populations inhabiting different natural biotopes and therefore the epidemic danger of studied foci.