Molecular identification of novel phlebovirus sequences in European ticks

Sequences Related to Chimay Rhabdovirus Are Widely Distributed in Ixodes ricinus Ticks across England and Wales

Ticks are the main arthropod vector of pathogens to humans and livestock in the British Isles. Despite their role as a vector of disease, many aspects of tick biology, ecology, and microbial association are poorly understood. To address this, we investigated the composition of the microbiome of adult and nymphal Ixodes ricinus ticks. The ticks were collected on a dairy farm in Southwest England and RNA extracted for whole genome sequencing. Sequences were detected from a range of microorganisms, particularly tick-associated viruses, bacteria, and nematodes. A majority of the viruses were attributed to phlebo-like and nairo-like virus groups, demonstrating a high degree of homology with the sequences present in I. ricinus from mainland Europe. A virus sharing a high sequence identity with Chimay rhabdovirus, previously identified in ticks from Belgium, was detected. Further investigations of I. ricinus ticks collected from additional sites in England and Wales also identified Chimay rhabdovirus viral RNA with varying prevalence in all tick populations. This suggests that Chimay rhabdovirus has a wide distribution and highlights the need for an extended exploration of the tick microbiome in the United Kingdom (UK).

Phlebovirus diversity in ticks from livestock in arid ecologies in Kenya

Phleboviruses are emerging pathogens of public health importance. However, their association with ticks is poorly described, particularly in Africa. Here, adult ticks infesting cattle, goats and sheep were collected in two dryland pastoralist ecosystems of Kenya (Baringo and Kajiado counties) and were screened for infection with phleboviruses. Ticks mainly belonged to the species Rhipicephalus appendiculatus, Hyalomma impeltatum, and Hyalomma rufipes. A fragment of the RNA-dependent RNA polymerase (RdRp) gene was identified in thirty of 671 tick pools, of which twenty-nine were from livestock sampled in Baringo county. Phylogenetic analyses revealed that twenty-five sequences were falling in three clades within the group of tick-associated phleboviruses. The sequences of the three clades showed nucleotide distances 8%, 19% and 22%, respectively, to previously known viruses suggesting that these sequence fragments may belong to three distinct viruses. Viruses of the group of tick-associated phleboviruses have been found in several countries and continents but so far have not been associated with disease in humans or animals. In addition, five sequences were found to group with the sandfly-associated phleboviruses Bogoria virus, Perkerra virus and Ntepes virus recently detected in the same region. Further studies are needed to investigate the transmission and maintenance cycles of these viruses, as well as to assess their potential to infect vertebrates.

Forest Diversity Reduces the Prevalence of Pathogens Transmitted by the Tick Ixodes ricinus

Tick-borne diseases represent the majority of vector-borne human diseases in Europe, with Ixodes ricinus, mostly present in forests, as the main vector. Studies show that vertebrate hosts diversification would decrease the prevalence of these pathogens. However, it is not well known whether habitat diversity can have similar impact on ticks and their infection rates. We measured the presence and abundance of different stages of I. ricinus, and the prevalence of associated pathogens in a large-scale forest experiment in which we manipulated tree diversity and moisture level. We showed that larval abundance was influenced by tree species identity, with larvae being more present in pine plots than in oak plots, while nymph abundance increased with canopy tree density. The proportion of Borrelia burgdorferi s.l.-infected nymphs decreased with increasing tree diversity. Our findings suggest that tree overstorey composition, structure and diversity, can affect tick abundance and pathogen prevalence. They support the idea that forest habitats may have “diluting” or “amplifying” effects on tick-borne diseases with direct relevance for human health.

Phlebovirus sequences detected in ticks collected in Russia: Novel phleboviruses, distinguishing criteria and high tick specificity

Phlebovirus is an abundant and rather heterogeneous genus within the Phenuiviridae family (order Bunyavirales). The genus Phlebovirus is divided into two antigenic complexes, which also correspond to the main vector: sandflies/mosquitoes and ticks. Previously, only sandfly/mosquito-borne phleboviruses were associated with human disease, such as Rift Valley fever virus, Toscana virus, Sicilian and Naples Sandfly fever viruses and others. Until recently, tick-borne phleboviruses were not considered as human pathogens. After the discovery of severe fever with thrombocytopenia syndrome, interest to tick-borne phleboviruses has increased dramatically. In the last decade, many novel phleboviruses have been reported in different regions. Despite this, the diversity, ecology and pathogenicity of these viruses still remain obscure. The aim of this work was to study the diversity of phleboviruses in ticks collected in several regions of Russia. We used pan-phlebovirus RT-PCR assays based on multiple degenerate primers targeting the polymerase gene fragment. Arthropod specimens were collected from 2005 to 2018. A total of 5901 Ixodidae ticks combined into 1116 pools were screened. A total of 160 specific amplicons were produced. In three cases RT-PCR assays amplified two distinct viruses from same tick pools. Direct sequencing of amplicons and subsequent phylogenetic analysis revealed twelve representatives of divergent phlebovirus groups. Based on the distribution of pairwise nucleotide sequence identity values, a cut-off (88%) was suggested to distinguish tick-borne phleboviruses. According to this provisional criterion, two viruses found here could be termed novel, while ten viruses have been described in previous studies. Detected phleboviruses demonstrated almost perfect specificity to a tick species or, at least, a genus. The same pattern was observed for tick-borne phleboviruses found in different studies around the world. Viruses that grouped together on a phylogenetic tree and differed less than this sequence identity threshold suggested above were hosted by ticks from the same genus.

