Assortative pairing in Ixodes ricinus (Acari: Ixodidae), the European vector of Lyme borreliosis

Aseptic Meningitis Linked to Borrelia afzelii Seroconversion in Northeastern Greece: An Emerging Infectious Disease Contested in the Region

Borreliosis (Lyme disease) is a zoonosis, mediated to humans and small mammals through specific vectors (ticks), with increasing global incidence. It is associated with a variety of clinical manifestations and can, if not promptly recognized and left untreated, lead to significant disability. In Europe, the main Borrelia species causing disease in humans are Borrelia burgdorferi s.s., Borrelia afzelii, Borrelia garinii, and Borrelia spielmanii. The Ixodes ricinus tick is their principal vector. Although Lyme disease is considered endemic in the Balkan region and Turkey, and all three main Lyme pathogens have been detected in ticks collected in these countries, autochthonous Lyme disease remains controversial in Greece. We report a case of aseptic meningitis associated with antibody seroconversion against Borrelia afzelii in a young female patient from the prefecture of Thasos without any relevant travel history. The patient presented with fever and severe headache, and the cerebrospinal fluid examination showed lymphocytic pleocytosis. Serum analysis was positive for specific IgG antibodies against Borrelia afzelii. In the absence of typical erythema migrans, serological evidence of infection is required for diagnosis. Although atypical in terms of clinical presentation, the seasonality and geographical location of potential disease transmission in the reported patient should raise awareness among clinicians for a still controversial and potentially underreported emerging infectious disease in Greece.

Male mating preference in an ixodid tick

Background

Mate choice is a fundamental element of sexual selection and has the potential to shape the evolution of traits. Mate choice based on body size has been shown to be a common trait in several arthropod species. In hard ticks, a taxon of medical and veterinary importance, engorgement weight is positively correlated with reproductive output but it is unknown whether adult males show mate choice. Here, we experimentally investigated whether males (i) use chemical cues to choose their mating partner, (ii) consistently choose for the same female individual and (iii) prefer females with highest weight after feeding.

Methods

We used two experimental setups which allowed chemical communication between ticks: (i) a horizontal tube preventing physical contact with the female and (ii) an arena where tactile cues were allowed. In total, we tested 62 different triads in 124 tests (66 tests in the horizontal tube and 58 in the arena) composed of one male that could choose between two engorged females. Specifically, we tested 42 triads in the tube and 46 in the arena; 24 triads were repeatedly tested in the tube while 38 triads were tested in both setups.

Results

We found no preference for individual or heavier females in either setup. However, in the horizontal tube setup, males significantly preferred females that were not visited by them in the previous test.

Conclusions

Our results suggest a lack of male mate choice despite heavier females having higher fecundity. However, future studies should take into account that males may recognize the potential mating partners they previously met.

Graphical Abstract

Comparative population genetic structure of two ixodid tick species (Acari:Ixodidae) (Ixodes ovatus and Haemaphysalis flava) in Niigata prefecture, Japan

Ixodid ticks (Acari:Ixodidae) are essential vectors of tick-borne diseases in Japan. In this study, we characterized the population genetic structure and inferred genetic divergence in two widespread and abundant ixodid species, Ixodes ovatus and Haemaphysalis flava. Our hypothesis was that genetic divergence would be high in I. ovatus because of the low mobility of their small rodent hosts of immature I. ovatus would limit their gene flow compared to more mobile avian hosts of immature H. flava. We collected 320 adult I. ovatus from 29 locations and 223 adult H. flava from 17 locations across Niigata Prefecture, Japan, and investigated their genetic structure using DNA sequences from fragments of two mitochondrial gene regions, cox1 and the 16S rRNA gene. For I. ovatus, pairwise F_ST and analysis of molecular variance (AMOVA) analyses of cox1 and 16S sequences indicated significant genetic variation among populations, whereas both markers showed non-significant genetic variation among locations for H. flava. A cox1 gene tree and haplotype network revealed three genetic groups of I. ovatus. One of these groups consisted of haplotypes distributed at lower altitudes (251-471 m.a.s.l.). The cox1 sequences of I. ovatus from Japan clustered separately from I. ovatus sequences reported from China, suggesting the potential for cryptic species in Japan. Our results support our hypothesis and suggest that the host preference of ticks at the immature stage may influence the genetic structure of the ticks. This information may be important for understanding the tick-host interactions in the field to better understand the tick-borne disease transmission and in designing an effective tick control program.

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.

