Skip to main content

Low MSP-1 haplotype diversity in the West Palearctic population of the avian malaria parasite Plasmodium relictum



Although avian Plasmodium species are widespread and common across the globe, limited data exist on how genetically variable their populations are. Here, the hypothesis that the avian blood parasite Plasmodium relictum exhibits very low genetic diversity in its Western Palearctic transmission area (from Morocco to Sweden in the north and Transcaucasia in the east) was tested.


The genetic diversity of Plasmodium relictum was investigated by sequencing a portion (block 14) of the fast-evolving merozoite surface protein 1 (MSP1) gene in 75 different P. relictum infections from 36 host species. Furthermore, the full-length MSP1 sequences representing the common block 14 allele was sequenced in order to investigate if additional variation could be found outside block 14.


The majority (72 of 75) of the sequenced infections shared the same MSP1 allele. This common allele has previously been found to be the dominant allele transmitted in Europe.


The results corroborate earlier findings derived from a limited dataset that the globally transmitted malaria parasite P. relictum exhibits very low genetic diversity in its Western Palearctic transmission area. This is likely the result of a recent introduction event or a selective sweep.


The amount of standing genetic variation in a pathogen population can affect its ability to adapt to changes in its host population or environment. It also reflects the pathogen’s population history: timing, trajectory and rate of past population fluctuations [1, 2]. If a genetic variant in a host–pathogen system is strongly favoured by selection, it is likely to rapidly increase in frequency, reducing the genetic diversity of that locus in the population. This, however, will have no effect on the genetic diversity of other pathogen populations where the novel variant may not exist or where selection does not favour it. Similarly, patterns of low genetic variation (across the genome) in a pathogen population might stem from an introduction event (founder effect) followed by rapid expansion in a host or vector population. While a founder effect is likely to result in low genetic diversity across the genome, selection for a particular allele would likely result in low diversity at a single locus and linked loci only. Thus, patterns of genetic variation in a pathogen population can elucidate whether selection has occurred, the rate of a new mutation becoming fixed, and whether there has been a recent introduction event into a new area or host population [3,4,5,6]. Furthermore, population genetic variation can be used to infer the geographical origin of pathogen populations [3, 5]

In contrast to the small number of malaria species infecting humans, avian Plasmodium parasites, commonly referred to as avian malaria parasites, are much more speciose [7,8,9]. Many decades ago, avian malaria was used as a model study system to better understand human malaria biology. However, the mouse malaria model soon replaced avian malaria as a model system [10, 11]. In recent years, avian malaria re-emerged as a model system for research into ecology and evolution of wildlife parasites [12]. While it previously has been difficult to sequence nuclear genes of these parasites, genome and transcriptome sequences are becoming increasingly available with the development of modern sequencing technologies [13,14,15,16,17,18]. These new data allow for investigations of intra- and interspecific genetic variation of avian malaria parasites at a higher resolution than was previously possible.

Plasmodium relictum is one of the most widespread avian malaria parasites [19,20,21] and is highly virulent in immunologically naïve avian species inhabiting remote islands and zoological gardens [22, 23]. The morphologically defined P. relictum consists of several mitochondrial cytochrome b (cyt b) lineages [19,20,21]. These lineages appear to be associated with different transmission areas. For example, the lineage SGS1 is primarily transmitted within tropical and temperate regions of the Old World whereas the lineage GRW4 occurs globally, although its transmission is limited to areas with warmer climates [21]. Furthermore, P. relictum cyt b lineages differ in evolutionary independence. SGS1 and GRW11, for example, share nuclear haplotypes in some populations likely suggesting a common recent evolutionary history, or that they belong to the same recombining population. In contrast, neither SGS1 nor GRW11 share nuclear haplotypes with GRW4, which points to the independent evolutionary history of GRW4 and SGS1/GRW11 [21].

