Tick microbial associations at the crossroad of horizontal and vertical transmission pathways
Parasites & Vectors volume 15, Article number: 380 (2022)
Microbial communities can affect disease risk by interfering with the transmission or maintenance of pathogens in blood-feeding arthropods. Here, we investigated whether bacterial communities vary between Ixodes ricinus nymphs which were or were not infected with horizontally transmitted human pathogens.
Ticks from eight forest sites were tested for the presence of Borrelia burgdorferi sensu lato, Babesia spp., Anaplasma phagocytophilum, and Neoehrlichia mikurensis by quantitative polymerase chain reaction (qPCR), and their microbiomes were determined by 16S rRNA amplicon sequencing. Tick bacterial communities clustered poorly by pathogen infection status but better by geography. As a second approach, we analysed variation in tick microorganism community structure (in terms of species co-infection) across space using hierarchical modelling of species communities. For that, we analysed almost 14,000 nymphs, which were tested for the presence of horizontally transmitted pathogens B. burgdorferi s.l., A. phagocytophilum, and N. mikurensis, and the vertically transmitted tick symbionts Rickettsia helvetica, Rickettsiella spp., Spiroplasma ixodetis, and Candidatus Midichloria mitochondrii.
With the exception of Rickettsiella spp., all microorganisms had either significant negative (R. helvetica and A. phagocytophilum) or positive (S. ixodetis, N. mikurensis, and B. burgdorferi s.l.) associations with M. mitochondrii. Two tick symbionts, R. helvetica and S. ixodetis, were negatively associated with each other. As expected, both B. burgdorferi s.l. and N. mikurensis had a significant positive association with each other and a negative association with A. phagocytophilum. Although these few specific associations do not appear to have a large effect on the entire microbiome composition, they can still be relevant for tick-borne pathogen dynamics.
Based on our results, we propose that M. mitochondrii alters the propensity of ticks to acquire or maintain horizontally acquired pathogens. The underlying mechanisms for some of these remarkable interactions are discussed herein and merit further investigation.
There is increasing evidence that members of the arthropod microbiome can decrease vector-borne disease risk [1,2,3]. Microbes can modulate the vectorial competence of an arthropod by decreasing their susceptibility to pathogens and reducing pathogen transmission, both of which are necessary for pathogen maintenance in enzootic cycles. For instance, an Enterobacter bacterium isolated from wild Anopheles gambiae mosquitoes in Zambia confers resistance to infection with the malaria parasite in 99% of mosquitoes within the population by interfering with its development prior to the invasion of the midgut epithelium . A microsporidian symbiont in another member of the An. gambiae complex is thought to play a similar role .
Like mosquitoes, ticks can transmit a plethora of pathogens causing disease in humans and domestic animals [5, 6]. Thus, an improved understanding of interactions between tick microorganisms is necessary for the development of novel strategies for controlling tick-borne diseases. A necessary first step is to understand whether and how tick microbiomes can impede or facilitate the transmission of pathogens. The microbiome of ticks consists of viruses, bacteria, and protozoa. Their relationships with ticks are often context-dependent and range from obligate mutualistic to exclusively parasitic [7,8,9]. In principle, all microorganisms (including tick-borne pathogens) relying on ticks for their survival are called symbionts. Here, we refer to microorganisms that have been shown to cause disease in humans and animals as pathogens.
In ticks, bacteria are the most commonly studied organisms in the microbiome-disease risk context. To date, studies have shown the mutually exclusive occurrence of pathogenic and non-pathogenic Rickettsia species in Amblyomma maculatum and Dermacentor andersoni ticks [10,11,12], as well as a negative association between the occurrence of Rickettsia bellii and Anaplasma marginale in D. andersoni ticks . Regarding facilitative interactions, Candidatus Midichloria mitochondrii (hereafter M. mitochondrii) has been shown to be a successful colonizing partner of pathogenic Rickettsia parkeri in A. maculatum ticks .
In Ixodes ricinus, the best-studied symbionts include obligate intracellular bacteria belonging to the Rickettsia, Midichloria, Rickettsiella, and Spiroplasma genera, which are predominantly transmitted vertically from a female tick to her eggs . Nevertheless, little is known about interactions between the tick microbiome members and pathogens such as Borrelia, Anaplasma, and Neoehrlichia. Although Borrelia burgdorferi sensu lato, responsible for Lyme disease, has been occasionally detected in larvae [16, 17], the aforementioned pathogens are predominantly transmitted horizontally [18,19,20]. In other words, ticks are born without and only acquire them as larvae or nymphs while feeding on vertebrates. The exception is Borrelia miyamotoi, a relapsing fever spirochete, in which vertical transmission is more efficient than horizontal .
The colonization of ticks by pathogens may be inhibited by vertically transmitted symbionts, which are already present in larvae. Evidence from other arthropod–pathogen systems suggests the involvement of several potential mechanisms, including direct killing, competition, and enhancement of host immune responses [22, 23]. Before describing these mechanisms, a crucial first step is to identify associations between bacterial species. For this purpose, co-infection analyses are frequently applied. However, until now, studies investigating co-infections of vertically and horizontally transmitted symbionts have reported contradictory results. For example, in Ixodes ticks, observed associations between Rickettsia and Borrelia have ranged from negative to positive [24,25,26].
At the same time, horizontally transmitted pathogens may engage in strategies to promote their replication in the presence of other microbes. Some studies have shown that Anaplasma and Borrelia may alter the tick gut microbiome to increase their colonization success by adjusting a tick antibacterial protein expression [27, 28]. However, existing evidence is limited to laboratory experiments, and it remains an outstanding question of how co-infection alters the host microbiome under natural conditions.
Here, we tested two hypotheses. The first hypothesis is that ticks infected with pathogenic bacteria such as B. burgdorferi s. l., possess a different bacterial microbiome than uninfected ticks or ticks infected with the pathogenic protozoan Babesia, bacterial B. miyamotoi, Anaplasma phagocytophilum, or N. mikurensis. To test this, we performed 16S rRNA amplicon sequencing on pools of I. ricinus nymphs infected with one of these pathogens and compared them with each other.
Importantly, I. ricinus is a three-host generalist tick, feeding once per life stage (larva, nymph, and adult) on various vertebrates, all of which may be amplification hosts of human pathogens . Because questing nymphs feed only once as larva, the nymphal microbiome predominantly consists of vertically transmitted symbionts and, potentially, of a horizontally acquired pathogen(s).
The second hypothesis is that vertically transmitted members of the tick microbiome either facilitate or impede the acquisition and maintenance of pathogens. To test this, we used data obtained in previous studies [30, 31], where approximately 14,000 individual questing nymphs from 19 locations scattered across the Netherlands were screened for the presence of the horizontally transmitted pathogens B. burgdorferi s.l., A. phagocytophilum, and N. mikurensis as well as the predominantly vertically transmitted symbionts M. mitochondrii, Rickettsia helvetica, Rickettsiella spp. and Spiroplasma ixodetis . Subsequently, we analysed microorganism co-infections within tick individuals using hierarchical modelling of species communities (HMSC), which tests whether microorganisms co-occur more (or less) often than by chance, given the local geographical prevalence.
We analysed questing I. ricinus from eight forest sites in the Netherlands (Fig. 1). The study sites were selected based on existing knowledge of Borrelia genospecies prevalence, the density of ticks, vegetation profile, and vertebrate composition obtained in a previous study . The full names of the sites and their coordinates are provided in Additional file 1: Table S1.
Sample collection, DNA extraction, and pathogen detection
For the microbiome analyses, we collected a total of 7874 questing I. ricinus nymphs in 2016 and 2017 (Additional file 1: Table S1). All ticks were washed three times in 70% ethanol, and the DNA of individual nymphs was extracted with ammonium hydroxide as described previously . The nymphs were analysed individually for the presence of tick-borne pathogens with multiplex real-time polymerase chain reaction (PCR), based on various target genes as described previously, including B. burgdorferi s.l. , B. miyamotoi , N. mikurensis , A. phagocytophilum , and Babesia spp.), which were designed to detect Babesia divergens, B. venatorum (formerly called EU1-3), B. capreoli, and B. odocoilei . A detailed description of the quantitative PCR (qPCR) protocol is provided in Additional file 1: Table S2. Samples positive for B. burgdorferi s.l. were subjected to conventional PCR followed by sequencing for genotype identification . After testing for tick-borne pathogens, ticks were pooled and re-extracted using the QIAGEN DNeasy Blood & Tissue Kit according to the manufacturer’s protocol (QIAGEN, Venlo, Netherlands). Pools consisted of DNA from five nymphs, which were positive for the same pathogen and negative for all the others. Sample metadata are provided in Additional file 1: Table S3, left table. Pools of nymphs negative for all pathogens were also included. For each location, a minimum of six and a maximum of 49 pools were processed (Additional file 1: Table S3, right table). Two negative controls, which underwent the same processes as tick samples including crushing, extraction, and amplification, and one company internal control were included. The DNA concentration in all samples was measured with the Qubit dsDNA HS Kit (Thermo Fisher Scientific, Ochten, Netherlands).
16S rRNA quantification and absolute bacterial density
In all samples, absolute bacterial density was quantified, and proportions were multiplied by a load to convert relative into absolute abundance. The quantification of total bacterial DNA load was determined by 16S rRNA qPCR as previously described [39,40,41]. The details on a positive control, primers, protocol, and in silico analysis can be found in Additional file 1: Table S4. It should be noted that the primers were not developed specifically for tick-associated microorganisms, and that in this study, the 16S rRNA qPCR was used in addition to other methods. The use of absolute bacterial density is a cost-effective and scalable solution for datasets of this size, since quantification methods through flow cytometry are not compatible with the sampling technique. Samples were normalized, and the normalized abundance values were scaled by the 16S rRNA qPCR loads (ng/µl) of each sample.
Microbial profiling and taxonomic clustering
Illumina MiSeq sequencing libraries targeting the V3-V4 region of the 16S rRNA amplicon were generated and sequenced by BaseClear (Leiden, Netherlands). In short, barcoded amplicons from the V3-V4 region of 16S rRNA genes were generated using a two-step PCR. Between 10 and 25 ng of genomic (g)DNA was used as a template for the first PCR with a total volume of 50 µl using the 341F (5′-CCTACGGGNGGCWGCAG-3′) and the 785R (5′-GACTACHVGGGTATCTAATCC-3′) primers appended with Illumina adaptor sequences. Control PCR reactions were performed alongside each separate amplification without the addition of a template. PCR products were purified, and the size of the PCR products was checked on a fragment analyser (Advanced Analytical) and quantified by fluorometric analysis. Purified PCR products were used for the second PCR in combination with sample-specific barcoded primers (Nextera XT index kit, Illumina). Subsequently, PCR products were purified, checked on the fragment analyser (Advanced Analytical) and quantified, followed by multiplexing, clustering, and sequencing on an Illumina MiSeq with the paired-end (2×) 300-base-pair (bp) protocol and indexing. The sequencing run was analysed with the Illumina CASAVA pipeline (v1.8.3) with demultiplexing based on sample-specific barcodes. The raw sequencing data produced were processed by removing the sequence reads of inadequate quality (only “passing filter” reads were selected) and discarding reads containing adaptor sequences or PhiX control using an in-house filtering protocol. A quality assessment on the remaining reads was performed using the FastQC version 0.10.0 quality control tool. Sequenced reads were imported to the QIAGEN CLC Genomics Workbench 10.0.1 supplemented with CLC Microbial Genomics Module 3.6.1 (www.clcbio.com). Overlapping pairs of raw reads were merged into single longer reads and trimmed with a quality score limit of 0.05 and two ambiguous nucleotides. At this stage, primer sequences were trimmed. Subsequently, reads were fixed-length trimmed (~ 400 bp). To identify operational taxonomic units (OTUs), reads were clustered using the SILVA 16S version 128 reference database with 97% identity as the clustering criterion. Chimeras were removed with a built-in tool in the CLC Genomics Workbench.
All analyses were carried out in R 3.6.3 . We used the R package vegan (version 2.5-6) for ordination . A principal coordinate analysis (PCoA) was carried out using Bray–Curtis dissimilarities. Results were visualized with ggplot2 (version 3.3.2; ).
Additionally, we carried out permutational multivariate analysis of variance (PERMANOVA)  with the adonis function from R package vegan to assess whether pathogen presence significantly affected community variation. We tested for multivariate dispersion through the betadisper function.
To investigate co-infections of tick microorganisms, we utilized the qPCR data on the infection prevalence of tick symbionts and pathogens generated in previous studies [30, 31, 46]. Briefly, questing nymphs of I. ricinus were collected from 19 forest sites in the Netherlands in 2013 and 2014 (Fig. 1). A total of 13,968 individual nymphs of I. ricinus were screened for S. ixodetis, R. helvetica, Rickettsiella spp., M. mitochondrii, A. phagocytophilum, N. mikurensis, and B. burgdorferi s.l. Borrelia miyamatoi and Babesia spp. data were not included in further investigations due to their low prevalence, which did not allow for performing statistical analyses with confidence.
In a recent study, we found geographical differences in tick microbial community structure . Therefore, to explore non-random residual associations (co-infections) of microorganisms within individual I. ricinus ticks while controlling for confounding factors associated with the between-individual and spatial variation, we used a recently developed statistical approach, HMSC . HMSC is a generalized framework based on a hierarchical joint species distribution approach (JSDM) . In the context of our study, JSDM considers information on the presence/absence of many species simultaneously, at the community level, and incorporates the effects of environmental factors and interspecific interactions on species incidence into a single model, which allows us to estimate species co-occurrence patterns (e.g. [49, 50]). An advantage of the HMSC methodology is that it allows the consideration of multiple hierarchical levels of the study design.
The R package Hmsc  was used to build the model. We included the tick’s identity as a sampling unit, which was nested as a random effect within the location where the tick was collected. We fitted a probit model to presence/absence data to ensure successful convergence. We sampled the posterior distribution with three Markov chain Monte Carlo (MCMC) chains, each with 1500 samples, thin 1000 and transient 750,000.
We then examined MCMC convergence by assessing the effective size of the posterior sample and Gelman and Rubin's  convergence diagnostic (potential scale reduction factor; PSRF) with the gelman.diag function. We performed variance portioning with the computeVariancePartitioning function to investigate which random effect was most responsible for the observed variation in the prevalence of individual microorganisms. Lastly, to compute the microorganism association matrices associated with each random level, we used the computeAssociations function. The corrplot function from the corrplot package was used to plot only those associations for which the posterior probability of being negative or positive was at least 0.95 .
Microbial profiling and taxonomic clustering
A total of 165 out of 166 processed samples successfully generated 6,590,988 raw reads on an Illumina MiSeq flow cell. One sample was excluded because it had an extremely low DNA load, which resulted in no amplification. A total of 4,454,814 reads were assigned taxonomy. A total of 184,682 unique reads were clustered into 8818 OTUs (Additional file 1: Table S5). However, 1966 OTUs, which accounted for 0.7% of all reads, could not be assigned to any known taxa and were discarded from further analyses. The top 10 most abundant genera included Borrelia, Midichloria, Neoehrlichia, Methylobacterium, Mycobacterium, Pseudomonas, Rickettsia, Rickettsiella, Sphingomonas, and Wolbachia (Additional file 2: Figure S1). Despite being the most abundant taxa in some samples, Wolbachia was not considered an I. ricinus tick symbiont, and we did not include it in the further qPCR analysis. Its origin in our tick samples was due to the presence of endoparasitoid Ixodiphagus hookeri eggs [54, 55].
Two types of negative controls were utilized: two extractions without template DNA and one only containing sequencing reagents. The PCoA showed that all three negative controls were highly similar, indicating that any errors introduced through reagent contamination were preserved (Additional file 2: Figure S2). The ordination also demonstrated that tick samples did not neatly cluster away from the negative controls. This lack of a distinction from the negative controls is likely attributable to low bacterial abundance in some samples and a high abundance of Pseudomonas, a known contaminant that was highly abundant in the negative controls.
The PCoA did not convincingly show that community composition was structured by pathogen-specific pools, or separated along the pools with or without pathogens (Fig. 2A). Instead, the PCoA suggested that samples were mostly separated by the geographical region, as indicated by longitude, of tick collection (Fig. 2B). Therefore, we tested with PERMANOVA whether communities were significantly different for ticks screened positively for different pathogens. Pathogen identity had a significant yet weak relationship with sample composition (R2 of 0.157, P = 0.001). However, significant overdispersion (P < 0.001) indicates that the P-values may be overinflated. Although the results suggest that microbiome composition was significantly different across samples with different pathogens, the R2 demonstrates that most variation could not be attributed to a pathogen status.
Co-infection of vertically and horizontally transmitted symbionts
Data on the prevalence of tick microorganisms in questing nymphs generated in previous studies were analysed for co-infections [30, 31, 46]. Briefly, 13,967 questing nymphs from 19 forest sites were tested with qPCR for the presence of S. ixodetis, R. helvetica, Rickettsiella spp., M. mitochondrii, A. phagocytophilum, N. mikurensis, and B. burgdorferi s.l. All microorganisms were present in all locations; however, infection rates varied strongly among the bacterial species and locations. The prevalence of the tick microorganisms per location and region is provided in Table S6, as well as in Krawczyk et al.  and Takumi et al. .
A total of 1732 (12.4%) tested questing nymphs were free of any of the seven studied microorganisms, and 12,225 (87.6%) nymphs were infected with one or more. The exact binomial test showed that infections with at least one, two, three, or four tick microorganisms were significantly less frequent than expected (P < 0.001 in all cases), and infections with at least five, six, or seven microorganisms occurred as often as expected (Additional file 1: Table S7). A single infection was detected in 4791 nymphs, a double infection in 4992 nymphs, a triple infection in 2177 nymphs, and a quadruple infection in 307 nymphs. A total of 27 nymphs were infected with five, one with six, and none with seven microorganisms, showing that co-infections are not randomly distributed in tick populations.
To investigate co-infections of tick microorganisms, we constructed an HMSC model in which we included the tick’s identity as a sampling unit, which was nested as a random effect within the location where the tick was collected. The MCMC convergence of the model was satisfactory: the potential scale reduction factors for the β-parameters were 1.056 on average (max = 1.211). Variance partitioning showed that the tick’s identity explained a substantial proportion of the observed variation in the prevalence of M. mitochondrii (80.0%), S. ixodetis (80.1%), N. mikurensis (70.6%), and B. burgdorferi s.l. (84.7%). Location explained variation in the prevalence of Rickettsiella spp. (99.6%), R. helvetica (84.3%), and A. phagocytophilum (61.2%).
Although all possible combinations of two microorganisms were found in nymphs, some combinations occurred significantly more or less often than by chance. The HMSC model showed that all microorganisms but Rickettsiella spp. had either significant negative (R. helvetica and A. phagocytophilum) or positive (S. ixodetis, N. mikurensis, and B. burgdorferi s.l.) associations with M. mitochondrii (Fig. 3). Borrelia burgdorferi s.l. and N. mikurensis had a significant positive association with each other and a negative association with A. phagocytophilum (Fig. 3). Neoehrlichia mikurensis also had a significant negative association with S. ixodetis and a positive association with R. helvetica (Fig. 3). Lastly, two tick symbionts, S. ixodetis and R. helvetica, had a significant negative association with each other (Fig. 3).
This study detected significant negative and positive associations between tick microorganisms known to utilize different transmission routes. Populations of ticks carrying M. mitochondrii were more likely to carry another tick symbiont, S. ixodetis, and the pathogens B. burgdorferi s.l. and N. mikurensis, than were ticks without M. mitochondrii (Fig. 3). Nevertheless, we did not observe differences in the bacterial microbiome of ticks infected with distinct horizontally acquired pathogens such as Borrelia afzelii, Borrelia garinii, Borrelia valaisiana, B. miyamotoi, Babesia, A. phagocytophilum, and N. mikurensis (Fig. 2).
Microbiome of I. ricinus
We did not detect differences between the bacterial microbiome of ticks with and without pathogens (Fig. 2). As the presence of a pathogen in nymphs indicates that it fed on a specific vertebrate host as a larva, our findings imply that the vertebrate host has little to no effect on the tick bacterial community. Our results are in line with previous findings by Hawlena et al. , but they contrast with a study by Swei and Kwan , who found differences in the microbiome of I. pacificus ticks feeding on lizards and mice. We are not aware of any of the nymphs studied here feeding on lizards, which are relatively uncommon tick hosts in our study areas . Nevertheless, we cannot rule out the potential influence of reptiles on the tick microbiome. Non-avian reptiles, as exothermic vertebrates, have less stringent requirements for regulating blood biochemical properties, and their blood may have vastly different osmotic pressure and pH compared with that of mammals and birds . While we could not investigate such effects in this study, this could lead to the absence of host effects compared with the work by Swei and Kwan .
A possible explanation for our finding could be that ticks are genetically equipped to resist opportunistic bacteria, such as host skin commensals. Hayes et al.  discovered that I. scapularis horizontally acquired an antimicrobial toxin gene from bacteria. The gene encodes lytic cell wall-degrading enzymes delivered to the host bite site via saliva, which is responsible for selectively killing the skin-associated bacteria of the vertebrate. Interestingly, the enzymes had no intrinsic lytic activity against B. burgdorferi . Therefore, pathogens may still affect the tick gut microbiome to facilitate their colonization, as previously observed in laboratory experiments [27, 28]. In future work, investigating individual organs deriving from field-collected ticks may extend our understanding of microbiome–pathogen interactions.
Beta diversity of the bacterial communities varied by geographical region, as indicated by longitude, which is in accord with previous studies [30, 60,61,62]. In our earlier study , we found a heterogeneous distribution of tick symbionts across space, with Rickettsia abundance being the most likely driver of the differences in observed bacterial communities in ticks.
Therefore, in our pairwise analysis exploring the co-infections between tick microorganisms in individual nymphs, we used a model controlling for location as well as tick identity.
Associations of vertically transmitted symbionts
We detected significant negative associations of R. helvetica with two other tick symbionts, M. mitochondrii and S. ixodetis, which were positively associated with each other (Fig. 3). Interestingly, in previous studies on the microbiome of I. ricinus, Aivelo et al.  and Lejal et al.  also detected a negative association between Rickettsia sp. and Spiroplasma. The former study did not investigate the association of Rickettsia sp. with M. mitochondrii, which was present in all samples, while the latter found a positive association.
It is possible that ticks cannot simultaneously maintain two symbiont species by vertical transmission, which may affect their distribution in tick populations. The phenomenon of interference between two vertically transmitted symbionts in hard ticks (Ixodidae) has been described in previous studies [10, 12, 65, 66]. Primary infection with one Rickettsia species has been suggested to block (by an unknown mechanism) transovarial transmission of a second Rickettsia species. However, to our knowledge, this so-called rickettsial interference has never been related to any other tick symbiont genus. If this interference between R. helvetica and S. ixodetis exists, to persist in a tick population, S. ixodetis would have to occasionally utilize horizontal transmission, which has been suggested but not undisputedly proven [67, 68].
The mutual exclusion among symbionts may also indicate that the roles which they play in ticks overlap and that a double infection comes with fitness costs without providing sufficient increases in fitness. For instance, in Acyrthosiphon pisum, both Serratia symbiotica and Hamiltonella defensa provide resistance to parasitoids. Although a double infection provided increased resistance in the laboratory, it was rarely observed under natural conditions .
The positive association between M. mitochondrii and S. ixodetis associations imply that these pairs of symbionts serve complementary functions and/or exhibit different tissue tropism in ticks. Different host tissues may constitute distinct microhabitats and be nutritionally favourable, immunotolerant, and easy to colonize for some but not for other microorganisms .
Associations of horizontally with vertically transmitted microorganisms
Our analyses showed that the most prevalent tick symbiont, M. mitochondrii, was positively associated with the pathogens B. burgdorferi s.l. and N. mikurensis, and negatively associated with the pathogen A. phagocytophilum (Fig. 3). Although these specific associations have not been reported previously, M. mitochondrii has been shown to be a successful colonizing partner of pathogenic R. parkeri in A. maculatum ticks .
It is possible that M. mitochondrii, which infects female ticks more often than males , enhances the ticks’ probability of becoming infected with zoonotic pathogens. Generally, immature female I. ricinus are larger and take lengthier and more extensive blood meals than males [72, 73], which might facilitate the acquisition of a pathogen from a host.
Additionally, the observed positive associations of the M. mitochondrii symbionts with B. burgdorferi s.l. and N. mikurensis could be indirect in nature and arise from this symbiont, providing fitness benefits to ticks. It is suggested that M. mitochondrii increases tick survival by supplying essential nutrients, which enhances reproductive fitness, in addition to benefiting energy production, maintenance of homeostasis and water balance, and antioxidant defence .
However, in our study, it appears that this facilitation occurs only for B. burgdorferi s.l. and N. mikurensis, which circulate among small mammals such as rodents and birds [75,76,77,78], and not for A. phagocytophilum. The acquisition of A. phagocytophilum occurs mostly while ticks feed on deer , a host fed upon by nymphal rather than larval tick stages . Therefore, it is possible that at the adult stage, where ticks had a chance to feed twice (once as a larva and once as a nymph), a positive association between A. phagocytophilum and M. mitochondrii could be detected.
Associations of horizontally transmitted pathogens
Different species of horizontally transmitted symbionts can be acquired from either a single host species or multiple species of vertebrate hosts [18, 77,78,79, 81], which is often reflected in their co-infection patterns in nymphs that fed only once as larvae. For instance, we observed a strong positive association between B. burgdorferi s.l. and N. mikurensis, which we expected given that both N. mikurensis and B. afzelii (the most prevalent genospecies of B. burgdorferi s.l. in I. ricinus in the Netherlands ) are amplified by rodents, which are often simultaneously infected with both microorganisms [68, 75, 76].
Anaplasma phagocytophilum may be acquired by I. ricinus from various vertebrate species, but it is most commonly acquired from roe deer . Therefore, the significant negative association of A. phagocytophilum with B. burgdorferi s.l. and N. mikurensis pathogens implies that nymphs that fed on deer had no chance to acquire rodent-borne pathogens .
Although here a few specific associations did not have a large effect on the composition of the entire microbiome, they can still be relevant for tick-borne pathogen dynamics. We observed that tick symbionts are heterogeneously distributed across the Dutch tick population and revealed a novel pattern of associations. Although the underlying mechanisms are unknown, they are important for understanding the role of vertically transmitted symbionts in the control of tick-borne diseases.
For instance, our findings imply that ticks carrying M. mitochondrii contribute more to the transmission cycle and the acarological risk of Lyme borreliosis and neoehrlichiosis than ticks lacking M. mitochondrii. Future studies on microorganism interference should investigate the dynamics of symbionts through the ontogeny of pre-infected ticks in comparison with non-infected ticks, applying bacterial quantification methods and identifying the microbiomes of ticks during different life events such as moulting, questing, and feeding.
Availability of data and materials
The datasets used and/or analysed during the current study are available in the supplementary data files and from the corresponding author on reasonable request.
Cirimotich CM, Ramirez JL, Dimopoulos G. Native microbiota shape insect vector competence for human pathogens. Cell Host Microbe. 2011;10:307–10.
Weiss B, Aksoy S. Microbiome influences on insect host vector competence. Trends Parasitol. 2011;27:514–22.
Wang GH, Du J, Chu CY, Madhav M, Hughes GL, Champer J. Symbionts and gene drive: two strategies to combat vector-borne disease. Trends Genet. 2022;38:708–23.
Herren JK, Mbaisi L, Mararo E, Makhulu EE, Mobegi VA, Butungi H, et al. A microsporidian impairs Plasmodium falciparum transmission in Anopheles arabiensis mosquitoes. Nat Commun. 2020;11:1–10.
de la Fuente J, Estrada-Pena A, Venzal JM, Kocan KM, Sonenshine DE. Overview: ticks as vectors of pathogens that cause disease in humans and animals. Front Biosci. 2008;13:6938–46.
Rochlin I, Toledo A. Emerging tick-borne pathogens of public health importance: a mini-review. J Med Microbiol. 2020;69:781–91.
Noda H, Munderloh UG, Kurtti TJ. Endosymbionts of ticks and their relationship to Wolbachia spp. and tick-borne pathogens of humans and animals. Appl Environ Microbiol. 1997;63:3926–32.
Sacchi L, Bigliardi E, Corona S, Beninati T, Lo N, Franceschi A. A symbiont of the tick Ixodes ricinus invades and consumes mitochondria in a mode similar to that of the parasitic bacterium Bdellovibrio bacteriovorus. Tissue Cell. 2004;36:43–53.
Scoles GA. Phylogenetic analysis of the Francisella-like endosymbionts of Dermacentor ticks. J Med Entomol. 2004;41:277–286.
Burgdorfer W, Hayes S, Thomas L, Lancaster Jr J. New spotted fever group Rickettsia from the lone star tick, Amblyomma americanum. Rickettsiae and rickettsial diseases/edited by W. Burgdorfer; RL Anacker. 1981.
Noden BH, Roselli MA, Loss SR. Rickettsia parkeri and Candidatus Rickettsia andeanae in Amblyomma maculatum group ticks. Emerg Infect Dis. 2020;26:371.
Paddock CD, Denison AM, Dryden MW, Noden BH, Lash RR, Abdelghani SS, et al. High prevalence of “Candidatus Rickettsia andeanae” and apparent exclusion of Rickettsia parkeri in adult Amblyomma maculatum (Acari: Ixodidae) from Kansas and Oklahoma. Ticks Tick Borne Dis. 2015;6:297–302.
Gall CA, Reif KE, Scoles GA, Mason KL, Mousel M, Noh SM, et al. The bacterial microbiome of Dermacentor andersoni ticks influences pathogen susceptibility. ISME J. 2016;10:1846–55.
Budachetri K, Kumar D, Crispell G, Beck C, Dasch G, Karim S. The tick endosymbiont Candidatus Midichloria mitochondrii and selenoproteins are essential for the growth of Rickettsia parkeri in the Gulf Coast tick vector. Microbiome. 2018;6:141.
Duron O, Binetruy F, Noel V, Cremaschi J, McCoy KD, Arnathau C, et al. Evolutionary changes in symbiont community structure in ticks. Mol Ecol. 2017;26:2905–21.
Kjelland V, Stuen S, Skarpaas T, Slettan A. Prevalence and genotypes of Borrelia burgdorferi sensu lato infection in Ixodes ricinus ticks in southern Norway. Scand J Infect Dis. 2010;42:579–85.
van Duijvendijk G, Coipan C, Wagemakers A, Fonville M, Ersöz J, Oei A, et al. Larvae of Ixodes ricinus transmit Borrelia afzelii and B. miyamotoi to vertebrate hosts. Parasit Vectors. 2016;9:1–7.
Burri C, Schumann O, Schumann C, Gern L. Are Apodemus spp. mice and Myodes glareolus reservoirs for Borrelia miyamotoi, Candidatus Neoehrlichia mikurensis, Rickettsia helvetica, R. monacensis and Anaplasma phagocytophilum? Ticks Tick Borne Dis. 2014;5:245–51.
van Duijvendijk G, Sprong H, Takken W. Multi-trophic interactions driving the transmission cycle of Borrelia afzelii between Ixodes ricinus and rodents: a review. Parasit Vectors. 2015;8:643.
Silaghi C, Beck R, Oteo JA, Pfeffer M, Sprong H. Neoehrlichiosis: an emerging tick-borne zoonosis caused by Candidatus Neoehrlichia mikurensis. Exp Appl Acarol. 2016;68:279–97.
Krause PJ, Fish D, Narasimhan S, Barbour AG. Borrelia miyamotoi infection in nature and in humans. Clin Microbiol Infect. 2015;21:631–9.
Pickard JM, Zeng MY, Caruso R, Núñez G. Gut microbiota: role in pathogen colonization, immune responses, and inflammatory disease. Immunol Rev. 2017;279:70–89.
Yordanova IA, Zakovic S, Rausch S, Costa G, Levashina E, Hartmann S. Micromanaging immunity in the murine host vs. the mosquito vector: microbiota-dependent immune responses to intestinal parasites. Front Cell Infect Microbiol. 2018;8:308.
Moutailler S, Valiente Moro C, Vaumourin E, Michelet L, Tran FH, Devillers E, et al. Co-infection of ticks: the rule rather than the exception. PLoS Negl Trop Dis. 2016;10:e0004539.
Raulf MK, Jordan D, Fingerle V, Strube C. Association of Borrelia and Rickettsia spp. and bacterial loads in Ixodes ricinus ticks. Ticks Tick Borne Dis. 2018;9:18–24.
Steiner FE, Pinger RR, Vann CN, Grindle N, Civitello D, Clay K, et al. Infection and co-infection rates of Anaplasma phagocytophilum variants, Babesia spp., Borrelia burgdorferi, and the rickettsial endosymbiont in Ixodes scapularis (Acari: Ixodidae) from sites in Indiana, Maine, Pennsylvania, and Wisconsin. J Med Entomol. 2014;45:289–97.
Abraham NM, Liu L, Jutras BL, Yadav AK, Narasimhan S, Gopalakrishnan V, et al. Pathogen-mediated manipulation of arthropod microbiota to promote infection. Proc Natl Acad Sci USA. 2017;114:E781–90.
Narasimhan S, Schuijt TJ, Abraham NM, Rajeevan N, Coumou J, Graham M, et al. Modulation of the tick gut milieu by a secreted tick protein favors Borrelia burgdorferi colonization. Nat Commun. 2017;8:184.
Hofmeester T, Coipan E, Van Wieren S, Prins H, Takken W, Sprong H. Few vertebrate species dominate the Borrelia burgdorferi sl life cycle. Environ Res Lett. 2016;11:043001.
Krawczyk AI, Röttjers L, Fonville M, Takumi K, Takken W, Faust K, et al. Quantitative microbial population study reveals geographical differences in bacterial symbionts of Ixodes ricinus. Microbiome. 2022;10:1–15.
Takumi K, Sprong H, Hofmeester TR. Impact of vertebrate communities on Ixodes ricinus-borne disease risk in forest areas. Parasit Vectors. 2019;12:1–12.
Wielinga PR, Gaasenbeek C, Fonville M, de Boer A, de Vries A, Dimmers W, et al. Longitudinal analysis of tick densities and Borrelia, Anaplasma, and Ehrlichia infections of Ixodes ricinus ticks in different habitat areas in The Netherlands. Appl Environ Microbiol. 2006;72:7594–601.
Heylen D, Adriaensen F, Van Dongen S, Sprong H, Matthysen E. Ecological factors that determine Ixodes ricinus tick burdens in the great tit (Parus major), an avian reservoir of Borrelia burgdorferi s.l. Int J Parasitol. 2013;43:603–11.
Hovius JW, de Wever B, Sohne M, Brouwer MC, Coumou J, Wagemakers A, et al. A case of meningoencephalitis by the relapsing fever spirochaete Borrelia miyamotoi in Europe. Lancet. 2013;382:658.
Jahfari S, Fonville M, Hengeveld P, Reusken C, Scholte EJ, Takken W, et al. Prevalence of Neoehrlichia mikurensis in ticks and rodents from North-west Europe. Parasit Vectors. 2012;5:74.
Courtney JW, Kostelnik LM, Zeidner NS, Massung RF. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J Clin Microbiol. 2004;42:3164–8.
Øines Ø, Radzijevskaja J, Paulauskas A, Rosef O. Prevalence and diversity of Babesia spp. in questing Ixodes ricinus ticks from Norway. Parasit Vectors. 2012;5:156.
Heylen D, Tijsse E, Fonville M, Matthysen E, Sprong H. Transmission dynamics of Borrelia burgdorferi s.l. in a bird tick community. Environ Microbiol. 2013;15:663–73.
Bogaert D, Keijser B, Huse S, Rossen J, Veenhoven R, Van Gils E, et al. Variability and diversity of nasopharyngeal microbiota in children: a metagenomic analysis. PLoS ONE. 2011;6:e17035.
Biesbroek G, Sanders EAM, Roeselers G, Wang X, Caspers MPM, Trzciński K, et al. Deep sequencing analyses of low density microbial communities: working at the boundary of accurate microbiota detection. PLoS ONE. 2021;7:e32942.
Hasrat R, Kool J, de Steenhuijsen Piters WA, Chu MLJN, Kuiling S, Groot JA, et al. Benchmarking laboratory processes to characterise low-biomass respiratory microbiota. Sci Rep. 2021;11:1–3.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2020.
Oksanen J, Kindt R, Legendre P, O’Hara B, Stevens MHH, Oksanen MJ, et al. The vegan package. Community Ecol Package. 2007;10:719.
Wickham H, Chang W, Henry L, Pedersen T, Takahashi K, Wilke C, et al. ggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics; R package version 3.3.0. 2020.
Anderson MJ. Permutational multivariate analysis of variance (PERMANOVA). Wiley statsref: statistics reference online. 2014;1–15.
Hofmeester TR, Sprong H, Jansen PA, Prins HHT, van Wieren SE. Deer presence rather than abundance determines the population density of the sheep tick, Ixodes ricinus, in Dutch forests. Parasit Vectors. 2017;10:433.
Ovaskainen O, Tikhonov G, Norberg A, Blanchet FG, Duan L, Dunson D, et al. How to make more out of community data? A conceptual framework and its implementation as models and software. Ecol Lett. 2017;20:561–76.
Pollock LJ, Tingley R, Morris WK, Golding N, O’Hara RB, Parris KM, et al. Understanding co-occurrence by modelling species simultaneously with a Joint Species Distribution Model (JSDM). Methods Ecol Evol. 2014;5:397–406.
Ovaskainen O, Abrego N, Halme P, Dunson D. Using latent variable models to identify large networks of species-to-species associations at different spatial scales. Methods Ecol Evol. 2016;7:549–55.
Dallas TA, Laine AL, Ovaskainen O. Detecting parasite associations within multi-species host and parasite communities. Proc Royal Soc B. 1912;2019:20191109.
Tikhonov G, Opedal ØH, Abrego N, Lehikoinen A, de Jonge MM, Oksanen J, et al. Joint species distribution modelling with the R-package Hmsc. Methods Ecol Evol. 2020;11:442–7.
Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;1:457–72.
Wei T, Simko V, Levy M, Xie Y, Jin Y, Zemla J. Package ‘corrplot.’ Statistician. 2017;56:e24.
Plantard O, Bouju-Albert A, Malard MA, Hermouet A, Capron G, Verheyden H. Detection of Wolbachia in the tick Ixodes ricinus is due to the presence of the hymenoptera endoparasitoid Ixodiphagus hookeri. PLoS ONE. 2012;7:e30692.
Tijsse-Klasen E, Braks M, Scholte EJ, Sprong H. Parasites of vectors–Ixodiphagus hookeri and its Wolbachia symbionts in ticks in The Netherlands. Parasit Vectors. 2011;4:228.
Hawlena H, Rynkiewicz E, Toh E, Alfred A, Durden LA, Hastriter MW, et al. The arthropod, but not the vertebrate host or its environment, dictates bacterial community composition of fleas and ticks. ISME J. 2013;7:221.
Swei A, Kwan JY. Tick microbiome and pathogen acquisition altered by host blood meal. ISME J. 2017;11:813–6.
Dessauer HC. Blood chemistry of reptiles: physiological and evolutionary aspects. Biol Reptil. 1970;3:1–72.
Hayes BM, Radkov AD, Yarza F, Flores S, Kim J, Zhao Z, et al. Ticks resist skin commensals with immune factor of bacterial origin. Cell. 2020;183:1562–71.
Carpi G, Cagnacci F, Wittekindt NE, Zhao F, Qi J, Tomsho LP, et al. Metagenomic profile of the bacterial communities associated with Ixodes ricinus ticks. PLoS ONE. 2011;6:e25604.
Clay K, Klyachko O, Grindle N, Civitello D, Oleske D, Fuqua C. Microbial communities and interactions in the lone star tick, Amblyomma americanum. Mol Ecol. 2008;17:4371–81.
Van Treuren W, Ponnusamy L, Brinkerhoff RJ, Gonzalez A, Parobek CM, Juliano JJ, et al. Variation in the microbiota of Ixodes ticks with regard to geography, species, and sex. Appl Environ Microbiol. 2015;81:6200–9.
Aivelo T, Norberg A, Tschirren B. Bacterial microbiota composition of Ixodes ricinus ticks: the role of environmental variation, tick characteristics and microbial interactions. PeerJ. 2019;7:e8217.
Lejal E, Chiquet J, Aubert J, Robin S, Estrada-Peña A, Rue O, et al. Temporal patterns in Ixodes ricinus microbial communities: an insight into tick-borne microbe interactions. Microbiome. 2021;9:1–20.
Macaluso KR, Sonenshine DE, Ceraul SM, Azad AF. Rickettsial infection in Dermacentor variabilis (Acari: Ixodidae) inhibits transovarial transmission of a second Rickettsia. J Med Entomol. 2002;39:809–13.
Sakai RK, Costa FB, Ueno TE, Ramirez DG, Soares JF, Fonseca AH, et al. Experimental infection with Rickettsia rickettsii in an Amblyomma dubitatum tick colony, naturally infected by Rickettsia bellii. Ticks Tick Borne Dis. 2014;5:917–23.
Binetruy F, Bailly X, Chevillon C, Martin OY, Bernasconi MV, Duron O. Phylogenetics of the Spiroplasma ixodetis endosymbiont reveals past transfers between ticks and other arthropods. Ticks Tick Borne Dis. 2019;10:575–84.
Krawczyk AI, Van Duijvendijk GL, Swart A, Heylen D, Jaarsma RI, Jacobs FH, et al. Effect of rodent density on tick and tick-borne pathogen populations: consequences for infectious disease risk. Parasit Vectors. 2020;13:1–17.
Oliver KM, Moran NA, Hunter MS. Costs and benefits of a superinfection of facultative symbionts in aphids. Proc Royal Soc B. 2006;273:1273–80.
Goto S, Anbutsu H, Fukatsu T. Asymmetrical interactions between Wolbachia and Spiroplasma endosymbionts coexisting in the same insect host. Appl Environ Microbiol. 2006;72:4805–10.
Lo N, Beninati T, Sassera D, Bouman EA, Santagati S, Gern L, et al. Widespread distribution and high prevalence of an alpha-proteobacterial symbiont in the tick Ixodes ricinus. Environ Microbiol. 2006;8:1280–7.
Daveu R, Laurence C, Bouju-Albert A, Sassera D, Plantard O. Symbiont dynamics during the blood meal of Ixodes ricinus nymphs differ according to their sex. Ticks Tick Borne Dis. 2021;12:101707.
Dusbabek F. Nymphal sexual dimorphism in the sheep tick Ixodes ricinus (Acari: Ixodidae). Folia Parasitol (Praha). 1996;43:75–9.
Olivieri E, Epis S, Castelli M, Boccazzi IV, Romeo C, Desirò A, et al. Tissue tropism and metabolic pathways of Midichloria mitochondrii suggest tissue-specific functions in the symbiosis with Ixodes ricinus. Ticks Tick Borne Dis. 2019;10:1070–7.
Andersson M, Bartkova S, Lindestad O, Raberg L. Co-infection with 'Candidatus Neoehrlichia Mikurensis’ and Borrelia afzelii in Ixodes ricinus ticks in southern Sweden. Vector Borne Zoonotic Dis. 2013;13:438–42.
Andersson M, Scherman K, Raberg L. Infection dynamics of the tick-borne pathogen “Candidatus Neoehrlichia mikurensis” and coinfections with Borrelia afzelii in bank voles in Southern Sweden. Appl Environ Microbiol. 2014;80:1645–9.
Hanincova K, Schafer SM, Etti S, Sewell HS, Taragelova V, Ziak D, et al. Association of Borrelia afzelii with rodents in Europe. Parasitology. 2003;126:11–20.
Hanincova K, Taragelova V, Koci J, Schafer SM, Hails R, Ullmann AJ, et al. Association of Borrelia garinii and B. valaisiana with songbirds in Slovakia. Appl Environ Microbiol. 2003;69:2825–30.
Jahfari S, Coipan EC, Fonville M, van Leeuwen AD, Hengeveld P, Heylen D, et al. Circulation of four Anaplasma phagocytophilum ecotypes in Europe. Parasit Vectors. 2014;7:365.
Pacilly FCA, Benning ME, Jacobs F, Leidekker J, Sprong H, Van Wieren SE, et al. Blood feeding on large grazers affects the transmission of Borrelia burgdorferi sensu lato by Ixodes ricinus. Ticks Tick Borne Dis. 2014;5:810–7.
Coipan CE, van Duijvendijk GLA, Hofmeester TR, Takumi K, Sprong H. The genetic diversity of Borrelia afzelii is not maintained by the diversity of the rodent hosts. Parasit Vectors. 2018;11:454.
Coipan EC, Jahfari S, Fonville M, Maassen CB, van der Giessen J, Takken W, et al. Spatiotemporal dynamics of emerging pathogens in questing Ixodes ricinus. Front Cell Infect Microbiol. 2013;3:36.
This study was financially supported by the Dutch Ministry of Health, Welfare and Sport (VWS), and by a grant from the European Interreg North Sea Region Programme, as part of the NorthTick project. The funders had no role in the study design and interpretation or the decision to submit the work for publication.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Details on studied forest sites. Table S2. The qPCR protocol for detection of tick symbionts and pathogens. Table S3. Bacterial loads of samples for the 16S rRNA sequencing; sample scheme. Table S4. Details on a positive control, primers, protocol, and in silico analysis used in the 16S rRNA quantification analyses. Table S5. OTU table. Table S6. Prevalence of veritcally and horizontally transmitted symbionts and pathogens per location. Table S7. Expected and observed co-infections in questing nymphs and results of the exact binomial test. Table S8. Tick microorganism species-to-species associations (co-infections) in ticks.
Abundance of the most abundant taxa, separated by quartiles of absolute bacterial density (16S rRNA content in ng/µL). Figure S2. Principal coordinate analysis of Bray–Curtis dissimilarities for tick microbiomes compared to blanks.
About this article
Cite this article
Krawczyk, A.I., Röttjers, S., Coimbra-Dores, M.J. et al. Tick microbial associations at the crossroad of horizontal and vertical transmission pathways. Parasites Vectors 15, 380 (2022). https://doi.org/10.1186/s13071-022-05519-w
- Ixodes ricinus
- Transmission dynamics
- Lyme borreliosis
- Tick-borne diseases