Tree-based delimitation of morphologically ambiguous taxa: a study of the lizard malaria parasites on the Caribbean island of Hispaniola

A new long-read mitochondrial-genome protocol (PacBio HiFi) for haemosporidian parasites: a tool for population and biodiversity studies

BACKGROUND

Studies on haemosporidian diversity, including origin of human malaria parasites, malaria's zoonotic dynamic, and regional biodiversity patterns, have used target gene approaches. However, current methods have a trade-off between scalability and data quality. Here, a long-read Next-Generation Sequencing protocol using PacBio HiFi is presented. The data processing is supported by a pipeline that uses machine-learning for analysing the reads.

METHODS

A set of primers was designed to target approximately 6 kb, almost the entire length of the haemosporidian mitochondrial genome. Amplicons from different samples were multiplexed in an SMRTbell® library preparation. A pipeline (HmtG-PacBio Pipeline) to process the reads is also provided; it integrates multiple sequence alignments, a machine-learning algorithm that uses modified variational autoencoders, and a clustering method to identify the mitochondrial haplotypes/species in a sample. Although 192 specimens could be studied simultaneously, a pilot experiment with 15 specimens is presented, including in silico experiments where multiple data combinations were tested.

RESULTS

The primers amplified various haemosporidian parasite genomes and yielded high-quality mt genome sequences. This new protocol allowed the detection and characterization of mixed infections and co-infections in the samples. The machine-learning approach converged into reproducible haplotypes with a low error rate, averaging 0.2% per read (minimum of 0.03% and maximum of 0.46%). The minimum recommended coverage per haplotype is 30X based on the detected error rates. The pipeline facilitates inspecting the data, including a local blast against a file of provided mitochondrial sequences that the researcher can customize.

CONCLUSIONS

This is not a diagnostic approach but a high-throughput method to study haemosporidian sequence assemblages and perform genotyping by targeting the mitochondrial genome. Accordingly, the methodology allowed for examining specimens with multiple infections and co-infections of different haemosporidian parasites. The pipeline enables data quality assessment and comparison of the haplotypes obtained to those from previous studies. Although a single locus approach, whole mitochondrial data provide high-quality information to characterize species pools of haemosporidian parasites.

Origin and diversity of malaria parasites and other Haemosporida

Symbionts, including parasites, are ubiquitous in all world ecosystems. Understanding the diversity of symbiont species addresses diverse questions, from the origin of infectious diseases to inferring processes shaping regional biotas. Here, we review the current approaches to studying Haemosporida's species diversity and evolutionary history. Despite the solid knowledge of species linked to diseases, such as the agents of human malaria, studies on haemosporidian phylogeny, diversity, ecology, and evolution are still limited. The available data, however, indicate that Haemosporida is an extraordinarily diverse and cosmopolitan clade of symbionts. Furthermore, this clade seems to have originated with their vertebrate hosts, particularly birds, as part of complex community level processes that we are still characterizing.

Plasmodium ouropretensis, n. sp., a new case of non-erythrocytic species within lizard malaria parasites

Delimiting and describing Plasmodium species in reptiles remains a pressing problem in Haemosporida taxonomy. The few morphological characters used can overlap, and the significance of some life-history traits is not fully understood. Morphologically identical lizard Plasmodium forms have been reported infecting different cell types (red and white blood cells) in the same host and have been considered the same species. An example is Plasmodium tropiduri tropiduri, a species known to infect erythrocytes, thrombocytes and lymphocyte-like cells. Here, both forms of P. t. tropiduri were analysed using light microscope-based morphological characteristics and phylogenetic inferences based on almost complete mitochondrial genomes of parasites naturally infecting lizards in southeastern Brazil. Although morphologically similar, two distinct phylogenetic lineages infecting erythrocytes and non-erythrocytic cells were found. The lineage found in the erythrocytes forms a monophyletic group with species from Colombia. However, the non-erythrocytic lineage shares a recent common ancestor with Plasmodium leucocytica, which infects leucocytes in lizards from the Caribbean islands. Here, Plasmodium ouropretensis n. sp. is described as a species that infects thrombocytes and lymphocyte-like cells.