The fast-evolving nuclear merozoite surface protein 1 (msp1) gene provides a more fine-scale genetic resolution of the parasite lineages than does the mitochondrial cyt b gene [24]. The msp1 gene is involved in the invasion of host’s red blood cells and a part of the gene codes for a protein segment that is exposed to the host immune system. This exposure is thought to lead to its relatively fast evolutionary rate likely due to selection to avoid detection. As a result, high msp1 allelic diversity has been found in natural populations of other Plasmodium species [25, 26]. However, studies of P. relictum found a remarkably low haplotype diversity of the msp1 gene in SGS1 and GRW11 lineages transmitted in Europe; nevertheless, msp1 diversity appeared higher in vectors engorged with passerine blood in Japan [21]. These findings raised the possibility that the P. relictum lineages colonized Europe from other areas, possibly Asia, where haplotype diversity of parasites appears higher, at least in vectors [21]. However, Hellgren et al. [21] analysed data from only five avian species, three of which were from the taxonomic family Passeridae. This limited host species sampling might have been the reason that low MSP1 diversity was found in Europe.

This study, significantly expands Hellgren et al. [21] dataset of P. relictum msp1 sequences by adding 27 host species sampled across the Western Palearctic in order to address: (1) whether the originally discovered low msp1 haplotype diversity of SGS1 and GRW11 in the Western Palearctic resulted from limited sampling or, alternatively, was in fact low; and, (2) if the diversity is high in the expanded dataset, whether there is an association between host species and parasite msp1 haplotypes, or if msp1 haplotypes are randomly distributed among host species. If the diversity of msp1 haplotypes is higher than originally estimated, the hypothesis is that selection imposed by different host immune systems should result in the non-random distribution of msp1 haplotypes among host species. This potential cryptic specialization of msp1 alleles may have been missed in studies using the slower evolving cyt b gene for lineage identification. Alternatively, the lack of msp1 haplotype structure among host species and its extremely low haplotype diversity would be consistent with previous work and suggest a recent population expansion of the parasites into the Western Palearctic.


Genetic diversity at the msp1 gene was analysed in two sample sets using two different sequencing methodologies. For the first, larger sample set, block 14 of the msp1 gene (249 bp out of 269 bp) was amplified and sequenced with Sanger sequencing. In order to investigate if there was variation in the msp1 gene that was missed by just sequencing block 14 (the complete msp1 gene is 4740 bp in length), the complete MSP1 gene was sequenced in a smaller sample set using a sequence capture method. The rationale behind this is that the paper by Hellgren et al. [24], where the complete msp1 gene of P. relictum was sequenced, relied on a PCR protocols, that are not nested as the block 14 protocol, using samples with extremely high parasitaemias obtained through infection experiments; such parasitaemia levels are almost never observed in the wild. Therefore, amplifying the full msp1 gene from samples obtained from natural infections using PCR and Sanger sequencing is more or less impossible. One way to circumvent this problem and to obtain full-length sequence data of nuclear genes is to use sequence capture methods [27].

Amplification and analysis of msp1 block 14

Block 14 of the msp1 gene was amplified from 27 avian blood samples infected with the P. relictum cyt b lineage GRW11 and 111 samples infected with cyt b lineage SGS1 using the primers and annealing temperatures described in Hellgren et al. [21] and sequenced with Sanger methodology (in-house). Note that the name of the internal primer pair, MSP1_3FN and MSP1_3RN, have been named incorrectly regarding their direction to the outside primers, MSP1_3F and MSP1_3R in the original article [21]. These samples were part of previous studies [28, 29], which identified single haemosporidian infections with one of the two Plasmodium cyt b lineages. All samples were collected under legal permits. These samples originated from 36 passerine species captured at 35 localities in Portugal, Morocco, Russia, and Armenia (Table 2). Blood was sampled by brachial venipuncture with a sterile needle and collected into a heparin-free glass capillary tube. Samples were immediately transferred into 2-ml vials with 96% ethanol and stored until DNA extraction as described previously [21].

