Phylogenetic Lineages and Postglacial Dispersal Dynamics Characterize the Genetic Structure of the Tick, Ixodes ricinus, in Northwest Europe

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.

A crucial nexus: Phylogenetic versus ecological support of the life-cycle of Ixodes ricinus (Ixodoidea: Ixodidae) and Borrelia spp. amplification

The tick Ixodes ricinus parasitizes a wide range of vertebrates. These hosts vary in the relative contribution to the feeding of the different tick life stages, and their interplay is pivotal in the transmission dynamics of tick-borne pathogens. We aimed to know if there is a phylogenetic signal in the feeding and propagation hosts of I. ricinus, independently of other traits, as well as in the amplification of Borrelia burgdorferi (s.l.) in feeding larvae. We used a compilation of 1127 published field surveys in Europe, providing data for 96,586 hosts, resulting in 265,124 larvae, 72,080 nymphs and 37,726 adults. The load of immature ticks on hosts showed a significant phylogenetic signal towards the genera Psammodromus, Podarcis, and Lacerta (nymphs only). We hypothesize that such signal is the background hallmark of the primitive hosts associations of I. ricinus, probably in the glaciation refugia. A secondary phylogenetic signal for tick immatures appeared for some genera of Rodentia and Eulipotyphla. Results suggest the notion that the tick gained these hosts after spread from glaciation refugia. Analyses support a phylogenetic signal in the tick adults, firmly linked to Cetartiodactyla, but not to Carnivora or Aves. This study provides the first demonstration of host preferences in the generalist tick I. ricinus. We further demonstrate that combinations of vertebrates contribute in different proportions supporting the tick life-cycle in biogeographical regions of the Western Palaearctic as each region has unique combinations of dominant hosts. Analysis of the amplification of B. burgdorferi (s.l.) demonstrated that each genospecies is better amplified by competent reservoirs with which a strong phylogenetic signal exists. These vertebrates are the same along the spatial range: environmental traits do not change the reservoirs along the large territory studied. The transmission of B. burgdorferi (s.l.) is amplified by a few species of vertebrates, that share biogeographical regions with the tick vector in variable proportions.

Vector biology of the cat flea Ctenocephalides felis

Ctenocephalides felis, the cat flea, is among the most prevalent and widely dispersed vectors worldwide. Unfortunately, research on C. felis and associated pathogens (Bartonella and Rickettsia spp.) lags behind that of other vectors and vector-borne pathogens. Therefore, we aimed to review fundamental aspects of C. felis as a vector (behavior, epidemiology, phylogenetics, immunology, and microbiome composition) with an emphasis on key techniques and research avenues employed in other vector species. Future laboratory C. felis experimental infections with Bartonella, Rickettsia, and Wolbachia species/strains should examine the vector-pathogen interface utilizing contemporary visualization, transcriptomic, and gene-editing techniques. Further environmental sampling will inform the range and prevalence of C. felis and associated pathogens, improving the accuracy of vector and pathogen modeling to improve infection/infestation risk assessment and diagnostic recommendations.

Genetic analysis challenges the presence of Ixodes inopinatus in Central Europe: development of a multiplex PCR to distinguish I. inopinatus from I. ricinus

BACKGROUND

Ixodes ricinus is an important vector of several pathogens, primarily in Europe. Recently, Ixodes inopinatus was described from Spain, Portugal, and North Africa and then reported from several European countries. In this study, a multiplex polymerase chain reaction (PCR) was developed to distinguish I. ricinus from I. inopinatus and used in the surveillance of I. inopinatus in Algeria (ALG) and three regions in the Czech Republic (CZ).

METHODS

A multiplex PCR on TROSPA and sequencing of several mitochondrial (16S rDNA, COI) and nuclear markers (TROSPA, ITS2, calreticulin) were used to differentiate these two species and for a subsequent phylogenetic analysis.

RESULTS