Novel Tick Phlebovirus Genotypes Lacking Evidence for Vertebrate Infections in Anatolia and Thrace, Turkey

We screened ticks and human clinical specimens to detect and characterize tick phleboviruses and pathogenicity in vertebrates. Ticks were collected at locations in Istanbul (Northwest Anatolia, Thrace), Edirne, Kırklareli, and Tekirdağ (Thrace), Mersin (Mediterranean Anatolia), Adiyaman and Şanlıurfa (Southeastern Anatolia) provinces from 2013-2018 and were analyzed following morphological identification and pooling. Specimens from individuals with febrile disease or meningoencephalitic symptoms of an unknown etiology were also evaluated. The pools were screened via generic tick phlebovirus amplification assays and products were sequenced. Selected pools were used for cell culture and suckling mice inoculations and next generation sequencing (NGS). A total of 7492 ticks were screened in 609 pools where 4.2% were positive. A phylogenetic sequence clustering according to tick species was observed. No human samples were positive. NGS provided near-complete viral replicase coding sequences in three pools. A comprehensive analysis revealed three distinct, monophyletic virus genotypes, comprised of previously-described viruses from Anatolia and the Balkans, with unique fingerprints in conserved amino acid motifs in viral replicase. A novel tick phlebovirus group was discovered circulating in the Balkans and Turkey, with at least three genotypes or species. No evidence for replication in vertebrates or infections in clinical cases could be demonstrated.

Field samplings of Ixodes ricinus ticks from a tick-borne encephalitis virus micro-focus in Northern Zealand, Denmark

In 2008-2009 a tick-borne encephalitis virus (TBEV) micro-focus was detected in Northern Zealand, Denmark. No new cases of TBE with an epidemiological link to Northern Zealand has been reported since. Here we undertook to investigate Ixodes ricinus ticks from this endemic micro-focus in 2016 and 2017. In addition to TBEV, I. ricinus ticks may host other pathogens that include Borrelia spp., Babesia spp., Rickettsia spp. and Neoehrlichia mikurensis, together with various endosymbiont microorganisms. To detect multiple organisms we used a metagenomics PanVirus microarray and next-generation sequencing to examine the persistence and evolution of other emerging viruses, bacteria and parasites. Here we report the rise and fall of the Danish TBEV micro-focus in Northern Zealand. However, we identify for the first time in Danish I. ricinus ticks the presence of Uukuniemi virus in addition to a tick-borne phlebovirus and a range of bacteria.

Geographic dispersal and genetic diversity of tick-borne phleboviruses (Phenuiviridae, Phlebovirus) as revealed by the analysis of L segment sequences

The large diversity of new tick-borne phleboviruses, and the negative impacts of the virulent viruses on human/animal health have led to a growing interest in their analysis. In this report, new insights are brought out into the diversity of putative phleboviruses circulating in Portugal (both the continental territory and the islands of São Miguel, in the Azores, and Madeira), as well as in the Spanish western regions of Extremadura and Castilla and León. Phlebovirus sequences were frequently detected (L-segment) from both questing and feeding ticks, but especially in Rhipicephalus sanguineus sensu lato (s.l.) specimens. These sequences were detected in adult ticks, as well as nymphs and eggs, supporting the hypothesis of viral maintenance by vertical transmission. Though multiple genetic groups could be identified in phylogenetic trees (AnLuc, KarMa, RiPar virus 1, and Spanish group 1 and 2), all the sequences from Portugal and Spain shared common ancestry with other viral sequence obtained from samples collected over a large geographic coverage. Spatiotemporal analysis placed Middle-East as the geographic origin of the most recent common ancestor (MRCA) of all phleboviruses analysed in the present study. More recent viral transitions might include migrations from Spain to continental Portugal, and from there to the Portuguese Islands. Our findings suggest that the time of the MRCA of phleboviruses was dated around 225 years ago [95% HPD: 124-387 year before the last sampling date].

Viral population analysis of the taiga tick, Ixodes persulcatus, by using Batch Learning Self-Organizing Maps and BLAST search