Molecular diversity and coalescent species delimitation of avian haemosporidian parasites in an endemic bird species of South America

Haemoproteus spp. and Plasmodium spp. are blood parasites that occur in birds worldwide. Identifying the species within this group is complex, especially in wild birds that present low parasitemia when captured, making morphological identification very difficult. Thus, the use of alternative tools to identify species may be useful in the elucidation of the distribution of parasites that circulate in bird populations. The objectives of this study were to determine the prevalence and parasitemia of the genera Plasmodium and Haemoproteus in Tachyphonus coronatus in the Atlantic Forest, Brazil, and to evaluate the molecular diversity, geographic distribution, and specificity of these parasites based on coalescent species delimitation methods. Microscopic analysis, PCR, cyt b gene sequencing, phylogenetic analysis and coalescent species delimitation using single-locus algorithms were performed (Poisson tree process (PTP) and multi-rate Poisson tree process (MPTP) methods). The analyses were performed in 117 avian host individuals. The prevalence was 55.5% for Plasmodium and 1.7% for Haemoproteus, with a mean parasitemia of 0.06%. Twenty-five Plasmodium and two Haemoproteus lineages were recovered. The MPTP method recovered seven different evolutionarily significant units (ESUs) of Plasmodium and one of Haemoproteus, whereas PTP presented fourteen ESUs of Plasmodium and one of Haemoproteus. The MPTP was more consistent with current taxonomy, while PTP overestimated the number of lineages. These ESUs are widely distributed and have already been found in 22 orders of birds that, all together, inhabit every continent, except Antarctica. The computational methods of species delimitation proved to be effective in cases where the classification of Haemosporida based just on morphology is insufficient.

Integrating coalescent species delimitation with analysis of host specificity reveals extensive cryptic diversity despite minimal mitochondrial divergence in the malaria parasite genus Leucocytozoon

BACKGROUND

Coalescent methods that use multi-locus sequence data are powerful tools for identifying putatively reproductively isolated lineages, though this approach has rarely been used for the study of microbial groups that are likely to harbor many unrecognized species. Among microbial symbionts, integrating genetic species delimitation methods with trait data that could indicate reproductive isolation, such as host specificity data, has rarely been used despite its potential to inform species limits. Here we test the ability of an integrative approach combining genetic and host specificity data to delimit species within the avian malaria parasite genus Leucocytozoon in central Alaska.

RESULTS

We sequenced seven nuclear loci for 69 Leucocytozoon samples and used multiple species delimitation methods (GMYC and BPP models), tested for differences in host infection patterns among putative species based on 406 individual infections, and characterized parasite morphology. We found that cryptic morphology has masked a highly diverse Leucocytozoon assemblage, with most species delimitation methods recovering support for at least 21 separate species that occur sympatrically and have divergent host infection patterns. Reproductive isolation among putative species appears to have evolved despite low mtDNA divergence, and in one instance two Leucocytozoon cytb haplotypes that differed by a single base pair (~ 0.2% divergence) were supported as separate species. However, there was no consistent association between mtDNA divergence and species limits. Among cytb haplotypes that differed by one to three base pairs we observed idiosyncratic patterns of nuclear and ecological divergence, with cytb haplotype pairs found to be either conspecific, reproductively isolated with no divergence in host specificity, or reproductively isolated with divergent patterns of host specialization.

CONCLUSION

Integrating multi-locus genetic species delimitation methods and non-traditional ecological data types such as host specificity provide a novel view of the diversity of avian malaria parasites that has been missed previously using morphology and mtDNA barcodes. Species delimitation methods show that Leucocytozoon is highly species-rich in Alaska, and the genus is likely to harbor extraordinary species-level diversity worldwide. Integrating genetic and ecological data will be an important approach for understanding the diversity and evolutionary history of microbial symbionts moving forward.

Primers targeting mitochondrial genes of avian haemosporidians: PCR detection and differential DNA amplification of parasites belonging to different genera

