Skip to main content

Along for the ride or missing it altogether: exploring the host specificity and diversity of haemogregarines in the Canary Islands



Host-parasite relationships are expected to be strongly shaped by host specificity, a crucial factor in parasite adaptability and diversification. Because whole host communities have to be considered to assess host specificity, oceanic islands are ideal study systems given their simplified biotic assemblages. Previous studies on insular parasites suggest host range broadening during colonization. Here, we investigate the association between one parasite group (haemogregarines) and multiple sympatric hosts (of three lizard genera: Gallotia, Chalcides and Tarentola) in the Canary Islands. Given haemogregarine characteristics and insular conditions, we hypothesized low host specificity and/or occurrence of host-switching events.


A total of 825 samples were collected from the three host taxa inhabiting the seven main islands of the Canarian Archipelago, including locations where the different lizards occurred in sympatry. Blood slides were screened to assess prevalence and parasitaemia, while parasite genetic diversity and phylogenetic relationships were inferred from 18S rRNA gene sequences.


Infection levels and diversity of haplotypes varied geographically and across host groups. Infections were found in all species of Gallotia across the seven islands, in Tarentola from Tenerife, La Gomera and La Palma, and in Chalcides from Tenerife, La Gomera and El Hierro. Gallotia lizards presented the highest parasite prevalence, parasitaemia and diversity (seven haplotypes), while the other two host groups (Chalcides and Tarentola) harbored one haplotype each, with low prevalence and parasitaemia levels, and very restricted geographical ranges. Host-sharing of the same haemogregarine haplotype was only detected twice, but these rare instances likely represent occasional cross-infections.


Our results suggest that: (i) Canarian haemogregarine haplotypes are highly host-specific, which might have restricted parasite host expansion; (ii) haemogregarines most probably reached the Canary Islands in three colonization events with each host genus; and (iii) the high number of parasite haplotypes infecting Gallotia hosts and their restricted geographical distribution suggest co-diversification. These findings contrast with our expectations derived from results on other insular parasites, highlighting how host specificity depends on parasite characteristics and evolutionary history.


Understanding what determines the ability of a parasite to colonize and remain associated with a host is central to parasite evolution [1,2,3,4]. One of the most decisive factors shaping parasite distribution is host specificity [5, 6]. This trait determines the range of hosts a parasite can successfully infect and its ability to establish in new environments. The degree of host specificity is influenced by two main filters: ‘encounter’ probability with potential hosts, and ‘compatibility’, i.e. if the host provides suitable resources and the parasite can overcome host defenses [7,8,9]. Moreover, host specificity is not fixed in time, and plays a key role in parasite diversification, for example through co-speciation or expansion to novel hosts [4, 10, 11]. As such, the range of hosts a parasite can infect depends on ecological, physiological, geographical and historical factors [5, 8], alongside parasite transmission dynamics and life-cycle complexity [12,13,14,15]. To study host specificity, ideally the whole community of potential hosts of a parasite has to be considered. In this context, oceanic islands present a naturally simplified version of ecological interactions, being characterized by a depauperate and dis-harmonic fauna with high levels of endemics, which resulted from the non-random arrival of a small subset of the mainland pool followed by diversification processes [16, 17]. Colonizers face different selective pressures, which promote loss of genetic variability, increase in densities, niche broadening, among others [18]. These features are collectively referred to as the Island Syndrome. In the case of parasites, studies suggest that in insular systems they exhibit reduced diversity, higher prevalence values and an enlargement of the ecological niche, which is reflected in more frequent host-switching and more generalist parasites, in comparison to the mainland [19, 20].

Haemogregarines (Apicomplexa, Coccidia, Adeleorina) are intraerythrocytic parasites that infect a wide variety of vertebrates [21]. They have a heteroxenous life-cycle involving a hematophagous invertebrate vector (definitive host) and one vertebrate host (intermediate host) [22]. Additionally, these parasites can be transmitted via predation between vertebrate hosts [23,24,25]. Haemogregarines include three families: Haemogregarinidae, Karyolysidae and Hepatozoidae [21]. According to several phylogenetic studies, the family Haemogregarinidae (represented by the genus Haemogregarina) forms a separate clade, while the other families might not be monophyletic [26]. The genera Hemolivia and Karyolysus (both from the family Karyolysidae) clustered separately inside the genus Hepatozoon (family Hepatozoidae), leading some authors to propose a rearrangement of their taxonomy [27]. Moreover, a growing number of studies have uncovered high levels of undescribed diversity within this group ([26] and references therein), further complicating the taxonomical status of these parasites. Haemogregarines have generally been regarded as having low host specificity toward their vertebrate hosts [22, 28, 29]. However, patterns remain partially obscure as a combination of different factors complicate this assessment, including local host availability and distribution, transmission dynamics, co-evolutionary history and taxonomical uncertainties [23, 29,30,31].

The Canary Islands have been the setting of several parasitological works (e.g. [32,33,34]), some of them reporting cryptic diversity in cestodes and nematodes [35, 36]. In a recent study on intestinal nematodes of the genus Spauligodon, Jorge et al. [37] found four lineages infecting Canarian lizards, three of which present host-switching events dating to the time of parasites’ colonization. However, this widening of host breadth seems to have been temporary, as contemporary host-parasite associations show a high level of host specificity, differing from the Island Syndrome expectations. Regarding haemogregarines, these parasites have been found in the archipelago infecting lizards of the genera Gallotia [38,39,40,41,42,43] and Tarentola [44]. However, in most cases detection was based on microscopy, and only one work suggests different haemogregarine morphological forms infecting Gallotia [37]. In the current study, we use microscopy and phylogenetic approaches to analyze the diversity, distribution and host specificity of haemogregarines infecting the Canarian lizard communities. We sampled the three native lizard groups (Tarentola, Chalcides and Gallotia) across the archipelago. These different lizard genera present multiple colonization events of the Canary Islands [45,46,47] (Fig. 1) and can occur on the same localities in close proximity. In fact, there are reports of Gallotia spp. consuming members of the other lizard genera, along with conspecifics and other vertebrates [48,49,50,51]. Given these factors, the regarded low host specificity of haemogregarines, their alternative modes of transmission, and what has been described for insular parasites [19, 20], all conditions seem to favor the occurrence of host-switches and the use of multiple hosts (i.e. low specificity). So, past-punctual or still-present expansions of the host range of Canarian haemogregarines are expected. Additionally, as the Canary Islands are a biodiversity hotspot with high levels of endemism [52], we also expect to find new haemogregarine haplotypes unique to the archipelago.

Fig. 1

Study area and host lizard colonization history. a Map of the Canary Islands. Beside each main island are the island emergence age and the lizard species that inhabit it (of the genus Gallotia, Chalcides and Tarentola, endangered species are indicated by an asterisk). b Estimated time and colonization routes for the genus Gallotia, with the phylogenetic relationships between species (according to Cox et al. [52]). c Estimated time and colonization routes for the genus Chalcides [53]. d Estimated time and colonization routes for the genus Tarentola [54]


Study area and hosts

The Canary Islands are a volcanic archipelago located 100 km off the northwest African coast (Fig. 1a). They have a well-known geological genesis, with a chronological east-west island emergence ranging from 20.6 to 1.12 Ma [52, 53]. Despite their proximity and considerable age, these islands were never connected to the mainland, and native fauna and flora originated from long-distance dispersals from Africa, Europe or other archipelagos [54]. The current reptile endemic fauna includes 16 native lizard species of three genera, each belonging to a different family. The three lizard groups can be found inhabiting the same locations, though they differ in their ecology: Gallotia lacertids are diurnal and ground-dwelling; Chalcides skinks are diurnal but semi-fossorial; and Tarentola geckos are crepuscular/nocturnal and saxicolous [55]. Gallotia ancestors colonized the archipelago shortly after the first island emerged [45]. This happened in a single event from the mainland (about 17 to 20 Ma), followed by gradual expansion and diversification through all the seven main islands (Fig. 1b). While eastern and central islands are inhabited by a single species each, the western islands harbor representatives of two other Gallotia lineages: one of giant lizards (all critically endangered or extinct), and another of smaller, common species. By contrast, the two other lizard groups arrived in multiple colonization events. The Chalcides skinks reached the archipelago from the mainland twice [46]: first to the western islands and Gran Canaria (around 7 Ma); and later to the eastern group (Fuerteventura and Lanzarote, 5 Ma) (Fig. 1c); and never colonized La Palma. The Tarentola geckos present the most complex history, with three independent colonizations [47] (Fig. 1d). The first established in the eastern islands (about 5.3–11.2 Ma), the second in Tenerife, La Gomera and La Palma (dated 4.1 to 8 Ma), and the third in Gran Canaria and El Hierro (5.3–6.7 Ma). It should however be noted that successive waves of colonization and ecological replacement may have taken place [56], with recent extinctions being known [57].

