Introduction
Parasites of the genus Plasmodium, the causative agent of malaria, are vector-borne Haemosporidians that greatly affect humans and wildlife (Sachs and Malaney, Reference Sachs and Malaney2002; Valkiūnas, Reference Valkiūnas2005). The avian malaria parasite Plasmodium shows a wide range of competent hosts belonging to different bird orders and families (Fallon et al., Reference Fallon, Bermingham and Ricklefs2005; Pérez-Tris et al., Reference Pérez-Tris, Hellgren, Križanauskienė, Waldenström, Secondi, Bonneaud, Fjeldså, Hasselquist and Bensch2007; Hellgren et al., Reference Hellgren, Pérez-Tris and Bensch2009). These parasites are transmitted by some competent mosquito species, in which they undergo a sexual reproduction phase (Valkiūnas, Reference Valkiūnas2005). Infective forms of the parasite migrate to the mosquito salivary glands and may be then transmitted to a new avian host. Consequently, mosquito–Plasmodium interactions may play an important role in the dynamics of parasite transmission (Kimura et al., Reference Kimura, Darbro and Harrington2010). As stated above, not all mosquito species are competent vectors of avian malaria parasites. The main vectors are mosquitoes of the genus Culex, although other genera such as Aedes or Anopheles could also be involved in their transmission (Santiago-Alarcon et al., Reference Santiago-Alarcon, Palinauskas and Schaefer2012). The finding of genetically related parasite lineages/species in different mosquito genera leads to the assumption of a generalist relationship between Plasmodium and mosquitoes (Kimura et al., Reference Kimura, Darbro and Harrington2010; Ferraguti et al., Reference Ferraguti, Martínez-de la Puente, Muñoz, Roiz, Ruiz, Soriguer and Figuerola2013; Schoener et al., Reference Schoener, Uebleis, Butter, Nawratil, Cuk, Flechl, Kothmayer, Obwaller, Zechmeister, Rubel, Lebl, Zittra and Fuehrer2017). However, most information regarding the mosquito species involved in avian Plasmodium transmission is based on the molecular identification of parasite DNA on mosquito pools, without a quantitative evaluation of the vector competence for different mosquito species (Kimura et al., Reference Kimura, Darbro and Harrington2010; Ferraguti et al., Reference Ferraguti, Martínez-de la Puente, Muñoz, Roiz, Ruiz, Soriguer and Figuerola2013; Schoener et al., Reference Schoener, Uebleis, Butter, Nawratil, Cuk, Flechl, Kothmayer, Obwaller, Zechmeister, Rubel, Lebl, Zittra and Fuehrer2017). Indeed, molecular detection of parasite DNA in insects' bodies does not imply these are competent vectors (Valkiūnas, Reference Valkiūnas2011) and interspecific differences in competence for the transmission of avian Plasmodium could be overlooked, as parasite DNA may also be isolated from the body of non-competent mosquito species (Beerntsen et al., Reference Beerntsen, James and Christensen2000; Ishtiaq et al., Reference Ishtiaq, Guillaumot, Clegg, Phillimore, Black, Owens, Mundy and Sheldon2008). Thus, it becomes crucial to evaluate the competence of different mosquito species for the transmission of different Plasmodium lineages to better understand the transmission of avian malaria parasites in the wild.
For the successful transmission of avian Plasmodium, vectors must survive long enough to allow parasites to complete their life cycle (between 8 and 13 days) (Valkiūnas, Reference Valkiūnas2005; LaPointe et al., Reference LaPointe, Goff and Atkinson2010). The development of Plasmodium in the mosquito may be affected by several environmental factors, such as temperature and humidity (Paaijmans et al., Reference Paaijmans, Imbahale, Thomas and Takken2010; Lefévre et al., Reference Lefévre, Vantaux, Dabire, Mouline and Cohuet2013). In addition, the particular species involved in the parasite–vector assemblage and the vertebrate host parasite load may further determine the success of development of parasites in mosquitoes (Cornet et al., Reference Cornet, Nicot, Rivero and Gandon2014). For example, mosquitoes feeding on birds with high parasite loads develop a high density of ookinetes (an initial non-infective phase of Plasmodium) in their abdomen, likely increasing parasite transmission success (Pigeault et al., Reference Pigeault, Vézilier, Cornet, Zélé, Nicot, Perret, Gandon and Rivero2015). However, Plasmodium development in the mosquito produces tissue damage, with potential negative consequences for mosquito survival. Previous studies on the impact of avian malaria parasites on vector survival have reported positive, negative or non-significant effects of parasite infection on mosquito longevity (Vézilier et al., Reference Vézilier, Nicot, Gandon and Rivero2012; Lalubin et al., Reference Lalubin, Deledevant, Glaizot and Christe2014; Delhaye et al., Reference Delhaye, Aletti, Glaizot and Christe2016; Pigeault and Villa, Reference Pigeault and Villa2018; Gutiérrez-López et al., Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Yan, Soriguer and Figuerola2019a). However, most of these studies focus on the interaction between Culex pipiens mosquitoes and Plasmodium relictum (lineage SGS1) (Cornet et al., Reference Cornet, Nicot, Rivero and Gandon2013; Pigeault et al., Reference Pigeault, Vézilier, Cornet, Zélé, Nicot, Perret, Gandon and Rivero2015; Martínez-de la Puente et al., Reference Martínez-de la Puente, Gutiérrez-López and Figuerola2018). Therefore, studies considering potential differences in virulence (i.e. the cost of the pathogen infections on their host) between parasite species/lineages on different mosquito species are necessary (Lachish et al., Reference Lachish, Knowles, Alves, Wood and Sheldon2011).
Here, we experimentally assessed the competence of two mosquito species, Cx. pipiens and Aedes (Ochlerotatus) caspius, for the transmission of four avian Plasmodium lineages. Both mosquito species are common in southern Spain, where they show different feeding patterns. While Cx. pipiens feed mainly on birds (Martínez-de la Puente et al., Reference Martínez-de la Puente, Ferraguti, Ruiz, Roiz, Soriguer and Figuerola2016), Ae. caspius prefers to bite mammals, although birds may represent up to 19% of their diet (Balenghien et al., Reference Balenghien, Fouque, Sabatier and Bicout2006; Muñoz et al., Reference Muñoz, Ruiz, Soriguer, Alcaide, Viana, Roiz, Vázquez and Figuerola2012; Gutiérrez-López et al., Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Soriguer and Figuerola2019b). Avian Plasmodium DNA has been isolated from both mosquito species (Ferraguti et al., Reference Ferraguti, Martínez-de la Puente, Muñoz, Roiz, Ruiz, Soriguer and Figuerola2013; Schoener et al., Reference Schoener, Uebleis, Butter, Nawratil, Cuk, Flechl, Kothmayer, Obwaller, Zechmeister, Rubel, Lebl, Zittra and Fuehrer2017), and the capacity of Cx. pipiens for the transmission of avian Plasmodium parasites has been previously demonstrated (Kazlauskiené et al., Reference Kazlauskienė, Bernotienė, Palinauskas, Iezhova and Valkiūnas2013; Gutiérrez-López et al., Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Yan, Soriguer and Figuerola2016; Palinauskas et al., Reference Palinauskas, Žiegytė, Iezhova, Ilgūnas, Bernotienė and Valkiūnas2016). In this study, mosquitoes were allowed to feed on birds naturally infected by Plasmodium to assess the effects of bird parasite load and parasite identity (i.e. different Plasmodium lineages grouped into main clades, see below) on the probability of mosquito infection and parasite transmission. We also analysed the impact of parasite development on mosquito survival. Finally, we estimated the impact of mosquito survival and parasite development on the risk of parasite transmission, based on the quantification of the relative basic reproductive number (R 0), modified from Ross (Reference Ross1911) and Macdonald (Reference Macdonald1955).
Materials and methods
Mosquito collection and rearing
Larvae of Cx. pipiens and Ae. caspius were collected from April to September in 2014 and 2016 in the natural reserve ‘La Cañada de los Pájaros’ (6°14′W, 36°57′N, Seville Province, Spain) and in marshlands of the Huelva Province (6°53′W, 37°17′N). Larvae were grown in plastic trays with fresh or brackish water, respectively, maintained following Gutiérrez-López et al. (Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Yan, Soriguer and Figuerola2019a) with food ad libitum. Adult female mosquitoes were identified to the species level following Schaffner et al. (Reference Schaffner, Angel, Geoffroy, Hervy, Rhaiem and Brunhes2001), placed in insect cages (BugDorm-43030F, 32.5 × 32.5 × 32.5 cm) and fed ad libitum with 1% sugar solution. One day prior to each experiment, 2–3-week-old female mosquitoes were deprived from sugar solution. Only F0 generation mosquitoes collected in the field were used in the experiment. Laboratory maintained colonies of mosquitoes were not used to minimize the potential effects of artificial selection on vector–host interactions (Franks et al., Reference Franks, Pratt and Tsutsui2011; Lagisz et al., Reference Lagisz, Port and Wolff2011).
Bird sampling and experimental procedure
A total of 60 wild house sparrows (Passer domesticus) were caught from May to September 2014 in ‘La Cañada de los Pájaros’ and from June to September 2016 in different localities from Huelva Province using mist nets. Birds were individually ringed and weighed. Blood samples (0.2 mL) were obtained by jugular venepuncture using sterile syringes. A drop of blood was smeared, air-dried, fixed with absolute methanol and stained with Giemsa for 45 min (Gering and Atkinson, Reference Gering and Atkinson2004). The rest of the blood was transferred to non-heparinized Eppendorf tubes to perform molecular detection of parasites (see below). A total of 4000–10 000 erythrocytes from each smear were scanned at high magnification (×1000). Plasmodium parasite load was estimated as the percentage of infected erythrocytes. Although the gametocytaemia (proportion of red blood cells infected by gametocytes, i.e. the sexual stage of the parasite that is transmitted to mosquitoes) may provide a more reliable quantitative measure of parasite infection than parasitaemia, both variables are strongly correlated (Pigeault et al., Reference Pigeault, Vézilier, Cornet, Zélé, Nicot, Perret, Gandon and Rivero2015).
Forty-five individual birds were enclosed for 30 min in an insect cage (BugDorm-43030F, 32.5 × 32.5 × 32.5 cm) containing either Cx. pipiens (mean ± s.d.: 68.48 ± 34.97, range 4–111) or Ae. caspius (57.29 ± 34.04, range 4–126) females. Although we tried to use a similar number of mosquitoes per box, this task was difficult to achieve due to the fact that mosquitoes were obtained from larvae in the field, which limited the number of mosquitoes of similar age/species available for each day/trial. Birds were immobilized to prevent defensive behaviours following Gutierrez-López et al. (Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Soriguer and Figuerola2019b). The feeding trials were undertaken from 7:30 to 12:00 h (GMT + 1 h). Because Gutiérrez-López et al. (Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Soriguer and Figuerola2019b) did not find any effect of house sparrow sex on Cx. pipiens and Ae. caspius biting rates, we grouped birds from both sexes in this study, including 18 females and 42 males. At the end of each trial, all birds were immediately released at the place of capture, with no apparent sign of damage. In addition, 15 of these 60 birds corresponded to control birds from the study of Yan et al. (Reference Yan, Martínez-de la Puente, Gangoso, Gutiérrez-López, Soriguer and Figuerola2018). These birds were injected with saline solution, maintained in captivity for 24 days before being exposed to mosquitoes in the context of Yan et al. (Reference Yan, Martínez-de la Puente, Gangoso, Gutiérrez-López, Soriguer and Figuerola2018) study. Twenty-four hours after this exposure, we re-exposed these 15 birds to mosquitoes in the context of this study following the same procedure described above.
After trials, engorged mosquitoes (see Table 1) were placed in insect cages, one for all mosquitoes fed on the same individual bird, and maintained under the same conditions detailed above. Mosquito survival was monitored every 12 h until 13 days post-exposure (dpe). At the end of this period, the saliva of surviving mosquitoes was obtained following the protocol detailed by Gutiérrez-López et al. (Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Yan, Soriguer and Figuerola2016). We chose the isolation of saliva over other conventional methods such as the analysis of mosquito salivary glands because the former allows the use of molecular methods for parasite detection and it has been widely used in studies on the competence of mosquitoes to transmit other pathogens such as viruses (Goddard et al., Reference Goddard, Roth, Reisen and Scott2002; Dubrulle et al., Reference Dubrulle, Mousson, Moutailler, Vazeille and Failloux2009; Gutiérrez-López et al., Reference Gutiérrez-López, Bialosuknia, Ciota, Montalvo, Martínez-de la Puente, Gangoso, Figuerola and Kramer2019c) and malarial parasites (Golenda et al., Reference Golenda, Burge and Schneider1992; Gutiérrez-López et al., Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Yan, Soriguer and Figuerola2019a). Subsequently, the head-thorax of each mosquito, which contains the salivary glands, was separated from the abdomen in a sterile Petri dish. Samples were kept at −80°C until further molecular analyses.
Table 1. Cx. pipiens and Ae. caspius engorged and analysed for the different Plasmodium lineages found in house sparrows
The number of individual house sparrows infected with each Plasmodium lineage is shown.
a Three mosquitoes fed on birds infected with P. relictum escaped.
Molecular analyses
Genomic DNA was extracted from bird blood samples and the head-thorax of mosquitoes using the MAXWELL® 16 LEV Blood DNA Kit while the Qiagen DNeasy® Kit Tissue and Blood (Qiagen, Hilden, Germany) was used to extract the DNA from mosquito saliva (Gutiérrez-López et al., Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Soriguer and Figuerola2015). The presence of Plasmodium/Haemoproteus parasites was determined following Hellgren et al. (Reference Hellgren, Waldenström and Bensch2004). Positive amplifications were sequenced in both directions using the BigDye technology (Applied Biosystems) or with the Macrogen Inc. (Amsterdam, The Netherlands) sequencing service. Sequences were edited using software Sequencher™ v4.9 (Gene Codes Corp., © 1991–2009, Ann Arbor, MI 48108, USA) and assigned to parasite lineages/morphospecies after comparison with public databases [GenBank and MalAvi (Bensch et al., Reference Bensch, Hellgren and Pérez-Tris2009; Kumar et al., Reference Kumar, Stecher and Tamura2016)]. The four Plasmodium lineages found in this study were grouped into two main clades (clade A and clade B, see ‘Results’ section) and uncorrected p-distances between lineages/clades were compared using MEGA7 software (Kumar et al., Reference Kumar, Stecher and Tamura2016). Five birds showed evidence of coinfection, as revealed by the existence of double peaks in the sequencing chromatogram, and were not included in this study to avoid potential confounding effects of multiple infections on parasite development and mosquito survival (Lover and Coker, Reference Lover and Coker2015). Ten house sparrows uninfected by avian Plasmodium (see ‘Results’ section) were used as control to study the survival of 80 mosquitoes that fed on them.
Statistical analyses
Statistical analyses were performed in R software 3.2.5 (R Core Development Team, Reference R Core2016). We fitted two similar generalized linear mixed models with binomial error and logit link function using the package lme4 (Bates et al., Reference Bates, Maechler, Bolker and Walker2015), to assess the effects of Plasmodium clade identity (fixed factor) and the bird parasite load (covariate) on the status of infection by Plasmodium (infected/uninfected) in the head-thorax or saliva, respectively. The variable bird parasite load was log-transformed to attain normality. In both models, bird identity was included as a random term. We fitted, by maximum likelihood using the package survival (Therneau and Lumley, Reference Therneau and Lumley2014), a mixed-effects Cox model to test the effects of bird parasite infection status (fixed factor, infected/uninfected birds) and parasite load (covariate) on mosquito survival (measured as the number of mosquitoes alive at each 12-h-period), and bird identity as a random factor, while controlling for the potential effect of mosquito age (2 or 3 weeks old). We also fitted a similar mixed-effects Cox model with parasite identity (clade A, clade B and uninfected) instead of infection status as explanatory variable. We restricted the analyses of survival to Cx. pipiens mosquitoes, as Plasmodium only developed successfully in this species (see ‘Results’ section).
Plasmodium transmission risk
We used a simplified equation of the R 0 model proposed by Macdonald (Reference Macdonald1955) to calculate relative R 0 values:
$$R_{0\comma \;{\rm rel}} = \;\displaystyle{c \over {\lpar -\ln P\rpar }}P^v$$where c represents the probability of a mosquito becoming infected after biting an infected host, P is the daily survival rate of mosquitoes measured as the probability that a mosquito survives for 1 day and v is the pathogen incubation period in the mosquito. In our study, c (hereafter, transmission rate) was considered as the probability of a mosquito carrying Plasmodium DNA in its saliva after feeding on an infected bird. In addition, we considered v as 13 days (Valkiūnas, Reference Valkiūnas2005; LaPointe et al., Reference LaPointe, Goff and Atkinson2010). The relative R 0 value was calculated considering the survival rate until 13 dpe and the proportion of mosquitoes with positive saliva samples infected with lineages of each Plasmodium clade. The relative R 0 provides an approach to quantify the impact of differences in the survival rate of mosquitoes or Plasmodium transmission to mosquito saliva on Plasmodium transmission.
Results
Avian malaria parasites in birds
From the 60 birds captured, 10 were uninfected (and the mosquitoes that fed on them were used as controls) and 50 were infected by at least one avian malaria lineage. Five of them showed coinfection by several lineages of haemosporidians and were removed from further analyses. Overall, we exposed 45 birds infected with four parasite lineages, including the P. relictum lineages SGS1 (N = 26) and GRW11 (N = 6), and the lineages COLL1 (N = 8) and PADOM01 (N = 5) (Table 1), to mosquitoes. The morphospecies for COLL1 and PADOM01 are unknown, but these lineages clustered with the lineage SEIAUR01, corresponding to Plasmodium cathemerium, which was already found in house sparrows in the same area (Ferraguti et al., Reference Ferraguti, Martínez-de la Puente, Bensch, Roiz, Ruiz, Viana, Soriguer and Figuerola2018) (Fig. 1). The uncorrected p-distance between lineages SGS1 and GRW1 and between COLL1 and PADOM01 was 0.002 (corresponding to a difference of a single base pair). By contrast, the uncorrected p-distance between lineages SGS1–GRW11 and COLL1–PADOM01 was 0.035. Thus, in further analyses these four lineages were grouped into two different clades: clade A corresponding to the P. relictum lineages SGS1 and GRW11, and clade B corresponding to the Plasmodium spp. lineages COLL1 and PADOM01, which are closely related to P. cathemerium. The parasite load in birds did not differ significantly between clades (mean ± s.d.: clade A: 1.39 ± 0.21, clade B: 1.15 ± 0.35, ANOVA; F 1,37 = 0.34, P = 0.56). In addition, the parasite load was similar between parasite lineages within clade A (ANOVA; F 1,27 = 0.03, P = 0.87) and clade B (ANOVA; F 1,8 = 0.06, P = 0.81).
Fig. 1. Bootstrap consensus tree inferred from 10 000 replications for the Plasmodium lineages found in house sparrows.
Parasite development in mosquitoes
Overall, 28 and 17 Plasmodium-infected house sparrows were exposed to 1713 Cx. pipiens and 974 Ae. caspius, respectively. Of these, 183 (10.7%) Cx. pipiens and 112 (11.5%) Ae. caspius fed on bird blood (Table 1). One house sparrow infected by the Plasmodium lineage PADOM01 was not bitten by any Cx. pipiens. Plasmodium infection status in the head-thorax was analysed for 126 Cx. pipiens fed on 27 infected birds (19 infected by parasites of clade A and 8 infected by parasites of clade B). Fifty-one (40.5%; N = 126) Cx. pipiens were positive for Plasmodium in the head-thorax. Eleven out of these 51 mosquitoes (21.6%) had Plasmodium DNA in their saliva. For Ae. caspius, Plasmodium infection status was analysed in the head-thorax of 45 mosquitoes fed on 12 infected birds (10 infected by parasites of clade A and 2 infected by parasites of clade B). None of the Ae. caspius that fed on three birds infected by Plasmodium SGS1, one bird infected by Plasmodium COLL1 and one bird infected by Plasmodium PADOM01 survived until 13 dpe. None of the 45 head-thoraxes of Ae. caspius analysed showed evidence of Plasmodium infection (Table 1).
Parasites were detected in the head-thorax of 37.6% (N = 93) and 48.8% (N = 33) Cx. pipiens fed on birds infected by Plasmodium lineages of clades A and B, respectively. Bird parasite load positively affected the prevalence of Plasmodium in the head-thorax of Cx. pipiens [estimate (est) = 0.86, Z = 2.99, P = 0.003], while parasite prevalence did not differ between clades (est = 0.42, Z = 0.98, P = 0.33; Fig. 2). By contrast, a higher prevalence of Plasmodium clade B (21.2%) than clade A (4.3%) was found in Cx. pipiens saliva (est = 1.81, Z = 2.68, P = 0.007; Fig. 2), while non-significant associations were found with parasite load in birds' blood (est = 0.52, Z = 1.40, P = 0.16). All Plasmodium lineages infecting house sparrows were detected in mosquito saliva and the same Plasmodium lineages were found in the head-thorax and saliva of each mosquito.
Fig. 2. Percentage of Cx. pipiens head-thoraxes and saliva with presence of Plasmodium DNA from parasites of clades A (red) and B (blue). Statistically significant differences are indicated with an asterisk (*). NS means non-significant differences.
Mosquito survival
We monitored the survival up to 13 dpe of 180 Cx. pipiens fed on infected birds and 80 Cx. pipiens fed on 10 uninfected control birds. Fifty-one mosquitoes died before 13 dpe. Of them, 33.3% fed on birds infected with Plasmodium parasites of clade A, while only 8.3% fed on birds infected with parasite lineages of clade B. From those mosquitoes fed on uninfected control birds, 32.5% died before 13 dpe. The survival of mosquitoes did not depend on the bird infection status (Z = −1.42, P = 0.16), the host infection intensity (Z = 1.31, P = 0.19) or the mosquito age (Z = 0.9; P = 0.4). However, when considering the identity of Plasmodium parasites instead of the bird infection status, Cx. pipiens fed on birds infected by parasites of clade B survived longer than those fed on birds infected by clade A (Z = 2.23; P = 0.03; Fig. 3), and on un-infected control birds (Z = 2.83; P = 0.005), while parasite load in the bird host (Z = 1.39, P = 0.16) and mosquito age (Z = 0.69; P = 0.49) did not affect mosquito survival.
Fig. 3. Proportion of Cx. pipiens that survived until 13 dpe to Plasmodium parasites of clade A (red line), clade B (blue line) or control (black line). The shaded areas comprise the standard errors.
Plasmodium transmission risk
Both parameters, i.e. transmission rate and vector survival, were affected by the parasite clade and consequently both clades differed in their risk of transmission. Mosquitoes fed on birds infected by Plasmodium lineages of clade B had a higher daily survival probability (P) than those fed on birds infected by Plasmodium lineages of clade A (daily survival probability = 0.99 and 0.97, respectively). Moreover, the transmission rate (c) of Plasmodium lineages of clade B by Cx. pipiens was higher than the transmission rate of Plasmodium lineages of clade A (transmission rate = 0.21 and 0.04, respectively). Consequently, Plasmodium parasites clade A showed a 20.3 times lower relative transmission risk number than those of clade B (R 0,rel = 18.48 and 0.9 for clades B and A, respectively).
Discussion
The successful transmission of mosquito-borne parasites largely depends on the availability of competent vector species in the area (Beerntsen et al., Reference Beerntsen, James and Christensen2000). Although mosquitoes of the genus Culex are considered the main vectors of avian Plasmodium, molecular detection of different Plasmodium lineages in other mosquito genera, such as Anopheles, Aedes and Lutzia (Kimura et al., Reference Kimura, Darbro and Harrington2010; Santiago-Alarcón et al., Reference Santiago-Alarcon, Palinauskas and Schaefer2012; Ferraguti et al., Reference Ferraguti, Martínez-de la Puente, Muñoz, Roiz, Ruiz, Soriguer and Figuerola2013; Schoener et al., Reference Schoener, Uebleis, Butter, Nawratil, Cuk, Flechl, Kothmayer, Obwaller, Zechmeister, Rubel, Lebl, Zittra and Fuehrer2017), suggests that avian Plasmodium spp. are not tightly coevolved with mosquito species. However, our results indicate that while Cx. pipiens being capable of transmitting the four Plasmodium lineages we isolated from birds, Ae. caspius does not, which suggests the existence of parasite/insect related mechanisms that prevents the successful development and subsequent transmission of the parasite by some mosquito species (Hardy et al., Reference Hardy, Houk, Kramer and Reeves1983; Black and Moore, Reference Black, Moore, Beaty and Marquardt1996). The insect midgut represents a strong barrier, being able to dramatically reduce the number of viable parasites from those initially ingested (Abraham and Jacobs-Lorena, Reference Abraham and Jacobs-Lorena2004; Siden-Kiamos et al., Reference Siden-Kiamos, Ecker, Nybäck, Louis, Sinden and Billker2006). In addition, Plasmodium penetrates the midgut intracellular epithelium by a complex mechanism involving numerous proteins of the membrane (Povelones et al., Reference Povelones, Waterhouse, Kafatos and Christophides2009). Thus, potential differences in the presence of these proteins between mosquito species could explain the inability of avian Plasmodium to develop in Ae. caspius mosquitoes. Furthermore, differences in the immune response against parasites or midgut microbiota between mosquito species may affect parasite development (Azambuja et al., Reference Azambuja, Garcia and Ratcliffe2005; Weiss and Aksoy, Reference Weiss and Aksoy2011). Our results strongly support that detection of Plasmodium DNA in a particular mosquito species body is not a good indicator of the species' vector competence, since Plasmodium DNA was previously detected, even at high prevalence, in Ae. caspius (Ferraguti et al., Reference Ferraguti, Martínez-de la Puente, Muñoz, Roiz, Ruiz, Soriguer and Figuerola2013; Schoener et al., Reference Schoener, Uebleis, Butter, Nawratil, Cuk, Flechl, Kothmayer, Obwaller, Zechmeister, Rubel, Lebl, Zittra and Fuehrer2017). Although mosquito head-thorax and saliva were sampled at 13 dpe, a period that exceeds the time needed for different Plasmodium spp. to develop in the mosquito salivary glands (Valkiūnas, Reference Valkiūnas2005; LaPointe et al., Reference LaPointe, Goff and Atkinson2010), there may exist mosquito–parasite combinations that require different pre-patent periods for parasite development. This could also explain the absence of Plasmodium DNA in Ae. caspius.
Despite Cx. pipiens is a competent vector of avian Plasmodium, we found that its competence for the parasite transmission differed between the two clades. We found a reduced impact on mosquito survival and a higher parasite prevalence in mosquito saliva of clade B lineages (COLL1 and PADOM01), while for clade A (SGS1 and GRW11), we found a higher impact on mosquito survival and lower prevalence in mosquito saliva. As a result, the vector competence of Cx. pipiens for clade B parasites was much higher than for those of clade A. However, these differences could be the result of unequal parasite development in the mosquitoes. It is possible that both clades differ in the time required to develop and reach the salivary glands, as has been reported between Plasmodium parasite species (LaPointe et al., Reference LaPointe, Goff and Atkinson2010; Palinauskas et al., Reference Palinauskas, Žiegytė, Iezhova, Ilgūnas, Bernotienė and Valkiūnas2016), with parasites of clade B producing sporozoites faster than those of clade A. This could be due to different pre-patent periods of parasite development in the mosquito. In fact, Maier (Reference Maier1973) found sporozoites of P. cathemerium in the salivary glands of Cx. pipiens mosquitoes at 7 dpe, and P. relictum sporozoites in the salivary glands of mosquitoes have been recorded as early as 4 and 5 dpe (Rosen and Reeves, Reference Rosen and Reeves1954; Work et al., Reference Work, Washino and Van Riper1990). However, Kazlauskienė et al. (Reference Kazlauskienė, Bernotienė, Palinauskas, Iezhova and Valkiūnas2013) did not find sporozoites of the P. relictum lineages SGS1 and GRW11 (clade A in this study) until 14 dpe in the salivary glands of Cx. pipiens mosquitoes (presence of sporozoites was analysed at 12 dpe, but not at 13 dpe). In addition to parasite identity, the parasite load of the vertebrate host may largely determine the success of parasite development in the insect vector and, potentially, its ability for being transmitted. In humans, Plasmodium gametocytaemia was positively associated with the mosquito infection rates (Bousema and Drakeley et al., Reference Bousema and Drakeley2011). In avian Plasmodium, however, non-significant associations between host parasitaemia and oocyst prevalence have been found (Pigeault and Villa, Reference Pigeault and Villa2018). However, the interaction of bird hosts and mosquitoes could potentially affect the within-host dynamics of Plasmodium and therefore, its transmission to mosquitoes. For example, it has been shown that daily variations in mosquito activity are correlated with bird parasitaemia and that mosquito bites increase P. relictum replication in the birds both in the acute and in the chronic phases of the infection (Pigeault et al., Reference Pigeault, Caudron, Nicot, Rivero and Gandon2018). Our results suggest that although Plasmodium load in the avian host facilitates the detection of parasite DNA in the mosquito head-thorax, the successful infection and transmission of Plasmodium into mosquito saliva was not directly related to infection load in the bird host. Consequently, the final development of parasites in mosquito saliva may be modulated by other factors, including specific mosquito–parasite assemblages.
The costs of Plasmodium infection for mosquito survival remain a subject of intense debate (Ferguson and Read, Reference Ferguson and Read2002; Martínez-de la Puente et al., Reference Martínez-de la Puente, Gutiérrez-López and Figuerola2018). Vézilier et al. (Reference Vézilier, Nicot, Gandon and Rivero2012) reported a decreased fecundity and an increased longevity of mosquitoes fed on infected birds, while Pigeault and Villa (Reference Pigeault and Villa2018) did not find any association between bird parasite load and mosquito survival. However, these studies focused on the interaction between Cx. pipiens and P. relictum. In our study, we found that mosquitoes fed on birds infected by Plasmodium clade B had higher survival than those fed on birds infected by clade A (P. relictum) or un-infected birds. Differences in the level of virulence between avian Plasmodium lineages on mosquitoes are currently unknown. Nonetheless, a differential impact on bird hosts of parasite lineages/morphospecies has been reported. For instance, Lachish et al. (Reference Lachish, Knowles, Alves, Wood and Sheldon2011) found that P. relictum had a lower virulence on birds than P. circumflexum. Our results suggest that this may also occur in mosquitoes, with differential cost [i.e. energetic cost (Hurd et al., Reference Hurd, Taylor, Adams, Underhill and Eggleston2005); survival, this study] imposed by different species/lineages of Plasmodium. Although birds infected by both clades (clades A and B) showed similar Plasmodium intensities of infection, mosquitos feeding on birds with high intensities of infection die sooner (Gutiérrez-López et al., Reference Gutiérrez-López, Martínez-de la Puente, Gangoso, Yan, Soriguer and Figuerola2019a), probably due to the cost of Plasmodium infection on mosquito survival.
We found that Plasmodium transmission risk differed between parasite clades mainly due to their differential impact on mosquito survival and transmission rate (i.e. presence in saliva). Lineages of P. relictum (clade A) had a higher virulence on mosquitoes and also showed a lower transmission rate than parasites of clade B. Consequently, transmission of Plasmodium was less effective when mosquitoes fed on birds infected by lineages of clade A than when infected by clade B. In addition to the variables measured here, the epidemiology of vector-borne parasites depends on a number of factors, such as host density (Gubbins et al., Reference Gubbins, Carpenter, Baylis, Wood and Mellor2008), host recovery rate (Macdonald, Reference Macdonald1955) and vector density (Hartemink et al., Reference Hartemink, Vanwambeke, Heesterbeek, Rogers, Morley, Pesson, Davies, Mahamdallie and Ready2011). In addition, our estimation of the parasite transmission risk should be considered with caution as we monitored the survival rate of mosquitoes until 13 dpe, while epidemiological models should consider the whole lifespan of individuals. In spite of these limitations, our results provide a first step towards the identification of the consequences of parasite identity for avian malaria epidemiology, a topic that has been traditionally neglected. Interestingly, P. relictum (clade A) is considered a generalist parasite infecting more than 300 species of birds belonging to 11 different orders worldwide, being transmitted by 20 different species of mosquitoes (Valkiūnas et al., Reference Valkiūnas, Ilgūnas, Bukauskaitė, Fragner, Weissenböck, Atkinson and Iezhova2018). However, recent studies have found that generalist parasites may also show a specialized behaviour in certain, phylogenetically related host species within their host range (Svensson-Coelho et al., Reference Svensson-Coelho, Loiselle, Blake and Ricklefs2016; Huang et al., Reference Huang, Ellis, Jönsson and Bensch2018). These findings suggest that the generalist/specialist character of a parasite should not be solely based on its host range, which may in turn vary among different host and parasite communities (Svensson-Coelho et al., Reference Svensson-Coelho, Loiselle, Blake and Ricklefs2016), and highlights that phylogenetic relationships among hosts must be considered when determining the specificity of a particular parasite (Huang et al., Reference Huang, Ellis, Jönsson and Bensch2018). However, all such studies refer to the generalist/specialist continuum from a vertebrate host perspective, while little information is available on virulence and specialization on the vectors. Our results showed a higher efficacy of transmission by Cx. pipiens of parasites of clade B as compared with those of clade A. This could be due to the generalist character of P. relictum, which may decrease its fitness when transmitted to Cx. pipiens, but not when infecting birds (Hellgren et al., Reference Hellgren, Pérez-Tris and Bensch2009). Further studies are necessary to understand the range of vectors successfully transmitting the different species of avian Plasmodium.
In conclusion, results from this study confirm the existence of inter- and intra-specific differences in the ability of Plasmodium lineages to develop in different mosquito species. While some mosquitoes such as Ae. caspius were refractory to parasite development, Cx. pipiens play a key role in the transmission of avian Plasmodium by regulating, for instance, its temporal (e.g. Lalubin et al., Reference Lalubin, Deledevant, Glaizot and Christe2014) and spatial dynamics (Martínez-de la Puente et al., Reference Martínez-de la Puente, Ferraguti, Ruiz, Roiz, Soriguer and Figuerola2016) in both natural and anthropized environments. Here, we add valuable information on the competence of Cx. pipiens for the transmission of four Plasmodium lineages with active circulation in Europe. Nevertheless, the identity of each vector–parasite assemblage may modulate the transmission success of Plasmodium lineages through differences in the parasite transmission rate in the mosquito and the costs of infection on mosquito survival. Consequently, Cx. pipiens was better at transmitting lineages related to P. cathemerium than P. relictum due to a lower impact on mosquito survival and higher ability to replicate and reach mosquito saliva of the former. Understanding how Plasmodium virulence in the host and the vector interact may be crucial to understand the maintenance and abundance of the different Plasmodium species and lineages in the wild.
Acknowledgements
Alberto Pastoriza, Alazne Díez, Esmeralda Pérez and Santiago Ruiz helped with bird and mosquito capture and identification. Isabel Martín and Laura Gómez helped in the laboratory. We would also like to thank the two anonymous reviewers for their constructive comments on the manuscript.
Author contributions
All authors designed the study. RGL and JMP conducted the experimental work. RGL and JMP performed the laboratory analyses. RGL, JMP, LG and JF performed the statistical analyses. All authors read and approved the final manuscript.
Financial support
This study was funded by projects CGL2015-65055-P and PGC2018-095704-B-100 from the Spanish Ministry of Science, Innovation and Universities and from the European Regional Development Fund (FEDER). RGL was supported by a FPI grant (BES-2013-065274). JMP was partially supported by a 2017 Leonardo Grant for Researchers and Cultural Creators, BBVA Foundation. The Foundation accepts no responsibility for the opinions, statements and contents included in the project and/or the results thereof, which are entirely the responsibility of the authors. LG was supported by a Marie Sklodowska-Curie Fellowship of the European Commission (grant number 747729 ‘EcoEvoClim’). Fieldwork facilities were provided by Doñana ICTS-RBD.
Conflict of interest
The authors declare that they have no competing interests.
Ethical standards
All experiments involving birds adhered to the requirements of the Spanish Legislative Decree ‘Real Decreto 53/2013 de 1 de Febrero’ on protection of animals used for experimentation and other scientific purposes, the European Community Council Directive no. 2010/63/UE on Laboratory Animal Protection. Regional Authorities and the CSIC Ethics Committee approved this project (ref. CEBA-EBD-12-40). Mosquito sampling was performed with all the necessary permissions from landowners and regional Department of Environment (Consejería de Medio Ambiente, Junta de Andalucía).