Haemosporida is a diverse group of vector-borne parasitic protozoa, ubiquitous in terrestrial vertebrates worldwide. The renewed interest in their diversity has been driven by the extensive use of molecular methods targeting mitochondrial genes. Unfortunately, most studies target a 478 bp fragment of the cytochrome b (cytb) gene, which often cannot be used to separate lineages from different genera found in mixed infections that are common in wildlife. In this investigation, an alignment constructed with 114 mitochondrial genome sequences belonging to four genera (Leucocytozoon, Haemoproteus, Plasmodium and Hepatocystis) was used to design two different sets of primers targeting the cytb gene as well as the other two mitochondrial DNA genes: cytochrome c oxidase subunit 1 and cytochrome c oxidase subunit 3. The design of each pair of primers required consideration of different criteria, including a set for detection and another for differential amplification of DNA from parasites belonging to different avian haemosporidians. All pairs of primers were tested in three laboratories to assess their sensitivity and specificity under diverse practices and across isolates from different genera including single and natural mixed infections as well as experimental mixed infections. Overall, these primers exhibited high sensitivity regardless of the differences in laboratory practices, parasite species, and parasitemias. Furthermore, those primers designed to separate parasite genera showed high specificity, as confirmed by sequencing. In the case of cytb, a nested multiplex (single tube PCR) test was designed and successfully tested to differentially detect lineages of Plasmodium and Haemoproteus parasites by yielding amplicons with different sizes detectable in a standard agarose gel. To our knowledge, the designed assay is the first test for detection and differentiation of species belonging to these two genera in a single PCR. The experiments across laboratories provided recommendations that can be of use to those researchers seeking to standardise these or other primers to the specific needs of their field investigations.

Plasmodium parasites in reptiles from the Colombia Orinoco-Amazon basin: a re-description of Plasmodium kentropyxi Lainson R, Landau I, Paperna I, 2001 and Plasmodium carmelinoi Lainson R, Franco CM, da Matta R, 2010

Colombia is a megadiverse country with about 600 species of reptiles; however, there are few studies on species of hemoparasites found in this taxonomic group. Here, we document the presence of Plasmodium spp. in four species of reptiles from the northern part of the Orinoco-Amazon region in Colombia. Individuals analyzed in this study were captured in localities between 200 and 500 m altitude, in the department of Guaviare. Each sample was screened for haemosporidian parasites by using morphology and a nested polymerase chain reaction (PCR) protocol that targets the mitochondrial cytochrome b (cytb) gene. Four morphotypes of the genus Plasmodium were found; two of these species are re-described using morphological and molecular data (cytb). For the other two morphotypes, it was not possible to assign a described species. Among those, Plasmodium screened one species was only detected by microscopy. Considering the potential species diversity, it is possible that commonly used primers may not detect all species, reinforcing the importance of using microscopy in haematozoa surveys. There was no correspondence between the morphological traits associated with the subgenera and the phylogenetic relationships that we found in our analyses. Additionally, we found an expansion in the geographical distribution of these two species, and a new host for P. kentropyxi, demonstrating that studies of tropical herpetofauna and their parasites deserve more attention.

Parasite spillover: indirect effects of invasive Burmese pythons

Identification of the origin of parasites of nonindigenous species (NIS) can be complex. NIS may introduce parasites from their native range and acquire parasites from within their invaded range. Determination of whether parasites are non-native or native can be complicated when parasite genera occur within both the NIS' native range and its introduced range. We explored potential for spillover and spillback of lung parasites infecting Burmese pythons (Python bivittatus) in their invasive range (Florida). We collected 498 indigenous snakes of 26 species and 805 Burmese pythons during 2004-2016 and examined them for lung parasites. We used morphology to identify three genera of pentastome parasites, Raillietiella, a cosmopolitan form, and Porocephalus and Kiricephalus, both New World forms. We sequenced these parasites at one mitochondrial and one nuclear locus and showed that each genus is represented by a single species, R. orientalis, P. crotali, and K. coarctatus. Pythons are host to R. orientalis and P. crotali, but not K. coarctatus; native snakes are host to all three species. Sequence data show that pythons introduced R. orientalis to North America, where this parasite now infects native snakes. Additionally, our data suggest that pythons are competent hosts to P. crotali, a widespread parasite native to North and South America that was previously hypothesized to infect only viperid snakes. Our results indicate invasive Burmese pythons have affected parasite-host dynamics of native snakes in ways that are consistent with parasite spillover and demonstrate the potential for indirect effects during invasions. Additionally, we show that pythons have acquired a parasite native to their introduced range, which is the initial condition necessary for parasite spillback.