Sample collection and parasite screening

Eleven of the sixteen Canarian lizard species were sampled for this study, as it was not possible to assess the endangered ones (Gallotia bravoana, Gallotia simonyi, Gallotia intermedia, Gallotia auaritae and Chalcides simonyi). Sampling was carried out under permits by the corresponding island authorities. Blood slides from 825 individuals of the three host genera (406 Gallotia, 266 Tarentola and 153 Chalcides) were collected from 46 locations across the seven main Canary Islands in several field trips between 2009 and 2014 (Fig. 2, Additional file 1: Table S1). In 27 of these locations, lizards of different genera occurred in sympatry (14 locations with Gallotia and Tarentola, two with Gallotia and Chalcides, two with Tarentola and Chalcides, and nine with the three genera). Blood was collected on slides for microscopy observation and on Whatman paper for genetic characterization, along with tail-tip host tissue samples which were used when the blood on paper was not available (for more details, see [58]). Additionally, 54 samples from the genus Psammodromus from North Africa and the Iberian Peninsula, the closest known relative of the genus Gallotia [59], were also screened.

Fig. 2

Prevalence variation across the sampling localities. For each locality, pie charts illustrate the proportion of infected (black) and uninfected (white) individuals for each genus and underneath are the corresponding number of sampled individuals. The numbers within squares represent the sampling locality codes (the same as used in the text and Additional file 1: Table S1, where further information is available). Abbreviations: G, Gallotia; T, Tarentola; C, Chalcides)

Blood slides were examined at 400× magnification, under an Olympus CX41 microscope with an in-built digital camera (SC30) (Olympus, Hamburg, Germany). All slides were initially screened for the presence of haemogregarines to assess their prevalence (i.e. the percentage of samples infected with haemogregarines). Each slide was examined until an infected cell was detected or for a minimum of 5 min (approximately ranging between 5–15 min) [41], after which the sample was considered as negative. Parasitaemia levels were then estimated for the infected samples, and calculated as the percentage of haemogregarine-infected cells per 2500 erythrocytes. For such, five random areas of each slide were photographed at 400× magnification with the cell^B software (Olympus, Münster, Germany), and 500 erythrocytes per area were counted using the Cell Counter plugin from the image processing software ImageJ v. 1.50b [60]. In 5.3% of the positive cases parasitaemia could not be estimated due to poor slide quality. Differences in overall prevalence between host genera and islands were estimated using a Generalized Linear Model (GLM) with a binomial distribution, with prevalence as the dependent variable and host genera and island as explanatory variables. Given that not all host lizard genera are present in all the islands, the interaction effect between the two explanatory variables was not considered. A second GLM analysis was performed for the Gallotia subset to confirm the differences in prevalence among islands. Analyses were performed in R version 3.3.1 [61].

Molecular identification and genetic analyses

Molecular characterization was performed on a subset of the infected samples. For Gallotia, up to five random representatives from each infected population were extracted (plus some extra individuals were sequenced to confirm parasite identification). Moreover, all infected Tarentola and Chalcides samples were also analyzed (although for one Chalcides individual, all PCR amplification attempts were unsuccessful). DNA from host blood or tissue was extracted using standard high-salt methods [62]. The PCR reactions were performed using primers specific for a 600 bp long region of the 18S rRNA gene, HepF300 and HepR900 [63]. For details on PCR conditions, see Harris et al. [64]. The 18S rRNA gene remains the most used genetic marker for the majority of haemogregarine clades [27, 65]. The pair of primers was chosen given its higher amplification success across different haemogregarine lineages infecting reptiles when compared to other primers (e.g. HEMO1 and HEMO2 [66], and EF and ER [67]), and because it provides comparable results to longer sequences in phylogenetic analysis [68]. Efforts to amplify other genetic markers available for haemogregarines [65], namely other fragments of the 18S rRNA gene [66, 67] and the ITS1 region [69], were unsuccessful. The amplified products were purified and sequenced by an external company (Beckman Coulter Genomics, UK).

Sequences were compared to the GenBank database to confirm the identity of the amplified products using the NCBI nucleotide BLAST. A total of 137 haemogregarine sequences were obtained from Gallotia hosts, 12 from Tarentola and four from Chalcides. Sequences were corrected and aligned in Geneious v5.6.7 [70], using the MAFFT algorithm [71]. Clean high-quality sequences were first used to identify the distinct haplotypes present, then the sequences of lower quality and with ambiguities were compared to these to ascertain their identity. Cases of sequences with double peaks were consistent with the previously identified haplotypes and are regarded as mixed infections. Uncorrected pairwise distances (p-distances) between haplotypes were calculated in MEGA7 [72], using a 569 bp alignment of the nine new haplotypes (Additional file 2: Table S2). New haplotype sequences are deposited in the GenBank database, under the accession numbers MG787243-MG787253.

For the phylogenetic analyses, 137 GenBank sequences of other haemogregarines were added (Additional file 3: Table S3). As putative outgroups of the Canarian haemogregarines, sequences of specimens infecting Chalcides, Tarentola and other lizard species from the African and Iberian mainland were included (previously assessed in [44, 68, 73, 74]), as well as the two new sequences from the host genus Psammodromus, the closest relative to Gallotia. In accordance with Barta et al. [29], Haemogregarina balli and Dactylosoma ranarum were used as outgroups for the overall phylogeny. In total, the final alignment matrix included 148 sequences and 583 nucleotide positions. The substitution model of evolution was chosen according to the BIC criterion selected by jModelTest 2 (model TIM1+I+G) [75]. Phylogenetic relationships were estimated using maximum likelihood (ML) and Bayesian inference (BI) methods. ML analysis was performed in PhyML 3.1 [76], with nodal support estimated using the bootstrap technique [77] with 1,000 replicates. For the BI analysis, MrBayes v.3.2.6 [78] was used, with parameters estimated as part of the analysis. The analysis was run for 1 × 107 generations, saving 1 tree each 1000 generations. The log-likelihood values of the sample point were plotted against the generation time and all the trees prior to reaching stationarity were discarded as ‘burn-in’ samples (25%). Remaining trees were combined in a 50% majority consensus tree, in which frequency of any particular clade represents the posterior probability. Additionally, phylogenetic networks for the clades containing Canarian haemogregarines were constructed using the statistical parsimony approach implemented in TCS [79] and displayed graphically using tcsBU [80].

A hierarchical analysis of molecular variance (AMOVA) was performed to test the hypotheses regarding genetic structure among haemogregarine haplotypes, using islands and host species as variables, and locations within these. Analysis was run both in the software Arlequin version 3.5.2 [81] and in the poppr package [82] in R [61], which implements the AMOVA tests from the packages ade4 [83] and pegas [84], with 1000 Monte Carlo permutations in all cases to assess statistical significance.


General assessment of haemogregarine prevalence and parasitaemia

n total, 301 specimens (36.4% of the individuals sampled) from 32 locations (69.6% of the localities) were infected with haemogregarines (see Fig. 2, Additional file 1: Table S1, Additional file 4: Table S4). The three reptile groups differed in total prevalence values (χ2 = 399.28, df = 2, P < 0.001), with the genus Gallotia being the most commonly infected (69.7%). By contrast, only 4.9% of Tarentola and 3.3% of Chalcides individuals were parasitized with haemogregarines. We also found a significant effect of the factor island (χ2 = 120.44, df = 6, P < 0.001). Differences of prevalence between islands were also observed when testing only the dataset of Gallotia spp. (χ2 = 31.44, df = 6, P < 0.001), with Gran Canaria showing the lowest infection levels (18.2%) while Tenerife had the highest (84.4%, Additional file 4: Table S4). As for parasitaemia (see Additional file 4: Table S4), Gallotia and Tarentola had similar mean values (1.07%), although the range for the former was wider (0.02–29.48 vs 0.02–4.36%, respectively). In Chalcides, parasitaemia was lower (mean = 0.05%, range = 0.02–0.16%).

Haemogregarine diversity and phylogenetic relationships