Sequencing of TROSPA, COI, and ITS2 separated these two species into two subclades, while 16S rDNA and calreticulin could not distinguish I. ricinus from I. inopinatus. Interestingly, 23 nucleotide positions in the TROSPA gene had consistently double peaks in a subset of ticks from CZ. Cloning of these PCR products led to a clear separation of I. ricinus and I. inopinatus indicating hybridization and introgression between these two tick taxa. Based on a multiplex PCR of TROSPA and analysis of sequences of TROSPA, COI, and ITS2, the majority of ticks in CZ were I. ricinus, no I. inopinatus ticks were found, and 10 specimens showed signs of hybridization. In contrast, most ticks in ALG were I. inopinatus, four ticks were I. ricinus, and no signs of hybridization and introgression were detected.

CONCLUSIONS

We developed a multiplex PCR method based on the TROSPA gene to differentiate I. ricinus and I. inopinatus. We demonstrate the lack of evidence for the presence of I. inopinatus in Central Europe and propose that previous studies be re-examined. Mitochondrial markers are not suitable for distinguishing I. inopinatus from I. ricinus. Furthermore, our data indicate that I. inopinatus and I. ricinus can hybridize, and the hybrids can survive in Europe.

Phylogenetic placement of Turkish populations of Ixodes ricinus and Ixodes inopinatus

Studies on phylogeography and population structure of Ixodes ricinus have been carried out in Europe for decades, but the number of specimens from the Middle East included in these analyses is relatively small, despite the wide distribution of the species in this area. This study aimed to clarify the phylogenetic positions of I. ricinus from Turkey as well as to investigate the presence of Ixodes inopinatus in Anatolia. For this purpose, one mitochondrial (mt 16S rDNA) and one nuclear gene (defensin) were used to generate molecular data from I. ricinus samples, which were collected from 17 locations across the species' distributional range in Turkey. Bayesian inference was used to investigate phylogenetic relationships. Globally, the mt 16S rDNA lineages correspond to the lineages revealed by defensin; I. ricinus and I. inopinatus sequences clustered separately. However, a discordant genetic pattern was observed between the phylogenetic position of turkish I. ricinus revealed by nuclear versus mitochondrial genes. All Turkish haplotypes of mt 16SrDNA clustered with I. ricinus samples from Europe, which might be the result of extensive gene flow between populations of Europe and the Middle East. On the other hand, a sample from Thrace Region grouped within I. inopinatus clade. Thus, the occurrence of I. inopinatus in Turkey was demonstrated for the first time using molecular data. Moreover, four individuals were found to be heterozygous for the defensin. The potential evolutionary processes that underlie this observed discrepancy between the phylogenetic trees of two genes have been discussed.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Assessing the importance of Molecular and Genetic perspectives in Prophesying the KFD transmission risk provinces in the Western Ghats, Kerala, INDIA in context with spatial distribution, Extensive genetic Diversity, and phylogeography

The deadly effects of KFD have been pointed in southern India; however, the infecting regions have been getting larger in recent epochs. People who live or work in regions where KFDV infected tick vectors are present are severely prone to procuring the infection. Being aware of tick vectors and infectious agents' geospatial location is vital to direct sustenance approaches to prevent and manage such infectious diseases as KFD. The present investigation has focussed on the spatial distribution, Extensive genetic Diversity, and phylogeography to forecast the probable KFD disease risk provinces in the Western Ghats. The statistical analysis for diversity indices and community comparison has been performed by using SPSS version 24.0.0 and R software version 3.4.2. The nucleotide sequences of the respective ticks and KFDV were retrieved from NCBI. The first strand of this investigation revealed that, around the world, the Indian province was found to exhibit a maximum range of diversity for tick vectors. The next strands prophesied the KFD transmission risk areas in the Western Ghats region, India, with computational spatial analysis and phylogeography. The final strand exposed the genetic diversity of the KFD virus and the tick vectors in terms of their spatial distribution worldwide.

What do we still need to know about Ixodes ricinus?