New msp1 block 14 sequences (GenBank; Pr10: MZ270634 and Pr11: MZ270635) were aligned to those from Hellgren et al. [21] using Geneious Prime 2019.2.3. ( The phylogenetic relationships among new and previously discovered [21] msp1 haplotypes were evaluated by Bayesian phylogenetic analysis, using MrBayes [30] as implemented in Geneious Prime 2019.2.3. Phylogenies was constructed using 1,100,000 iterations and sampling every 200th tree under a general time reversal + invariable site model (GTR + I) allowing for six different substitution rates. After discarding 10% of the sampled trees as a burn-in, the remaining trees were used to construct a majority consensus tree. The phylogeny contained all known msp1 haplotypes originating from cyt b SGS1 and GRW11 infections. As outgroups (for outgroup sequences see [21]) msp1 haplotypes Pr4 (KJ850282) and Pr5 (KJ850281) originating from the cyt b lineage GRW4 was used. For visualization of the relationship and frequency of obtained haplotypes a haplotype network which was inferred using PopArt version 1.7.2 [31], with the TCS inference method [32] was constructed.

Sequence capture of msp1

For an additional 23 samples, collected under legal permits, with SGS1 (n = 20) or GRW11 (n = 3) infections confirmed using the molecular screening described in Hellgren et al. [33], the complete msp1 gene was sequenced following the sequence capture protocol described in Huang et al. [18], but using a probe set designed to target P. relictum instead (available upon request) after the presence of SGS1 (n = 20) or GRW11 (n = 3) infection was confirmed using the molecular screening protocol described in Hellgren et al. [33]. Capture probes targeting 1,036 genes from the P. relictum genome [15] was constructed and used in a SureSelectXT Target Enrichment kit (Agilent Technologies) as previously described [18], with a few modifications. The protocol started with 200 ng of sheared DNA from each sample and followed the manufacturer’s protocol for sequence capture, including hybridizing probes with the samples for 24 h. For the PCR, 18 cycles of was used in the post-capture amplification with indexed primers (instead of the recommended 10 to 16 cycles) to ensure sufficient amount of DNA for sequencing. Captured fragments were sequenced on an Illumina MiSeq instrument at the Lund University DNA Sequencing Facility.

Raw sequencing reads was processed by removing adapter sequences and low-quality base calls with a sliding window using Trimmomatic v0.39 [34]. Read quality was assessed with FastQC v0.11.9. The reads were mapped to the P. relictum genome (, release 39 [35]) using Nextgenmap v0.5.5 [36] with default options. The resulting sam files were compressed, sorted and indexed using samtools v1.10 [37]. MarkDuplicates v2.20.8 in Picard tools (Broad Institute, 2019) was used to identify duplicate reads. To ensure that the infections were indeed single infections and matched SGS1 or GRW11 cyt b sequences, the consensus sequence was first checked against the SGS1 and GRW11 reference sequences [9]. Then double infections were identified by visually detecting SNPs in the mapped reads of the cyt b gene. For individuals that showed no evidence of double infections at the cyt b gene, the majority consensus sequences of the msp1 gene were generated for sites with a depth of coverage > 5 as implemented in Geneious Prime 2019.2.3. Discarding double infections left 11 samples, of which five had incomplete coverage; these 11 samples were used for further analysis at different degrees. The full-length consensus sequences were first aligned against each other using the MUSCLE alignment algorithm with eight iterations in Geneious Prime 2019.2.3. This allowed the identification of SNPs outside the msp1 block 14 section (targeted in the larger sample set). Next, the full-length msp1 genes were mapped against the shorter reference sequences of msp1 block 14, found in Hellgren et al. [21], to assess the haplotype identity of the complete msp1 sequences. Those previously identified msp1haplotypes of P. relictum have been named Pr1-Pr9 [21] and will be referred to as such thereafter; novel haplotypes were given numbers > Pr9, numbered in the order in which they were discovered.


msp1 block 14

Block 14 of the msp1 gene was successfully amplified and sequenced in 11 samples of GRW11 originating from nine avian species and in 64 samples of SGS1 from 25 avian species sampled in Russia, Armenia, Morocco, and Portugal (GenBank accession numbers MZ270578-MZ270624). Haplotype Pr2 was discovered in 72 of these 75 infections. Samples that were not sequenced failed in amplification of the nuclear gene, a common problem when trying to sequence nuclear genes of avian malaria due to low parasitaemia in wild hosts [27]. A single adult common whitethroat (Sylvia communis) infected with GRW11 from Armenia had the Pr3 haplotype, which previously had been found only in vectors caught and sampled in the Eastern Palearctic, i.e., Japan [21]. The remaining two haplotypes were new and both found in SGS1 infections. Pr10 was found in an adult rock-thrush (Monticola saxatilis) from Armenia and Pr11 was found in a juvenile Cetti’s warbler (Cettia cetti) from Portugal (Table 2). The majority of the Pr2 haplotypes linked to both GRW11 and SGS1 infections were found in resident bird species confirming their local transmission in the sampling areas.

The phylogenetic reconstruction and the haplotype network placed the Pr11 haplotype together with other haplotypes transmitted either in the Western Palearctic (Pr2) or in the Eastern Palearctic (Fig. 1). In fact, the Pr11 haplotype only differed by a single synonymous substitution from Pr2, the dominant lineage transmitted in the Western Palearctic. The other new haplotype, Pr10, found in an adult Paleotropical migrant (Monticola saxatilis) was a sister to the SGS1/GRW11 clade and differed by a single non-synonymous substitution from Pr1, which has Afrotropical transmission (Fig. 1).

Fig. 1

A Haplotype network of MSP1 haplotypes found in Plasmodium relictum (cytochrome b lineages SGS1 and GRW11). Each circle represents a haplotype with the size proportional to its sampling frequency. Colours represent the area where the infected bird was sampled, not taking into account whether it was a migratory species or if it was a juvenile before its first migration or an adult. Black dots on lines represent the number of additional mutational steps between two haplotypes. B Phylogenetic reconstruction of the MSP1 haplotypes found in the SGS1 and GRW11 cyt b lineages of the Plasmodium relictum. Transmission area is determined if the haplotype have been recovered from juvenile or resident birds in an area. Scale bar represents nucleotide substitutions per site. 1 Pr2 has also been found in two cases in a non-migratory bulbul species (Pycnonotus capensis)endemic to a region in South-Africa with a Mediterranean climate. Pr10 and Pr11 are new haplotypes identified in this study. 2 Pr10 has only been found in a single adult African-migratory bird in Armenia and is therefore likely transmitted in the Afrotropical region

msp1 sequence capture

After discarding samples that showed evidence of double infections at the cyt b gene, 11 samples remained and are referred to in the subsequent analysis. The complete msp1 gene, 4740 bp, was sequenced with a coverage depth > 5 in two GRW11 and three SGS1 infected samples. For the remaining six samples, 0.2–43.4% of the nucleotide positions had a coverage depth < 5 × and were excluded from the analysis (Table 3). All five complete samples were unambiguously assigned to the haplotype Pr2 at the MSP1 block 14 and all consensus sequences where monomorphic across the rest of the gene, showing no genetic variation across these samples (represented by GenBank entry MW385258). However, all of them (GenBank accession numbers MZ270625-MZ270635) differed by a single substitution from the reference sequence (KC969175) on GenBank but were identical to the msp1 gene in the genome of the SGS1 parasite available at PlasmoDB [15, 35]. This discrepancy may point to a potential sequencing error in the reference sequence KC969175 [24].


The results confirm previous work that documented very low genetic diversity of the avian blood parasite P. relictum in its Western Palearctic transmission area [21]. The lack of diversity suggests that a single variant of the pathogen may have gone through a bottleneck or a selective sweep in the Western Palearctic. Furthermore, this extremely low diversity was found at the usually highly diverse merozoite surface protein 1 (msp1) gene. Earlier work found a single msp1 haplotype transmitted in Europe [21] in a very limited number of host species, i.e., house martins (Delichon urbicum), common crossbills (Loxia Curvirostra) and house sparrows (Passer domesticus) (Table 1). In this study the msp1 haplotypes was identified in 75 P. relictum infections (11 GRW11 and 64 SGS1) from 36 different host species. The vast majority of these infections (72 of 75) had the same msp1 haplotype (Pr2). Of the three remaining haplotypes found, −Pr11 was found in a non-migratory bird in Portugal that suggests its transmission in Europe. These results strengthen the picture that there is one major msp1 haplotype that is transmitted in Europe and neighbouring temperate regions: in this study there were 72 cases (92%) of a single haplotype (Pr2), and only a single case of a haplotype (Pr11) other than the dominant one found in individual hosts that either belonged to resident bird species or juveniles caught before they had migrated to another area (Table 2).

Table 1 Summary of MSP1 (block 14) diversity in Plasmodium relictum cyt b lineages SGS1 and GRW11 found in previous studies
Table 2 Host and sampling location of birds infected with Plasmodium relictum cyt b lineages SGS1 and GRW11 and the msp1 haplotype of the infection

The other haplotypes that were found were Pr10, a haplotype closely related to those previously found being transmitted in Afrotropical avian species [21] or identified in adult migratory birds that have been captured in Europe while returning from tropical Africa. This haplotype was found in an adult rock-thrush in Armenia. Rock-thrushes migrate to tropical Africa and it is therefore likely that this haplotype’s origin is in tropical Africa since it has not been found in juveniles or resident bird species in the Western Palearctic. If Pr10 is transmitted in Africa only as suggested (and Pr11 is transmitted in Europe only), the phylogenetic tree of the haplotypes suggests a single jump to temperate zone transmission by the ancestor of haplotypes Pr2, Pr3, Pr11 and Pr7 (Fig. 1b).

In this study the haplotype Pr3 was found in an adult common whitethroat. This is a haplotype that sits in the temperate transmission cluster in the phylogeny (Fig. 1b) with only one non-synonymous nucleotide difference from the dominant Pr2 (Fig. 1a). The haplotype was sequenced from an infection of the cyt b lineage GRW11, a lineage that is thought to be transmitted exclusively in temperate regions [21]. The haplotype Pr3 has previously only been found in vectors that have taken blood meals from birds residing in Eastern Asia (the jungle crow Corvus macrorhynchos and Asian house martin Delichon dasypus) [21]. In no case has the Pr3 been detected in tropical wintering areas of the common whitethroats which suggests that the transmission in this case had likely taken place at the breeding grounds in Armenia.

Although the diversity of host species samples used for sequence capture was not high, the lack of genetic diversity across the entire msp1 gene showed that the msp1 diversity in this and previous studies by investigating only block 14 (269 bp) of the almost 5 kb long gene was not underestimated (Table 3).

Table 3 Host species and sampling location of infections of Plasmodium relictum cyt b lineages SGS1 and GRW11 used for sequencing the full msp1 gene


Plasmodium relictum is one of the most common avian malaria parasites across the globe [9, 21] and is a morphologically defined species consisting of a clade of cyt b lineages which show marked differences in distribution and transmission areas [21]. In order to understand how this parasite has come to inhabit almost every corner of the Earth and seems to have fine-scale genetic differentiation that is linked to differences in ecological limitations [9, 21], variation in virulence [38, 39], and transmission [21, 40], there is a need to accurately describe its population genetics and the relationships among genetic lineages. For the P. relictum parasites transmitted in Europe, this study confirms that there is little genetic variation at msp1 throughout the parasite population. One possible explanation for the lack of genetic variation could be a recent introduction event into the temperate region, with reduced genetic variation as a result of a rapid population expansion. To fully understand the evolutionary history of this parasite there is a need to generate more multi-locus or genomic data, a task that has recently begun [13, 27], to document whether the low genetic diversity at msp1 is reflective of the entire genome or perhaps a result of selection. There is also a need for more extensive sampling in tropical areas in order to contrast the genetic variation found in Europe with that found in other transmission areas. Comparisons of the population genetics of multiple parasite lineages that differ in ecological traits such as transmission area or host preference will ultimately provide the opportunity to uncover the genetic basis of such traits and potentially reveal the historical circumstances that led to lineage differentiation and pathogen dispersal.

Availability of data and materials

The datasets analysed in this study are available from the corresponding author on reasonable request.


  1. 1.

    Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–92.

    CAS  Article  Google Scholar 

  2. 2.

    Heled J, Drummond AJ. Bayesian inference of population size history from multiple loci. Bmc Evol Biol. 2008;8:289.

    Article  Google Scholar 

  3. 3.

    Miao M, Yang Z, Patch H, Huang Y, Escalante AA, Cui L. Plasmodium vivax populations revisited: mitochondrial genomes of temperate strains in Asia suggest ancient population expansion. BMC Evol Biol. 2012;12:22.

    Article  Google Scholar 

  4. 4.

    Rachowicz LJ, Hero J, Alford RA, Taylor JW, Morgan JAT, Vredenburg VT, et al. The novel and endemic pathogen hypotheses: competing explanations for the origin of emerging infectious diseases of wildlife. Conserv Biol. 2005;19:1441–8.

    Article  Google Scholar 

  5. 5.

    Moya A, Holmes EC, González-Candelas F. The population genetics and evolutionary epidemiology of RNA viruses. Nat Rev Microbiol. 2004;2:279–88.

    CAS  Article  Google Scholar 

  6. 6.

    Morgan JAT, Vredenburg VT, Rachowicz LJ, Knapp RA, Stice MJ, Tunstall T, et al. Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis. Proc Natl Acad Sci USA. 2007;104:13845–50.

    CAS  Article  Google Scholar 

  7. 7.

    Bensch S, Perez-Tris J, Waldenstrom J, Hellgren O. Linkage between nuclear and mitochondrial DNA sequences in avian malaria parasites: multiple cases of cryptic speciation? Evolution. 2004;58:1617–21.

    CAS  Article  Google Scholar 

  8. 8.

    Valkiunas G. Avian malaria parasites and other haemosporidia. New York: CRC; 2005.

    Google Scholar 

  9. 9.

    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:1353–8.

    Article  Google Scholar 

  10. 10.

    Stephens R, Culleton RL, Lamb TJ. The contribution of Plasmodium chabaudi to our understanding of malaria. Trends Parasitol. 2012;28:73–82.

    Article  Google Scholar 

  11. 11.

    Rivero A, Gandon S. Evolutionary ecology of avian malaria: past to present. Trends Parasitol. 2018;34:712–26.

    Article  Google Scholar 

  12. 12.

    Pigeault R, Vézilier J, Cornet S, Zélé F, Nicot A, Perret P, et al. Avian malaria: a new lease of life for an old experimental model to study the evolutionary ecology of Plasmodium. Philos Trans R Soc Lond B Biol Sci. 2015;370:20140300.

    Article  Google Scholar 

  13. 13.

    Videvall E. Genomic advances in avian malaria research. Trends Parasitol. 2019;35:254–66.

    CAS  Article  Google Scholar 

  14. 14.

    Videvall E, Cornwallis CK, Ahren D, Palinauskas V, Valkiūnas G, Hellgren O. The transcriptome of the avian malaria parasite Plasmodium ashfordi displays host-specific gene expression. Mol Ecol. 2017;26:2939–58.

    CAS  Article  Google Scholar 

  15. 15.

    Böhme U, Otto TD, Cotton JA, Steinbiss S, Sanders M, Oyola SO, et al. Complete avian malaria parasite genomes reveal features associated with lineage-specific evolution in birds and mammals. Genome Res. 2018;28:547–60.

    Article  Google Scholar 

  16. 16.

    Garcia-Longoria L, Palinauskas V, Ilgūnas M, Valkiūnas G, Hellgren O. Differential gene expression of Plasmodium homocircumflexum (lineage pCOLL4) across two experimentally infected passerine bird species. Genomics. 2020;112:2857–65.

    CAS  Article  Google Scholar 

  17. 17.

    Barrow LN, Allen JM, Huang X, Bensch S, Witt CC. Genomic sequence capture of haemosporidian parasites: methods and prospects for enhanced study of host-parasite evolution. Mol Ecol Resour. 2018;19:400–10.

    Article  Google Scholar 

  18. 18.

    Huang X, Hansson R, Palinauskas V, Valkiūnas G, Hellgren O, Bensch S. The success of sequence capture in relation to phylogenetic distance from a reference genome: a case study of avian haemosporidian parasites. Int J Parasitol. 2018;48:947–54.

    CAS  Article  Google Scholar 

  19. 19.

    Valkiūnas G, Ilgunas M, Bukauskaite D, Fragner K, Weissenböck H, Atkinson CT, et al. Characterization of Plasmodium relictum, a cosmopolitan agent of avian malaria. Malar J. 2018;17:184.

    Article  Google Scholar 

  20. 20.

    Beadell JS, Ishtiaq F, Covas R, Melo M, Warren BH, Atkinson CT, et al. Global phylogeographic limits of Hawaii’s avian malaria. Proc Biol Sci. 2006;273:2935–44.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Hellgren O, Atkinson CT, Bensch S, Albayrak T, Dimitrov D, Ewen JG, et al. Global phylogeography of the avian malaria pathogen Plasmodium relictum based on MSP1 allelic diversity. Ecography. 2015;38:842–50.

    Article  Google Scholar 

  22. 22.

    Vanriper C, Vanriper SG, Goff ML, Laird M. The epizootiology and ecological significance of malaria in Hawaiian land birds. Ecol Monogr. 1986;56:327–44.

    Article  Google Scholar 

  23. 23.

    Bueno MG, Lopez RPG, Menezes RMT, Costa-Nascimento M, Lima G, Araújo RA, et al. Identification of Plasmodium relictum causing mortality in penguins (Spheniscus magellanicus) from São Paulo Zoo Brazil. Vet Parasitol. 2010;173:123–7.

    Article  Google Scholar 

  24. 24.

    Hellgren O, Kutzer M, Bensch S, Valkiūnas G, Palinauskas V. Identification and characterization of the merozoite surface protein 1 (msp1) gene in a host-generalist avian malaria parasite, Plasmodium relictum (lineages SGS1 and GRW4) with the use of blood transcriptome. Malar J. 2013;12:381.

    Article  Google Scholar 

  25. 25.

    Kadekoppala M, Holder AA. Merozoite surface proteins of the malaria parasite: the MSP1 complex and the MSP7 family. Int J Parasitol. 2010;40:1155–61.

    CAS  Article  Google Scholar 

  26. 26.

    Tanabe K, Sakihama N, Rooth I, Björkman A, Färnert A. High frequency of recombination-driven allelic diversity and temporal variation of Plasmodium falciparum msp1 in Tanzania. Am J Trop Med Hyg. 2007;76:1037–45.

    CAS  Article  Google Scholar 

  27. 27.

    Bensch S, Hellgren O. The Use of Molecular Methods in Studies of Avian Haemosporidians. In: Santiago-Alarcon D, Marzal R, Alfonso AG, editors. Avian malaria in the tropics. In avian malaria and related parasites in the tropics, ecology, evolution and systematics. Switzerland: Springer; 2020. p. 113–35.

    Chapter  Google Scholar 

  28. 28.

    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:711–9.

    Article  Google Scholar 

  29. 29.

    Drovetski SV, Aghayan SA, Mata VA, Lopes RJ, Mode NA, Harvey JA, et al. Does the niche breadth or trade-off hypothesis explain the abundance–occupancy relationship in avian Haemosporidia? Mol Ecol. 2014;23:3322–9.

    Article  Google Scholar 

  30. 30.

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

    CAS  Article  Google Scholar 

  31. 31.

    Leigh JW, Bryant D. popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.

    Article  Google Scholar 

  32. 32.

    Clement M, Snell Q, Walke P, Posada D, Crandall K. TCS: estimating gene genealogies. In: Proceedings of 16th international parallel distributed process symposium. 2002, p 7.

  33. 33.

    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:797–802.

    CAS  Article  Google Scholar 

  34. 34.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    CAS  Article  Google Scholar 

  35. 35.

    Aurrecoechea C, Brestelli J, Brunk BP, Dommer J, Fischer S, Gajria B, et al. PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res. 2009;37:D539–43.

    CAS  Article  Google Scholar 

  36. 36.

    Sedlazeck FJ, Rescheneder P, von Haeseler A. NextGenMap: fast and accurate read mapping in highly polymorphic genomes. Bioinformatics. 2013;29:2790–1.

    CAS  Article  Google Scholar 

  37. 37.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    Article  Google Scholar 

  38. 38.

    Videvall E, Palinauskas V, Valkiūnas G, Hellgren O. Host transcriptional responses to high- and low-virulent avian malaria parasites. Am Nat. 2020;195:1070–84.

    Article  Google Scholar 

  39. 39.

    Dimitrov D, Palinauskas V, Iezhova TA, Bernotienė R, Ilgunas M, Bukauskaite D, et al. Plasmodium spp.: an experimental study on vertebrate host susceptibility to avian malaria. Exp Parasitol. 2015;148:1–16.

    CAS  Article  Google Scholar 

  40. 40.

    Garcia-Longoria L, Hellgren O, Bensch S, Lope FD, Marzal A. Detecting transmission areas of malaria parasites in a migratory bird species. Parasitology. 2015;142:1215–20.

    CAS  Article  Google Scholar 

Download references


Tomas Johansson sheared DNA samples, ran quality assessments, and performed sequencing on the sequence capture samples. Staffan Bensch and Tomas Johansson gave critical advice on the sequence capture protocol. Jane Jönsson screened samples and provided invaluable assistance in all aspects of the laboratory work. Paulius Rapševičius performed sequence capture laboratory work with VAE. We further thank Vaidas Palinauskas and Helena Westerdahl for kindly providing samples for sequence capture. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the US Government.


Open access funding provided by Lund University. The study was funded by the Swedish research foundation to OH (VR 2016–03419). VAE was supported by a postdoctoral fellowship from the Carl Tryggers Foundation. This study was partially supported by FEDER funds through the COMPETE programme, POPH/ QREN/FSE funds to SVD, by the Fundação para a Ciência e Tecnologia/MEC to SVD (FCOMP-01–0124-FEDER-008941; PTDC/BIA-BEC/103435/2008) and by the National Geographic Society to SVD. RJL was supported by Fundação para a Ciência e Tecnologia (FCT)—Transitory Norm contract DL57/2016/CP1440/CT0006.

Author information




OH designed the study and conducted molecular screening of MSP1. OH and SVD interpreted the results and wrote the first draft. VK was responsible for the bioinformatical work, VAE together with AC, MD and XH developed and conducted all the sequence capture laboratory work. SVD, RJL, VAM were responsible for the field sampling of the majority of the used samples. VAM and SAA conducted initial screening of large sample sets to identify those with single SGS1 or GRW11 infections from Portugal/Morocco and Armenia/Russia, respectively. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Olof Hellgren.

Ethics declarations

Ethics approval and consent to participate

All samples were collected under legal permits.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Hellgren, O., Kelbskopf, V., Ellis, V.A. et al. Low MSP-1 haplotype diversity in the West Palearctic population of the avian malaria parasite Plasmodium relictum. Malar J 20, 265 (2021).

Download citation


  • Avian malaria
  • Haemosporida
  • Host specificity
  • Merozoite surface protein 1
  • Parasite population structure
  • Plasmodium relictum