Genetic diversity of Ixodes ricinus (Ixodida: Ixodidae) ticks in sympatric and allopatric zones in Baltic countries

Ixodes ricinus (Linnaeus 1758) and Ixodes persulcatus (Schulze 1930) ticks are involved in the transmission of a wide variety of pathogens with considerable impact on human and animal health. The co-distribution zone of these two tick species is situated in the Baltic countries, which provides a special setting for the population studies. In the present study, genetic variability of I. ricinus ticks collected in allopatric and sympatric locations in the Baltic countries has been investigated using a sequence analysis of the mitochondrial DNA control region, 16S rRNA and cytb genes. There were 32 haplotypes (Hd: 0.8551) and 27 haplotypes (Hd:0.8213) of control region sequences from ticks in allopatric and sympatric zones detected, respectively. Out of 47 16S rRNA gene haplotypes, 32 haplotypes (Hd: 0.7213) were found in the allopatric zone and 27 (Hd:0.9572) in the sympatric zone. The Cytb gene was very conserved and monomorphic in ticks from the allopatric zone, whereas three unique haplotypes were observed in the sympatric zone. The higher number of unique haplotypes of the control region was detected in the allopatric zone. Median joining network and Fst analysis did not reveal a clear separation between ticks from the two zones.

Host, habitat and climate preferences of Ixodes angustus (Acari: Ixodidae) and infection with Borrelia burgdorferi and Anaplasma phagocytophilum in California, USA

The Holarctic tick Ixodes angustus is a competent vector for Borrelia burgdorferi, the etiologic agent of Lyme disease, and possibly Anaplasma phagocytophilum, the etiologic agent of granulocytic anaplasmosis, as well. From 2005 to 2013, we collected host-feeding I. angustus individuals from live-trapped small mammals and by flagging vegetation from 12 study sites in northern and central California, and tested for B. burgdorferi sensu lato, A. phagocytophilum, and Rickettsia spp. DNA by real-time PCR. Among 261 I. angustus collected (259 from hosts and two by flagging), the most common hosts were tree squirrels (20 % of ticks) and chipmunks (37 %). The PCR-prevalence for A. phagocytophilum and B. burgdorferi in ticks was 2 % and zero, respectively. The minimum infection prevalence on pooled DNA samples was 10 % for Rickettsia spp. DNA sequencing of the ompA gene identified this rickettsia as Candidatus Rickettsia angustus, a putative endosymbiont. A zero-inflated negative binomial mixed effects model was used to evaluate geographical and climatological predictors of I. angustus burden. When host species within study site and season within year were included in the model as nested random effects, all significant variables revealed that I. angustus burden increased as temperature decreased. Together with published data, these findings suggest that I. angustus is a host generalist, has a broad geographic distribution, is more abundant in areas with lower temperature within it's range, and is rarely infected with the pathogens A. phagocytophilum and B. burgdorferi.

Lyme disease/borreliosis as a systemic disease

Lyme disease/borreliosis (LD) is a well-known arthropod-transmitted entity in the northern hemisphere. The incidence of LD is reportedly rising throughout the world, although better diagnostic facilities may be contributory. The disease distribution is expanding in Europe, with its presence being now documented at higher altitudes and latitudes. Borrelia burgdorferi sensu lato is the most important genospecies leading to LD, although newer ones continue to be discovered. The variations in clinical spectrum with genospecies involved are an interesting feature. The alteration in gene expression while the organism cycles between two very different hosts is intriguing and has been described. The disease presents in three stages-namely, the early localized, early disseminated, and late stage. Erythema chronicum migrans is the pathognomic early lesion, and its diagnosis is purely clinical; however, laboratory diagnosis is essential for later manifestations. Two-tier serologic testing using an enzyme-linked immunosorbent assay (ELISA) as the first tier and immunoglobulin M (IgM) and IgG immunoblot as the second, if ELISA is positive or equivocal, is the mainstay of diagnosis. Doxycycline is the cornerstone of treatment, whereas parenteral therapy, mainly with ceftriaxone, is indicated in a few specific scenarios.

A seventeen-year epidemiological surveillance study of Borrelia burgdorferi infections in two provinces of northern Spain