Ticks transmit a wide range of viral, bacterial, and protozoal pathogens, which are often zoonotic. Several novel tick-borne viral pathogens have been reported during the past few years. The aim of this study was to investigate a diversity of tick viral populations, which may contain as-yet unidentified viruses, using a combination of high throughput pyrosequencing and Batch Learning Self-Organizing Map (BLSOM) program, which enables phylogenetic estimation based on the similarity of oligonucleotide frequencies. DNA/cDNA prepared from virus-enriched fractions obtained from Ixodes persulcatus ticks was pyrosequenced. After de novo assembly, contigs were cataloged by the BLSOM program. In total 41 different viral families and order including those previously associated with human and animal diseases such as Bunyavirales, Flaviviridae, and Reoviridae, were detected. Therefore, our strategy is applicable for viral population analysis of other arthropods of medical and veterinary importance, such as mosquitos and lice. The results lead to the contribution to the prediction of emerging tick-borne viral diseases. A sufficient understanding of tick viral populations will also empower to analyze and understand tick biology including vector competency and interactions with other pathogens.

Infection of newly identified phleboviruses in ticks and wild animals in Hokkaido, Japan indicating tick-borne life cycles

Recent discoveries of tick-borne pathogens have raised public health concerns on tick-borne infectious diseases and emphasize the need to assess potential risks of unrecognized tick-borne pathogens. First, to determine the existence of tick-borne phleboviruses (TBPVs), genetic surveillance of phleboviruses in ticks was conducted mainly in Hokkaido, the northernmost island in Japan from 2013 to 2015. Genes of two TBPVs, previously reported as Mukawa virus (MKWV) and a newly identified relative of MKWV, Kuriyama virus (KURV), were detected and the viruses were isolated from Ixodes persulcatus collected in Hokkaido, but not in I. persulcatus collected from other areas of Japan. These viruses were phylogenetically and antigenically similar to each other. Next, to investigate the infection of MKWV in mammals, serum samples from wildlife captured in Hokkaido from 2007 to 2011 were used for serological screening. Neutralizing antibodies against MKWV were detected in both Yezo-deer (Cervus nippon yesoensis) (2/50) and raccoons (Procyon lotor) (16/64). However, no infectious MKWV was recovered from laboratory mice in experimental infections, though viral RNAs were detected in their tissues. Thus, MKWV and KURV may maintain tick-mammalian life cycles in Hokkaido, suggesting their potential as causative agents of tick-borne diseases in mammals.

Habitat properties are key drivers of Borrelia burgdorferi (s.l.) prevalence in Ixodes ricinus populations of deciduous forest fragments

BACKGROUND

The tick Ixodes ricinus has considerable impact on the health of humans and other terrestrial animals because it transmits several tick-borne pathogens (TBPs) such as B. burgdorferi (sensu lato), which causes Lyme borreliosis (LB). Small forest patches of agricultural landscapes provide many ecosystem services and also the disservice of LB risk. Biotic interactions and environmental filtering shape tick host communities distinctively between specific regions of Europe, which makes evaluating the dilution effect hypothesis and its influence across various scales challenging. Latitude, macroclimate, landscape and habitat properties drive both hosts and ticks and are comparable metrics across Europe. Therefore, we instead assess these environmental drivers as indicators and determine their respective roles for the prevalence of B. burgdorferi in I. ricinus.

METHODS

We sampled I. ricinus and measured environmental properties of macroclimate, landscape and habitat quality of forest patches in agricultural landscapes along a European macroclimatic gradient. We used linear mixed models to determine significant drivers and their relative importance for nymphal and adult B. burgdorferi prevalence. We suggest a new prevalence index, which is pool-size independent.

RESULTS

During summer months, our prevalence index varied between 0 and 0.4 per forest patch, indicating a low to moderate disservice. Habitat properties exerted a fourfold larger influence on B. burgdorferi prevalence than macroclimate and landscape properties combined. Increasingly available ecotone habitat of focal forest patches diluted and edge density at landscape scale amplified B. burgdorferi prevalence. Indicators of habitat attractiveness for tick hosts (food resources and shelter) were the most important predictors within habitat patches. More diverse and abundant macro- and microhabitat had a diluting effect, as it presumably diversifies the niches for tick-hosts and decreases the probability of contact between ticks and their hosts and hence the transmission likelihood.

CONCLUSIONS

Diluting effects of more diverse habitat patches would pose another reason to maintain or restore high biodiversity in forest patches of rural landscapes. We suggest classifying habitat patches by their regulating services as dilution and amplification habitat, which predominantly either decrease or increase B. burgdorferi prevalence at local and landscape scale and hence LB risk. Particular emphasis on promoting LB-diluting properties should be put on the management of those habitats that are frequently used by humans. In the light of these findings, climate change may be of little concern for LB risk at local scales, but this should be evaluated further.