Skip to main content

Prevalence of co-infection and genetic diversity of avian haemosporidian parasites in two rehabilitation facilities in Iran: implications for the conservation of captive raptors

Abstract

Background

Various haemosporidian parasites infect raptors, especially captive hosts who may be more exposed. Diagnosis of threatening factors such as infectious diseases indirectly has a significant role in protecting endangered or threatened species that may boost the mortality or extinction resulting from declined reproduction. Few investigations have been performed in captive hosts to detect the prevalence of haemosporidian parasites and define genetic diversity in west Asia. For the first time, the current study was designed to determine the prevalence and genetic diversity of haemosporidian parasites in captive raptors by molecular methods in two rehabilitation facilities in North and North-east Iran and to define phylogenetic relationships of detected lineages circulating in raptors.

Results

Molecular characterization of the haemosporidian parasite was accomplished by PCR-based method and DNA sequencing in 62 captive raptors. The overall prevalence was ~ 36% with higher infection of Haemoproteus spp. than Leucocytozoon spp. Plasmodium infection was not detected in any host. Results showed that 22 individuals (of 10 species) were infected with unique lineages. Genus Haemoproteus was detected in 26.66% of examined individuals (of eight species) and Leucocytozoon was found in 10% of individuals (of four species). The molecular analysis could detect ten lineages (nine Haemoproteus spp. and one Leucocytozoon spp.) which were categorizes as new and six lineages which have been previously detected in the other investigations.

Conclusions

The Bayesian phylogenetic analysis derived from obtained data in the present study and published lineages in previous investigations indicated the probable host specificity of Haemoproteus and Leucocytozoon parasites in several sub-clades at hosts’ order and genus level. As monitoring the parasite loads of captive birds when admitted reduce the risk of infecting hosts in captivity at those locations, we designed this study to determine infection prevalence and genetic diversity of blood parasites in raptors examined in Iran. These results allow mapping of haemosporidian distribution and shed light on the depth of their diversity in Iran to protect species by identification of risk in rehabilitation facilities.

Peer Review reports

Background

Haemosporidian parasites (Apicomplexa: Haemosporida); Haemoproteus, Plasmodium, and Leucocytozoon species are cosmopolitan vector-borne organisms that are transmitted by biting hematophagous arthropods [1,2,3,4,5,6,7]. Haemoproteosis, avian malaria, and Leucocytozoonosis resulting from blood parasites may have great influences on their hosts’ life [2, 8,9,10,11,12,13,14,15]. Molecular characterization of haemosporidian parasites by using PCR-based methods has been extensively applied as quick, specific, and sensitive approaches in wildlife parasite screening and estimation of genetic diversity [16, 17]. Some genetic markers, such as the mitochondrial cytb gene have been barcoded more than 4100 distinct lineages in MalAvi [18]. These compiled datasets from various bird hosts around the world [2, 18,19,20] assign about 20% of molecular lineages to morphospecies [21]. Traditional microscopy method have described more than 250 haemosporidian morphospecies which can be applied as a complementary technique for screening parasites. The infection of each parasite is influenced by parasite strain, host species, and host-parasite-vector interactions in various geographical zones [22, 23]. Symptoms of haemosporidian infection may vary from no obvious clinical signs (asymptomatic) to severe signs; anemia, breath problems, thinness, lameness, poor appetite, and death [3, 14, 22, 24,25,26,27]. Sub-lethal influences of blood parasites in wild birds have not been completely understood. These parasites appeared harmless but some of these diverse intracellular parasitic species are the causative agents for serious avian diseases may promote clinical cases and mortality in non-adapted avian hosts [11, 28,29,30]. Consequently, these parasites pose a health concern, particularly in captive animals maintained at high densities in rehabilitation facilities, zoos, and gardens with imperfect resistance and increased stress due to human handling [31,32,33,34].

In this case, bird species that are under conservation concerns should be prioritized [35,36,37,38,39]. Although, birds have been transferred among zoological aviaries, rehabilitation centers, or zoos, there is no document on their health status or previous exposure to disease.. Thus, it is paramount to screening of blood parasites in captive hosts at admission to any rehabilitation facility, zoo, or gardens. Diagnosis of threatening factors such as infectious diseases indirectly plays a vital role in protecting endangered or threatened species by identification of risk factors that may boost the mortality or extinction caused by the reduced success of reproduction [40]. Monitoring the parasite loads of captive birds when admitted reduces the risk of infecting hosts in captivity at those locations [41,42,43,44,45,46]. No published investigation is available for the potential distribution of insect vectors of avian haemosporidian parasites in Iran. However, a low abundance of competent vectors, sampling time, and sample size may have influenced the low prevalence of detected parasites in the current study results. The small sample size was a restriction to defining host specificity. Bayesian analysis provide comparative results with published lineages in raptors that demonstrated the probable host specificity of several sub-clades at the order, family, and genus level. Haemoproteus and Leucocytozoon species clades displayed host specificity to some degree compared with Plasmodium spp.

Birds of prey are important in ecosystems and many studies have examined the genetic diversity, prevalence variability, host specificity, and description of new parasite lineages and species [40, 47,48,49,50,51]. Few investigations have targeted haemosporidian parasites in wild birds and captive hosts [12, 15, 52,53,54,55,56,57] to define the distribution and genetic diversity of parasites and vectors [12]. Because the identification of haemosporidia in the erythrocytic phase is required for conservation measures and protection of birds during rehabilitation programs [58], the present investigation aimed to determine the molecular prevalence of haemosporidian in captive raptors in two rehabilitation facilities and to assess the lineages diversity, and perceive the phylogenetic relationships using nested PCR amplification and DNA sequencing.

Methods

Blood samples collection and DNA extraction

This investigation was carried out by collecting blood samples from captive birds in two rehabilitation facilities in Khorasan Razavi and Golestan provinces, north and north-east of Iran, April to July 2015 and 2016 (Table 2). 50–100 µl of whole blood were taken from the brachial vein of sixty-two adult captive birds, which belong to three orders; Accipitriformes (~ 66%), Falconiformes (~ 26%), and Strigiformes (~ 8%) (Table 2), and to 13 species; Accipiter badius (n = 8), Aquila chrysaetos (n = 2), Aquila rapax (n = 5), Buteo buteo (n = 21), Circus aeruginosus (n = 2), Haliaeetus albicilla (n = 1), Milvus migrans (n = 1), Neophron percnopterus (n = 1), Falco subbuteo (n = 1), Falco tinnunculus (n = 15), Athene noctua (n = 3), Bubo bubo (n = 1), and Otus scops (n = 1). The examined species consist of 11 least concern species, one vulnerable (Aquila rapax) and one endangered species (Neophron percnopterus). Eight species are resident in Iran and the remaining are categorized as winter resident in some parts of country, migrants or summer breeders, shown in Table 2 [59]. Blood was preserved in Queen’s buffer for molecular experiments [60]. Genomic DNA was extracted using PrimePrep Genomic DNA Isolation Kit (GENETBIO Inc. Daejeon, South Korea) according to the manufacturer's instructions. The isolated DNA concentration and quality were estimated using a spectrophotometer (DeNovix Inc. USA).

Nested PCR assay for parasites screening