In total, nine haplotypes (A, B1, B2, C, D1, D2, E, F and T) were found infecting the Canarian lizards. The phylogenetic analysis showed congruent patterns for the ML and BI analyses. Haplotypes grouped into two well-supported clades (Fig. 3a). The first includes only Haplotype T, which clusters within a recently discovered clade of haemogregarines infecting geckos from Morocco, Algeria and Oman. Haplotype T was only found in three localities from La Palma and La Gomera, with an overall prevalence of 19.4%. Haplotype T diverged by 0.058 (uncorrected p-distance) from the other Canarian haplotypes. The other eight haplotypes (haplotypes A to F) are part of a second ‘Karyolysus’ clade (sensu Karadjian et al. [27]) and are further distributed in two subclades (Fig. 3a). Haplotype A is placed in a subclade together with haemogregarines infecting lizards from the Iberian Peninsula and North-western Africa, while the second subclade comprises the remaining haplotypes (p-distance = 0.0018 between haplotype A and haplotypes B1 to F). This second subclade comprises haemogregarines from a wide range of reptile hosts, including lacertid lizards, snakes and skinks from Europe, Morocco and Oman. Unfortunately, the relationships within this subclade are not highly supported (Fig. 3a). The highest p-distance within this group was 0.011 between B2 and D1, and the lowest was 0.002 (1 mutation) between B1 and B2, and also D1 and D2 (see Additional file 2: Table S2). Given the sequence similarity and association between them, we treat these sister haplotypes as a single unit each and hereafter they will be called haplogroups B and D when referred to collectively. Regarding the parasites from Gallotia’s sister taxa, two haemogregarine haplotypes were discovered in Psammodromus algirus from Morocco (identified as PSA 1 and 2 in Fig. 3a and Additional file 5: Figure S1a). One grouped within the same subclade as haplotype A, while the other one clustered with the remaining haplotypes found in Gallotia spp., although with low support. Haplotype networks showing the relationships among haplotypes are detailed in Additional file 5: Figure S1.

Fig. 3

Phylogenetic analysis (a) and distribution (b) of the haemogregarine haplotypes infecting the three lizards studied. Evolutionary relationships were estimated using the 18S rRNA gene. Presented is the tree derived from the BI analysis. Bayesian posterior probabilities are given above nodes and bootstrap values for ML below them (only values above 75 for the former and 70 for the latter are shown). For better reading, some branches of the phylogenetic tree have been collapsed or shortened (Additional file 3: Table S3 contains all the GenBank codes and extra information on the sequences used). Sequence KU680457, identified with an asterisk in the tree, corresponds to the infection in T. angustimentalis previously reported by Tomé et al. [46]. The nine Canarian haplotypes are differentiated by colors in the tree and maps (here the striped pattern identifies mixed infections). Pie charts represent the proportion of each haplotype per location, and beside them is the number of sequenced infections

Geographical distribution of haplotypes and host specificity

All recovered haplotypes were restricted to a specific lizard host genus. Seven of them (haplotypes A, B1, B2, C, D1, D2 and E) infected Gallotia lizards, while haplotype F was exclusive to Chalcides and haplotype T to Tarentola (Fig. 3b, Additional file 1: Table S1). The exception was haplotype B1 which, although mostly found in Gallotia, was also detected in two instances infecting Chalcides (one individual, locality 2) and Tarentola (one individual, locality 12). Haplotype F was found in a single site in Tenerife (location 2), infecting three Chalcides individuals. Haplotype T was found in Tarentola from two islands (La Palma and La Gomera). In Gallotia, haplotypes A and E were found in the two easternmost islands, haplotype C in the central island of Gran Canaria, and haplotypes B1, B2, D1 and D2 in the western islands. AMOVA tests were performed on the set of haplotypes with a plausible common origin: B1, B2, C, D1, D2 and E. Results were the same for the three different methodologies and showed a high correlation of genetic distances to both islands (69.32% of molecular variance) and Gallotia host species (72.82%) (Table 1).

Table 1 AMOVA tests of genetic structure for haplotypes B1 to E

It should also be noted that Gallotia individuals presented mixed infections (in 13 of the 32 localities with infected Gallotia, see Additional file 1: Table S1 and Fig. 3b). From the island perspective, the number of haplotypes present varied (Fig. 3b). Gran Canaria and La Gomera had only one haplotype each (C and D1, respectively). Two haplotypes were present in Fuerteventura, Lanzarote (A and E) and El Hierro (B and D1). Four haplotypes were found in Tenerife: B1, B2, D1, together with F infecting Chalcides. Finally, La Palma harbored the highest number of haplotypes, with haplotypes B1, B2, D1 and D2 infecting Gallotia and T infecting Tarentola.

Differences in prevalence and parasitaemia

Because only a subset of parasite specimens was sequenced for Gallotia, we could not thoroughly assess the prevalence for the different haplotypes. However, as sequenced samples were randomly selected, the relative frequency between sympatric Gallotia haplotypes gives a reasonable approximation of their prevalence (Additional file 1: Table S1 and Fig. 3b). Within the eastern islands (Lanzarote and Fuerteventura) haplotype A was more frequent than haplotype E. In the Central island of Gran Canaria, only haplotype C was present, in a very low frequency (18.2%). In both Tenerife and La Palma, haplotype B was more frequent than haplotype D, while in El Hierro the reverse occurred. Lastly, haplotype D was the only one present in La Gomera, in a lower frequency (52.2%) than the overall mean prevalence in the genus Gallotia (69.7%). Regarding parasitaemia, haplotype A had the highest mean (1.15%) and E the lowest (0.13%), though it should be noted that sample size was quite disparate between haplotypes (Additional file 4: Table S4). Haplotype B was the only one for which cross-infections were detected, having a mean parasitaemia of 1.38% in the genus Gallotia, while for the only two cases of Chalcides and Tarentola infected with this haplotype, values were 0.02% and 0.04%, respectively.


Haemogregarine distribution and transmission

This is the first comprehensive overview on the haemogregarines infecting the endemic lizards of the Canary Islands. This archipelago is considered a biodiversity hotspot due to its high levels of endemism [52] and our results further support this, as we identified a high diversity of haemogregarine haplotypes previously unknown, none matching sequences from other geographical areas. Moreover, and although sampling included host sympatric localities, haplotypes were mostly restricted to a particular lizard genus, suggesting that these parasites are highly specific to their vertebrate hosts. Seven of the nine haplotypes were found infecting a single host species (haplotype C in Gallotia stehlini, E in Gallotia atlantica, and F in Chalcides viridanus) or species with a recent common ancestry (haplotype D in Gallotia galloti and Gallotia caesaris, and haplotype T in Tarentola gomerensis and Tarentola delalandii). In contrast, haplotype B1, which mostly infected western Gallotia, was also found in Tarentola and Chalcides. Similarly, haplotype A, in this study frequently found in G. atlantica (but not geckos), was previously reported infecting one Tarentola angustimentalis from Fuerteventura [44]. Given the low prevalence of these two haplotypes in Tarentola and Chalcides hosts, we consider these cases as occasional cross-infections (although these lizards might still be competent hosts, as peripheral blood stages were observed). This also highlights that, to correctly assess host range and specificity of a parasite, studies should extend to the community of potential sympatric hosts [9, 15, 85]. Regarding mixed infections, interestingly the few cases confirmed were only detected in Gallotia hosts (Fig. 3 and Additional file 1: Table S1). These cases involved only Gallotia specific haplotypes that co-occurred in the same island, although some of them were distantly related (A and E).

It has been suggested that haemogregarines are more host-specific to their invertebrate hosts than to their vertebrate hosts [22, 29]. This is relevant as the level of specialization of the vector is expected to condition the opportunity of a parasite to encounter a given vertebrate host [86, 87]. Mites of the genus Ophionyssus are the main suspected vector and definitive host of the “Karyolysus” group ([21, 88, 89], and references within), although this has not been confirmed as general. In the Canary Islands, three Ophionyssus species have been described parasitizing Gallotia spp. [90, 91], and their role as vector for haemogregarines has already been suggested [36]. Their distribution seems congruent with the observed patterns for haemogregarine haplotypes from Gallotia, with a species reported for eastern, Gran Canaria, and western islands. Additionally, Bannert et al. [92] examined the three lizard groups in Tenerife, and found neither geckos nor skinks infected by Ophionyssus spp. (while these were highly prevalent in Gallotia lizards of the same location). On the other hand, Tarentola from the Canary Islands are reported to be infected by mites of the genus Geckobia [93], a specialist mite of geckos [94], although again, their role as vectors of haemogregarines has not been confirmed. Given the high dependency of mites on the lizard host for their life-cycle [92, 94], it seems reasonable to assume that they colonized the archipelago while carried by their lizard hosts, possibly along with the haemogregarines. Surely, the colonization and expansion success of their invertebrate hosts must have also conditioned the distribution of the haemogregarines. Therefore, future work should identify and assess the role of these vectors in the colonization success and adaptation to new environments, as well as in the host specificity of parasites with heteroxenous life-cycles. On the other hand, Gallotia spp. are known to consume other vertebrates, including conspecifics [48,49,50,51], which might provide an additional infection pathway [95]. Indeed, transmission of haemogregarines via predation has been shown in other hosts [24, 25]. Also worth noting, in Gallotia lizards, Sarcocytis gallotiae has adapted to insularity by using the same host species as both definitive and intermediate host by exploiting its cannibalistic behavior [96, 97].

Phylogenetic relationships and biogeography

