Silvia-Diana Borşan1*Angela Monica Ionică1,2Clémence Galon3
Andra Toma-Naic1
Cosmin Peştean4Attila D. Sándor1,5Sara Moutailler3Andrei Daniel Mihalca1- 1Department of Parasitology and Parasitic Diseases, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, Cluj-Napoca, Romania
- 2“Regele Mihai I al României” Life Sciences Institute, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, Cluj-Napoca, Romania
- 3UMR BIPAR, Animal Health Laboratory, ANSES, INRAE, Ecole Nationale Vétérinaire d’Alfort, Paris-Est Sup, Maisons-Alfort, France
- 4Department of Surgery, Anesthesiology and Intensive Therapy, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, Cluj-Napoca, Romania
- 5Department of Parasitology and Zoology, University of Veterinary Medicine, Budapest, Hungary
Despite the increasingly recognized eco-epidemiological importance of ticks as vectors for numerous zoonotic pathogens in urban areas, data regarding the pathogen diversity and co-infection rates in ticks and wildlife hosts in urban and peri-urban Romania are scanty. We aimed to establish the risk of human exposure to co-infected ticks in Cluj-Napoca, a major city in Romania. DNA was isolated from 151 questing ticks: Ixodes ricinus (n = 95), Haemaphysalis punctata (n = 53), Dermacentor reticulatus (n = 2), and Dermacentor marginatus (n = 1); 222 engorged ticks: I. ricinus (n = 164), I. hexagonus (n = 36), H. punctata (n = 16), H. concinna (n = 6), and 70 tissue samples collected from wildlife hosts during 2018 in five urban, and two peri-urban sites. Using a pre-designed Fluidigm real-time PCR dynamic array, all DNA samples were individually screened for the presence of 44 vector-borne pathogens. Subsequently, conventional PCRs were performed for a selection of samples to allow validation and sequencing. In total, 15 pathogens were identified to species and 6 to genus level. In questing ticks, single infections were more common than co-infections. Seven Borrelia spp. were detected in questing I. ricinus, and three in H. punctata ticks. An overall high prevalence 26.35% (95% CI: 19.46–34.22) and diversity of Borrelia burgdorferi sensu lato was seen in urban questing ticks. Other pathogens of the order Rickettsiales were present with variable prevalence. Co-infections occurred in 27.4% (95% CI: 18.72-37.48) of all infected questing ticks. In engorged ticks the overall Bo. burgdorferi sensu lato prevalence was 35.6% (95% CI: 29.29–42.27), with five species present. Pathogens of the order Rickettsiales were also frequently detected. We report for the first time in Romania the presence of Rickettsia aeschlimannii and Rickettsia felis. Overall, from the infected engorged ticks, 69.2% showcased co-infections. In Ixodes spp., dual co-infections, namely Borrelia spp. and Anaplasma phagocytophilum, and Rickettsia helvetica and A. phagocytophilum were the most prevalent. Given the outcome, we underline the need to establish proper tick-surveillance programs in cities and include co-infections in the management plan of tick-borne diseases in Romania.
Introduction
Ticks are arthropods that can transmit pathogenic microorganisms including protozoa, bacteria, and viruses. In Europe, the majority of human and animal arthropod-borne diseases are vectored by ticks (Jongejan and Uilenberg, 2004; Colwell et al., 2011). Both humans and pets face a significantly higher risk of contracting tick-borne pathogens (TBPs) due to the emergence of ticks in urban areas (Rizzoli et al., 2014).
Ixodes ricinus is the predominant tick species reported in Europe (Rizzoli et al., 2014) and the most widespread questing tick species collected in Romania’s natural (Mihalca et al., 2012a) and urban habitats (Borşan et al., 2020). Moreover, it is also the most prevalent tick species reported to bite humans in Romania (Briciu et al., 2014; Andersson et al., 2018).
The habitat range of I. ricinus includes both natural and urban environments such as recreational areas, parks, and gardens, which can ensure the abiotic and biotic requirements for optimal development of the off-host stages (Rizzoli et al., 2014).
To date, a considerable number of studies describe the pathogens vectored by I. ricinus worldwide (Keesing et al., 2010; Rizzoli et al., 2014; Strnad et al., 2017). The TBPs which pose the greatest risk for the public health are the spirochetes of the Borrelia burgdorferi sensu lato complex, the causative agents of human Lyme borreliosis (LB), and the European tick-borne encephalitis virus, which can lead to tick-borne encephalitis (TBE) (Gritsun et al., 2003; Rizzoli et al., 2011). In Romania, various molecular approaches have been used to detect the prevalence of infection with the Borrelia spp. (including the relapsing fever spirochete B. miyamotoi) (Kalmár et al., 2016) in questing ticks (Kalmár et al., 2013; Raileanu et al., 2017), ticks collected from animal hosts (Gherman et al., 2012; Dumitrache et al., 2015; Kalmár et al., 2019), or humans (Andersson et al., 2014; Briciu et al., 2014; Andersson et al., 2018; Kalmár et al., 2020). During 2018 a total of 532 human Lyme disease cases were confirmed by serology in Romania (NCSCC, 2018), placing the country at the inferior position of the European incidence (Rizzoli et al., 2011). Other pathogens vectored by I. ricinus that are of rising importance for medical and veterinary health are bacteria of the order Rickettsiales. Despite the wide distribution of Anaplasma phagocytophilum across Romania (Matei et al., 2015), and of the presence of this pathogen in I. ricinus ticks collected from humans (Matei et al., 2017), no clinical human cases were reported so far in the country. Yet to be described from humans in Romania, nonetheless detected in questing (Kalmár et al., 2016) and engorged I. ricinus ticks collected from humans (Andersson et al., 2014; Kalmár et al., 2020), Neoehrlichia mikurensis is an emerging pathogen which can either lead to severe febrile illness in immunocompromised patients (Grankvist et al., 2014) or fever in clinically healthy humans (Li et al., 2012). The spotted fever group rickettsiae (SFG) cause rickettsioses in humans (Parola et al., 2013). In Romania, R. conorii (the Mediterranean Spotted Fever), R. massiliae, and R. slovaca and R. raoultii (SENLAT-scalp eschar and neck lymphadenopathy after a tick bite syndrome) were reported from human patients (Serban et al., 2009; Zaharia et al., 2016). While neither of these bacterial species are vectored by I. ricinus, several other Rickettsia spp. were also identified in ticks collected from wildlife hosts (Mărcuţan et al., 2016; Sándor et al., 2017) or the environment (Ioniţa et al., 2013). The role of I. ricinus is also suspected in the human transmission of Bartonella species such as B. quintana and B. henselae (Socolovschi et al., 2009; Vu Hai et al., 2014). This tick species is also a vector for zoonotic apicomplexans of the genus Babesia, such as B. divergens and B. microti (Gray et al., 2010). To date, B. microti and B. venatorum are reported in ticks collected from humans in Romania (Kalmár et al., 2020).
Due to the generalist feeding behavior of I. ricinus, co-infections with several micro-organisms are frequent in this tick (Reis et al., 2011). Multiple strains of bacteria, parasites, and viruses can be acquired by ticks either from a host with multiple infections, through feeding on subsequent hosts (along with the individual development), or through co-feeding mechanisms (Piesman and Happ, 2001). Transstadial, or in the case of some TBPs (i.e., Borrelia spp., Rickettsia spp., and TBE-complex virus), transovarial transmission in ticks can also contribute to the ecology of such pathogens (Sprong et al., 2009; Rizzoli et al., 2011; Karbowiak and Biernat, 2016). It is noteworthy that the transmission of pathogens from co-infected ticks is likely to alter the severity of clinical signs in humans or animals (Cutler et al., 2020), sometimes causing delays or errors in diagnosis as reported for concurrent babesiosis and Lyme disease (Grunwaldt et al., 1983; Golightly et al., 1989). Due to the potential implications of co-infections in urban ticks and the likelihood of co-transmission of TBPs it is vital to identify local enzootic cycles, especially in recreational areas.
Co-infection prevalence in ticks in European countries ranges from 3.2% to 45% (Reye et al., 2010; Reis et al., 2011; Lommano et al., 2012; Moutailler et al., 2016; Klitgaard et al., 2019; Nebbak et al., 2019; Kalmár et al., 2020). Little is known about the co-infection rates in questing ticks and wildlife hosts in urban and peri-urban Romania (Raileanu et al., 2017, 2018). Effective tick-based surveillance is essential for monitoring human and animal disease emergence. Therefore, by using a powerful broad-spectrum high-throughput approach, our study aimed to: 1. detect the TBPs in questing, engorged ticks, and in tissue samples from wildlife fauna collected in seven locations in Cluj-Napoca, Romania; 2. determine the co- infection rates; and 3. perform a comparative statistical analysis of infection rates and pathogen diversity in ticks from urban and peri-urban habitats.
Materials and Methods
Sampling Protocol
The characteristics of the urban and peri-urban locations assessed in this study, the method of collection of questing ticks, the sampling and trapping protocols for urban wildlife hosts and their associated tick fauna, followed by species-specific identification of all organisms as well as the research and ethical permits are described in detail elsewhere (Borşan et al., 2020). The five urban locations assessed during 2018 consisted of two parks: “Iuliu Haţieganu Park” and the campus of the “University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca (USAMV Campus); two gardens: “Alexandru Borza Botanical Garden” and a centrally located private garden; and “Mănăştur Cemetery”. The peri-urban sites were represented by “Hoia” and “Făget” forest. All locations (except the private garden) are open to the public year-round. The sampling activities were performed from March until November 2018 and included flagging (bimestrial) to collect questing ticks and collecting of wildlife using standardized methods (i.e., rodent trapping with “snap-traps”; bird sampling using ornithological mist nets, and hedgehog sampling by “torch-based” searches). All wildlife hosts were searched for ticks, blood samples were collected from birds and hedgehogs (if feasible), followed by release, while the trapped rodents underwent necropsy.
Questing Ticks
From the 3383 total questing ticks collected during 2018, we randomly selected individual ticks for DNA isolation as follows: 10% (or 5 individual ticks if the numbers were too low to meet the 10% criteria) from each tick species and developmental stage (questing larvae were excluded) per location during each month of sampling. If the tick number (from the same species, stage, location, and month) was less than 5, all sampled ticks were included. A total of (n = 95) Ixodes ricinus (n = 77 ticks collected in urban, and n = 18 in peri-urban sites) and (n = 53) Haemaphysalis punctata ticks (n = 32 ticks collected in urban, and n = 21 in peri-urban sites) were selected. In addition, we included two individuals of Dermacentor reticulatus and one individual of D. marginatus accidentally collected in peri-urban locations during the flagging campaigns. Overall, 151 questing ticks were used for DNA isolation (Supplementary File 1).
Urban Wildlife Hosts
All the engorged ticks collected from rodents and birds were individually tested, while in the case of ticks collected from hedgehogs (Erinaceus roumanicus), samples were selected using a similar algorithm as for questing ticks (including larvae). Therefore, we selected for DNA isolation 222 engorged ticks (n = 20 ticks from rodents; n = 22 ticks from birds; n = 180 ticks from hedgehogs) as follows: a total of 215 ticks consisting of I. ricinus (n = 157), I. hexagonus (n = 36), H. punctata (n = 16) and H. concinna (n = 6) collected from wildlife hosts found in urban sites, and 7 I. ricinus ticks found on hosts from peri-urban areas) (Table 1 and Supplementary File 1).
Table 1. Wildlife and associated tick species collected in Cluj-Napoca during 2018.
Following the visual inspection, all micromammals except hedgehogs underwent necropsy. During the necropsy, the heart, liver tissue, and two skin biopsies (one from the interscapular region and the second from the ear pavilion that were pooled together for DNA isolation) were collected from each animal. All the tissue samples were labeled accordingly and stored at −20°C. During the anesthesia -protocol described in Borşan et al. (2020), 0.5 ml of blood were collected from the jugular or saphenous vein of eight hedgehogs. Blood sampling was unsuccessful in the case of three hedgehogs. Regarding the birds (Borşan et al., 2020), a 150 μl blood sample was collected from the brachial vein using a microcapillary tube. All blood samples were stored in a 3.2% citrate tube at −20°C. Overall, DNA was individually isolated from the heart, liver, and skin biopsy tissues from 29 mammal hosts, 33 blood samples from birds, and eight blood samples from hedgehogs (Supplementary File 1).
DNA Isolation
The genomic DNA isolation was performed individually for all tissues and tick samples using the Isolate II Genomic DNA Kit (Bioline, London, United Kingdom), according to the manufacturer’s instructions. Each tick was dried, cut into halves using a sterile scalpel blade, and crushed with a sterile pestle. For tissue samples, up to 25 mg of tissue was cut into small pieces and crushed with a sterile pestle, before the lysis. To ensure proper lysis, overnight digestion was performed for both ticks and tissue samples. The blood samples were processed using the same kit. A quantity of 200 μl of blood was used from the hedgehog samples, and 100 μl for birds. The DNA was stored at −20°C until further processing.
Detection of Tick-Borne Pathogens
DNA Pre-amplification
The DNA pre-amplification steps were followed as described in (Michelet et al., 2014).
High-Throughput Real-Time PCR System
The BioMark real-time polymerase chain reaction (PCR) system (Fluidigm, San Francisco, CA, United States), was used for high-throughput microfluidic real-time PCR amplification using the 48.48 dynamic arrays. The chips dispensed 48 PCR assays and 48 samples into individual wells, after which on-chip microfluidics assemble PCR reactions in individual chambers before thermal cycling resulting in 2304 individual reactions (Michelet et al., 2014).
Subsequent pre-amplification, Real-Time PCR were performed using FAM- and black hole quencher (BHQ1)-labeled TaqMan probes with PerfeCTa qPCR ToughMix, Low ROX (QuantaBio, Beverly, MA, United States) following the protocol by Michelet et al. (2014). Thermal cycling conditions were as follows: 50°C for 2 min, 95°C for 10 min, 40 cycles of 2-step amplification at 95°C for 15 s, and 60°C for 1 min. Data were acquired on the BioMark Real-Time PCR system and processed using the Fluidigm Real-Time PCR Analysis software to obtain a cut-off (Ct) value (Michelet et al., 2014; Gondard et al., 2020).
The BioMark real-time PCR system (Fluidigm, San Francisco, CA, United States) was used for high-throughput microfluidic real-time PCR for the most common bacterial and parasitic TBP species known to circulate or recently emerging in Europe. The real-time PCR system developed for the screening of known and potential TBPs in Romanian ticks included 47 sets of primers and probes (Michelet et al., 2014; Sprong et al., 2019; Gondard et al., 2020). Among them, 37 primers were used for the detection of pathogens to species level (n = 30 bacterial and n = 7 apicomplexan) and 8 primers to genus level (n = 5 bacterial and n = 3 apicomplexan). Three sets of primers and probes were used for the molecular identification of two tick species found in Romania: I. ricinus and D. reticulatus. Lastly, a primer targeting a conserved region of the 16S rRNA gene in ticks, called “Tick spp.” was used as a control for the DNA extraction. To determine if factors present in the sample could inhibit the PCR, the Escherichia coli strain EDL933 DNA was added to each sample as an internal inhibition control (Nielsen and Andersen, 2003) (Supplementary File 2).
Validation of the BioMark Real-Time PCR Results
Conventional or nested PCRs using primers that targeted different genes or regions than those of the BioMark system were performed on several samples presenting low Ct values or co-infections with multiple pathogen species with individual low Ct values. Each reaction was carried out in a 25 μl reaction volume containing 12.5 μl of 2x Green PCR Mastermix (Rovalab, GmbH, Teltow, Germany), 1 μM of each primer, and 4 μl of DNA sample. The amplification reactions were carried out in C1000 Thermal Cyclers (Bio-Rad, CA, United States), using previously published primers and protocols (Michelet et al., 2014) (Table 2).
Table 2. Primer sets used for pathogen DNA amplification by conventional PCR/nested PCR.
Amplicons were further processed by sequencing (performed by Macrogen Europe B.V., Amsterdam, Netherlands), for the final confirmation of pathogen species. Thus, the results obtained by the real-time microfluidic PCR assay for 89 Borrelia spp., 63 Anaplasma spp., 89 Rickettsia spp., and 36 Bartonella spp. samples were re-tested by conventional or nested PCRs. Following the molecular analysis, 80 Borrelia spp., 53 A. phagocytophilum, and 44 Rickettsia spp. samples were sequenced. Identity percentages of the sequences obtained with reference sequences available in GenBank (NCBI) are presented (Table 3).
Table 3. Homology between obtained sequences and reference sequences in GenBank.
The sequences were compared to other GenBank entries by BLAST (Basic Local Alignment Search Tool) analysis and further submitted to the GenBank under the following accession numbers: Bo. afzelii (MW272725, MW272726, MW272727, MW272728, MW272729, MW272730, MW272731, MW272732, MW272733, MW272734); Bo. garinii (MW272735, MW272736, MW272737, MW272738, MW272739, MW272740); Bo. lusitaniae (MW272741, MW272742); Bo. spielmanii (MW272743, MW272744); Bo. valaisiana (MW272745, MW272746, MW272747); Bo. bavariensis (MW272749, MW272750); Bo. miyamotoi (MW272748); A. phagocytophilum (MW272751, MW272752); R. helvetica (MW272753, MW272755, MW272756, MW272757); R. monacensis (MW272758); R. aeschlimannii (MW272754). Five Borrelia spp. samples showcased equal identity percentages with both Bo. bavariensis and Bo. garinii reference sequences. Since this title formulation is not accepted by GenBank the two individual sequences were submitted as Bo. bavariensis.
Statistical Analysis
Statistical calculations were performed using EpiInfo 7 software (CDC, Atlanta, GE, United States). The prevalence was established and differences between various groups were assessed using chi-square tests. All differences were considered statistically significant for p < 0.05.
Results
By using the microfluidic PCR assay 443 samples were analyzed: 151 questing ticks, 222 engorged ticks, and 70 tissue and blood samples collected from urban wildlife. Overall, by considering the results of both methods used (microfluidic PCR, conventional, and nested PCRs), of the targeted pathogens, 15 were detected to species level, and 6 to genus level in the seven locations in Cluj-Napoca (Supplementary File 1).
Questing Ticks
Borrelia spp. in Ixodes ricinus
Borrelia burgdorferi s.l DNA was detected in 36/95 (37.9%) of questing I. ricinus ticks collected from all sampling sites. Infection rates were not significantly different between females 43.8% (14/32); males 36.3% (12/33), and nymphs 33.3% (10/30) (χ2 = 1.65; d.f. = 2; p = 0.437), or between ticks collected in urban areas compared to peri-urban sites (χ2 = 0.5; d.f. = 1; p = 0.475). Seven species of Borrelia were identified: Bo. afzelii (12.6%), Bo. lusitaniae (9.5%), Bo. garinii (7.4%), Bo. spielmanii (3.2%), Bo. burgdorferi s.s. (2.1%), Bo. valaisiana (2.1%), Bo. bavariensis/Bo. garinii (1.05%), and Bo. miyamotoi (1.05%), with statistically significant differences among locations only in the case of I. ricinus ticks from Mănăştur Cemetery and Făget forest which were more frequently infected with Bo. lusitaniae (χ2 = 21.10; d.f. = 6; p = 0.001). All aforementioned Borrelia spp. were detected in ticks in the urban locations, while Bo. afzelii and Bo. lusitaniae were the only two species present in peri-urban ticks (Figure 1 and Supplementary File 3).
Figure 1. The distribution of Borrelia spp. in questing and engorged ticks in the seven locations assessed in Cluj-Napoca. (1) USAMV Campus; (2) Mănăştur Cemetery; (3) Iuliu Haţieganu Park; (4) Alexandru Borza Botanical Garden; (5) Hoia forest; (6) Făget forest; (7) Private garden.
Other Tick-Borne Pathogens in Ixodes ricinus
Anaplasma phagocytophilum DNA had a prevalence of 24.2% (23/95) in questing I. ricinus ticks: females 21.9% (7/32), males 33.3% (11/33), nymphs 16.7% (5/30), and was found in ticks collected from all the seven locations, without statistically significant differences among locations (χ2 = 4.38; d.f. = 6; p = 0.624), or stages (χ2 = 2.52; d.f. = 2; p = 0.283), except for Iuliu Haţieganu Park, which recorded a significantly higher prevalence of infection in adult ticks: 40% in females (95% CI: 5.27–85.34) and 50% in males (95% CI: 1.26–98.74) (χ2 = 6.21; d.f. = 2; p = 0.044) compared to the other locations.
Two Rickettsia species were found in I. ricinus ticks. Rickettsia helvetica had an overall prevalence of 22.1% (21/95): females 18.75% (6/32), males 18.2% (6/33), nymphs 30% (9/30), without statistical differences among life stages (χ2 = 1.59; d.f. = 2; p = 0.451). This species was detected only in urban I. ricinus ticks (χ2 = 4.81; d.f. = 1; p = 0.028). Rickettsia monacensis had a prevalence of 12.6% (12/95): females 18.75% (6/32), males 12.1% (4/33), nymphs 6.7% (2/30), without statistical differences among life stages (χ2 = 2.06; d.f. = 2; p = 0.356), but with a statistically significant difference of the prevalence rate among locations (χ2 = 13.16; d.f. = 6; p = 0.040), respectively higher for the peri-urban environment (16.7%) (95% CI: 3.58–41.42) (Figure 2 and Supplementary File 3).
Figure 2. The distribution of Rickettsia spp. in questing and engorged ticks in the seven locations assessed in Cluj-Napoca. (1) USAMV Campus; (2) Mănăştur Cemetery; (3) Iuliu Haţieganu Park; (4) Alexandru Borza Botanical Garden; (5) Hoia forest; (6) Făget forest; (7) Private garden.
Tick-Borne Pathogens in Other Tick Species
Three Borrelia spp. were found in 3/53 (5.7%) of the Haemaphysalis punctata ticks collected in three locations (Iuliu Haţieganu Park, Mănăştur Cemetery, and Hoia forest) as follows: Bo. afzelii (1.9%), Bo. garinii (1.9%), and Bo. lusitaniae (1.9%) (Supplementary File 3). Between locations, the prevalence of infection with Bo. burgdorferi s.l was statistically higher in H. punctata nymphs than adults (7.14%) (95% CI: 0.88–23.5) (χ2 = 15.43; d.f. = 4; p = 0.003).
The prevalence of infection with A. phagocytophilum in H. punctata ticks from Făget forest and Mănăştur Cemetery was 18.9% (10/53): females 9.09% (1/11); males 28.6% (4/14); and nymphs 17.9% (5/28). No statistical differences in prevalence among locations or life stages were recorded.
Three Rickettsia spp. were detected in H. punctata ticks: R. monacensis 15.09% (8/53): females 27.3% (3/11); males 0%; nymphs 17.9% (5/28); R. helvetica 13.2% (7/53): females 27.3% (3/11); males 21.4% (3/14); nymphs 3.6% (1/28); and R. conorii 1.9%: females 9.09% (1/11). No statistical differences were recorded for the prevalence of infection with Rickettsiales among locations or life stages in H. punctata ticks.
The only D. marginatus specimen included tested positive for A. phagocytophilum DNA. Also, both D. reticulatus individuals were positive for Rickettsia spp. DNA.
Overall, among the tick species analyzed, the prevalence of Bo. burgdorferi s.l. was statistically higher in questing ticks collected in urban areas (30.28%; 95% CI: 21.84–39.81) than peri-urban ones (15.38%; 95% CI: 5.86–30.53) (χ2 = 10.82; d.f. = 1; p < 0.001).
Urban Wildlife Hosts
Borrelia spp. in Engorged Ticks
Borrelia burgdorferi s.l. DNA was present in engorged ticks collected from urban wildlife as follows: I. ricinus 36.6% (60/164), of which 53.3% (8/15) in females, 46.7% (7/15) in males, 42.6% (26/61) in nymphs, and 26% in larvae (19/73); I. hexagonus 60% (21/35) of which 66.7% (18/27) in females and 37.5% (3/8) in males; H. punctata 6.25% (1/16) in nymphs. The overall prevalence of infection was 24.8% for Bo. afzelii (with statistically significant differences between locations (χ2 = 11.21; d.f. = 3; p = 0.01), 4.05% for B. garinii, 2.25% for Bo. spielmanii, 2.25% for Bo. valaisiana, 2.25% for Bo. bavariensis/Bo. garinii, and 1.8% for Bo. miyamotoi. The prevalence of Bo. afzelii was significantly higher in urban I. hexagonus (44.4%; 95% CI: 27.94–61.9) compared to I. ricinus (24.2; 95% CI: 17.73-31.67) (χ2 = 4.99; d.f. = 1; p = 0.025). There were statistically significant differences in the urban areas among the prevalence of Bo. burgdorferi s.l. between tick species (χ2 = 13.37; d.f. = 2; p = 0.001), with a significantly higher prevalence in I. hexagonus (58.3%; 95% CI: 40.76–74.49). Borrelia miyamotoi DNA was found in two co-feeding I. ricinus ticks (1 larva and 1 nymph) collected from the same E. roumanicus individual. Nevertheless, the blood sample collected from the respective hedgehog was negative for Bo. miyamotoi (Figure 1 and Supplementary File 3).
Other Tick-Borne Pathogens in Engorged Ticks
Anaplasma phagocytophilum DNA had a prevalence of 60.4% (99/164) in I. ricinus ticks, of which 66.7% (10/15) in females, 80% (12/15) in males, 68.9% (42/61) in nymphs, and 48% (35/73) in larvae, with significant differences in prevalence between life stages (χ2 = 9.2; d.f. = 3; p = 0.026). Also, 92.6% (25/27) of females and 100% (8/8) of males of I. hexagonus, and 25% (4/16; 2 nymphs and 2 larvae) of H. punctata tested positive for the presence of A. phagocytophilum DNA. The overall prevalence of A. phagocytophilum infection was significantly higher in I. hexagonus (91.7%; 95% CI: 77.54–98.25), compared to I. ricinus (60.4%; 95% CI: 52.44–67.91) (χ2 = 11.53; d.f. = 1; p = 0.0006). A percentage of 98.4% (134/136) of the A. phagocytophilum-positive engorged ticks were collected from urban hedgehogs. The remaining two ticks were collected from birds (Phylloscopus collybita and Turdus merula).
Rickettsia helvetica DNA was found in 30.5% (50/164) of I. ricinus ticks (53.3%, 8/15 females; 53.3% 8/15 males; 29.5%, 18/61 nymphs, and 21.9% 16/73 larvae) and 52.8% (19/36) of I. hexagonus ticks (44.4% 12/27 females, and 87.5% 7/8 males) with statistical differences regarding the prevalence of infection between the two species (χ2 = 5.54; d.f. = 1; p = 0.018), and among locations (χ2 = 29.93; d.f. = 3; p = 0). R. monacensis was the second most prevalent Rickettsia spp., detected in 6.1% of I. ricinus, 8.3% of I. hexagonus, 12.5% of H. punctata, and 16.7% of H. concinna ticks. Also, in USAMV Campus, one I. ricinus nymph tested positive for the presence of R. felis and one H. concinna nymph for R. aeschlimannii DNA (Figure 2). Hepatozoon spp. had a 0.9% prevalence in the engorged ticks analyzed, in which we also detected Theileria spp., with a 2.25% prevalence (Supplementary File 3).
Urban Wildlife Tissue Samples
Rodents
Rodents were considered positive if pathogenic DNA was detected in any of the tissue samples collected. Therefore, three Borrelia spp. (Bo. afzelii, Bo. spielmanii, and Bo. miyamotoi) were found in 13.8% (4/29) of rodents, as follows: the skin biopsy and heart tissue of one Arvicola terrestris from USAMV Campus, and the skin biopsy of one Apodemus agrarius from Iuliu Haţieganu Park were positive for Bo. afzelii, while both Bo. miyamotoi, and Bo. spielmanii DNA were detected in the skin biopsy of two different A. agrarius (one pathogen/rodent) from Iuliu Haţieganu Park. Anaplasma phagocytophilum DNA was detected in the heart tissue of one A. flavicollis. Also, the prevalence of R. monacensis was 17.2% (5/29), detected individually in the skin biopsies of three A. agrarius and one Apodemus sylvaticus from USAMV Campus, and one Mus musculus from Hoia forest, while N. mikurensis had a prevalence of 2.9% (2/29) in the skin biopsy and liver tissue of one A. terrestris from USAMV Campus, and in the heart and liver tissue of one A. agrarius from Iuliu Haţieganu Park (Supplementary File 1).
Birds
Borrelia afzelii was detected in the blood sample of one urban Parus major 3% (1/33). Also, 30.3% (10/33) of urban birds [Corvus frugilegus (n = 1), Erithacus rubecula (n = 1), Garrulus glandarius (n = 1), Parus major (n = 3), Sturnus vulgaris (n = 1), and Turdus merula (n = 3)] harbored A. phagocytophilum DNA, 12.1% (4/33) R. helvetica DNA [Corvus frugilegus (n = 1), and Parus major (n = 3)], and 3% (1/33) R. monacensis DNA [Turdus merula (n = 1)] (Supplementary File 1).
Hedgehogs
Of the eight blood samples tested, two were positive for A. phagocytophilum DNA, and one for R. helvetica DNA.
Co-infections Between Tick-Borne Pathogens
Questing Ticks
Ixodes ricinus
Co-infections occurred in 34.3% (23/67) of all I. ricinus infected ticks. Co-infection prevalence was 36% (9/25) in females, 29.2% (7/24) in males and 38.9% (7/18) in nymphs. A statistically significant difference was recorded regarding the prevalence of co-infections among life stages in the peri-urban sites: 0% in females; 25% (95% CI: 0.63–80.59) in males; and 100% (95% CI: 15.81–100) in nymphs (χ2 = 7.21; d.f. = 2; p = 0.027). The most frequent dual co-infections were between Rickettsia spp. and Borrelia spp., followed by Rickettsia spp. and A. phagocytophilum, and A. phagocytophilum and Borrelia spp. (Figure 3).
Figure 3. Single and multiple infections (percent of total collected ticks) detected in questing Ixodes ricinus (n = 95) and H. punctata (n = 53) ticks collected in Cluj-Napoca.
Co-infections with three pathogens were less common and consisted of combinations between R. helvetica, A. phagocytophilum, and Borrelia. spp. (Table 4).
Table 4. Co-infections with tick-borne pathogens in questing ticks in Cluj-Napoca.
Haemaphysalis punctata
Co-infections were present in 20% (6/30) of all infected H. punctata ticks. Co-infection prevalence was 50% (4/8) in females, and 13.3% (2/15) in nymphs, with a statistically higher prevalence of infection in females 57.14% (95% CI: 18.41–90.1), compared to nymphs (0% prevalence) in urban areas (χ2 = 6.85; d.f. = 2; p = 0.032). The most prevalent co-infections were between Rickettsia spp. and Borrelia spp., followed by Rickettsia spp. and A. phagocytophilum (Table 4 and Figure 3). One D. reticulatus from Făget forest was co-infected with Borrelia spp. and Rickettsia spp.
Globally, there were no statistically significant correlations between the co-infection rates in questing ticks and the environment (urban/peri-urban) (χ2 = 0.49; d.f. = 1; p = 0.483). Despite the lack of other significant correlations, the urban areas showed a more diverse array of pathogen species compared to the peri-urban sites (Figures 1, 2).
Engorged Ticks
Of the engorged infected ticks, 69.2% (108/156) were co-infected with various TBPs.
Ixodes spp.
From the total infected I. ricinus ticks, 69% (79/115) were co-infected with multiple pathogens. Of these 72.7% (8/11) were females, 91.7% (11/12) males, 66% (33/50) nymphs and 65.1% (27/42) were larvae. Of the I. hexagonus analyzed, 79.4% (27/34) were co-infected: 76.9% (20/26) were females, and 87.5% (7/8) males (Figure 4). The most frequent co-infections in Ixodes spp. were dual co-infections equally prevalent with A. phagocytophilum and Borrelia spp., and A. phagocytophilum and R. helvetica, followed by infections with three pathogens namely R. helvetica, A. phagocytophilum, and Borrelia spp. (Supplementary File 4).
Figure 4. Co-infections detected in engorged Ixodes ricinus (n = 164), I. hexagonus (n = 36), Haemaphysalis punctata (n = 16), and H. concinna (n = 6) ticks (percent of total collected ticks) collected from wildlife hosts in Cluj-Napoca.
Haemaphysalis spp.
One H. punctata showed a co-infection with A. phagocytophilum and Bo. afzelii, while another one was co-infected with A. phagocytophilum and R. monacensis.
All the co-infected engorged ticks were collected from hosts found in urban locations (Figure 4 and Supplementary File 4). Overall, there were significant differences among the engorged tick species between locations in terms of co-infection rates (χ2 = 8.6; d.f. = 2; p = 0.013).
Urban Wildlife Tissue Samples
Co-infections were also detected in two tissue samples collected from rodents. One A. agrarius from Iuliu Haţieganu Park tested positive for Bo. miyamotoi and N. mikurensis, and one A. terrestris from USAMV Campus was co-infected with Bo. afzelii and N. mikurensis. Also, three samples collected from birds (one Parus major and one Turdus merula from Alexandru Borza Botanical Graden, and one Corvus frugilegus from Iuliu Haţieganu Park) presented dual co-infections with A. phagocytophilum and Rickettsia spp., while one sample (Parus major from Alexandru Borza Botanical Garden) showed a triple co-infection with A. phagocytophilum, R. helvetica, and Bo. afzelii (Table 5).
Table 5. Co-infections with tick-borne pathogens in tissue samples collected from urban wildlife hosts in Cluj-Napoca.
Discussion
Using the BioMark system we performed a comprehensive survey of the various TBPs that co-circulate in tick-host cycles in five urban, and two peri-urban locations in Cluj-Napoca, a major city in Romania. Even though local prevalence studies have limited value in terms of epidemiological risk assessment, the prevalence of Bo. burgdorferi s.l. spirochetes in questing ticks has been considered an essential element of risk assessment for Lyme borreliosis (LB) (Rauter and Hartung, 2005). Therefore, collectively, the results of this study may have important implications in terms of public health, especially for the urban areas, since until recently LB risk was considered to be correlated with residency in rural areas (Rizzoli et al., 2011).
Nowadays, a higher LB incidence is registered in urban environments (3.2%/100.000 inhabitants) than in rural settlements (2.5%/100.000 inhabitants) in Romania (NCSCC, 2018). Increased access to information and better accessibility and addressability to medical services of the urban population could explain the more common reference of patients to hospitals for diagnostic purposes, hence, the higher incidence (NCSCC, 2018).
The overall Borrelia burgdorferi s.l. prevalence (37.9%) in questing I. ricinus ticks across all locations assessed in this study was much higher compared to previous data on the prevalence of Lyme spirochetes reported in questing ticks from Romania by conventional PCR studies (3.8–18%) (Coipan and Vladimirescu, 2011; Kalmár et al., 2013), but similar to the prevalence reported in Iaşi county (25.8%) by microfluidic real-time PCR (Raileanu et al., 2017).
Since some clinical manifestations are specific to particular Borrelia species, their prevalence in a certain area is important for risk assessment (Strnad et al., 2017). As in our study, a Europe-wide meta-analysis of Bo. burgdorferi s.l. species in questing ticks (Estrada-Peña et al., 2018) and previous reports in Romania (Kalmár et al., 2013) showed that the most prevalent Borrelia spp. are Bo. afzelii, and Bo. garinii.
Accounting for most of the LB human cases in Europe, Bo. afzelii is mostly isolated from medium-sized and small rodents (Coipan et al., 2018). Borrelia garinii is commonly hosted by birds (Dubska et al., 2009), particularly species that can reach high densities in urban sites (Taragel’ová et al., 2008). Nevertheless, a distinct and highly pathogenic ecotype of Bo. garinii, now confirmed to species status, Bo. bavariensis (formerly known as OspA type 4) uses rodents as reservoir hosts (Huegli et al., 2002; Margos et al., 2013). The relatively high prevalence of these Borrelia spp. could be linked to the high diversity and abundance of rodent species in Romania (Mihalca et al., 2012b). Commonly associated with birds (Rizzoli et al., 2011), Bo. valaisiana had a lower prevalence in questing ticks from Cluj-Napoca than in Iaşi county (Raileanu et al., 2017) and was only present in urban locations.
We report higher Bo. lusitaniae infection rates in questing ticks in Cluj-Napoca compared to previous studies in Romania (Kalmár et al., 2013; Raileanu et al., 2017). Borrelia lusitaniae is mainly associated with lizards (Rizzoli et al., 2011). Heltai et al. (Heltai et al., 2015) reported that cemeteries contribute significantly to the habitats of lizards in urban areas due to the presence of stony habitats, their size, heterogeneity, and reduced levels of human disturbance. Thus, the significantly higher prevalence of infection with Bo. lusitaniae of questing ticks from Mănăştur Cemetery is most likely linked to the presence and abundance of lizards (confirmed through visual inspection - data not shown) in the respective site.
Borrelia miyamotoi, the only relapsing fever agent transmitted by Ixodes species in Europe (Cutler et al., 2019), was initially reported in questing ticks in central Romania (Kalmár et al., 2016). Recently, Raileanu et al. confirmed the low infection rate in ticks from eastern Romania (Raileanu et al., 2017). We also confirm the presence of Bo. miyamotoi at a low prevalence in ticks and wildlife hosts in recreational areas of Cluj-Napoca.
Borrelia spp. infection and co-infection rates were not significantly different between adults and nymphs of questing I. ricinus ticks. The prevalence of Bo. burgdorferi s.l. in Europe is higher in adults than in nymphs (Strnad et al., 2017). Nonetheless, nymphs are mainly responsible for transmitting Borrelia spp. to humans (Rizzoli et al., 2011) and can be encountered in suburban and urban environments (Pejchalová et al., 2007), and even roadsides (Haemig et al., 2008).
The public health relevance of H. punctata ticks is considered to be rather limited (Briciu et al., 2014). Moreover, its vectorial role for Borrelia spp. has not been clearly demonstrated. As in our study, LB spirochetes have formerly been reported in questing H. punctata (Tälleklint, 1996), but at lower prevalence compared to I. ricinus.
The presence of A. phagocytophilum has been investigated in detail in Romania (Matei et al., 2015; Kalmár et al., 2016; Raileanu et al., 2017; Raileanu et al., 2018). The higher prevalence of A. phagocytophilum infection in adult I. ricinus compared to nymphs could be linked to the greater number of bloodmeals, since transstadial transmission of A. phagocytophilum is improbable (Raileanu et al., 2018). Therefore, despite no reports of human infection with A. phagocytophilum in Romania so far, the risk of acquiring this pathogen following tick bites in recreational areas is possible.
In the current study, the dominant SFG Rickettsia species in questing ticks from both urban and peri-urban areas were R. helvetica and R. monacensis, in concordance with former mentions in urban sites in Romania (Raileanu et al., 2018), and Europe (Rizzoli et al., 2014; Kowalec et al., 2019). Here we report for the first time the presence of R. aeschlimannii in Romania, identified in H. concinna. Rickettsia aeschlimannii is an emerging human and animal pathogen, reported from various ticks in Europe and Africa, including several Hyalomma spp. ticks collected from migrant bird species (Parola et al., 2013; Chisu et al., 2016). Another SFG Rickettsia present in this study is R. felis which is known to be transmitted via cat flea bites (Brown and Macaluso, 2016). Despite seldom reports regarding the presence and prevalence of this TBP, several other studies have also identified R. felis in questing I. ricinus ticks in Europe (Vayssier-Taussat et al., 2013; Lejal et al., 2019). A previous study in Italy (Ciervo et al., 2006) also reports the presence of the human pathogen R. conorii in H. punctata questing ticks. However, further research is required to assess the vectorial competence of Haemaphysalis spp. and I. ricinus and their implication in the transmission of these SFG Rickettsia spp.
Nowadays, thanks to more sensitive and efficient detection tools, co-infections with different TBPs are more frequently reported in ticks (Michelet et al., 2014; Raileanu et al., 2017; Lejal et al., 2019; Gondard et al., 2020). The pathogenesis and aftermath of co-infections in humans is a complex process that still needs further research (Baneth, 2014). Pathogens can synergistically colonize more favorably their hosts through processes initiated by co-transmission and entering of multiple pathogens inside the respective host’s organism (Baneth, 2014).
We report a high prevalence of co-infections in both questing and engorged ticks. Contrary to previous reports which mentioned Bo. afzelii and Bo. garinii as the most prevalent co-infection in ticks in Romania (Raileanu et al., 2017), or France (Moutailler et al., 2016), we hereby detected co-infections among Rickettsia spp. and Borrelia spp., followed by Rickettsia spp. and A. phagocytophilum, and A. phagocytophilum and Borrelia spp., as most prevalent. A high prevalence of co-infection with R. helvetica and A. phagocytophilum in I. ricinus ticks was also reported by Lejal et al. (Lejal et al., 2019) who suggested that the superior acclimatization of these two pathogens in ticks might portend them as stronger competitors than other pathogen species. Our findings also uncover new risks for urban inhabitants since the co-infection with A. phagocytophilum and Bo. burgdorferi s.l. was shown to enhance the colonization ability of Bo. burgdorferi s.l. (Grab et al., 2007).
Despite our statistical analysis revealing a significantly higher prevalence of pathogens in urban ticks compared to the peri-urban ones, the sample sizes assessed were uneven, as the majority of samples were collected from urban sites. This was not related to bias in the selection algorithm but because of the availability and abundance of wildlife hosts and questing ticks when the collection was performed. A previous study we conducted in these seven recreational locations in Cluj-Napoca showed a higher abundance of ticks in the urban versus the peri-urban locations, linked to the abundance and diversity of local wildlife species, particularly hedgehogs (E. roumanicus) (Borşan et al., 2020). These findings may explain the more diverse TBPs community we detected in urban sites, further highlighting the importance of urban dwellers such as hedgehogs, rodents, and birds in the ecology of tick-borne diseases. Given the results of the present study, E. roumanicus could facilitate pathogen exchange among infected and uninfected ticks without displaying a systemic infection (through co-feeding mechanisms, pathogens stationed in tissues rather than in the bloodstream) (Randolph, 2011; Voordouw, 2015), and therefore can be considered an amplifier host and an epidemiologically important wildlife species for the urban environment (Jahfari et al., 2017).
Since I. ricinus ticks have a high affinity for biting humans and the level of co-infections detected in this tick species in Cluj-Napoca is high, co-transmission and enhanced disease severity in humans are possible scenarios for the city inhabitants (Moutailler et al., 2016).
Conclusion
The most noteworthy outcomes of this study are (1) the detection of a high prevalence of Bo. burgdorferi s.l. in urban questing ticks; (2) the overall great diversity and prevalence of TBPs in engorged ticks collected from urban sites (3) co-infections were frequent in both questing and engorged ticks.
Therefore, additional tick-surveillance and awareness programs should be implemented, especially in recreational areas, since the TBPs detected in ticks in Cluj-Napoca pose a significant risk to human health.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.
Ethics Statement
The flagging, sampling, and trapping campaigns were performed with the consent of Cluj-Napoca City Hall (Technical Department: Decision No. 13345/442/18.01.2018 and Urban Ecology and Green Spaces Department: Decision No. 13.351/10.01.2018) and the support of the Romanian Ornithological Society (SOR), Babeş Bolyai University (Decision No. 210/090/2018), the University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca and the owners of the private garden. All the activities were performed according to ethical permits and national legislation.
Author Contributions
AM and S-DB conceptualized the study. S-DB, AT-N, and AS performed the field work. CP executed the hedgehog anesthesia and sampling. S-DB performed the necropsy examination of the rodents, morphological identification of the ticks, and wrote the first draft of the manuscript. S-DB and AI executed the laboratory work and data analysis. AI implemented the statistical analysis of results and designed the maps. CG and SM designed the microfluidic PCR protocol, processed the samples, and interpreted the results. AM reviewed the manuscript for important intellectual content. All authors have read and approved the final manuscript.
Funding
This study was performed under the framework of the VectExel project “Multidisciplinary One Health excellence research platform for neglected and emerging vector-borne diseases,” project number 57 PCCDI/2018, grant agency The Executive Unit for Funding Higher Education and University Scientific Research (UEFISCSU) Romania.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Many thanks to DVM Mircea Coroian for its help during laboratory work.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.645002/full#supplementary-material
References
Alberti, A., Addis, M. F., Sparagano, O., Zobba, R., Chessa, B., Cubeddu, T., et al. (2005). Anaplasma phagocytophilum, Sardinia, Italy. Emerg. Infect. Dis. 11, 1322–1324. doi: 10.3201/eid1108.050085
Andersson, M., Zaghdoudi-Allan, N., Tamba, P., Stefanache, M., and Chitimia, L. (2014). Co-infection with ‘Candidatus Neoehrlichia mikurensis’ and Borrelia afzelii in an Ixodes ricinus tick that has bitten a human in Romania. Ticks Tick Borne Dis. 5, 706–708. doi: 10.1016/j.ttbdis.2014.05.013
Andersson, M. O., Marga, G., Banu, T., Dobler, G., and Chitimia-Dobler, L. (2018). Tick-borne pathogens in tick species infesting humans in Sibiu County, central Romania. Parasitol. Res. 117, 1591–1597. doi: 10.1007/s00436-018-5848-0
Baneth, G. (2014). Tick-borne infections of animals and humans: a common ground. Int. J. Parasitol. 44, 591–596. doi: 10.1016/j.ijpara.2014.03.011
Barbour, A. G., Carter, C. J., Bundoc, V., and Hinnebusch, J. (1996). The nucleotide sequence of a linear plasmid of Borrelia burgdorferi reveals similarities to those of circular plasmids of other prokaryotes. J. Bacteriol. 178, 6635–6639. doi: 10.1128/jb.178.22.6635-6639.1996
Borşan, S. D., Toma-Naic, A., Péter, Á, Sándor, A. D., Peştean, C., and Mihalca, A. D. (2020). Impact of abiotic factors, habitat type and urban wildlife on the ecology of hard ticks (Acari: Ixodidae) in urban and peri-urban habitats. Parasit. Vectors 13:476. doi: 10.1186/s13071-020-04352-3
Briciu, V. T., Meyer, F., Sebah, D., Ţăţulescu, D. F., Coroiu, G., Lupşe, M., et al. (2014). Real-time PCR-based identification of Borrelia burgdorferi sensu lato species in ticks collected from humans in Romania. Ticks Tick Borne Dis. 5, 575–581. doi: 10.1016/j.ttbdis.2014.04.007
Brown, L. D., and Macaluso, K. R. (2016). Rickettsia felis, an emerging flea-borne rickettsiosis. Curr. Trop. Med. Rep. 3, 27–39. doi: 10.1007/s40475-016-0070-6
Chisu, V., Zobba, R., Foxi, C., Pisu, D., Masala, G., and Alberti, A. (2016). Molecular detection and groEL typing of Rickettsia aeschlimannii in Sardinian ticks. Parasitol. Res. 115, 3323–3328. doi: 10.1007/s00436-016-5091-5
Ciervo, A., Simeoni, J., Khoury, C., Mancini, F., and Ciceroni, L. (2006). Spotted fever group Rickettsiae in Ixodes ricinus and Haemaphysalis punctata ticks in Italy. Eur. J. Inflamm. 4, 177–180. doi: 10.1177/1721727X0600400306
Coipan, C. E., van Duijvendijk, G. L., Hofmeester, T. R., Takumi, K., and Sprong, H. (2018). The genetic diversity of Borrelia afzelii is not maintained by the diversity of the rodent hosts. Parasit. Vectors 11:454. doi: 10.1186/s13071-018-3006-2
Coipan, E. C., and Vladimirescu, A. F. (2011). Ixodes ricinus ticks (Acari: Ixodidae): vectors for Lyme disease spirochetes in Romania. Exp. Appl. Acarol. 54, 293–300. doi: 10.1007/s10493-011-9438-4
Colwell, D. D., Dantas-Torres, F., and Otranto, D. (2011). Vector-borne parasitic zoonoses: emerging scenarios and new perspectives. Vet. Parasitol. 182, 14–21. doi: 10.1016/j.vetpar.2011.07.012
Cutler, S., Vayssier-Taussat, M., Estrada-Peña, A., Potkonjak, A., Mihalca, A. D., and Zeller, H. (2019). A new Borrelia on the block: Borrelia miyamotoi–a human health risk? Euro. Surveill. 24:1800170. doi: 10.2807/1560-7917.ES.2019.24.18.1800170
Cutler, S. J., Vayssier-Taussat, M., Estrada-Peña, A., Potkonjak, A., Mihalca, A. D., and Zeller, H. (2020). Tick-borne diseases and co-infection: current considerations. Tiks Tick Borne Dis. 12:101607. doi: 10.1016/j.ttbdis.2020.101607
Dubska, L., Literak, I., Kocianova, E., Taragelova, V., and Sychra, O. (2009). Differential role of passerine birds in distribution of Borrelia spirochetes, based on data from ticks collected from birds during the postbreeding migration period in Central Europe. Appl. Environ. Microbiol. 75, 596–602. doi: 10.1128/AEM.01674-08
Dumitrache, M. O., Matei, I. A., Ionică, A. M., Kalmár, Z., D’Amico, G., Sikó-Barabási, S., et al. (2015). Molecular detection of Anaplasma phagocytophilum and Borrelia burgdorferi sensu lato genospecies in red foxes (Vulpes vulpes) from Romania. Parasit. Vectors 8:514. doi: 10.1186/s13071-015-1130-9
Estrada-Peña, A., Cutler, S., Potkonjak, A., Vassier-Tussaut, M., Van Bortel, W., Zeller, H., et al. (2018). An updated meta-analysis of the distribution and prevalence of Borrelia burgdorferi sl in ticks in Europe. Int. J. Health Geogr. 17:41. doi: 10.1186/s12942-018-0163-7
Gherman, C. M., Sándor, A. D., Kalmár, Z., Marinov, M., and Mihalca, A. D. (2012). First report of Borrelia burgdorferi sensu lato in two threatened carnivores: the Marbled polecat, Vormela peregusna and the European mink, Mustela lutreola (Mammalia: Mustelidae). BMC Vet. Res. 8:137. doi: 10.1186/1746-6148-8-137
Golightly, L. M., Hirschhorn, L. R., and Weller, P. F. (1989). Infectious disease rounds: fever and headache in a splenectomized woman. Rev. Infect. Dis. 11, 629–637. doi: 10.1093/clinids/11.4.629
Gondard, M., Delannoy, S., Pinarello, V., Aprelon, R., Devillers, E., Galon, C., et al. (2020). Upscaling the surveillance of tick-borne pathogens in the French Caribbean islands. Pathogens 9:176. doi: 10.3390/pathogens9030176
Grab, D. J., Nyarko, E., Barat, N. C., Nikolskaia, O. V., and Dumler, J. S. (2007). Anaplasma phagocytophilum-Borrelia burgdorferi coinfection enhances chemokine, cytokine, and matrix metalloprotease expression by human brain microvascular endothelial cells. Clin. Vaccine Immunol. 14, 1420–1424. doi: 10.1128/CVI.00308-07
Grankvist, A., Andersson, P. O., Mattsson, M., Sender, M., Vaht, K., Höper, L., et al. (2014). Infections with the tick-borne bacterium “Candidatus Neoehrlichia mikurensis” mimic noninfectious conditions in patients with B cell malignancies or autoimmune diseases. Clin. Infect. Dis. 58, 1716–1722. doi: 10.1093/cid/ciu189
Gray, J., Zintl, A., Hildebrandt, A., Hunfeld, K. P., and Weiss, L. (2010). Zoonotic babesiosis: overview of the disease and novel aspects of pathogen identity. Ticks Tick Borne Dis. 1, 3–10. doi: 10.1016/j.ttbdis.2009.11.003
Gritsun, T. S., Lashkevich, V. A., and Gould, E. A. (2003). Tick-borne encephalitis. Antiviral Res. 57, 129–146. doi: 10.1016/S0166-3542(02)00206-1
Grunwaldt, E., Barbour, A. G., and Benach, J. L. (1983). Simultaneous occurrence of babesiosis and Lyme disease. N. Engl. J. Med. 308:1166. doi: 10.1056/NEJM198305123081919
Haemig, P. D., Waldenstrom, J., and Olsen, B. (2008). Roadside ecology and epidemiology of tick-borne diseases. Scand. J. Infect. Dis. 40, 11–12. doi: 10.1080/00365540802270003
Heltai, B., Sály, P., Kovács, D., and Kiss, I. (2015). Niche segregation of sand lizard (Lacerta agilis) and green lizard (Lacerta viridis) in an urban semi-natural habitat. Amphib. Reptil. 36, 389–399. doi: 10.1163/15685381-00003018
Huegli, D., Hu, C. M., Humair, P. F., Wilske, B., and Gern, L. (2002). Apodemus species mice are reservoir hosts of Borrelia garinii OspA serotype 4 in Switzerland. J. Clin. Microbiol. 40, 4735–4737. doi: 10.1128/JCM.40.12.4735-4737.2002
Ioniţa, M., Mitrea, I. L., Pfister, K., Hamel, D., and Silaghi, C. (2013). Molecular evidence for bacterial and protozoan pathogens in hard ticks from Romania. Vet. Parasitol. 196, 71–76. doi: 10.1016/j.vetpar.2013.01.016
Jahfari, S., Ruyts, S. C., Frazer-Mendelewska, E., Jaarsma, R., Verheyen, K., and Sprong, H. (2017). Melting pot of tick-borne zoonoses: the European hedgehog contributes to the maintenance of various tick-borne diseases in natural cycles urban and suburban areas. Parasit. Vectors 10:134. doi: 10.1186/s13071-017-2065-0
Jongejan, F., and Uilenberg, G. (2004). The global importance of ticks. Parasitol Cambridge 129, S3–S14. doi: 10.1017/s0031182004005967
Kalmár, Z., Dumitrache, M. O., D’Amico, G., Matei, I. A., Ionică, A. M., Gherman, C. M., et al. (2020). Multiple Tick-borne pathogens in Ixodes ricinus ticks collected from humans in Romania. Pathogens 9:390. doi: 10.3390/pathogens9050390
Kalmár, Z., Mihalca, A. D., Dumitrache, M. O., Gherman, C. M., Magdaş, C., Mircean, V., et al. (2013). Geographical distribution and prevalence of Borrelia burgdorferi genospecies in questing Ixodes ricinus from Romania: a countrywide study. Ticks Tick Borne Dis. 4, 403–408. doi: 10.1016/j.ttbdis.2013.04.007
Kalmár, Z., Sándor, A. D., Matei, I. A., Ionică, A. M., D’Amico, G., and Gherman, C. M. (2019). Borrelia spp. in small mammals in Romania. Parasit. Vectors 12:461. doi: 10.1186/s13071-019-3713-3
Kalmár, Z., Sprong, H., Mihalca, A. D., Gherman, C. M., Dumitrache, M. O., Coipan, E. C., et al. (2016). Borrelia miyamotoi and Candidatus Neoehrlichia mikurensis in Ixodes ricinus ticks, Romania. Emerg. Infect. Dis. 22, 550–551. doi: 10.3201/eid2203.150140
Karbowiak, G., and Biernat, B. (2016). The role of particular tick developmental stages in the circulation of tick-borne pathogens affecting humans in Central Europe. 2. Tick-borne encephalitis virus. Ann. Parasitol. 62, 3–9. doi: 10.17420/ap6201.25
Keesing, F., Belden, L. K., Daszak, P., Dobson, A., Harvell, C. D., Holt, R. D., et al. (2010). Impacts of biodiversity on the emergence and transmission of infectious diseases. Nature 468, 647–952. doi: 10.1038/nature09575
Klitgaard, K., Kjær, L. J., Isbrand, A., Hansen, M. F., and Bødker, R. (2019). Multiple infections in questing nymphs and adult female Ixodes ricinus ticks collected in a recreational forest in Denmark. Ticks Tick Borne Dis. 10, 1060–1065. doi: 10.1016/j.ttbdis.2019.05.016
Kowalec, M., Szewczyk, T., Welc-Falêciak, R., Siński, E., Karbowiak, G., and Bajer, A. (2019). Rickettsiales Occurrence and Co-occurrence in Ixodes ricinus Ticks in Natural and Urban Areas. Microb. Ecol. 77, 809–904. doi: 10.1007/s00248-018-1269-y
Lejal, E., Moutailler, S., Šimo, L., Vayssier-Taussat, M., and Pollet, T. (2019). Tick-borne pathogen detection in midgut and salivary glands of adult Ixodes ricinus. Parasit. Vectors 12:152. doi: 10.1186/s13071-019-3418-7
Li, H., Jiang, J. F., Liu, W., Zheng, Y. C., Huo, Q. B., Tang, K., et al. (2012). Human infection with Candidatus Neoehrlichia mikurensis, China. Emerg. Infect. Dis. 18, 1636–1639. doi: 10.3201/eid1810.120594
Loh, S. M., Gillett, A., Ryan, U., Irwin, P., and Oskam, C. (2017). Molecular characterization of ‘Candidatus Borrelia tachyglossi’(family Spirochaetaceae) in echidna ticks, Bothriocroton concolor. Int. J. Syst. Evol. Microbiol. 67, 1075–1080. doi: 10.1099/ijsem.0.001929
Lommano, E., Bertaiola, L., Dupasquier, C., and Gern, L. (2012). Infections and coinfections of questing Ixodes ricinus ticks by emerging zoonotic pathogens in Western Switzerland. Appl. Environ. Microbiol. 78, 4606–4612. doi: 10.1128/AEM.07961-11
Mărcuţan, I. D., Kalmár, Z., Ionică, A. M., D’Amico, G., Mihalca, A. D., Vasile, C., et al. (2016). Spotted fever group rickettsiae in ticks of migratory birds in Romania. Parasit. Vectors 9:294. doi: 10.1186/s13071-016-1565-7
Margos, G., Wilske, B., Sing, A., Hizo-Teufel, C., Cao, W. C., Chu, C., et al. (2013). Borrelia bavariensis sp. nov. is widely distributed in Europe and Asia. Int. J. Syst. Evol. Microbiol. 63, 4284–4288. doi: 10.1099/ijs.0.052001-0
Matei, I. A., Kalmár, Z., Lupşe, M., D’Amico, G., Ionică, A. M., Dumitrache, M. O., et al. (2017). The risk of exposure to rickettsial infections and human granulocytic anaplasmosis associated with Ixodes ricinus tick bites in humans in Romania: a multiannual study. Ticks Tick Borne Dis. 8, 375–378. doi: 10.1016/j.ttbdis.2016.12.013
Matei, I. A., Kalmár, Z., Magda, C., Magdaş, V., Toriay, H., Dumitrache, M. O., et al. (2015). Anaplasma phagocytophilum in questing Ixodes ricinus ticks from Romania. Ticks Tick borne Dis. 6, 408–413. doi: 10.1016/j.ttbdis.2015.03.010
Michelet, L., Delannoy, S., Devillers, E., Umhang, G., Aspan, A., Juremalm, M., et al. (2014). High-throughput screening of tick-borne pathogens in Europe. Front. Cell. Infect. Microbiol. 4:103. doi: 10.3389/fcimb.2014.00103
Mihalca, A. D., Dumitrache, M. O., Sándor, A. D., Magdaş, C., Oltean, M., Györke, A., et al. (2012b). Tick parasites of rodents in Romania: host preferences, community structure and geographical distribution. Parasit. Vectors 5:266. doi: 10.1186/1756-3305-5-266
Mihalca, A. D., Gherman, C. M., Magdaş, C., Dumitrache, M. O., Györke, A., Sándor, A. D., et al. (2012a). Ixodes ricinus is the dominant questing tick in forest habitats in Romania: the results from a countrywide dragging campaign. Exp. Appl. Acarol. 58, 175–182. doi: 10.1007/s10493-012-9568-3
Moutailler, S., Valiente Moro, C., Vaumourin, E., Michelet, L., Tran, F. H., Devillers, E., et al. (2016). Co-infection of ticks: the rule rather than the exception. PLoS Negl. Trop. Dis. 10: e0004539. doi: 10.1371/journal.pntd.0004539
NCSCC (2018). Evolution Analysis of Communicable Diseases Under Surveillance in Romania. Report for 2018. Romania: National Center for Surveillance and Control of Communicable Diseases-NCSCC. In Romanian.
Nebbak, A., Dahmana, H., Almeras, L., Raoult, D., Boulanger, N., Jaulhac, B., et al. (2019). Co-infection of bacteria and protozoan parasites in Ixodes ricinus nymphs collected in the Alsace region, France. Ticks Tick Borne Dis. 10:101241. doi: 10.1016/j.ttbdis.2019.06.001
Nielsen, E. M., and Andersen, M. T. (2003). Detection and characterization of verocytotoxin-producing Escherichia coli by automated 5’ nuclease PCR assay. J. Clin. Microbiol. 41, 2884–2893. doi: 10.1128/JCM.41.7.2884-2893.2003
Norman, A. F., Regnery, R., Jameson, P., Greene, C., and Krause, D. C. (1995). Differentiation of Bartonella-like isolates at the species level by PCR-restriction fragment length polymorphism in the citrate synthase gene. J Clin Microbiol 33, 1797–1803. doi: 10.1128/JCM.33.7.1797-1803.1995
Parola, P., Paddock, C. D., Socolovschi, C., Labruna, M. B., Mediannikov, O., Kernif, T., et al. (2013). Update on tick-borne rickettsioses around the world: a geographic approach. Clin. Microbiol. Rev. 26, 657–702. doi: 10.1128/CMR.00032-13
Pejchalová, K., Zakovska, A., Mejzlíková, M., Halouzka, J., and Dendis, M. (2007). Isolation, cultivation and identification of Borrelia burgdorferi genospecies from Ixodes ricinus from the city of Brno, Czech Republic. Ann. Agric. Environ. Med. 14, 75–79.
Piesman, J., and Happ, C. M. (2001). The efficacy of co-feeding as a means of maintaining Borrelia burgdorferi: a North American model system. J. Vector Ecol. 26, 216–220.
Raileanu, C., Moutailler, S., Pavel, I., Porea, D., Mihalca, A. D., Savuta, G., et al. (2017). Borrelia diversity and co-infection with other tick borne pathogens in ticks. Front. Cell. Infect. Microbiol. 7:36. doi: 10.3389/fcimb.2017.00036
Raileanu, C., Moutailler, S., Porea, D., Oslobanu, L., Anita, D., Anita, A., et al. (2018). Molecular evidence of Rickettsia spp., Anaplasma phagocytophilum, and “Candidatus Neoehrlichia mikurensis” in Ticks from Natural and Urban Habitats in Eastern Romania. Vector Borne Zoonotic Dis. 18, 343–349. doi: 10.1089/vbz.2017.2221
Randolph, S. E. (2011). Transmission of tick-borne pathogens between co-feeding ticks: Milan Labuda’s enduring paradigm. Ticks Tick Borne Dis. 2, 179–182. doi: 10.1016/j.ttbdis.2011.07.004
Rauter, C., and Hartung, T. (2005). Prevalence of Borrelia burgdorferi sensu lato genospecies in Ixodes ricinus ticks in Europe: a metaanalysis. Appl. Environ. Microbiol. 71, 7203–7216. doi: 10.1128/AEM.71.11.7203-7216.2005
Regnery, R. L., Spruill, C. L., and Plikaytis, B. D. (1991). Genotypic identification of rickettsiae and estimation of intraspecies sequence divergence for portions of two rickettsial genes. J. Bacteriol. 173, 1576–1589. doi: 10.1128/jb.173.5.1576-1589.1991
Reis, C., Cote, M., Paul, R. E., and Bonnet, S. (2011). Questing ticks in suburban forest are infected by at least six tick-borne pathogens. Vector Borne Zoonotic Dis. 11, 907–916. doi: 10.1089/vbz.2010.0103
Reye, A. L., Hübschen, J. M., Sausy, A., and Muller, C. P. (2010). Prevalence and seasonality of tick-borne pathogens in questing Ixodes ricinus ticks from Luxembourg. Appl. Environ. Microbiol. 76, 2923–2931. doi: 10.1128/AEM.03061-09
Rizzoli, A., Hauffe, H., Carpi, G., Vourc, H. G., Neteler, M., and Rosa, R. (2011). Lyme borreliosis in Europe. Euro. Surveill. 16:19906.
Rizzoli, A., Silaghi, C., Obiegala, A., Rudolf, I., Hubálek, Z., Földvári, G., et al. (2014). Ixodes ricinus and its transmitted pathogens in urban and peri-urban areas in Europe: new hazards and relevance for public health. Front. Public Health 2:251. doi: 10.3389/fpubh.2014.00251
Sándor, A. D., Kalmár, Z., Matei, I., Ionică, A. M., and Mărcuţan, I. D. (2017). Urban breeding corvids as disseminators of ticks and emerging tick-borne pathogens. Vector Borne Zoonotic Dis. 17, 152–154. doi: 10.1089/vbz.2016.2054
Serban, R., Pistol, A., Neguţ, M., and Cucuiu, R. (2009). Rickettsia conorii infection in Romania, 2000-2008. Bacteriol. Virusol. Parazitol. Epidemiol. 54, 177–183.
Socolovschi, C., Mediannikov, O., Raoult, D., and Parola, P. (2009). Update on tick-borne bacterial diseases in Europe. Parasite 16, 259–273. doi: 10.1051/parasite/2009164259
Sprong, H., Fonville, M., van Leeuwen, A. D., Devillers, E., Ibanez-Justicia, A., Stroo, A., et al. (2019). Detection of pathogens in Dermacentor reticulatus in northwestern Europe: evaluation of a high-throughput array. Heliyon 5:e01270. doi: 10.1016/j.heliyon.2019.e01270
Sprong, H., Wielinga, P. R., Fonville, M., Reusken, C., Brandenburg, A. H., Borgsteede, F., et al. (2009). Ixodes ricinus ticks are reservoir hosts for Rickettsia helvetica and potentially carry flea-borne Rickettsia species. Parasit. Vectors 2:41. doi: 10.1186/1756-3305-2-41
Strnad, M., Hönig, V., Růžek, D., Grubhoffer, L., and Rego, R. O. (2017). Europe-wide meta-analysis of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks. Appl. Environ. Microbiol. 83:e00609–17. doi: 10.1128/AEM.00609-17
Tälleklint, L. (1996). Lyme borreliosis spirochetes in Ixodes ricinus and Haemaphysalis punctata ticks (Acari: Ixodidae) on three islands in the Baltic Sea. Exp. Appl. Acarol. 20, 467–476. doi: 10.1007/BF00053310
Taragel’ová, V., Koči, J., Hanincová, K., Kurtenbach, K., Derdáková, M., and Ogden, N. H. (2008). Blackbirds and song thrushes constitute a key reservoir of Borrelia garinii, the causative agent of borreliosis in Central Europe. Appl. Environ. Microbiol. 74, 1289–1293. doi: 10.1128/AEM.01060-07
Vayssier-Taussat, M., Moutailler, S., Michelet, L., Devillers, E., Bonnet, S., and Cheval, J. (2013). Next generation sequencing uncovers unexpected bacterial pathogens in ticks in western Europe. PloS one 8:e81439. doi: 10.1371/journal.pone.0081439
Voordouw, M. J. (2015). Co-feeding transmission in Lyme disease pathogens. Parasitology 142, 290–302. doi: 10.1017/S0031182014001486
Vu Hai, V., Almeras, L., and Socolovschi, C. (2014). Monitoring human tick-borne disease risk and tick bite exposure in Europe: available tools and promising future methods. Ticks Tick Borne Dis. 5, 607–619. doi: 10.1016/j.ttbdis.2014.07.022
Keywords: urban, hard ticks, wildlife hosts, tick-borne pathogens, co-infections
Citation: Borşan S-D, Ionică AM, Galon C, Toma-Naic A, Peştean C, Sándor AD, Moutailler S and Mihalca AD (2021) High Diversity, Prevalence, and Co-infection Rates of Tick-Borne Pathogens in Ticks and Wildlife Hosts in an Urban Area in Romania. Front. Microbiol. 12:645002. doi: 10.3389/fmicb.2021.645002
Received: 22 December 2020; Accepted: 05 February 2021;
Published: 09 March 2021.
Edited by:
François JMA Meurens, UMR INRAE-Oniris 1300 Oniris - Nantes Atlantic National College of Veterinary Medicine, FranceReviewed by:
Snorre Stuen, Norwegian School of Veterinary Science, NorwayRenata Welc-Falêciak, University of Warsaw, Poland
Copyright © 2021 Borşan, Ionică, Galon, Toma-Naic, Peştean, Sándor, Moutailler and Mihalca. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Silvia-Diana Borşan, c2lsdmlhLmJvcnNhbkB1c2FtdmNsdWoucm8=