This paper reports a 17-year seroepidemiological surveillance study of Borrelia burgdorferi infection, performed with the aim of improving our knowledge of the epidemiology of this pathogen. Serum samples (1,179) from patients (623, stratified with respect to age, sex, season, area of residence and occupation) bitten by ticks in two regions of northern Spain were IFA-tested for B. burgdorferi antibodies. Positive results were confirmed by western blotting. Antibodies specific for B. burgdorferi were found in 13.3% of the patients; 7.8% were IgM positive, 9.6% were IgG positive, and 4.33% were both IgM and IgG positive. Five species of ticks were identified in the seropositive patients: Dermacentor marginatus (41.17% of such patients) Dermacentor reticulatus (11.76%), Rhiphicephalus sanguineus (17.64%), Rhiphicephalus turanicus (5.88%) and Ixodes ricinus (23.52%). B. burgdorferi DNA was sought by PCR in ticks when available. One tick, a D. reticulatus male, was found carrying the pathogen. The seroprevalence found was similar to the previously demonstrated in similar studies in Spain and other European countries.

Development of genomic resources for the tick Ixodes ricinus: isolation and characterization of single nucleotide polymorphisms

Assessing the genetic variability of the tick Ixodes ricinus-an important vector of pathogens in Europe-is an essential step for setting up antitick control methods. Here, we report the first identification of a set of SNPs isolated from the genome of I. ricinus, by applying a reduction in genomic complexity, pyrosequencing and new bioinformatics tools. Almost 1.4 million of reads (average length: 528 nt) were generated with a full Roche 454 GS FLX run on two reduced representation libraries of I. ricinus. A newly developed bioinformatics tool (DiscoSnp), which isolates SNPs without requiring any reference genome, was used to obtain 321 088 putative SNPs. Stringent selection criteria were applied in a bioinformatics pipeline to select 1768 SNPs for the development of specific primers. Among 384 randomly SNPs tested by Fluidigm genotyping technology on 464 individuals ticks, 368 SNPs loci (96%) exhibited the presence of the two expected alleles. Hardy-Weinberg equilibrium tests conducted on six natural populations of ticks have shown that from 26 to 46 of the 384 loci exhibited significant heterozygote deficiency.

Isolation and characterization of ten polymorphic microsatellite loci in Ixodes arboricola, and cross-amplification in three other Ixodes species

We characterized ten polymorphic microsatellite loci from the tree-hole tick, Ixodes arboricola. Loci were screened in 11-18 individuals from three Belgian populations and five to ten alleles were found at each locus. Seven loci did not show deviations from Hardy-Weinberg equilibrium conditions and there were no indications for null alleles at these loci. The three other loci showed significant heterozygote deficiencies in at least one population, and a high potential for the occurrence of null alleles. We observed no effect of potential host DNA on the scoring of the microsatellites. Cross-amplification of the microsatellites was tested in eight specimens of three congeneric species: I. ricinus, I. hexagonus and I. frontalis. Depending on the species, six or seven of the loci were amplified in ≥ 4 of the 8 specimens and were polymorphic in each of these species (except for Ixaf 11 in I. frontalis and I. ricinus). These loci thus provide a tool for population genetic analysis of I. arboricola. The suitability of these markers needs to be further investigated in its congeners.

Host specialization in ticks and transmission of tick-borne diseases: a review

Determining patterns of host use, and the frequency at which these patterns change, are of key importance if we are to understand tick population dynamics, the evolution of tick biodiversity, and the circulation and evolution of associated pathogens. The question of whether ticks are typically host specialists or host generalists has been subject to much debate over the last half-century. Indeed, early research proposed that morphological diversity in ticks was linked to host specific adaptations and that most ticks were specialists. Later work disputed this idea and suggested that ticks are largely limited by biogeographic conditions and tend to use all locally available host species. The work presented in this review suggests that the actual answer likely lies somewhere between these two extremes. Although recent observational studies support the view that phylogenetically diverse host species share ticks when found on similar ecological ranges, theory on host range evolution predicts that host specialization should evolve in ticks given their life history characteristics. Contemporary work employing population genetic tools to examine host-associated population structure in several tick systems support this prediction and show that simple species records are not enough to determine whether a parasite is a true host generalist; host specialization does evolve in ticks at local scales, but may not always lead to speciation. Ticks therefore seem to follow a pattern of being global generalists, local specialists. Given this, the notion of host range needs to be modified from an evolutionary perspective, where one simply counts the number of hosts used across the geographic distribution, to a more ecological view, where one considers host use at a local scale, if we are to better understand the circulation of tick-borne pathogens and exposure risks for humans and livestock.

Changing distributions of ticks: causes and consequences