In spite of many decades of intensive research on Ixodes ricinus, the castor bean tick of Europe, several important aspects of its basic biology remain elusive, such as the factors determining seasonal development, tick abundance and host specificity, and the importance of water management. Additionally, there are more recent questions about the geographical diversity of tick genotypes and phenotypes, the role of migratory birds in the ecoepidemiology of I. ricinus, the importance of protective immune responses against I. ricinus, particularly in the context of vaccination, and the role of the microbiome in pathogen transmission. Without more detailed knowledge of these issues, it is difficult to assess the likely effects of changes in climate and biodiversity on tick distribution and activity, to predict potential risks arising from new and established tick populations and I. ricinus-borne pathogens, and to improve prevention and control measures. This review aims to discuss the most important outstanding questions against the backdrop of the current state of knowledge of this important tick species.

Perpetuation of Borreliae

With one exception (epidemic relapsing fever), borreliae are obligately maintained in nature by ticks. Although some Borrelia spp. may be vertically transmitted to subsequent generations of ticks, most require amplification by a vertebrate host because inheritance is not stable. Enzootic cycles of borreliae have been found globally; those receiving the most attention from researchers are those whose vectors have some degree of anthropophily and, thus, cause zoonoses such as Lyme disease or relapsing fever. To some extent, our views on the synecology of the borreliae has been dominated by an applied focus, viz., analyses that seek to understand the elements of human risk for borreliosis. But, the elements of borrelial perpetuation do not necessarily bear upon risk, nor do our concepts of risk provide the best structure for analyzing perpetuation. We identify the major global themes for the perpetuation of borreliae, and summarize local variations on those themes, focusing on key literature to outline the factors that serve as the basis for the distribution and abundance of borreliae.

Mitochondrial DNA in Ixodus ricinus (Acari: Ixodidae) on birds reflects ticks' transportation routes to Lista, Norway

Ticks are important pathogen vectors, and large mammals and birds have the greatest potential for dispersing them. To study tick dispersal by migrating birds, we have analysed genetic variations in mitochondrial DNA control region from Ixodes ricinus from northward migrating blackbird, Turdus merula, and (European) robin, Erithacus rubecula, at the Lista Bird Observatory in southwestern Norway. We compared their genetic structure with that of resident tick populations from areas covering their expected last stop (i.e. Great Britain and Jutland, Denmark) before taking off for southern Norway, and the resident tick population at Lista. The statistical analysis showed that the I. ricinus found on blackbirds differed significantly from those found on robins, which is consistent with the birds' differential migration routes. I. ricinus from robins did not differ genetically from those flagged at Jutland, suggesting that the former mainly originate in continental Europe. Bayesian analysis indicated that most of the blackbirds caught early in the spring (i.e. before or on the 1st of April) carried ticks of a mixed origin from both Great Britain and continental Europe, while blackbirds caught later in the season carried an increasing amount of ticks acquired locally.

Strong genetic structure among populations of the tick Ixodes ricinus across its range

Ixodes ricinus is the most common and widely distributed tick species in Europe, responsible for several zoonotic diseases, including Lyme borreliosis. Population genetics of disease vectors is a useful tool for understanding the spread of pathogens and infection risks. Despite the threat to the public health due to the climate-driven distribution changes of I. ricinus, the genetic structure of tick populations, though essential for understanding epidemiology, remains unclear. Previous studies have demonstrated weak to no apparent spatial pattern of genetic differentiation between European populations. Here, we analysed the population genetic structure of 497 individuals from 28 tick populations sampled from 20 countries across Europe, the Middle-East, and northern Africa. We analysed 125 SNPs loci after quality control. We ran Bayesian and multivariate hierarchical clustering analyses to identify and describe clusters of genetically related individuals. Both clustering methods support the identification of three spatially-structured clusters. Individuals from the south and north-western parts of Eurasia form a separated cluster from northern European populations, while central European populations are a mix between the two groups. Our findings have important implications for understanding the dispersal processes that shape the spread of zoonotic diseases under anthropogenic global changes.

Geographical distribution and prevalence of tick-borne encephalitis virus in questing Ixodes ricinus ticks and phylogeographic structure of the Ixodes ricinus vector in Norway