Canarian lizards share common ancestors with lizards from North Africa and the Iberian Peninsula [45,46,47], so it is likely that their parasites share this origin and our results support this. In fact, the haplotypes detected here clustered in haemogregarine clades found in European and African reptiles (Fig. 3a). Particularly, some of them were related to haemogregarines from the hosts’ closest relatives, such as the genus Psammodromus and other Tarentola species from the continent. The pattern is, however, not so clear for the haemogregarines infecting skinks (haplotype F). Nonetheless, further conclusions are difficult due to the low support of the phylogenetic relationships within this group. Additionally, our results provide preliminary insights on how these haemogregarines colonized the Canary Archipelago. We note, however, that some inferences remain preliminary as additional data are necessary, namely the identity of haemogregarine definitive hosts, the haplotypes infecting the endangered Canarian lizard species, and phylogenetic information from additional faster-evolving molecular markers.

Our results suggest that haemogregarines arrived to the Canary Islands along with the ancestors of the three lizard genera. The haemogregarine with the clearest colonization is haplotype T. The ancestors of one lineage of Tarentola arrived from the African mainland to the Canary Islands, reaching first Tenerife and then dispersing to La Palma and La Gomera (Fig. 1d) [47]. The haemogregarine ancestors of haplotype T might have followed the same path, and its absence in Tenerife could be explained by extinction events or low detectability as a consequence of low prevalence values. Regarding haemogregarines infecting Gallotia, haplotypes A and E were found in the two easternmost islands. These were the first islands of the archipelago colonized by Gallotia’s ancestors (Fig. 1b) [45]. Although haplotypes A and E belong to two distinct clades, they could have “shared the boat”, colonizing these islands in the same event aboard their Gallotia ancestor hosts (and/or vector hosts). The absence of haplotype A (or any close relatives) on the remaining islands suggests that it “missed the boat” when Gallotia further expanded through the archipelago (or went extinct during these new colonizations). Contrarily, haplotype E has closely related haplotypes (B, C and D) in the other islands, supporting further colonizations of this haplotype ancestor to other islands together with its Gallotia host. Unfortunately, phylogenetic relationships within this clade are not fully resolved, preventing us from drawing more conclusions. The results from the AMOVA show a high level of genetic structure concordant with both Gallotia host species and geography, supporting the hypothesis that these haplotypes followed the colonization paths of their Gallotia hosts. Finally, regarding haplotype F, which was detected in skinks from Tenerife, the simplest scenario would be that it arrived there along with the Chalcides ancestors that colonized the western islands (Fig. 1c) [46], and then either did not further expand or later retracted (as it was only found in one location). Alternatively, host-switching events might have taken place, from the skinks to the Gallotia or vice-versa, although with the current data it is not possible to infer if this is the case. The presence of paralogs is a characteristic of the 18S rRNA gene that might cloud the interpretation of phylogenetic results, and it has been confirmed for some apicomplexans (e.g. [98,99,100], suspected for haemogregarines in [101]). However, the haemogregarine haplotypes detected here present a non-random distribution (as shown by the aforementioned geographical structuration) and, although we did find a few cases with more than one haplotypic form within the same host individual, these were instead more suggestive of mixed infections (as they were rare and correspond to the haplotypes retrieved in single infections). Given these reasons, we do not suspect paralogy is an issue in the current dataset, though this can only be ruled out completely when other sequence data are available for these parasites, for example from mitochondrial genes.

Factors shaping haemogregarine diversity and host specificity

Parasite prevalence, parasitaemia and genetic diversity were higher in Gallotia hosts than in Tarentola and Chalcides (values in Gallotia spp. were consistent with previous studies [39,40,41,42,43]). Interestingly, a similar trend was also observed in nematodes of the genus Spauligodon infecting the same hosts [37]. Such congruence for two radically different parasites might suggest that host ecological characteristics play a role in parasite abundance, diversity distribution and transmission dynamics [31, 102], alongside parasite adaptation or host immunology. The three lizard genera can be found inhabiting the same localities but differ in daily activity and microhabitat use [55]. This segregation and differences in host densities might lead to heterogeneity in the encounter opportunity, affecting the chances of parasite transmission and expansion [103,104,105]. For example, Martín et al. [106] have conjectured that a fossorial lifestyle (such as in Chalcides [46]) could lead to a lower exposure to vectors, and consequently to haemogregarines. The lizard groups also differ in their body temperatures, with Gallotia lizards thermoregulating at higher temperatures than Tarentola geckos [107, 108]. Host internal temperature might be also relevant since the infection levels of haemogregarines seem to be affected by this factor [109]. Additionally, distinct lizard populations are exposed to different environmental stressors (either natural climatic or human related), which might locally influence prevalence and parasitaemia, e.g. [110]. Nonetheless, nematodes and haemogregarines did not show the same patterns of diversity distribution and potential colonization pathways. While both Spauligodon nematodes and haemogregarines show high levels of diversity and host specificity, past host-switch events were detected for the former, but not for haemogregarines. The two parasites differ in many characteristics; most notably Spauligodon nematodes have a direct life-cycle while haemogregarines have a multiple-host one. This might have relevant repercussions, as a heteroxenous life-cycle complicates the maintenance of the transmission dynamics, and increases the constraints for a successful establishment ([95, 111] but see [15]).

There is a marked contrast in the number of haemogregarine haplotypes between eastern (Gran Canaria, Fuerteventura and Lanzarote) and western islands (El Hierro, La Palma, La Gomera and Tenerife). This is partly due to the presence of the haplotypes infecting skinks and geckos but, considering only the haplotypes infecting Gallotia lizards, this pattern is still maintained. There is a lower Gallotia species richness in the East (only two species), compared to the radiation that took place in the West (in total seven species of two distinct lineages [45]). Here, co-speciation with the host [112] and/or isolation in the different western islands might have fomented parasite diversification (as also hypothesized for nematodes by Jorge et al. [113]). On the other hand, an impoverished biodiversity in the eastern islands relative to the rest of the archipelago is a general pattern [114, 115], which has also been observed in the helminth fauna infecting G. atlantica [116], and relates to the nonlinear relation between the age of an oceanic island and its biodiversity [117, 118]. Also, Fuerteventura and Lanzarote were connected several times during the Pleistocene [119, 120], which would dilute the effect of their previous isolation. Within-island geological events could explain the heterogeneity in prevalence and in the relative abundance of the parasite haplotypes. For example, volcanic eruptions and landslides have led to repetitive isolations and secondary contacts of host populations [121,122,123]. Another possible phenomenon is parasite spillover from now-extinct or introduced hosts [124], as insular systems are particularly prone to extinction-recolonization events and to introductions [17, 125]. However, we found no evidence of this in the current results, given the congruence between haemogregarine and host distributions, and the fact that we did not identify parasite haplotype exchange between islands with known lizard translocations or introductions [126].

Our results showed that haemogregarines infecting the lizards of the Canary Island did not have a low host specificity or represent switching events between host groups. This contrasts with our initial predictions, as it does not seem to support the hypothesized parasite niche broadening expected in insular systems [19]. Although such might have still occurred, it is not evident based on the currently available data. While the mechanism behind the high host specificity pattern of these haemogregarines is not elucidated, we list some of the factors that might be involved. Segregation between lizard groups due to their ecology and differences in transmission (such as their vectors) seem major suspects, but molecular factors may (also or instead) play a role, including host immunity response or the ability of the parasite to exploit host resources [127]. It should however be noted that occasional infections show that the different lizard groups might be competent hosts, at least for some of the Canarian haemogregarine haplotypes. Thus, even if these haemogregarines can have a lower specificity regarding their potential host range, their realized niche seems to actually be restrictive in the hosts used. This is relevant as haemogregarines are generally considered as having low host specificity for their vertebrate hosts [22, 29, 128]. This has been supported by experimental studies (e.g. [69], references within [129]), but might not reflect the reality in the wild [129]. Additionally, it should be noted that host specificity and host-switch potential might depend on parasite lineage [31, 37], so the patterns might differ between haemogregarine groups and explain incongruences across studies.


This is the first comprehensive study on haemogregarines of the Canary Islands, specifically those infecting lizards. Our results seem to suggest a high level of host specificity for haemogregarine parasites infecting lizards of the Canary Islands, with no evidence of host-switching events. Moreover, for haemogregarines infecting Gallotia hosts, co-diversification apparently took place, as a high diversity of geographically-structured parasite haplotypes is present. The possible factors shaping the host specificity and biogeography of these parasites include differences in host ecology, vector identity and parasite transmission. The observed patterns contrast with our initial expectations, and similarly to Spauligodon nematodes, no long-lasting host range broadening was observed. This highlights how host specificity depends on parasite characteristics and evolutionary history, and questions what is generally assumed for insular parasites and haemogregarines. Lastly, we note that the nine parasite haplotypes found in the Canary Islands are distinct from those detected until this date on the mainland. Although more markers and samples are needed to fully confirm our results, the evolutionary history inferred with the 18S rRNA gene suggests that the Canarian haemogregarines underwent a strong diversification process, representing one more example of the high diversity present in this archipelago.



Analysis of molecular variance