A partial amplification of mt-DNA cytb gene (479 base pair) of the haemosporidian parasites was performed by nested PCR [16, 61]. The final volume of PCR reactions included 50 ng/µl of whole genomic DNA, 12.5 μl Ampliqon PCR master mix (AMPLIQON, Denmark), 0.6 mM of each primer and nuclease free water (up to 25 μl). To achieve the concise prevalence of haemosporidian parasites in raptors, we used two sets of primers. For the first PCR set, three parasitic genera were amplified by using general standard primers HaemNF1/HaemNR3 in the first reactions and HaemF/HaemR2 and HaemFL/HaemR2L in the second reactions (Table 1) [16, 61]. PCR cycling was performed as primary denaturation at 95 °C (5 min), annealing at 50 °C (30 s), extension at 72 °C (45 s), and finally followed by a final extension at 72 °C (10 min) which was run for 20 cycles in 1st nested PCRs and 35 cycles in 2nd nested PCRs. Positive amplicons of the previous study and ultra-pure ddH2O were utilized as positive and negative controls for each PCR reaction set. PCR products were visualized on 2% agarose gel. Purification and sequencing were carried out using 20 μl of PCR products (479 bp) by Macrogene Co. (Seoul, South Korea). In addition, specific primers Plas1/HaemNR3 and 3760F/HaemJR4 for the detection of these parasites in raptors were used according to suggested procedures (Table 1) [47]. Purification and bi-directional sequencing of PCR products (524 bp) was performed by Codon Genetic Group (Tehran, Iran). Two set of primers were used for all host samples. The result showed that specific primers didn’t amplify all lineages which were detected by standard primers and vice versa; some of detected lineages were new.

Table 1 Primers used for nested PCR

DNA sequencing and phylogenetic analysis

The novelty of amplified sequences’ was confirmed by “Nucleotide BLAST” in National Center for Biotechnology Information (NCBI) (Basic Local Alignment Search Tool) search [24] and MalAvi blast hits [18]. Mitochondrial lineage with one or more substitutions was considered a new lineage [16], [16]. The amplified sequences were named based on MalAvi nomenclature [18]. Raw sequences were edited and aligned using Bioedit [63] and online MAFFT programs [64]. Three sequences from two individuals of Buteo buteo and two Aquila rapax, and one Athene noctua who had mature Haemoproteus species gametocytes in their blood smears, were excluded from the analysis because of having double peaks in the electropherograms.

To perceive the phylogenetic relationship of haemosporidian parasites with previously described lineages circulating in raptors, a dataset containing our sequences and those present on the online database MalAvi was constructed to insert our sequences in a wider context (sequences of similar hosts’ haemosporidian sequences were retrieved) (Fig. 1). A table containing the general information of the sequences belonging to the MalAvi dataset is given in Additional File 1. The sequences divergence of blood parasite lineages was measured with the Kimura 2-parameter (K2P) distance matrix, implemented in MEGA6.0 using the pairwise comparison of all sequences [65] and interspecific pairwise distances were calculated by converter program, Excalibar v.1.1, is shown in Fig. 2 and Table 4 [66].

Fig. 1
figure 1figure 1

Bayesian analysis of haemosporidian cytb lineages in raptors data deposited in MalAvi, GenBank and present study sequences. Posterior probability of Bayesian analysis are specified for each branch (> 80). Detected lineages of the current investigation are in bold. Host species of Accipitriformes, Falconiformes, and Strigiformes are specified with red, blue, and green colors. Schematic image of entire tree which is separated to each genus (a), Haemoproteus spp. (b), Plasmodium spp. (c), Leucocytozoon spp. (d)

Fig. 2
figure 2

Comparison of interspecific pairwise Kimura 2-parameter (K2P) distances for haemosporidian parasitic genus detected in raptors. Haemoproteus (red), Plasmodium (green), and Leucocytozoon (blue) are specified