A Phylogeographic Assessment of the Malagasy Giant Chameleons (Furcifer verrucosus and Furcifer oustaleti)

The Malagasy giant chameleons (Furcifer oustaleti and Furcifer verrucosus) are sister species that are both broadly distributed in Madagascar, and also endemic to the island. These species are also morphologically similar and, because of this, have been frequently misidentified in the field. Previous studies have suggested that cryptic species are nested within this chameleon group, and two subspecies have been described in F. verrucosus. In this study, we utilized a phylogeographic approach to assess genetic diversification within these chameleons. This was accomplished by (1) identifying clades within each species supported by both mitochondrial and nuclear DNA, (2) assessing divergence times between clades, and (3) testing for niche divergence or conservatism. We found that both F. oustaleti and F. verrucosus could be readily identified based on genetic data, and within each species, there are two well-supported clades. However, divergence times are not contemporary and spatial patterns are not congruent. Diversification within F. verrucosus occurred during the Plio-Pleistocene, and there is evidence for niche divergence between a southwestern and southeastern clade, in a region of Madagascar that shows no obvious landscape barriers to dispersal. Diversification in F. oustaleti occurred earlier in the Pliocene or Miocene, and niche conservatism is supported with two genetically distinct clades separated at the Sofia River in northwestern Madagascar. Divergence within F. verrucosus is most consistent with patterns expected from ecologically mediated speciation, whereas divergence in F. oustaleti most strongly matches the patterns expected from the riverine barrier hypothesis.

Clonal reproduction shapes evolution in the lizard malaria parasite Plasmodium floridense

The preponderant clonal evolution hypothesis (PCE) predicts that frequent clonal reproduction (sex between two clones) in many pathogens capable of sexual recombination results in strong linkage disequilibrium and the presence of discrete genetic subdivisions characterized by occasional gene flow. We expand on the PCE and predict that higher rates of clonal reproduction will result in: (1) morphologically cryptic species that exhibit (2) low within-species variation and (3) recent between-species divergence. We tested these predictions in the Caribbean lizard malaria parasite Plasmodium floridense using 63 single-infection samples in lizards collected from across the parasite's range, and sequenced them at two mitochondrial, one apicoplast, and five nuclear genes. We identified 11 provisionally cryptic species within P. floridense, each of which exhibits low intraspecific variation and recent divergence times between species (some diverged approximately 110,000 years ago). Our results are consistent with the hypothesis that clonal reproduction can profoundly affect diversification of species capable of sexual recombination, and suggest that clonal reproduction may have led to a large number of unrecognized pathogen species. The factors that may influence the rates of clonal reproduction among pathogens are unclear, and we discuss how prevalence and virulence may relate to clonal reproduction.

Avian haemosporidians from Neotropical highlands: Evidence from morphological and molecular data