Bayesian inference


Bayesian information criterion


Deoxyribonucleic acid


Generalized linear model


Internal transcribed spacer 1


Maximum likelihood


Polymerase chain reaction


Ribosomal ribonucleic acid


  1. 1.

    Dunn RR, Harris NC, Colwell RK, Koh LP, Sodhi NS. The sixth mass coextinction: are most endangered species parasites and mutualists? Proc Biol Sci. 2009;276:3037–45.

    PubMed Central  Article  PubMed  Google Scholar 

  2. 2.

    Hatcher MJ, Dick JTA, Dunn AM. Diverse effects of parasites in ecosystems: Linking interdependent processes. Front Ecol Environ. 2012;10:186–94.

    Article  Google Scholar 

  3. 3.

    Gonzalez-Quevedo C, Davies RG, Richardson DS. Predictors of malaria infection in a wild bird population: Landscape-level analyses reveal climatic and anthropogenic factors. J Anim Ecol. 2014;83:1091–102.

    Article  PubMed  Google Scholar 

  4. 4.

    Hoberg EP, Brooks DR. Evolution in action: climate change, biodiversity dynamics and emerging infectious disease. Philos Trans R Soc B Biol Sci. 2015;370:20130553.

    Article  Google Scholar 

  5. 5.

    Poulin R. Are there general laws in parasite ecology? Parasitology. 2007;134:763–76.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Falk BG, Perkins SL. Host specificity shapes population structure of pinworm parasites in Caribbean reptiles. Mol Ecol. 2013;22:4576–90.

    Article  PubMed  Google Scholar 

  7. 7.

    Remold S. Understanding specialism when the Jack of all trades can be the master of all. Proc Biol Sci. 2012;279:4861–9.

    PubMed Central  Article  PubMed  Google Scholar 

  8. 8.

    Combes C. Parasitism: the ecology and evolution of intimate interactions. Chicago: University of Chicago Press; 2001.

    Google Scholar 

  9. 9.

    Pinheiro RBP, Félix GMF, Chaves AV, Lacorte GA, Santos FR, Braga ÉM, et al. Trade-offs and resource breadth processes as drivers of performance and specificity in a host-parasite system: a new integrative hypothesis. Int J Parasitol. 2016;46:115–21.

  10. 10.

    Araujo SBL, Braga MP, Brooks DR, Agosta SJ, Hoberg EP, Von Hartenthal FW, et al. Understanding host-switching by ecological fitting. PLoS One. 2015;10:e0139225.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  11. 11.

    Agosta SJ, Janz N, Brooks DR. How specialists can be generalists: resolving the “parasite paradox” and implications for emerging infectious disease. Zool (Curitiba, Impresso). 2010;27:151–62.

    Article  Google Scholar 

  12. 12.

    Poulin R, Keeney DB. Host specificity under molecular and experimental scrutiny. Trends Parasitol. 2008;24:24–8.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Dick CW, Patterson BD. Against all odds: Explaining high host specificity in dispersal-prone parasites. Int J Parasitol. 2007;37:871–6.

    Article  PubMed  Google Scholar 

  14. 14.

    Harbison CW, Clayton DH. Community interactions govern host-switching with implications for host-parasite coevolutionary history. Proc Natl Acad Sci USA. 2011;108:9525–9.

  15. 15.

    Malcicka M, Agosta SJ, Harvey JA. Multi level ecological fitting: Indirect life cycles are not a barrier to host switching and invasion. Glob Chang Biol. 2015;21:3210–8.

    Article  PubMed  Google Scholar 

  16. 16.

    MacArthur RH, Wilson EO. The theory of island biogeography. Princeton University Press: Princeton; 1967.

    Google Scholar 

  17. 17.

    Whittaker RJ, Fernández-Palacios JM. Island biogeography: ecology, evolution, and conservation. Oxford: Oxford University Press; 2007.

    Google Scholar 

  18. 18.

    Adler GH, Levins R. The island syndrome in rodent populations. Q Rev Biol. 1994;69:473–90.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Nieberding C, Morand S, Libois R, Michaux JR. Parasites and the island syndrome: The colonization of the western Mediterranean islands by Heligmosomoides polygyrus (Dujardin, 1845). J Biogeogr. 2006;33:1212–22.

    Article  Google Scholar 

  20. 20.

    Pérez-Rodríguez A, Ramírez Á, Richardson DS, Pérez-Tris J. Evolution of parasite island syndromes without long-term host population isolation: parasite dynamics in Macaronesian blackcaps Sylvia atricapilla. Glob Ecol Biogeogr. 2013;22:1272–81.

  21. 21.

    Telford SR Jr. Hemoparasites of the Reptilia: Color Atlas and Text. Boca Raton: CRC Press; 2009.

    Google Scholar 

  22. 22.

    Smith TG. The genus Hepatozoon (Apicomplexa: Adeleina). J Parasitol. 1996:565–85.

  23. 23.

    Tomé B, Maia JPMC, Harris DJ. Hepatozoon infection prevalence in four snake genera: influence of diet, prey parasitemia levels, or parasite type? J Parasitol. 2012;98:913–7.

    Article  PubMed  Google Scholar 

  24. 24.

    Sloboda M, Kamler M, Bulantová J, Votýpka J, Modrý D. Rodents as intermediate hosts of Hepatozoon ayorgbor (Apicomplexa: Adeleina: Hepatozoidae) from the African ball python, Python regius? Folia Parasitol. 2008;55:13–6.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Johnson EM, Panciera RJ, Allen KE, Sheets ME, Beal JD, Ewing SA, et al. Alternate pathway of infection with Hepatozoon americanum and the epidemiologic importance of predation. J Vet Intern Med. 2009;23:1315–8.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Maia JP, Carranza S, Harris DJ. Comments on the systematic revision of adeleid haemogregarines: is more data needed? J Parasitol. 2016;102:549–52.

    Article  PubMed  Google Scholar 

  27. 27.

    Karadjian G, Chavatte J-M, Landau I. Systematic revision of the adeleid haemogregarines, with creation of Bartazoon n. g., reassignment of Hepatozoon argantis Garnham, 1954 to Hemolivia, and molecular data on Hemolivia stellata. Parasite. 2015;22:31.

    PubMed Central  Article  PubMed  Google Scholar 

  28. 28.

    Hoff G. Diseases of amphibians and reptiles. Springer Science & Business Media; 2012.

  29. 29.

    Barta JR, Ogedengbe JD, Martin DS, Smith TG. Phylogenetic position of the adeleorinid coccidia (Myzozoa, Apicomplexa, Coccidia, Eucoccidiorida, Adeleorina) inferred using 18S rDNA sequences. J Eukaryot Microbiol. 2012;59:171–80.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Maia JP, Álvares F, Boratyński Z, Brito JC, Leite JV, Harris DJ. Molecular assessment of Hepatozoon (Apicomplexa: Adeleorina) infections in wild canids and rodents from North Africa, with implications for transmission dynamics across taxonomic groups. J Wildl Dis. 2014;50:837–48.

    Article  PubMed  Google Scholar 

  31. 31.

    Maia JP, Harris DJ, Carranza S, Goméz-Díaz E. Assessing the diversity, host-specificity and infection patterns of apicomplexan parasites in reptiles from Oman, Arabia. Parasitology. 2016;143:1730–47.

    Article  PubMed  Google Scholar 

  32. 32.

    Roca V, Carretero M, Llorente G, Montori A, Martin J. Helminth communities of two lizard populations (Lacertidae) from Canary Islands (Spain): Host diet-parasite relationships. Amphib Reptil. 2005;26:535–42.

    Article  Google Scholar 

  33. 33.

    Feliu C, López M, Gómez MS, Torres J, Sánchez S, Miquel J, et al. Parasite fauna of rodents (Murinae) from El Hierro (Canary Islands, Spain): a multidisciplinary approach. Acta Parasitol. 2012;57:171–8.

    Article  PubMed  Google Scholar 

  34. 34.

    Foronda PR, Figueruelo EO, Ortego AR, Abreu NA, Casanova JC. Parasites (viruses, coccidia and helminths) of the wild rabbit (Oryctolagus cuniculus) introduced to Canary Islands from Iberian Peninsula. Acta Parasitol. 2005;50:80–4.

    Google Scholar 

  35. 35.

    Foronda P, López-González M, Hernández M, Haukisalmi V, Feliu C. Distribution and genetic variation of hymenolepidid cestodes in murid rodents on the Canary Islands (Spain). Parasit Vectors. 2011;4:185.

    PubMed Central  Article  PubMed  Google Scholar 

  36. 36.

    Jorge F, Perera A, Carretero MA, Harris DJ, Roca V. Cryptic species unveiled: The case of the nematode Spauligodon atlanticus. J Zool Syst Evol Res. 2013;51:187–202.

    Article  Google Scholar 

  37. 37.

    Jorge F, Perera A, Poulin R, Roca V, Carretero MA. Getting there and around: Host range oscillations during colonization of the Canary Islands by the parasitic nematode Spauligodon. Mol Ecol. 2018;00:1–17.

  38. 38.

    Bannert B, Lux E, Sedlaczek J. Studies on endo- and ectoparasites of Canarian Lizards. Sci Herpetol. 1995;1995:293–6.

    Google Scholar 

  39. 39.

    Oppliger A, Vernet R, Baez M. Parasite local maladaptation in the Canarian lizard Gallotia galloti (Reptilia: Lacertidae) parasitized by haemogregarian blood parasite. J Evol Biol. 1999;12(5):951.

    Article  Google Scholar 

  40. 40.

    Martínez-Silvestre A, Mateo JA, Silveira L, Bannert B. Presencia de protozoos intraeritrocitarios en el lagarto gigante de La Gomera. Bol Asoc Herpetol Esp. 2001;12:90–2.

    Google Scholar 

  41. 41.

    García-Ramírez A, Delgado-García JD, Foronda-Rodríguez P, Abreu-Acosta N. Haematozoans, mites and body condition in the oceanic island lizard Gallotia atlantica (Peters and Doria, 1882) (Reptilia: Lacertidae). J Nat Hist. 2005;39:1299–305.

    Article  Google Scholar 

  42. 42.

    Foronda P, Santana-Morales MA, Orós J, Abreu-Acosta N, Ortega-Rivas A, Lorenzo-Morales J, et al. Clinical efficacy of antiparasite treatments against intestinal helminths and haematic protozoa in Gallotia caesaris (lizards). Exp Parasitol. 2007;116:361–5.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Megía-Palma R, Martínez J, Merino S. A structural colour ornament correlates positively with parasite load and body condition in an insular lizard species. Sci Nat. 2016;103:52.

    Article  CAS  Google Scholar 

  44. 44.

    Tomé B, Rato C, Perera A, Harris DJ. High diversity of Hepatozoon spp. in geckos of the genus Tarentola. J Parasitol. 2016;102:476–80.

    Article  PubMed  Google Scholar 

  45. 45.

    Cox SC, Carranza S, Brown RP. Divergence times and colonization of the Canary Islands by Gallotia lizards. Mol. Phylogenet Evol. 2010;56:747–57.

    Article  PubMed  Google Scholar 

  46. 46.

    Carranza S, Arnold EN, Geniez P, Roca J, Mateo JA. Radiation, multiple dispersal and parallelism in the skinks, Chalcides and Sphenops (Squamata: Scincidae), with comments on Scincus and Scincopus and the age of the Sahara Desert. Mol Phylogenet Evol. 2008;46:1071–1094.

  47. 47.

    Carranza S, Arnold EN, Mateo JA, Geniez P. Relationships and evolution of the North African geckos, Geckonia and Tarentola (Reptilia: Gekkonidae), based on mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2002;23:244–56.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Rodríguez A, Nogales M, Rumeu B, Rodríguez B. Temporal and spatial variation in the diet of the endemic lizard Gallotia galloti in an insular Mediterranean scrubland. J Herpetol. 2008;42:213–22.

    Article  Google Scholar 

  49. 49.

    Martínez-Silvestre A, Mateo JA, Jones M, Gallego S, Soler J. Comportamiento alimentario inusual en Gallotia atlantica. Bol Asoc Herpetol Esp. 2012;23:23–7.

    Google Scholar 

  50. 50.

    Fariña B, Martín A. Ornitofagia de Gallotia intermedia en Guaza (Tenerife, Islas Canarias). Zoology. 2013;105:203–14.

    Google Scholar 

  51. 51.

    Mateo JA, Pleguezuelos JM. Cannibalism of an endemic island lizard (genus Gallotia). Zool Anz. 2015;259:131–4.

    Article  Google Scholar 

  52. 52.

    Juan C, Emerson BC, Oromí P, Hewitt GM. Colonization and diversification: Towards a phylogeographic synthesis for the Canary Islands. Trends Ecol Evol. 2000;15:104–9.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Ancochea E, Hernán F, Huertas MJ, Brändle JL, Herrera R. A new chronostratigraphical and evolutionary model for La Gomera: implications for the overall evolution of the Canarian Archipelago. J Volcanol Geotherm Res. 2006;157:271–93.

  54. 54.

    Illera JC, Spurgin LG, Rodriguez-Exposito E, Nogales M, Rando JC. What are we learning about speciation and extinction from the Canary Islands? Ardeola. 2016;63:5–23.

    Article  Google Scholar 

  55. 55.

    Mateo JA, Afonso OM, Geniez P. Los reptiles de Canarias, una nueva sinopsis puesta al día. Bol Asoc Herpetol Esp. 2007;12:2–10.

    Google Scholar 

  56. 56.

    López-Jurado LF. Synopsis of the Canarian herpetofauna. Rev Esp Herp. 1991;6:107–18.

    Google Scholar 

  57. 57.

    Maca-Meyer N, Carranza S, Rando JC, Arnold EN, Cabrera VM. Status and relationships of the extinct giant Canary Island lizard Gallotia goliath (Reptilia: Lacertidae), assessed using ancient mtDNA from its mummified remains. Biol J Linn Soc. 2003;80:659–70.

    Article  Google Scholar 

  58. 58.

    Maia JP, Harris DJ, Carranza S, Gómez-Díaz E. A comparison of multiple methods for estimating parasitemia of hemogregarine hemoparasites (Apicomplexa: Adeleorina) and its application for studying infection in natural populations. PLoS One. 2014;9:e95010.

    PubMed Central  Article  PubMed  Google Scholar 

  59. 59.

    Mendes J, Harris DJ, Carranza S, Salvi D. Evaluating the phylogenetic signal limit from mitogenomes, slow evolving nuclear genes, and the concatenation approach. New insights into the Lacertini radiation using fast evolving nuclear genes and species trees. Mol Phylogenet Evol. 2016;100:254–67.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Abràmoff MD, Magalhães PJ, Ram SJ. Image processing with imageJ. Biophotonics Int. 2004;11:36–41.

    Google Scholar 

  61. 61.

    R Core Team. R: A language and environment for statistical computing. Version 3.3.1. 2016. Accessed 1 Sept 2016.

  62. 62.

    Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: A laboratory manual. New York: Cold Spring Harbor Laboratory Press; 1989.

    Google Scholar 

  63. 63.

    Ujvari B, Madsen T, Olsson M. High prevalence of Hepatozoon spp. (Apicomplexa, Hepatozoidae) infection in water pythons (Liasis fuscus) from tropical Australia. J. Parasitol. 2004;90:670–2.

    Article  PubMed  Google Scholar 

  64. 64.

    Harris DJ, Maia JPMC, Perera A. Molecular characterization of Hepatozoon species in reptiles from the Seychelles. J Parasitol. 2011;97:106–10.

    Article  PubMed  Google Scholar 

  65. 65.

    Modrý D, Beck R, Hrazdilová K, Baneth G. A review of methods for detection of Hepatozoon infection in carnivores and arthropod vectors. Vector Borne Zoonotic Dis. 2017;17:66–72.

    Article  PubMed  Google Scholar 

  66. 66.

    Perkins SL, Keller AK. Phylogeny of nuclear small subunit rRNA genes of hemogregarines amplified with specific primers. J Parasitol. 2001;87:870–6.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Kvičerová J, Pakandl M, Hypša V. Phylogenetic relationships among Eimeria spp. (Apicomplexa, Eimeriidae) infecting rabbits: Evolutionary significance of biological and morphological features. Parasitology. 2008;135:443–52.

    Article  CAS  PubMed  Google Scholar 

  68. 68.

    Maia JPMC, Perera A, Harris DJ. Molecular survey and microscopic examination of Hepatozoon Miller, 1908 (Apicomplexa: Adeleorina) in lacertid lizards from the western Mediterranean. Folia Parasitol. 2012;59:241–8.

    Article  PubMed  Google Scholar 

  69. 69.

    Kim B, Smith TG, Desser SS. The life history and host specificity of Hepatozoon clamatae (Apicomplexa: Adeleorina) and ITS-1 nucleotide sequence variation of Hepatozoon species of frogs and mosquitoes from Ontario. J Parasitol. 1998;84:789–97.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

  71. 71.

    Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

  72. 72.

    Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Maia JPMC, Harris DJ, Perera A. Molecular survey of Hepatozoon species in lizards from North Africa. J Parasitol. 2011;97:513–7.

    Article  PubMed  Google Scholar 

  74. 74.

    Tomé B, Maia J, Salvi D, Brito J. Patterns of genetic diversity in Hepatozoon spp. infecting snakes from North Africa and the Mediterranean Basin. Syst Parasitol. 2014;87:249–58.

    Article  PubMed  Google Scholar 

  75. 75.

    Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  76. 76.

    Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.

    Article  PubMed  Google Scholar 

  77. 77.

    Felsenstein J. Use of Bootstrap in phylogenetic calculations. Evolution. 1985;39:783–91.

    Article  PubMed  Google Scholar 

  78. 78.

    Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: Efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.

    PubMed Central  Article  PubMed  Google Scholar 

  79. 79.

    Clement M, Posada D, Crandall KA. TCS: A computer program to estimate gene genealogies. Mol Ecol. 2000;9:1657–9.

    CAS  Article  PubMed  Google Scholar 

  80. 80.

    Múrias Dos Santos A, Cabezas MP, Tavares AI, Xavier R, Branco M. TcsBU: A tool to extend TCS network layout and visualization. Bioinformatics. 2015;32:627–8.

    Article  CAS  PubMed  Google Scholar 

  81. 81.

    Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.

  82. 82.

    Kamvar ZN, Tabima JF, Grünwald NJ. Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ. 2014;2:e281.

    PubMed Central  Article  PubMed  Google Scholar 

  83. 83.

    Dray S, Dufour AB. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22:1–20.

    Article  Google Scholar 

  84. 84.

    Paradis E. pegas: an R package for population genetics with an integrated-modular approach. Bioinformatics. 2010;26:419–20.

    CAS  Article  PubMed  Google Scholar 

  85. 85.

    Norton DA, Carpenter MA. Mistletoes as parasites: host specificity and speciation. Trends Ecol Evol. 1998;13:101–5.

  86. 86.

    Whiteman NK, Matson KD, Bollmer JL, Parker PG. Disease ecology in the Galápagos Hawk (Buteo galapagoensis): host genetic diversity, parasite load and natural antibodies. Proc Biol Sci. 2006;273:797–804.

    Article  PubMed  Google Scholar 

  87. 87.

    Medeiros MCI, Hamer GL, Ricklefs RE. Host compatibility rather than vector-host-encounter rate determines the host range of avian Plasmodium parasites. Proc R Soc B. 2013;280:20122947.

    PubMed Central  Article  PubMed  Google Scholar 

  88. 88.

    Svahn K. Blood parasites of the genus Karyolysus (Coccidia, Adeleidae) in Scandinavian lizards. Description and life cycle. Nor J Zool. Universitetsforlaget. 1975;23:277–95.

    Google Scholar 

  89. 89.

    Haklová-Kočíková B, Hižňanová A, Majláth I, Račka K, Harris DJ, Földvári G, et al. Morphological and molecular characterization of Karyolysus - a neglected but common parasite infecting some European lizards. Parasit Vectors. 2014;7:555.

  90. 90.

    Fain A, Bannert B. Two new species of Ophionyssus Megnin (Acari: Macronyssidae) parasitic on lizards of the genus Gallotia Boulenger (Reptilia: Lacertidae) from the Canary Islands. Internat J Acarol. 2000;26:41–50.

    Article  Google Scholar 

  91. 91.

    Fain A, Bannert B. New observations on species of the genus Ophionyssus Mégnin (Acari: Macronyssidae) parasitic on lizards of the genus Gallotia Boulenger (Reptilia: Lacertidae) from the Canary Islands, Spain with description of a new species. Int J Acarol. 2002;28:361–6.

    Article  Google Scholar 

  92. 92.

    Bannert B, Karaca HY, Wohltmann A. Life cycle and parasitic interaction of the lizard-parasitizing mite Ophionyssus galloticolus (Acari: Gamasida: Macronyssidae), with remarks about the evolutionary consequences of parasitism in mites. Exp Appl Acarol. 2000;24:597–613.

    CAS  Article  PubMed  Google Scholar 

  93. 93.

    Zapatero-Ramos LM, Gonzalez-Santiago PM, Solera-Puertas MA, Carvajal-Gallardo MM. Estudio de nuevas especies de Pterigosomidae (Acari: Antinedida) sobre Gekkónidos de las Islas Canarias. Descripción de Geckobia canariensis n. sp. y de Geckobia tinerfensis n. sp. Rev Iber Parasitol. 1989;49:51–64.

    Google Scholar 

  94. 94.

    Bertrand M, Pfliegler WP, Sciberras A. Does the African native host explain the African origin of the parasite? The Maltese Geckobia estherae n. sp. parasitic on Tarentola mauritanica (Acari: Raphignathoidea: Pterygosomatidae). Acarologia. 2012;52:353–66.

    Article  Google Scholar 

  95. 95.

    Rudolf VHW, Antonovics J. Disease transmission by cannibalism: rare event or common occurrence? Proc Biol Sci. 2007;274:1205–10.

    PubMed Central  Article  PubMed  Google Scholar 

  96. 96.

    Matuschka FR, Bannert B. Cannibalism and autotomy as predator-prey relationship for monoxenous Sarcosporidia. Parasitol Res. 1987;74:88–93.

    CAS  Article  PubMed  Google Scholar 

  97. 97.

    Matuschka FR, Bannert B. Recognition of cyclic transmission of Sarcocystis stehlinii n. sp. in the Gran Canarian giant lizard. J Parasitol. 1989;75:383–7.

    CAS  Article  PubMed  Google Scholar 

  98. 98.

    Nishimoto Y, Arisue N, Kawai S, Escalante AA, Horii T, Tanabe K, Hashimoto T. Evolution and phylogeny of the heterogeneous cytosolic SSU rRNA genes in the genus Plasmodium. Mol Phylogenet Evol. 2008;47:45–53.

  99. 99.

    El-Sherry S, Ogedengbe ME, Hafeez MA, Barta JR. Divergent nuclear 18S rDNA paralogs in a turkey coccidium, Eimeria meleagrimitis, complicate molecular systematics and identification. Int J Parasitol. 2013;43:679–85.

    CAS  Article  PubMed  Google Scholar 

  100. 100.

    Stenger BL, Clark ME, Kváč M, Khan E, Giddings CW, Dyer NW, et al. Highly divergent 18S rRNA gene paralogs in a Cryptosporidium genotype from eastern chipmunks (Tamias striatus). Infect Genet Evol. 2015;32:113–23.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  101. 101.

    Hamšíková Z, Silaghi C, Rudolf I, Venclíková K, Mahríková L, Slovák M, et al. Molecular detection and phylogenetic analysis of Hepatozoon spp. in questing Ixodes ricinus ticks and rodents from Slovakia and Czech Republic. Parasitol Res. 2016;115:3897–904.

    Article  PubMed  Google Scholar 

  102. 102.

    Roca V, Jorge F, Carretero MA. Synopsis of the helminth communities of the lacertid lizards from the Balearic and Canary Islands. Basic Appl Herpetol. 2012;26:107–16.

    Google Scholar 

  103. 103.

    Loot G, Park YS, Lek S, Brosse S. Encounter rate between local populations shapes host selection in complex parasite life cycle. Biol J Linn Soc. 2006;89:99–106.

    Article  Google Scholar 

  104. 104.

    Patterson JEH, Ruckstuhl KE. Parasite infection and host group size: a meta-analytical review. Parasitology. 2013;140:803–13.

    PubMed Central  Article  PubMed  Google Scholar 

  105. 105.

    Gourbière S, Morand S, Waxman D. Fundamental factors determining the nature of parasite aggregation in hosts. PLoS One. 2015;10:e0116893.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  106. 106.

    Martín J, Garrido M, Ortega J, García-Roa R, Ibáñez A, Marzal A. Absence of haemoparasite infection in the fossorial amphisbaenian Trogonophis wiegmanni. Parasitology. 2016;4:1–4.

    Google Scholar 

  107. 107.

    Brown RP. Thermal biology of the gecko Tarentola boettgeri: comparisons among populations from different elevations within Gran Canaria. Herpetologica. 1996;52:396–405.

    Google Scholar 

  108. 108.

    Márquez R, Cejudo D, Pérez-Mellado V. Selected body temperatures of four lacertid lizards from the Canary Islands. Herpetol J. 1997;7:122–4.

    Google Scholar 

  109. 109.

    Oppliger A, Celerier ML, Clobert J. Physiological and behaviour changes in common lizards parasitized by haemogregarines. Parasitology. 1996;113:433–8.

    Article  Google Scholar 

  110. 110.

    Lazić MM, Carretero MA, Živkovi U, Crnobrnja-Isailović J. City life has fitness costs: lower body condition and increased parasite intensity in urban common wall lizards Podarcis muralis. Salamandra. 2017;53:10–7.

    Google Scholar 

  111. 111.

    Choisy M, Brown SP, Lafferty KD, Thomas F. Evolution of trophic transmission in parasites: why add intermediate hosts? Am Nat. 2003;162:172–81.

    Article  PubMed  Google Scholar 

  112. 112.

    Yoder JB, Nuismer SL. When does coevolution promote diversification? Am Nat. 2010;176:802–17.

    Article  PubMed  Google Scholar 

  113. 113.

    Jorge F, Roca V, Perera A, Harris DJ, Carretero MA. A phylogenetic assessment of the colonisation patterns in Spauligodon atlanticus Astasio-Arbiza et al., 1987 (Nematoda: Oxyurida: Pharyngodonidae), a parasite of lizards of the genus Gallotia Boulenger: no simple answers. Syst Parasitol. 2011;80:53–66.

  114. 114.

    Reyes-Betancort JA, Santos Guerra A, Guma IR, Humphries CJ, Carine MA. Diversity, rarity and the evolution and conservation of the Canary Islands endemic flora. An Jard bot Madr. 2008;65:25–45.

    Article  Google Scholar 

  115. 115.

    Sanmartín I, Van Der Mark P, Ronquist F. Inferring dispersal: A Bayesian approach to phylogeny-based island biogeography, with special reference to the Canary Islands. J Biogeogr. 2008;35:428–49.

    Article  Google Scholar 

  116. 116.

    Carretero MA, Jorge F, Llorente GA, Roca V. Relationships between helminth communities and diet in Canarian lizards: the evidence from Gallotia atlantica (Squamata: Lacertidae). J Nat Hist. 2014;48:1199–216.

    Article  Google Scholar 

  117. 117.

    Borregaard MK, Amorim IR, Borges PA, Cabral JS, Fernández-Palacios JM, Field R, et al. Oceanic island biogeography through the lens of the general dynamic model: assessment and prospect. Biol Rev Camb Philos Soc. 2016;92:830–53.

  118. 118.

    Whittaker RJ, Ladle RJ, Araújo MB, Fernández-Palacios JM, Delgado JD, Arévalo JR. The island immaturity - speciation pulse model of island evolution: an alternative to the “diversity begets diversity” model. Ecography. 2007;30:321–7.

  119. 119.

    Rijsdijk KF, Hengl T, Norder SJ, Otto R, Emerson BC, Ávila SP, et al. Quantifying surface-area changes of volcanic islands driven by Pleistocene sea-level cycles: biogeographical implications for the Macaronesian archipelagos. J Biogeogr. 2014;41:1242–54.

  120. 120.

    Weigelt P, Steinbauer MJ, Cabral JS, Kreft H. Late Quaternary climate change shapes island biodiversity. Nature. 2016;532:99–102.

    CAS  Article  PubMed  Google Scholar 

  121. 121.

    Gubitz T, Thorpe RS, Malhotra A. Phylogeography and natural selection in the Tenerife gecko Tarentola delalandii: testing historical and adaptive hypotheses. Mol Ecol. 2000;9:1213–21.

  122. 122.

    Brown RP, Hoskisson PA, Welton JH, Báez M. Geological history and within-island diversity: a debris avalanche and the Tenerife lizard Gallotia galloti. Mol Ecol. 2006;15:3631–40.

  123. 123.

    Bloor P, Kemp SJ, Brown RP. Recent volcanism and mitochondrial DNA structuring in the lizard Gallotia atlantica from the island of Lanzarote. Mol Ecol. 2008;17:854–66.

    CAS  Article  PubMed  Google Scholar 

  124. 124.

    Roy HE, Lawson Handley LJ. Networking: a community approach to invaders and their parasites. Funct Ecol. 2012;26:1238–48.

    Article  Google Scholar 

  125. 125.

    Loope LL, Mueller-Dombois D. Characteristics of invaded islands, with special reference to Hawaii. In: Drake JA, Mooney HA, di Castri F, Groves RH, Kruger FJ, Rejmánek M, Williamson M, editors. Biological invasions: a global perspective. New York: Wiley; 1989. p. 257–280.

  126. 126.

    Mateo JA, Ayres C, Lopez-Jurado LF. Los anfibios y reptiles naturalizados en España. Historia y evolución de una problemática creciente. Bol Asoc Herpetol. Esp. 2011;22:2–42.

  127. 127.

    Dickson CM, Ogbuah CT, Smith TG. The role of gamont entry into erythrocytes in the specificity of Hepatozoon species (Apicomplexa: Adeleida) for their frog hosts. J Parasitol. 2013;99:1028–33.

    Article  PubMed  Google Scholar 

  128. 128.

    Dvořáková N, Kvičerová J, Hostovský M, Široký P. Haemogregarines of freshwater turtles from Southeast Asia with a description of Haemogregarina sacaliae sp. n. and a redescription of Haemogregarina pellegrini Laveran and Pettit, 1910. Parasitology. 2015;142:816–26.

    Article  PubMed  Google Scholar 

  129. 129.

    Telford SR Jr, Wozniak EJ, Butler JF. Haemogregarine specificity in two communities of Florida snakes, with descriptions of six new species of Hepatozoon (Apicomplexa: Hepatozoidae) and a possible species of Haemogregarina (Apicomplexa: Haemogregarinidae). J Parasitol. 2001;87:890–905.

    Article  PubMed  Google Scholar 