Barcoding gap reconstruction for each parasite genus was performed by ABGD analysis [67] run online (https://bioinfo.mnhn.fr/abi/public/abgd/abgdweb.html). This method finds a range of previous intraspecific distances (calculated from Pmin to Pmax, with P Steps). A proxy for the minimum relative gap width is shown by the X value (relative width). The selected distance calculation method was Kimura (K80) TS/TV and relative width = 1 for Haemoproteus and Leucocytozoon and relative width = 0.75 for Plasmodium. The other variables were based on default settings (Additional file 2 and Additional file 3).

Phylogenetic analyses were performed by Bayesian inference executed in MrBayes v3.2 software [68] with two concurrent MCMC simulations for 10,000,000 generations and sampled every 1000 generations. The posterior probabilities were measured at the end of the analysis when the burn-in period of 50% was set and the chains reached stationary status. ModelTest v.7 software was used for the best evolutionary model selection (GTR + I + G) [69]. The resultant phylogenetic tree is visualized by FigTree v1.4 [70]. Moreover, the analysis of Poisson Tree Processes (PTP) was performed to reveal the species delimitation [71] on the online server https://species.h-its.org/. PTP algorithm uses the number of mutations (inferred from branch lengths) for speciation modeling on branches. The resultant tree showed the position of lineages at the genus level (Additional file 4).

Results

Prevalence of haemosporidian parasites and genetic diversity of identified lineages

The overall prevalence of detected haemosporidian parasites by using nested PCR assay was 36.66% (95% CI 24.11, 49.22) in comparison with the reported prevalence of 28.33% (95% CI 16.59, 40.07) by microscopic examinations in our former study (Table 2). The molecular method revealed that 22 individuals of each ten species were infected with one haemosporidian lineage. Compared results of prevalence for genus and identification method are presented as follows: morphological method; Haemoproteus spp. 25% (95% CI 13.72, 36.28), Leucocytozoon spp. 3.33% (95% CI 1.32, 8.01) and molecular method; Haemoproteus spp. 26.66% (95% CI 15.15, 38.19), Leucocytozoon spp. 10% (95% CI 2.18, 17.82). Microscopic results have shown that 15 individuals (24.1%) were infected by Haemoproteus spp. and two individuals were infected (3.22%) by Leucocytozoon spp. The obtained result was reported at the genus level. The compared results of microscopic [52] and molecular methods presented in this study revealed that the molecular method could detect more infections rather than morphological examination but failed to amplify five morphologically positive parasites which may be co-infected. Regardless of amplification during PCR, the blast hit showed no results, and the presence of double peaks seen in sequences enabled us to omit them from analysis and prevalence calculation.

Table 2 Haemoproteus and Leucocytozoon lineages isolated from raptors in rehabilitation facilities in Iran

All bird hosts, except for the golden eagle (Aquila chrysaetos), tawny eagle (Aquila rapax), and Egyptian vulture (Neophron percnopterus) were parasitized by haemosporidian lineages (Table 2). 16 of 62 individuals were positive in the nested-PCR test for Haemoproteus spp., and six samples were positive targeting Leucocytozoon parasite infections. One Eurasian buzzard individual (Buteo buteo) was co-infected by Haemoproteus and Leucocytozoon species. Plasmodium parasite infection was not observed in morphological and molecular examinations. 14 identified infected individuals sampled from a rehabilitation facility in Khorasan Razavi (one co-infected sample) and eight individuals from Golestan province centers. The distribution of previously detected lineages in similar hosts examined in this study and other countries is summarized in Table 3. The compared data showed that most infections in other countries were from Leucocytozoon spp. Buteo buteo and Milvus migrans have shown the most variety of infecting lineages. Similar to our investigation, no haemosporidian lineage has been discovered in Neophron percnopterus and Aquila rapax around the world (Table 3).

Table 3 Contribution of known haemosporidian lineages in raptors

Detected Haemoproteus spp. lineages include hACCBAD02, hBUTBUT04, hBUTBUT05, hBUTBUT06, hCIAE04, hMILMIG01, hFALTIN01, hLK03, hFALSUB02, hTYTAL04, and hOTSCO07. In addition to, the identified Leucocytozoon spp. lineages consist of lATN01, lMULVIS01, lBUBT03, lCIAE02, and lBUBO05. Haemoproteus spp. lineages were detected from eight bird species; Accipiter badius, Buteo buteo, Circus aeruginosus, Milvus migrans, Falco tinnunculus, Falco subbuteo, Athene noctua, and Otus scops, and Leucocytozoon spp. lineages were discovered in four raptors species; Buteo buteo, Haliaeetus albicilla, Athene noctua, and Bubo bubo. New lineages were detected in Shikra (hACCBAD02), Eurasian buzzards (hBUTBUT04, hBUTBUT05, hBUTBUT06), Western marsh harrier (hCIAE04), Black kite (hMILMIG01), Common kestrel (hFALTIN01), Eurasian hobby (hFALSUB02), Eurasian scops owl (hOTSCO07), and Eurasian eagle owl (lBUBO05). Common kestrel, Little owl, Eurasian buzzard, Little owl, Eurasian buzzard, and White-tailed eagle parasitized by lineages hLK03, hTYTAL04, lATN01, lBUBT03, lMULVIS01, and lCIAE02, and were introduced as new host records examined in Iran (Scientific names of infected hosts are given in Table 2). Lineage hBUTBUT04 was found in two individuals of the same host species, and lATN01 (Accipitriformes and Strigiformes) and hLK03 (Falconiformes and Strigiformes) were discovered in two different orders (Table 2). Furthermore, using molecular approach revealed the co-infection of one Buteo buteo by two lineages; hBUTBUT04 and lBUBT03 (Table 4).

Table 4 Genetic information for haemosporidian parasite’ sequences detected in raptors and the retrieved data from MalAvi

Interspecific distances and phylogeny of known haemosporidian lineages in raptors

The Bayesian phylogenetic tree was reconstructed by the present data detected in this investigation and deposited sequences of all lineages that infect raptors in MalAvi (n = 249). The results demonstrated the placement on three distinctive clades; A (Leucocytozoon spp.), D (Plasmodium spp.), and C (Haemoproteus spp.) (Fig. 1). The available data of infection in raptors around the world included sequences retrieved from wild (87%) and captive individuals (13%) in Asia (49%), Europe (22.3%), South America (8.7%), North America (15.3%), Africa (3.9%), and unknown places (0.9%). Average K2P distance was calculated to show whether there is a sequence variability for each parasitic genus. These results showed different genetic variations for each genus (Fig. 2). The highest and lowest interspecific distances range belonged to Leucocytozoon spp. (0–40.4%) and Haemoproteus spp. (0–16.7%). The maximum and the minimum mean of K2P genetic distances (%) was calculated to Leucocytozoon spp. (17.42%) and Haemoproteus spp. (6.45%). The number of conserved, variable, parsimony, and singleton sites is mentioned in Table 3. The highest number of conserved sites (73.4%) belonged to Haemoproteus and Leucocytozoon had the highest proportion of variable sites (45.81%), parsimony informative sites (39.74%), and GC content (28.4%). The maximum number of singleton sites belonging genus Plasmodium (6.48%) (Table 3). In the Bayesian phylogenetic tree, some Leucocytozoon spp. sequences are placed within host order-specific sub-clades. Except for some lineages, A2 is Accipitriformes-specific and A3 is Strigiformes-specific. Within A5, all but one raptor parasite fall into Accipiter-specific and Buteo-specific monophyletic sub-clades. Returning to sub-clade A6, all lineages are parasites of Buteo except MILVUS01 detected in Milvus milvus and M. migrans from Europe. Disregarding the placement of BUTBUT03 and BUBT3 detected in Buteo buteo, the remaining raptor parasite lineages fall into Accipiter sub-clade A7. Within A1, excluding several lineages, remaining lineages in A3 are positioned as Strigiformes-specific sub-clade while the other sub-clade (A4) includes parasites from three raptors orders. Plasmodium spp. clade is sister taxon to Haemoproteus spp. clade (clade C) and contains parasites associated with three raptors orders. Except for pMYCAME02 and pRTSR1, the remaining Plasmodium lineages are discovered in one order (clade D). Within Haemoproteus spp., raptor parasites are monophyletic and consisted of two sub-clades, C1 (C2; Falconiformes-specific and C3; Accipitriformes-specific) and C4 (Strigiformes-specific, except for several lineages). Most of the detected lineages of this study are grouped within the C2 sub-clade (Fig. 1a–d). Sub-clade C4 and A1 comprise exclusively Strigiformes parasites (with some exceptions) while clade A2 and C3 contain parasites hosted by Accipitriformes and A6 consist of Buteo parasites from North American and Accipiter. In addition, the sub-clade of A1 (Leucocytozoon spp.) and clade D (Plasmodium spp.) are included more generalist lineages (Fig. 1a–d).

The position of lineages in the maximum likelihood (ML) tree showed the similarity of sub-clades and also some differences at the genus level. In comparison with the Bayesian tree where Haemoproteus and Plasmodium lineages were positioned as sister clades and Leucocytozoon was situated in a basal out-group, the location of lineages in the ML tree was different and three parasite genera were placed in a paraphyletic manner (Additional file 4). So, for the precise interpretation of species delimitation, using an extensive dataset is essential and the analysis should be concluded more cautiously by different algorithms.

Discussion

In the current investigation, the molecular characterization of haemosporidian parasites infecting ten species of captive birds was established in two rehabilitation facilities in Iran. Several raptors were new hosts from this region required such investigations to clarify the diversity of blood parasite species and their hosts through conservation programs. Numerous studies have reported the high prevalence of haemosporidian parasite infections in captive birds leading to immense mortality than in wild hosts [72, 73]. Investigating birds in captivity provides an opportunity for researchers and students to explore the maximum host range of parasite lineages which are very hard to sample in nature [74]. Birds in captivity are not in a natural environment and may have less resistance to parasitic agents in a situation imposing increased stress. Such investigations may provide more information about health status, deterioration factors, and death in captive birds referable wild hosts [75, 76]. Various haemosporidian parasites infect raptors, especially captive hosts who may be more exposed. The captivity situations may lead to infection outbreaks in domestic and wild hosts [32, 77]. Improved and periodic screening, molecular epidemiologic investigations, and control of possible vectors in managerial programs from a conservation viewpoint are recommended as essential procedures to protect valuable raptor hosts. Lack of suitable foods and stressful conditions can significantly increase infection prevalence and decrease the immune responses [73, 77] that highlight the requirement of protecting wild hosts that may have interacted with captive birds or be infected by potential vectors in adjacent locations. These unusual interactions may cause the transmission of parasites into new wild or migrant birds, such as the mortality that occurred in captive psittacines parrots, Cyanoramphus auriceps, Bolborhynchus lineola, and Melopsittacus undulates in two different aviaries from Switzerland and Germany which showed Haemoproteus infections [78, 79]. The mortality of Egyptian geese (Alopochen aegyptiaca) in São Paulo zoo and European parrots were attributed to P. nucleophilum and H. minutus infections, respectively [77, 80]. Histopathological and molecular approaches documented the presence of Plasmodium spp., as a causative agent for mortality in raptor species at a safari park in Italy [73]. A recorded prevalence of 12.6% in 677 hosts sampled in the São Paulo zoo, showed infection to 14 lineages of Plasmodium spp. and 2 Haemoproteus spp. of which 8 were novel [62]. Wild perching birds have been considered as potentially infected hosts for transmission of infection to captive birds by vectors’ biting [81]. The high rate of mortality in captive Magellanic penguins (Spheniscus magellanicus) in a zoo in Southern Brazil was confirmed by Plasmodium parasite infections (pTURALB01), a previously detected lineage in a passerine host, Turdus albicollis [82]. Moreover, migratory species are acknowledged frequently as active carriers of various haemosporidian lineages to infect captive and/or non-migrant species [83].

The infection prevalence and genetic diversity of haemoparasites of captive raptors examined in Iran have not been studied in detail. In present study, ten new lineages and six previously known Haemoproteus spp. and Leucocytozoon spp. lineages are detected in examined hosts from this region. In the current study, no captured bird was infected by Plasmodium spp. Frequent elements may be related to infection prevalence of haemosporidian in birds, encompassing host-related features; species, sex, host-vector-parasite relationships, immune system, ecological factors; study site, habitats, sampling time, detection procedures; and vector’s related characteristics; host specificity and vector competency [2, 20, 50, 72, 84]. PCR-based identification method revealed Haemoproteus and Leucocytozoon of ~ 36% of examined individuals. The highest infection of Haemoproteus spp. was found in Accipitriformes, and Buteo buteo had the highest number of captured individuals in both rehabilitation facilities. These results are consistent with previous investigations in Iran that have shown a high diversity of Haemoproteus spp. lineages in detected hemoparasites [52,53,54,55,56]. The higher prevalence of Haemoproteus spp. (78%) in comparison with Leucocytozoon spp. (35%) in 55 captive raptors in another part of the world was consistent with our investigation [85]. Both a very high overall prevalence of 97.6% (genus Plasmodium) and a very low prevalence of 2.4% (genus Haemoproteus) have been recorded in captive birds in Brazil [62]. Application of PCR-based methods for detection of blood haemosporidian parasites in 167 individuals of owls belonging to 12 species from Thailand showed an overall prevalence of 34.1% and reported 17 new detected lineages [58]. Their findings may indicate high transmission of Haemoproteus and Plasmodium species in owls in comparison with previous studies in Asia, Korea (62.8%) [86], and Japan (57.1%) [87]. The low prevalence of blood parasites (16.7%) in a collection of 324 Falco eleonorae nestling and adult individuals highlighted the possible effect of ecological factors such as occupied habitats, host-vector-parasite relationship, and hosts immune system efficacy [50].

Among Accipitriformes, Accipiter badius has not been previously identified as a Haemoproteus novel host, except for an unpublished record in MalAvi in which pACCBAD01 was recorded from Thailand (Salakij et al. unpubl, recorded in MalAvi). Haemoproteus spp. lineage hACCBAD02 isolated from Accipiter badius, displayed 99% similarity as opposed to the nearest lineage hARBRU02. Circus aeruginosus as a new host examined in Iran for Haemoproteus spp. join records of this species parasitized by pCIAE01 and lCIAE02 in Germany [88]. Novel lineage hCIAE04 detected in Circus aeruginosus demonstrated 99% similarity with hPLONEL01 recorded in the MalAvi database. In addition, Haliaeetus albicilla is a new host for lCIAE02 detected in birds examined in Iran. Lineage hLK03 as the reference lineage for hFALTIN01 and hFALSUB02 has been detected in Falco naumanni and F. tinnunculus from Spain and Germany [88, 89]. Although there have been limited studies on blood parasites of raptors, the lineage of hLK03 has not yet been found in other avian families. Our results presented the infection to hLK03 in F. tinnunculus (Falconiformes) and Athene noctua, as a new order member infection (Strigiformes). Lineage hTYTAL04 detected in Athene noctua, is reported as another new host from this region. We isolated Haemoproteus spp. hBUTBUT04 from two individuals of Buteo buteo; this lineage shows 99% similarity with hSISKIN1 has been reported in various species of Passeriformes [90, 91] and Piciformes [56]. The newly discovered lineage hBUTBUT06 was 99% identical to the reference sequence hPIRUB01 has been previously recorded in Pitangus sulphuratus from Brazil [92]. The other new Haemoproteus spp. lineage hBUTBUT05 discovered from Buteo buteo, showed 99% similarity in comparison with hARBRU02 was recovered in Buarremon brunneinucha from South America [93, 94]. Haemoproteus spp. lineage hCELEC01 was the nearest lineage to newly isolated hOTSCO07 in Otus scops (99% similarity) had been registered on the MalAvi online database. Haemoproteus spp. lineage hMILMIG01 found in Milvus migrans in the present study was similar to hMILANS03 (99% similarity) was isolated from the same host species in Spain [47]. The shared lineage lATN01 was discovered in Athene noctua and Buteo buteo belonging to two different orders examined in Iran. This lineage was previously recovered in Athene noctua on other continents from Morocco and Portugal [95]. However, lineages lMULVIS01 and lBUBT03 detected in Buteo buteo join records of Leucocytozoon spp. in this species.

From a conservation perspective, the discovery of co-infections is related to virulence [96] and there are some limitations in the identification of Plasmodium and Haemoproteus species by PCR-based methods, therefore microscopic method can serve as a diagnostic method [97]. In the current study, we focused on molecular results, and our morphological examination results have confirmed the parasite infection at the genus level [52].However, we used two sets of primers to identify general and specific lineages for raptors ambiguous sequences with double peaks, due to un-successful amplification of parasite or co-infections within the same host, necessitate the re-examination by other diagnostic tools in further investigation. Molecular diagnostic methods have been established comparatively as applicable as traditional techniques using light microscopy for haemosporidian parasites detection [16, 78, 96, 98]. On the other side, morphological procedures may even undervalue low parasitemia, abortive infections (unsuccessful development of the parasite in which gametocytes are not in smears), or cryptic species, when compared with molecular methods are known for underestimation co-infections [97, 99]. As the detection of simultaneous infections is frequently common in avian hosts [100], routine molecular methods not able to validate co-infections. The microscopic method can identify co-infections. Therefore, a combination of molecular and microscopic methodologies increases the probability of discovery of haemosporidian infections at the species level in wild birds [101, 102]. Besides, the parallel application of diverse primers [96, 103], multiplex PCR [9], RT-PCR [10, 33], high-resolution melt (HRM) analyses [23], and next-generation sequencing (NGS) are recommended to increase the sensitivity of parasite diversity in haemosporidian detection [7]. Nonetheless, these approaches are expensive and time-consuming and may not be applicable in every geographical region in surveillance investigations.

The significant finding of the current investigation was the molecular screening of avian haemosporidian blood parasites in captive raptors maintained in two rehabilitation facilities and documentation of novel information on the genetic diversity of haemosporidian lineages in an unexamined geographical region in west Asia. As the monitoring parasite loads of captive birds when admitted reduce the risk of infecting hosts in captivity at those locations, we designed this study to determine blood parasites prevalence and genetic diversity in raptors examined in Iran. The results allow mapping of haemosporidian distribution and shed light on the depth of their diversity in Iran to protect species by identification of risk in rehabilitation facilities.

Availability of data and materials

Amplified sequences were submitted in the GenBank under accession numbers MN224217-MN224231 and MW209702-MW209713.

References

  1. Santiago-Alarcon D, Havelka P, Schaefer HM, Segelbacher G. Bloodmeal analysis reveals avian Plasmodium infections and broad host preferences of Culicoides (Diptera: Ceratopogonidae) vectors. PLoS ONE. 2012;7(2): e31098.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Valkiunas G. Avian malaria parasites and other haemosporidia. CRC Press; 2004.

    Book  Google Scholar 

  3. Atkinson CT, Dusek RJ, Woods KL, Iko WM. Pathogenicity of avian malaria in experimentally-infected Hawaii Amakihi. J Wildl Dis. 2000;36(2):197–201.

    Article  CAS  PubMed  Google Scholar 

  4. Remple JD. Intracellular hematozoa of raptors: a review and update. J Avian Med Surg. 2004;18(2):75–88.

    Article  Google Scholar 

  5. Tanigawa M, Sato Y, Ejiri H, Imura T, Chiba R, Yamamoto H, Kawaguchi M, Tsuda Y, Murata K, Yukawa M. Molecular identification of avian haemosporidia in wild birds and mosquitoes on Tsushima Island, Japan. J Vet Med Sci 2012:12–0359.

  6. Perkins SL. Malaria’s many mates: past, present, and future of the systematics of the order Haemosporida. J Parasitol. 2014;100(1):11–26.

    Article  PubMed  Google Scholar 

  7. Jarvi S, Farias M, Lapointe D, Belcaid M, Atkinson C. Next-generation sequencing reveals cryptic mtDNA diversity of Plasmodium relictum in the Hawaiian Islands. Parasitology. 2013;140(14):1741–50.

    Article  CAS  PubMed  Google Scholar 

  8. Asghar M, Hasselquist D, Bensch S. Are chronic avian haemosporidian infections costly in wild birds? J Avian Biol. 2011;42(6):530–7.

    Article  Google Scholar 

  9. Ciloglu A, Ellis VA, Bernotienė R, Valkiūnas G, Bensch S. A new one-step multiplex PCR assay for simultaneous detection and identification of avian haemosporidian parasites. Parasitol Res. 2019;118(1):191–201.

    Article  PubMed  Google Scholar 

  10. Sijbranda D, Gartrell B, Grange Z, Howe L. Use of a real-time PCR to explore the intensity of Plasmodium spp. infections in native, endemic and introduced New Zealand birds. Parasitology. 2017;144(13):1743–51.

    Article  CAS  PubMed  Google Scholar 

  11. Marzal A, De Lope F, Navarro C, Møller AP. Malarial parasites decrease reproductive success: an experimental study in a passerine bird. Oecologia. 2005;142(4):541–5.

    Article  PubMed  Google Scholar 

  12. Nourani L, Zakeri S, Djadid ND. Dynamics of prevalence and distribution pattern of avian Plasmodium species and its vectors in diverse zoogeographical areas-a review. Infect Genet Evol 2020:104244.

  13. Synek P, Albrecht T, Vinkler M, Schnitzer J, Votýpka J, Munclinger P. Haemosporidian parasites of a European passerine wintering in South Asia: diversity, mixed infections and effect on host condition. Parasitol Res. 2013;112(4):1667–77.

    Article  CAS  PubMed  Google Scholar 

  14. Merino S, Moreno J, Sanz JJ, Arriero E. Are avian blood parasites pathogenic in the wild? A medication experiment in blue tits (Parus caeruleus). Proc Royal Soc London B Biol Sci. 2000;267(1461):2507–10.

    Article  CAS  Google Scholar 

  15. Nourani L, Dinparast Djadid N. Molecular identification of avian malaria in wild birds from Iran. In: 2nd International Congress of Climate Change & vector borne Diseases and 4th National Congress of Medical; Shiraz, Iran. 2019.

  16. Bensch S, Stjernman M, Hasselquist D, Örjan Ö, Hannson B, Westerdahl H, Pinheiro RT. Host specificity in avian blood parasites: a study of Plasmodium and Haemoproteus mitochondrial DNA amplified from birds. Proc Royal Soc London B Biol Sci. 2000;267(1452):1583–9.

    Article  CAS  Google Scholar 

  17. Ricklefs RE, Fallon SM, Bermingham E. Evolutionary relationships, cospeciation, and host switching in avian malaria parasites. Syst Biol. 2004;53(1):111–9.

    Article  PubMed  Google Scholar 

  18. Bensch S, Hellgren O, Pérez-Tris J. MalAvi: a public database of malaria parasites and related haemosporidians in avian hosts based on mitochondrial cytochrome b lineages. Mol Ecol Resour. 2009;9(5):1353–8.

    Article  PubMed  Google Scholar 

  19. Clark NJ, Clegg SM, Lima MR. A review of global diversity in avian haemosporidians (Plasmodium and Haemoproteus: Haemosporida): new insights from molecular data. Int J Parasitol. 2014;44(5):329–38.

    Article  PubMed  Google Scholar 

  20. Valkiūnas G, Iezhova TA. Keys to the avian malaria parasites. Malar J. 2018;17(1):212.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Clark NJ, Adlard RD, Clegg SM. Molecular and morphological characterization of Haemoproteus (ParaHaemoproteus) ptilotis, a parasite infecting Australian honeyeaters (Meliphagidae), with remarks on prevalence and potential cryptic speciation. Parasitol Res. 2015;114(5):1921–8.

    Article  PubMed  Google Scholar 

  22. Knowles SC, Wood MJ, Alves R, Wilkin TA, Bensch S, Sheldon BC. Molecular epidemiology of malaria prevalence and parasitaemia in a wild bird population. Mol Ecol. 2011;20(5):1062–76.

    Article  PubMed  Google Scholar 

  23. Schoener E, Hunter S, Howe L. Development of a rapid HRM qPCR for the diagnosis of the four most prevalent Plasmodium lineages in New Zealand. Parasitol Res. 2017;116(7):1831–41.

    Article  CAS  PubMed  Google Scholar 

  24. Palinauskas V, Markovets MY, Kosarev VV, Efremov VD, Sokolov LV, Valkiûnas G. Occurrence of avian haematozoa in Ekaterinburg and Irkutsk districts of Russia. Ekologija. 2005;4:8–12.

    Google Scholar 

  25. Palinauskas V, Valkiūnas G, Bolshakov CV, Bensch S. Plasmodium relictum (lineage P-SGS1): effects on experimentally infected passerine birds. Exp Parasitol. 2008;120(4):372–80.

    Article  PubMed  Google Scholar 

  26. Palinauskas V, Valkiūnas G, Križanauskienė A, Bensch S, Bolshakov CV. Plasmodium relictum (lineage P-SGS1): further observation of effects on experimentally infected passeriform birds, with remarks on treatment with Malarone™. Exp Parasitol. 2009;123(2):134–9.

    Article  CAS  PubMed  Google Scholar 

  27. Palinauskas V, Valkiūnas G, Bolshakov CV, Bensch S. Plasmodium relictum (lineage SGS1) and Plasmodium ashfordi (lineage GRW2): the effects of the co-infection on experimentally infected passerine birds. Exp Parasitol. 2011;127(2):527–33.

    Article  PubMed  Google Scholar 

  28. Atkinson CT, Thomas NJ, Hunter DB. Parasitic diseases of wild birds. John Wiley & Sons; 2009.

    Google Scholar 

  29. Valkiūnas G, Palinauskas V, Ilgūnas M, Bukauskaitė D, Dimitrov D, Bernotienė R, Zehtindjiev P, Ilieva M, Iezhova TA. Molecular characterization of five widespread avian haemosporidian parasites (Haemosporida), with perspectives on the PCR-based detection of haemosporidians in wildlife. Parasitol Res. 2014;113(6):2251–63.

    Article  PubMed  Google Scholar 

  30. Braga ÉM, Silveira P, Belo NO, Valkiūnas G. Recent advances in the study of avian malaria: an overview with an emphasis on the distribution of Plasmodium spp in Brazil. Mem Inst Oswaldo Cruz. 2011;106:3–11.

    Article  PubMed  Google Scholar 

  31. Peirce M. The significance of avian haematozoa in conservation strategies. In: Disease and threatened birds Proceedings of a Symposium held at the XIX World Conference of the International Council for Bird Preservation, June 1986, Queens University, Kingston, Ontario, Canada: 1989. International Council for Bird Preservation: 69–76.

  32. Alley M, Fairley R, Martin D, Howe L, Atkinson T. An outbreak of avian malaria in captive yellowheads/mohua (Mohoua ochrocephala). N Z Vet J. 2008;56(5):247–51.

    Article  CAS  PubMed  Google Scholar 

  33. Bell JA, Weckstein JD, Fecchio A, Tkach VV. A new real-time PCR protocol for detection of avian haemosporidians. Parasit Vectors. 2015;8(1):383.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Brossy J. Malaria in wild and captive Jackass Penguins Spheniscus demersus along the southern African coast. Ostrich. 1992;63(1):10–2.

    Article  Google Scholar 

  35. Motta ROC, Marques MVR, Junior FCF, de Assis Andery D, Horta RS, Peixoto RB, Lacorte GA, de Abreu Moreira P, Leme FdOP, Melo MM. Does haemosporidian infection affect hematological and biochemical profiles of the endangered Black-fronted piping-guan (Aburria jacutinga)? PeerJ. 2013;1:e45.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Castle MD, Christensen BM. Hematozoa of wild turkeys from the midwestern United States: translocation of wild turkeys and its potential role in the introduction of Plasmodium kempi. J Wildl Dis. 1990;26(2):180–5.

    Article  CAS  PubMed  Google Scholar 

  37. Moore R, Nekaris K. Compassionate conservation, rehabilitation and translocation of Indonesian slow lorises. Endangered Species Res. 2014;26(2):93–102.

    Article  Google Scholar 

  38. Pacheco MA, Escalante AA, Garner MM, Bradley GA, Aguilar RF. Haemosporidian infection in captive masked bobwhite quail (Colinus virginianus ridgwayi), an endangered subspecies of the northern bobwhite quail. Vet Parasitol. 2011;182(2–4):113–20.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Zabransky CJ, Webb SL, Fedynich AM, Bryan AL Jr. Blood parasites in wood storks (Mycteria americana) from the southeastern United States. J Parasitol. 2008;94(5):1178–9.

    Article  PubMed  Google Scholar 

  40. Ishak HD, Dumbacher JP, Anderson NL, Keane JJ, Valkiūnas G, Haig SM, Tell LA, Sehgal RN. Blood parasites in owls with conservation implications for the spotted owl (Strix occidentalis). PLoS ONE. 2008;3(5): e2304.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Schettler E, Fickel J, Hotzel H, Sachse K, Streich WJ, Wittstatt U, Frandölich K. Newcastle disease virus and Chlamydia psittaci in free-living raptors from eastern Germany. J Wildl Dis. 2003;39(1):57–63.

    Article  PubMed  Google Scholar 

  42. Forbes NA. Raptors: parasitic disease. In: BSAVA Manual of Raptors, Pigeons and Passerine Birds. BSAVA Library; 2008: 202–211.

  43. Joseph V. Aspergillosis in raptors. In: Seminars in avian and exotic pet medicine: 2000. 66–74.

  44. Fowler ME, Schulz T, Ardans A, Reynolds B, Behymer D. Chlamydiosis in captive raptors. Avian Dis 1990:657–662.

  45. Morishita TY, Lowenstine LJ, Hirsh DC, Brooks DL. Pasteurella multocida in raptors: prevalence and characterization. Avian Dis 1996:908–918.

  46. McClure CJ, Westrip JR, Johnson JA, Schulwitz SE, Virani MZ, Davies R, Symes A, Wheatley H, Thorstrom R, Amar A. State of the world’s raptors: distributions, threats, and conservation recommendations. Biol Cons. 2018;227:390–402.

    Article  Google Scholar 

  47. Pérez-Rodríguez A, de la Puente J, Onrubia A, Pérez-Tris J. Molecular characterization of haemosporidian parasites from kites of the genus Milvus (Aves: Accipitridae). Int J Parasitol. 2013;43(5):381–7.

    Article  PubMed  Google Scholar 

  48. Sehgal RN, Hull AC, Anderson NL, Valkiunas G, Markovets MJ, Kawamura S, Tell LA. Evidence for cryptic speciation of Leucocytozoon spp. (Haemosporida, Leucocytozoidae) in diurnal raptors. J Parasitol. 2006;92(2):375–9.

    Article  PubMed  Google Scholar 

  49. Ishak HD, Loiseau C, Hull AC, Sehgal RN. Prevalence of blood parasites in migrating and wintering California hawks. J Raptor Res. 2010;44(3):215–23.

    Article  Google Scholar 

  50. Gutierrez-Lopez R, Gangoso L, Martinez-de la Puente J, Fric J, Lopez-Lopez P, Mailleux M, Munoz J, Touati L, Samraoui B, Figuerola J. Low prevalence of blood parasites in a long-distance migratory raptor: the importance of host habitat. Parasit Vectors. 2015;8:189.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Walther E, Valkiūnas G, Wommack EA, Bowie RC, Iezhova TA, Sehgal RN. Description and molecular characterization of a new Leucocytozoon parasite (Haemosporida: Leucocytozoidae), Leucocytozoon californicus sp. nov., found in American kestrels (Falco sparverius sparverius). Parasitol Res. 2016;115(5):1853–62.

    Article  PubMed  Google Scholar 

  52. Nourani L, Aliabadian M, Amini H, Mirshamsi O. Prevalence of haemosporidian parasite infections in raptors of Iran. J Raptor Res 2020; 54.

  53. Nourani L, Aliabadian M, Dinparast-Djadid N, Mirshamsi O. New Host Records for Haemoproteus spp. (Apicomplexa: Haemosporidiasina) in Passeriformes from North-West of Iran. J Arthropod-borne Dis. 2017;11(2):236.

    PubMed  PubMed Central  Google Scholar 

  54. Nourani L, Aliabadian M, Djadid ND, Mirshamsi O. Occurrence of Haemoproteus spp. (Haemosporida: Haemoproteidae) in New Host Records of Passerine Birds from the East of Iran. Iran J Parasitol. 2018;13(2):267–74.

    PubMed  PubMed Central  Google Scholar 

  55. Nourani L, Aliabadian M, Mirshamsi O, Djadid ND. Detection of the avian blood parasites from northeast of Iran. In: 1st International Congress on Vector-Borne Diseases and Climate Change: 2017.

  56. Nourani L, Aliabadian M, Mirshamsi O, Djadid ND. Molecular detection and genetic diversity of avian haemosporidian parasites in Iran. PLoS ONE. 2018;13(11): e0206638.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Nourani L, Djadid ND, Rabiee K, Mezerji MS, Shakiba M, Bakhshi H, Shokrollahi B, Farahani RK. Detection of haemosporidian parasites in wild and domestic birds in northern and central provinces of Iran: introduction of new lineages and hosts. Int J Parasitol 2020.

  58. Pornpanom P, Chagas CRF, Lertwatcharasarakul P, Kasorndorkbua C, Valkiūnas G, Salakij C. Molecular prevalence and phylogenetic relationship of Haemoproteus and Plasmodium parasites of owls in Thailand: data from a rehabilitation centre. Int J Parasitol 2019.

  59. Kaboli M, Aliabadian M, Tohidifar M, Hashemi A, Musavi SB, Roselaar CC. Atlas of birds of Iran: Jahad Daneshgahi, Karazmi Branch; 2016.

  60. Seutin G, White BN, Boag PT. Preservation of avian blood and tissue samples for DNA analyses. Can J Zool. 1991;69(1):82–90.

    Article  CAS  Google Scholar 

  61. Hellgren O, Waldenström J, Bensch S. A new PCR assay for simultaneous studies of Leucocytozoon, Plasmodium, and Haemoproteus from avian blood. J Parasitol. 2004;90(4):797–802.

    Article  CAS  PubMed  Google Scholar 

  62. Chagas CRF, Valkiūnas G, de Oliveira GL, Monteiro EF, Guida FJV, Simões RF, Rodrigues PT, de Albuquerque Luna EJ, Kirchgatter K. Diversity and distribution of avian malaria and related haemosporidian parasites in captive birds from a Brazilian megalopolis. Malar J. 2017;16(1):83.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In: Nucleic acids symposium series: 1999. [London]: Information Retrieval Ltd., c1979-c2000.: 95–98.

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Aliabadian M, Nijman V, Mahmoudi A, Naderi M, Vonk R, Vences M. ExcaliBAR: a simple and fast software utility to calculate intra-and interspecific distances from DNA barcodes. Contrib Zool. 2014;83(1):79–84d.

    Article  Google Scholar 

  67. Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol Ecol. 2012;21(8):1864–77.

    Article  CAS  PubMed  Google Scholar 

  68. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572–4.

    Article  CAS  PubMed  Google Scholar 

  69. Posada D, Crandall KA. Modeltest: testing the model of DNA substitution. Bioinformatics (Oxford, England). 1998;14(9):817–8.

    Article  CAS  Google Scholar 

  70. Rambaut A. FigTree v1. 4. Molecular evolution, phylogenetics and epidemiology. Edinburgh: University of Edinburgh, Institute of Evolutionary Biology 2012.

  71. Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29(22):2869–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Belo N, Passos L, Júnior L, Goulart C, Sherlock T, Braga E. Avian malaria in captive psittacine birds: detection by microscopy and 18S rRNA gene amplification. Prev Vet Med. 2009;88(3):220–4.

    Article  CAS  PubMed  Google Scholar 

  73. Scaglione F, Cannizzo F, Chiappino L, Sereno A, Ripepi M, Salamida S, Manuali E, Bollo E. Plasmodium spp. In a captive raptor collection of a safaripark in northwest Italy. Res Vet Sci. 2016;104:123–5.

    Article  CAS  PubMed  Google Scholar 

  74. Stauber E. The value of wildlife rehabilitation-opportunities for medical training, research, education, conservation. Jpn J Zoo Wildl Med. 2002;7(1):1–4.

    Article  Google Scholar 

  75. Deem SL, Karesh WB, Weisman W. Putting theory into practice: wildlife health in conservation. Conserv Biol. 2001;15(5):1224–33.

    Article  Google Scholar 

  76. Deem SL, Ladwig E, Cray C, Karesh WB, Amato G. Health assessment of the ex situ population of St Vincent parrots (Amazona guildingii) in St Vincent and the Grenadines. J Avian Med Surg. 2008;22:114–22.

    Article  PubMed  Google Scholar 

  77. Olias P, Wegelin M, Zenker W, Freter S, Gruber AD, Klopfleisch R. Avian malaria deaths in parrots, Europe. Emerg Infect Dis. 2011;17(5):950.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Ricklefs RE, Fallon SM. Diversification and host switching in avian malaria parasites. Proc Royal Soc London B Biol Sci. 2002;269(1494):885–92.

    Article  Google Scholar 

  79. Krizanauskiene A, Hellgren O, Kosarev V, Sokolov L, Bensch S, Valkiunas G. Variation in host specificity between species of avian hemosporidian parasites: evidence from parasite morphology and cytochrome B gene sequences. J Parasitol. 2006;92(6):1319–24.

    Article  CAS  PubMed  Google Scholar 

  80. Palinauskas V, Iezhova TA, Krizanauskiene A, Markovets MY, Bensch S, Valkiunas G. Molecular characterization and distribution of Haemoproteus minutus (Haemosporida, Haemoproteidae): a pathogenic avian parasite. Parasitol Int. 2013;62(4):358–63.

    Article  CAS  PubMed  Google Scholar 

  81. Dinhopl N, Mostegl MM, Richter B, Nedorost N, Maderner A, Fragner K, Weissenbock H. Application of in-situ hybridization for the detection and identification of avian malaria parasites in paraffin wax-embedded tissues from captive penguins. Avian Pathol. 2011;40(3):315–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Taunde PA, Bianchi MV, Perles L, da Silva FS, Guim TN, Stadler RA, André MR, Driemeier D, Pavarini SP. Pathological and molecular characterization of avian malaria in captive Magellanic penguins (Spheniscus magellanicus) in South America. Parasitol Res. 2019;118(2):599–606.

    Article  PubMed  Google Scholar 

  83. Waldenstrom J, Bensch S, Kiboi S, Hasselquist D, Ottosson U. Cross-species infection of blood parasites between resident and migratory songbirds in Africa. Mol Ecol. 2002;11(8):1545–54.

    Article  CAS  PubMed  Google Scholar 

  84. Kimura M, Darbro J, Harrington L. Avian malaria parasites share congeneric mosquito vectors. J Parasitol. 2010;96(1):144–51.

    Article  CAS  PubMed  Google Scholar 

  85. Ziman M, Colagross-Schouten A, Griffey S, Stedman B. Haemoproteus spp. and Leukocytozoon spp. in a captive raptor population. J Wildl Dis. 2004;40(1):137–40.

    Article  PubMed  Google Scholar 

  86. Rhim H, Bae J, Kim H, Han J-I. Prevalence and phylogenetic analysis of avian haemosporidia in wild birds in the republic of Korea. J Wildl Dis. 2018;54(4):772–81.

    Article  CAS  PubMed  Google Scholar 

  87. Inumaru M, Murata K, Sato Y. Prevalence of avian haemosporidia among injured wild birds in Tokyo and environs, Japan. Int J Parasitolo. 2017;6(3):299–309.

    Google Scholar 

  88. Krone O, Waldenström J, Valkiūnas G, Lessow O, Müller K, Iezhova T, Fickel J, Bensch S. Haemosporidian blood parasites in European birds of prey and owls. J Parasitol. 2008;94(3):709–15.

    Article  CAS  PubMed  Google Scholar 

  89. Ortego J, Cordero PJ, Aparicio JM, Calabuig G. No relationship between individual genetic diversity and prevalence of avian malaria in a migratory kestrel. Mol Ecol. 2007;16(22):4858–66.

    Article  CAS  PubMed  Google Scholar 

  90. Marzal A, Ricklefs RE, Valkiūnas G, Albayrak T, Arriero E, Bonneaud C, Czirják GA, Ewen J, Hellgren O, Hořáková D. Diversity, loss, and gain of malaria parasites in a globally invasive bird. PLoS ONE. 2011;6(7): e21905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Palinauskas V, Iezhova TA, Križanauskienė A, Markovets MY, Bensch S, Valkiūnas G. Molecular characterization and distribution of Haemoproteus minutus (Haemosporida, Haemoproteidae): a pathogenic avian parasite. Parasitol Int. 2013;62(4):358–63.

    Article  CAS  PubMed  Google Scholar 

  92. Lacorte GA, Felix GM, Pinheiro RR, Chaves AV, Almeida-Neto G, Neves FS, Leite LO, Santos FR, Braga EM. Exploring the diversity and distribution of neotropical avian malaria parasites–a molecular survey from Southeast Brazil. PLoS ONE. 2013;8(3): e57770.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Mantilla JS, González AD, Lotta IA, Moens M, Pacheco MA, Escalante AA, Valkiūnas G, Moncada LI, Pérez-Tris J, Matta NE. Haemoproteus erythrogravidus n. sp. (Haemosporida, Haemoproteidae): description and molecular characterization of a widespread blood parasite of birds in South America. Acta Trop. 2016;159:83–94.

    Article  PubMed  Google Scholar 

  94. Fecchio A, Wells K, Bell JA, Tkach VV, Lutz HL, Weckstein JD, Clegg SM, Clark NJ. Climate variation influences host specificity in avian malaria parasites. Ecol Lett. 2019;22(3):547–57.

    Article  PubMed  Google Scholar 

  95. Mata VA, da Silva LP, Lopes RJ, Drovetski SV. The Strait of Gibraltar poses an effective barrier to host-specialised but not to host-generalised lineages of avian Haemosporidia. Int J Parasitol. 2015;45(11):711–9.

    Article  PubMed  Google Scholar 

  96. Bernotienė R, Palinauskas V, Iezhova T, Murauskaitė D, Valkiūnas G. Avian haemosporidian parasites (Haemosporida): a comparative analysis of different polymerase chain reaction assays in detection of mixed infections. Exp Parasitol. 2016;163:31–7.

    Article  PubMed  Google Scholar 

  97. Valkiūnas G, Bensch S, Iezhova TA, Križanauskienė A, Hellgren O, Bolshakov CV. Nested cytochrome b polymerase chain reaction diagnostics underestimate mixed infections of avian blood haemosporidian parasites: microscopy is still essential. J Parasitol. 2006;92(2):418–22.

    Article  PubMed  Google Scholar 

  98. Chagas CR, Binkienė R, Ilgūnas M, Iezhova T, Valkiūnas G. The buffy coat method: a tool for detection of blood parasites without staining procedures. Parasit Vectors. 2020;13(1):1–12.

    Article  Google Scholar 

  99. Palinauskas V, Ziegyte R, Ilgunas M, Iezhova TA, Bernotiene R, Bolshakov C, Valkiunas G. Description of the first cryptic avian malaria parasite, Plasmodium homocircumflexum n. sp., with experimental data on its virulence and development in avian hosts and mosquitoes. Int J Parasitol. 2015;45(1):51–62.

    Article  PubMed  Google Scholar 

  100. Pérez-Tris J, Bensch S. Dispersal increases local transmission of avian malarial parasites. Ecol Lett. 2005;8(8):838–45.

    Article  Google Scholar 

  101. Dimitrov D, Palinauskas V, Iezhova TA, Bernotiene R, Ilgunas M, Bukauskaite D, Zehtindjiev P, Ilieva M, Shapoval AP, Bolshakov CV, et al. Plasmodium spp.: an experimental study on vertebrate host susceptibility to avian malaria. Exp Parasitol. 2015;148:1–16.

    Article  CAS  PubMed  Google Scholar 

  102. Dimitrov D, Valkiunas G, Zehtindjiev P, Ilieva M, Bensch S. Molecular characterization of haemosporidian parasites (Haemosporida) in yellow wagtail (Motacillaflava), with description of in vitro ookinetes of Haemoproteus motacillae. Zootaxa. 2013;3666:369–81.

    Article  PubMed  Google Scholar 

  103. Schmid S, Dinkel A, Mackenstedt U, Tantely ML, Randrianambinintsoa FJ, Boyer S, Woog F. Avian malaria on Madagascar: bird hosts and putative vector mosquitoes of different Plasmodium lineages. Parasit Vectors. 2017;10(1):6.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to announce their truthful gratitude to who supported in the sampling especially, Dr. Shirzoor, Mr. Khadar, Mr. Kheirabadi, and Mr. Poorghaz. We also appreciate Andrew Purches for his valuable comments in English edition. We express sincere gratitude to M.R. Nourani for his help during sampling and S. Yaghoubi for her kind support throughout our investigation.

Funding

The Office of Research Affairs, Ferdowsi University of Mashhad, Iran, funded this investigation. (L.N., Ph.D. project number 3/32568). M.A was supported by the Iranian National Science Foundation (INSF: 99025967).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: LN, MA, OM, Formal analysis: LN, MA, Software and writhing first draft of manuscript: LN, Sampling: LN, Project administration: LN, MA, and all authors contributed to writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mansour Aliabadian.

Ethics declarations

Ethics approval and consent to participate

All the experimental protocols were approved by the ethics committee for the care and use of animals for scientific purposes, Ferdowsi University of Mashhad, Iran (No: IR.MUM.FUM.REC.1397.035). All efforts were performed to reduce birds suffering during sampling. All the experiments were performed in accordance with relevant guide‑lines and regulations of ethics committee of the Ferdowsi University of Mashhad, Iran. All methods are reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments. Field sampling was authorized by the Department of the Environment of Iran, Tehran (No: 93/53552).

Consent for publication

Not applicable.

Competing interests

We have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1.

The general information of the sequences retrieved from the MalAvi dataset. This table includes the lineages’ names, genus, host family, host species, host statues, continent, country, and references.

Additional file 2.

Barcoding gap reconstruction for each parasite genus is done by ABGD analysis. The resultant histograms are shown for Haemoproteus, Leucocytozoon, and Plasmodium. These graphs illustrate the genetic distances (X) and the number ofinterspecific pairwise comparisons (Y).

Additional file 3.

More detailed information about the number sequences for each genus analysis, the number of partitions, and the frequency of groups dedicated to each partition are given.

Additional file 4.

The Poisson Tree Processes (PTP) was performed to reconstruct maximum likelihood (ML) tree analysis on the online server https://species.h-its.org/. The resultant tree showed the position of lineages at the genus level.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nourani, L., Aliabadian, M., Mirshamsi, O. et al. Prevalence of co-infection and genetic diversity of avian haemosporidian parasites in two rehabilitation facilities in Iran: implications for the conservation of captive raptors. BMC Ecol Evo 22, 114 (2022). https://doi.org/10.1186/s12862-022-02068-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-022-02068-9

Keywords