Today, we are witnessing changes in the spatial distribution and abundance of many species, including ticks and their associated pathogens. Evidence that these changes are primarily due to climate change, habitat modifications, and the globalisation of human activities are accumulating. Changes in the distribution of ticks and their invasion into new regions can have numerous consequences including modifications in their ecological characteristics and those of endemic species, impacts on the dynamics of local host populations and the emergence of human and livestock disease. Here, we review the principal causes for distributional shifts in tick populations and their consequences in terms of the ecological attributes of the species in question (i.e. phenotypic and genetic responses), pathogen transmission and disease epidemiology. We also describe different methodological approaches currently used to assess and predict such changes and their consequences. We finish with a discussion of new research avenues to develop in order to improve our understanding of these host-vector-pathogen interactions in the context of a changing world.

Testing local-scale panmixia provides insights into the cryptic ecology, evolution, and epidemiology of metazoan animal parasites

When every individual has an equal chance of mating with other individuals, the population is classified as panmictic. Amongst metazoan parasites of animals, local-scale panmixia can be disrupted due to not only non-random mating, but also non-random transmission among individual hosts of a single host population or non-random transmission among sympatric host species. Population genetics theory and analyses can be used to test the null hypothesis of panmixia and thus, allow one to draw inferences about parasite population dynamics that are difficult to observe directly. We provide an outline that addresses 3 tiered questions when testing parasite panmixia on local scales: is there greater than 1 parasite population/species, is there genetic subdivision amongst infrapopulations within a host population, and is there asexual reproduction or a non-random mating system? In this review, we highlight the evolutionary significance of non-panmixia on local scales and the genetic patterns that have been used to identify the different factors that may cause or explain deviations from panmixia on a local scale. We also discuss how tests of local-scale panmixia can provide a means to infer parasite population dynamics and epidemiology of medically relevant parasites.

Cryptic vector divergence masks vector-specific patterns of infection: an example from the marine cycle of Lyme borreliosis

Vector organisms are implicated in the transmission of close to a third of all infectious diseases. In many cases, multiple vectors (species or populations) can participate in transmission but may contribute differently to disease ecology and evolution. The presence of cryptic vector populations can be particularly problematic as differences in infection can be difficult to evaluate and may lead to erroneous evolutionary and epidemiological inferences. Here, we combine site-occupancy modeling and molecular assays to evaluate patterns of infection in the marine cycle of Lyme borreliosis, involving colonial seabirds, the tick Ixodes uriae, and bacteria of the Borrelia burgdorferi s.l. complex. In this cycle, the tick vector consists of multiple, cryptic (phenotypically undistinguishable but genetically distinct) host races that are frequently found in sympatry. Our results show that bacterial detection varies strongly among tick races leading to vector-specific biases if raw counts are used to calculate Borrelia prevalence. These differences are largely explained by differences in infection intensity among tick races. After accounting for detection probabilities, we found that overall prevalence in this system is higher than previously suspected and that certain vector-host combinations likely contribute more than others to the local dynamics and large-scale dispersal of Borrelia spirochetes. These results highlight the importance of evaluating vector population structure and accounting for detection probability when trying to understand the evolutionary ecology of vector-borne diseases.

Wahlund effects and sex-biased dispersal in Ixodes ricinus, the European vector of Lyme borreliosis: new tools for old data

Population genetics can help us better understand species microevolution and population biology, but inferences made from the genetic polymorphisms of field-collected organisms critically rely on sampling design. The population structure of the tick Ixodes ricinus L. (Acari, Ixodidae), a commonly encountered ectoparasite and the principal vector of human Lyme disease in Europe, has been the focus of some study, but many ecological aspects of this species remain poorly understood. Here, we apply a Bayesian clustering approach to observed and simulated data to examine within-population structure in I. ricinus, and to re-analyse patterns of sex-biased dispersal based on this substructure. We found between 18 and 27 distinct clusters within each of the 12 subsamples examined with a significant drop of heterozygote deficits. Parallel analyses on a comparable species, the seabird tick Ixodes uriae, indicated that these clusters can reflect important ecological features of the species (i.e., local host-associations). Analyses that considered the within-population clustering pattern of I. ricinus showed reversed patterns of sex-biased dispersal as compared to raw data (i.e., female-biased instead of male-biased dispersal). Simulated data supported the hypothesis that these scale-dependant patterns could be due to a combination of sex-specific dispersal and mortality. These different results raise new questions on the dispersal and host use strategies of I. ricinus and the potential importance of these ecological features for disease transmission. Furthermore, this work underlines the importance of taking into account patterns of genetic substructure when investigating sex-biased dispersal in natural populations.