Download references


We thank the Cabildos Insulares (Island Authorities) from Lanzarote, Fuerteventura, Gran Canaria, Tenerife, La Palma, La Gomera and El Hierro of Spain for the lizard collection permits; the Centro Insular para la Cría en Cautividad del lagarto gigante de La Gomera (Gallotia bravoana) and the Centro de Recuperación del Lagarto Gigante de El Hierro. Special thanks to M. López-Darias, B. Fariña, Z. Tobias, A. Martínez-Silvestre, A. Kaliontzopoulou, and Applied Phylogenetics CIBIO research group for their contribution to this work. Lastly, we acknowledge the contribution of the anonymous reviewers, whose comments helped to improve this manuscript.


BT is funded by a Fundação para a Ciência e Tecnologia (FCT) PhD grant (PD/BD/52601/2014), and this study is part of her thesis dissertation. AP2 (Ana Perera) and DJH are funded through FCT contracts (IF/01257/2012 and IF/01627/2014) under the Programa Operacional Potencial Humano - Quadro de Referência Estratégico Nacional from the European Social Fund and Portuguese Ministério da Educação e Ciência. MAC is supported by project NORTE-01-0145-FEDER-000007. This study was funded by the IF exploratory project IF01257/2012/CP0159/CT0005, and by PTDC/BIA-BDE/67678/2006 and PTDC/BIA-BEC/105327/2008 both by FCT and FCOMP-01-0124-FEDER-007062 and FCOMP-01-0124-FEDER-008970 COMPETE program, respectively.