The tick-borne encephalitis virus (TBEV), a zoonotic flaviviral infection, is endemic in large parts of Norway and Eurasia. Humans are mainly infected with TBEV via bites from infected ticks. In Norway, the main geographical distribution of ticks is along the Norwegian coastline from southeast (~59°N) and up to the southern parts of Nordland County (~65°N). In this study, we collected ticks by flagging along the coast from Østfold County to Nordland County. By whole-genome sequencing of the mitochondrial genome of Ixodes ricinus, the phylogenetic tree suggests that there is limited phylogeographic structure both in Norway and in Europe. The overall TBEV prevalence is 0.3% for nymphs and 4.3% for adults. The highest estimated TBEV prevalence in adult ticks was detected in Rogaland and Vestfold County, while for nymphs it is highest in Vestfold, Vest-Agder and Rogaland. The present work is one of the largest studies on distribution and prevalence of TBEV in ticks in Scandinavia, showing that the virus is wider distributed in Norway than previously anticipated.

How general are generalist parasites? The small mammal part of the Lyme disease transmission cycle in two ecosystems in northern Europe

The pathogens causing Lyme disease are all vectored by generalist tick species found on a wide range of vertebrates, but spatial and annual variation in host use has rarely been quantified. We here compare the load of Ixodes ricinus (the vector) on small mammals and investigate the infection prevalence of Borrelia burgdorferi s.l. (the pathogen) involved in the enzootic transmission cycle of Lyme disease in two contrasting ecosystems in Norway from 2014 to 2016. The most common larval tick host in the eastern region was the bank vole, while the common shrew dominated in the western region of Norway. However, the wood mouse and the bank vole had consistently higher larval tick loads than the common shrew in both ecosystems. Hence, the evidence indicated that species are differently suitable as hosts, regardless of their abundances. The pathogen infection prevalence was similar among small mammal species, but markedly higher in the region with larger small mammal populations and higher tick loads, while the seasonal and annual variation was less marked. Our study indicated that the generalist I. ricinus shows consistent patterns of load on species of small vertebrate hosts, while B. burgdorferi s.l. (B. afzelii) was a true generalist. The similar roles of host species across regions suggest that disease dynamics can be predicted from host community composition, but predicting the role of host community composition for disease dynamics requires a detailed understanding of the different species population limitations under global change.

Multi-locus sequence typing of Ixodes ricinus and its symbiont Candidatus Midichloria mitochondrii across Europe reveals evidence of local co-cladogenesis in Scotland

Ticks have relatively complex microbiomes, but only a small proportion of the bacterial symbionts recorded from ticks are vertically transmitted. Moreover, co-cladogenesis between ticks and their symbionts, indicating an intimate relationship over evolutionary history driven by a mutualistic association, is the exception rather than the rule. One of the most widespread tick symbionts is Candidatus Midichloria, which has been detected in all of the major tick genera of medical and veterinary importance. In some species of Ixodes, such as the sheep tick Ixodes ricinus (infected with Candidatus Midichloria mitochondrii), the symbiont is fixed in wild adult female ticks, suggesting an obligate mutualism. However, almost no information is available on genetic variation in Candidatus M. mitochondrii or possible co-cladogenesis with its host across its geographic range. Here, we report the first survey of Candidatus M. mitochondrii in I. ricinus in Great Britain and a multi-locus sequence typing (MLST) analysis of tick and symbiont between British ticks and those collected in continental Europe. We show that while the prevalence of the symbiont in nymphs collected in England is similar to that reported from the continent, a higher prevalence in nymphs and adult males is apparent in Wales. In general, Candidatus M. mitochondrii exhibits very low levels of sequence diversity, although a consistent signal of host-symbiont coevolution was apparent in Scotland. Moreover, the tick MLST scheme revealed that Scottish specimens form a clade that is partially separated from other British ticks, with almost no contribution of continental sequence types in this north-westerly border of the tick's natural range. The low diversity of Candidatus M. mitochondrii, in contrast with previously reported high rates of polymorphism in I. ricinus mitogenomes, suggests that the symbiont may have swept across Europe recently via a horizontal, rather than vertical, transmission route.