Avian haemosporidian parasites have been scarcely studied in the Neotropical highlands despite the high avian diversity reported and the uniqueness of these ecosystems. The aims of this study were to examine Haemoproteus and Plasmodium diversity based on morphological and molecular data, as well as to explore the concordance between these two approaches, when identifying species. We sampled 1487 birds belonging to 166 species, in localities of the Colombian Andean region at elevations ranging from 2100 to 4000 m above sea level. Here, we report twelve morphological parasite species, of which five are undescribed. Thirty parasite cytochrome b lineages are reported, 17 of which for the first time. We provide morphological information and illustrations, as well as, cytochrome b lineages for six morphospecies: Haemoproteus columbae, Haemoproteus witti, Haemoproteus coatneyi, Haemoproteus vireonis, Plasmodium lutzi, and Plasmodium unalis. This is the first report to provide a linkage between morphology and a molecular lineage for H. witti. Cytochrome b gene proved to be useful for species determination as DNA barcoding. Differences in parasite composition between lowlands and highlands in Colombia suggest a replacement of avian Plasmodium fauna. Parasite lineages restricted to either Colombian resident or Nearctic migratory birds were found; but a single lineage common in both has not been recorded in Nearctic non-migratory birds. We generated valuable information by using both morphological and molecular data representing competent host-parasite relationships which are based on observation of gametocytes in circulation; and increased the taxon sampling of avian haemosporidian.

Malaria parasite genetics: doing something useful

Genetics has informed almost every aspect of the study of malaria parasites, and remains a key component of much of the research that aims to reduce the burden of the disease they cause. We describe the history of genetic studies of malaria parasites and give an overview of the utility of the discipline to malariology.

Species limits in avian malaria parasites (Haemosporida): how to move forward in the molecular era

Delimiting species of malaria parasites (Haemosporida) has become increasingly problematic as new lineages of parasites are identified solely by molecular information, particularly mitochondrial cytochrome b sequence data. In this review, we highlight some of the issues, both historical and contemporary, that have hindered the development of objective criteria to diagnose, delimit and define species of haemosporidians. Defining species is not the focal interest of most researchers, most of whom merely wish to determine whether lineages identified in their samples match those of other researchers, and if so, where and in which host species. Rather than revisiting all the issues with respect to delimiting and naming species, we instead focus on finding a practical near-term resolution to the 'species problem' that utilizes the community's largest resource: mitochondrial cytochrome b DNA sequences. We recommend a standardized procedure to 'tag' these sequences, based on per cent sequence similarity, that will allow researchers to directly assess the novelty, known hosts and geographic distribution of avian malaria parasite lineages.

Malaria's many mates: past, present, and future of the systematics of the order Haemosporida

Malaria has been one of the most important diseases of humans throughout history and continues to be a major public health concern. The 5 species of Plasmodium that cause the disease in humans are part of the order Haemosporida, a diverse group of parasites that all have heteroxenous life cycles, alternating between a vertebrate host and a free-flying, blood-feeding dipteran vector. Traditionally, the identification and taxonomy of these parasites relied heavily on life-history characteristics, basic morphological features, and the host species infected. However, molecular approaches to resolving the phylogeny of the group have sometimes challenged many of these traditional hypotheses. One of the greatest debates has concerned the origin of the most virulent of the human-infecting parasites, Plasmodium falciparum, with early results suggesting a close relationship with an avian parasite. Subsequent phylogenetic studies placed it firmly within the mammalian clade instead, but the avian origin hypothesis has been revived with recent genome-based analyses. The rooting of the tree of Haemosporida has also been inconsistent, and the various topologies that result certainly affect our interpretation of the history of the group. There is clearly a pressing need to obtain a much more complete degree of taxon sampling of haemosporidians, as well as a greater number of characters before confidence can be placed in any hypothesis regarding the evolutionary history of the order. There are numerous challenges moving forward, particularly for generating complete genome sequences of avian and saurian parasites.

A new species of Haemoproteus from a tortoise (Testudo graeca) in Turkey, with remarks on molecular phylogenetic and morphological analysis

Haemoproteus anatolicum n. sp. was identified in the tortoise Testudo graeca. The new species is described based on the morphology of its blood stages and a segment of the mitochondrial cytochrome b gene, which can be used for molecular identification and diagnosis. Chelonian haemoproteids recorded in the past were defined solely on the basis of their morphological characteristics. The chelonian haemoproteid we describe as a new species has a close genetic relationship to lizard haemoproteids, i.e., Haemoproteus ptyodactylii and Haemoproteus kopki. The new species description provides significant new information for little-known chelonian haemoproteids.