Availability of data and materials

Summary data are available in the supplementary material and genetic sequences are deposited in the GenBank database under the accession numbers MG787243-MG787253). Any additional analysis or data included in this study are available upon request from the authors.

Author information




BT, FJ, DJH and AP2 (Ana Perera) conceived and designed the study; MAC and AP2 provided the funding; FJ, MAC and AP2 collected the samples; BT, AP1 (Ana Pereira) and AP2 did the lab work (blood slide screening and PCR amplifications); BT performed the phylogenetic analysis; BT and AP2 did the statistical analysis and main manuscript writing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Beatriz Tomé.

Ethics declarations

Ethics approval and consent to participate

Regarding ethical consideration on animal sampling and handling, all samples were collected under permits issued by the environmental authorities of the Canary Islands (Cabildos Insulares of Lanzarote, Fuerteventura, Gran Canaria, Tenerife, La Palma, La Gomera and El Hierro), and of Morocco (Haut Commissariat aux Eaux et Forêts et à la Lutte contre la Désertification).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Sampling locations and prevalence values. (DOCX 31 kb)

Additional file 2:

Table S2. Uncorrected p-distances between the nine haemogregarine haplotypes discovered in the Canarian lizards. (DOCX 13 kb)

Additional file 3:

Table S3. GenBank codes and additional information on the genetic sequences used for the phylogenetic analyses. (DOCX 26 kb)

Additional file 4:

Table S4 Statistical summary of the variation of haemogregarine prevalence and parasitaemia from Canarian lizards. (DOCX 18 kb)

Additional file 5:

Figure S1. Networks of the haemogregarine haplotypes infecting Gallotia, Chalcides (a) and Tarentola lizards (b). Represented is the network estimations performed with TCS. Black nodes represent mutations, while coloured ones correspond to the sequences of Canarian haemogregarines (colours match the ones used in Fig. 3), and the white ones correspond to the remaining samples. The network of a includes the samples of the haemogregarine clade from lacertid, snake, skink and varanid hosts (as identified in Additional file 4: Table S4) and b includes the samples of the clade from gecko hosts. (PDF 911 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tomé, B., Pereira, A., Jorge, F. et al. Along for the ride or missing it altogether: exploring the host specificity and diversity of haemogregarines in the Canary Islands. Parasites Vectors 11, 190 (2018).

Download citation


  • Apicomplexa
  • Hemoparasites
  • Lizards
  • Gallotia
  • Tarentola
  • Chalcides
  • Insularity
  • Phylogeography
  • 18S rRNA