Plasmodium parasites of birds have the most AT-rich genes of eukaryotes

Lost in translation: conserved amino acid usage despite extreme codon bias in foraminifera

Analyses of codon usage in eukaryotes suggest that amino acid usage responds to GC pressure so AT-biased substitutions drive higher usage of amino acids with AT-ending codons. Here, we combine single-cell transcriptomics and phylogenomics to explore codon usage patterns in foraminifera, a diverse and ancient clade of predominantly uncultivable microeukaryotes. We curate data from 1,044 gene families in 49 individuals representing 28 genera, generating perhaps the largest existing dataset of data from a predominantly uncultivable clade of protists, to analyze compositional bias and codon usage. We find extreme variation in composition, with a median GC content at fourfold degenerate silent sites below 3% in some species and above 75% in others. The most AT-biased species are distributed among diverse non-monophyletic lineages. Surprisingly, despite the extreme variation in compositional bias, amino acid usage is highly conserved across all foraminifera. By analyzing nucleotide, codon, and amino acid composition within this diverse clade of amoeboid eukaryotes, we expand our knowledge of patterns of genome evolution across the eukaryotic tree of life.IMPORTANCEPatterns of molecular evolution in protein-coding genes reflect trade-offs between substitution biases and selection on both codon and amino acid usage. Most analyses of these factors in microbial eukaryotes focus on model species such as Acanthamoeba, Plasmodium, and yeast, where substitution bias is a primary contributor to patterns of amino acid usage. Foraminifera, an ancient clade of single-celled eukaryotes, present a conundrum, as we find highly conserved amino acid usage underlain by divergent nucleotide composition, including extreme AT-bias at silent sites among multiple non-sister lineages. We speculate that these paradoxical patterns are enabled by the dynamic genome structure of foraminifera, whose life cycles can include genome endoreplication and chromatin extrusion.

High subtelomeric GC content in the genome of a zoonotic Cryptosporidium species

Cryptosporidium canis is a zoonotic species causing cryptosporidiosis in humans in addition to its natural hosts dogs and other fur animals. To understand the genetic basis for host adaptation, we sequenced the genomes of C. canis from dogs, minks, and foxes and conducted a comparative genomics analysis. While the genomes of C. canis have similar gene contents and organisations, they (~41.0 %) and C. felis (39.6 %) have GC content much higher than other Cryptosporidium spp. (24.3-32.9 %) sequenced to date. The high GC content is mostly restricted to subtelomeric regions of the eight chromosomes. Most of these GC-balanced genes encode Cryptosporidium-specific proteins that have intrinsically disordered regions and are involved in host-parasite interactions. Natural selection appears to play a more important role in the evolution of codon usage in GC-balanced C. canis, and most of the GC-balanced genes have undergone positive selection. While the identity in whole genome sequences between the mink- and dog-derived isolates is 99.9 % (9365 SNVs), it is only 96.0 % (362 894 SNVs) between them and the fox-derived isolate. In agreement with this, the fox-derived isolate possesses more subtelomeric genes encoding invasion-related protein families. Therefore, the change in subtelomeric GC content appears to be responsible for the more GC-balanced C. canis genomes, and the fox-derived isolate could represent a new Cryptosporidium species.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

Avian Malaria in Penguins: Diagnostics and Future Direction in the Context of Climate Change

Simple Summary

Avian malaria is caused by infection with protozoa of the genus Plasmodium. This vector-borne parasite is spread by mosquitoes and has a variable significance depending on environmental, host, mosquito and parasite factors. Captive penguins in non-native environments are exposed to the protozoa without having coevolved with them and are especially sensitive to infection. The most common presentation of the disease in affected penguins is acute death. Infection of wild penguins is reported and a greater understanding of the significance of such infections is required. Global warming and related surges in vector availability present an increasing threat to conservation in captive environments and targeted research into the early diagnosis of disease is required. Current diagnostic methods predominantly rely upon direct microscopy and/or molecular testing on tissues obtained from penguin postmortem examinations, and frequently fail to identify the causative agent at a species level. There are several barriers to the development of a rapid method to detect infection and the causative species; however, this information would further our understanding of this disease, and development of such a method is a valuable undertaking. This paper provides a summary of current diagnostic methods, identifies the likely future impacts of avian malaria in penguins, and highlights the need to improve both the speed and scope of available diagnostics.

Abstract

Avian malaria is caused by infection with haemoprotozoa of the genus Plasmodium. Infection is endemic in large parts of the world and is typically subclinical in birds that are native to these regions. Several penguin species have evolved in non-endemic regions without the selective pressure that these parasites exert and are highly susceptible to infection when transplanted to endemic regions, for example, in the context of zoological collections or rehabilitation centers. Avian malaria in penguins typically causes acute mortality without premonitory signs, or less commonly, nonspecific signs of morbidity, followed by mortality. Additionally, infection is reported in wild penguins, though the significance of these infections remains equivocal. As global temperatures continue to increase, avian malaria is likely to pose a continued and further threat to conservation efforts in captive environments. Intra vitam diagnosis currently relies on the evaluation of blood smears and molecular methods. The former is unreliable in penguins, as the acute clinical course typically does not allow the development of parasitemia. This absence of parasitemia also makes speciation challenging. Current molecular methods typically target the Cytochrome B or 18s subunit and have proven variably sensitive and specific. Reliable intra vitam diagnosis of avian malaria and further information about the causative agents at a species level would be very valuable in understanding the epidemiology and likely future course of avian malaria infection in penguins, and in particular, the implications avian malaria may have for conservation efforts. This paper provides an overview of malaria in penguins, discusses its changing impact on management and conservation, offers a summary of current diagnostics, and suggests future direction for the development of diagnostic tests. The latter will be key in understanding and managing this disease.

Cryptosporidium felis differs from other Cryptosporidium spp. in codon usage

Cryptosporidium spp. are important enteric pathogens in a wide range of vertebrates including humans. Previous comparative analysis revealed conservation in genome composition, gene content, and gene organization among Cryptosporidium spp., with a progressive reductive evolution in metabolic pathways and invasion-related proteins. In this study, we sequenced the genome of zoonotic pathogen Cryptosporidium felis and conducted a comparative genomic analysis. While most intestinal Cryptosporidium species have similar genomic characteristics and almost complete genome synteny, fewer protein-coding genes and some sequence inversions and translocations were found in the C. felis genome. The C. felis genome exhibits much higher GC content (39.6 %) than other Cryptosporidium species (24.3–32.9 %), especially at the third codon position (GC3) of protein-coding genes. Thus, C. felis has a different codon usage, which increases the use of less energy costly amino acids (Gly and Ala) encoded by GC-rich codons. While the tRNA usage is conserved among Cryptosporidium species, consistent with its higher GC content, C. felis uses a unique tRNA for GTG for valine instead of GTA in other Cryptosporidium species. Both mutational pressures and natural selection are associated with the evolution of the codon usage in Cryptosporidium spp., while natural selection seems to drive the codon usage in C. felis. Other unique features of the C. felis genome include the loss of the entire traditional and alternative electron transport systems and several invasion-related proteins. Thus, the preference for the use of some less energy costly amino acids in C. felis may lead to a more harmonious parasite–host interaction, and the strengthened host-adaptation is reflected by the further reductive evolution of metabolism and host invasion-related proteins.

Fifty Generations of Amitosis: Tracing Asymmetric Allele Segregation in Polyploid Cells with Single-Cell DNA Sequencing

Amitosis is a widespread form of unbalanced nuclear division whose biomedical and evolutionary significance remain unclear. Traditionally, insights into the genetics of amitosis have been gleaned by assessing the rate of phenotypic assortment. Though powerful, this experimental approach relies on the availability of phenotypic markers. Leveraging Paramecium tetraurelia, a unicellular eukaryote with nuclear dualism and a highly polyploid somatic nucleus, we probe the limits of single-cell whole-genome sequencing to study the consequences of amitosis. To this end, we first evaluate the suitability of single-cell sequencing to study the AT-rich genome of P. tetraurelia, focusing on common sources of genome representation bias. We then asked: can alternative rearrangements of a given locus eventually assort after a number of amitotic divisions? To address this question, we track somatic assortment of developmentally acquired Internal Eliminated Sequences (IESs) up to 50 amitotic divisions post self-fertilization. To further strengthen our observations, we contrast empirical estimates of IES retention levels with in silico predictions obtained through mathematical modeling. In agreement with theoretical expectations, our empirical findings are consistent with a mild increase in variation of IES retention levels across successive amitotic divisions of the macronucleus. The modest levels of somatic assortment in P. tetraurelia suggest that IESs retention levels are largely sculpted at the time of macronuclear development, and remain fairly stable during vegetative growth. In forgoing the requirement for phenotypic assortment, our approach can be applied to a wide variety of amitotic species and could facilitate the identification of environmental and genetic factors affecting amitosis.

Phylogenomics of a new fungal phylum reveals multiple waves of reductive evolution across Holomycota

Compared to multicellular fungi and unicellular yeasts, unicellular fungi with free-living flagellated stages (zoospores) remain poorly known and their phylogenetic position is often unresolved. Recently, rRNA gene phylogenetic analyses of two atypical parasitic fungi with amoeboid zoospores and long kinetosomes, the sanchytrids Amoeboradix gromovi and Sanchytrium tribonematis, showed that they formed a monophyletic group without close affinity with known fungal clades. Here, we sequence single-cell genomes for both species to assess their phylogenetic position and evolution. Phylogenomic analyses using different protein datasets and a comprehensive taxon sampling result in an almost fully-resolved fungal tree, with Chytridiomycota as sister to all other fungi, and sanchytrids forming a well-supported, fast-evolving clade sister to Blastocladiomycota. Comparative genomic analyses across fungi and their allies (Holomycota) reveal an atypically reduced metabolic repertoire for sanchytrids. We infer three main independent flagellum losses from the distribution of over 60 flagellum-specific proteins across Holomycota. Based on sanchytrids' phylogenetic position and unique traits, we propose the designation of a novel phylum, Sanchytriomycota. In addition, our results indicate that most of the hyphal morphogenesis gene repertoire of multicellular fungi had already evolved in early holomycotan lineages.

Evolutionary pressures and codon bias in low complexity regions of plasmodia

The biological meaning of low complexity regions in the proteins of Plasmodium species is a topic of discussion in evolutionary biology. There is a debate between selectionists and neutralists, who either attribute or do not attribute an effect of low-complexity regions on the fitness of these parasites, respectively. In this work, we comparatively study 22 Plasmodium species to understand whether their low complexity regions undergo a neutral or, rather, a selective and species-dependent evolution. The focus is on the connection between the codon repertoire of the genetic coding sequences and the occurrence of low complexity regions in the corresponding proteins. The first part of the work concerns the correlation between the length of plasmodial proteins and their propensity at embedding low complexity regions. Relative synonymous codon usage, entropy, and other indicators reveal that the incidence of low complexity regions and their codon bias is species-specific and subject to selective evolutionary pressure. We also observed that protein length, a relaxed selective pressure, and a broad repertoire of codons in proteins, are strongly correlated with the occurrence of low complexity regions. Overall, it seems plausible that the codon bias of low-complexity regions contributes to functional innovation and codon bias enhancement of proteins on which Plasmodium species rest as successful evolutionary parasites.

Identification of long regulatory elements in the genome of Plasmodium falciparum and other eukaryotes

Long regulatory elements (LREs), such as CpG islands, polydA:dT tracts or AU-rich elements, are thought to play key roles in gene regulation but, as opposed to conventional binding sites of transcription factors, few methods have been proposed to formally and automatically characterize them. We present here a computational approach named DExTER (Domain Exploration To Explain gene Regulation) dedicated to the identification of candidate LREs (cLREs) and apply it to the analysis of the genomes of P. falciparum and other eukaryotes. Our analyses show that all tested genomes contain several cLREs that are somewhat conserved along evolution, and that gene expression can be predicted with surprising accuracy on the basis of these long regions only. Regulation by cLREs exhibits very different behaviours depending on species and conditions. In P. falciparum and other Apicomplexan organisms as well as in Dictyostelium discoideum, the process appears highly dynamic, with different cLREs involved at different phases of the life cycle. For multicellular organisms, the same cLREs are involved in all tissues, but a dynamic behavior is observed along embryonic development stages. In P. falciparum, whose genome is known to be strongly depleted of transcription factors, cLREs are predictive of expression with an accuracy above 70%, and our analyses show that they are associated with both transcriptional and post-transcriptional regulation signals. Moreover, we assessed the biological relevance of one LRE discovered by DExTER in P. falciparum using an in vivo reporter assay. The source code (python) of DExTER is available at https://gite.lirmm.fr/menichelli/DExTER.

Genomic and proteomic biases inform metabolic engineering strategies for anaerobic fungi

Anaerobic fungi (Neocallimastigomycota) are emerging non-model hosts for biotechnology due to their wealth of biomass-degrading enzymes, yet tools to engineer these fungi have not yet been established. Here, we show that the anaerobic gut fungi have the most GC depleted genomes among 443 sequenced organisms in the fungal kingdom, which has ramifications for heterologous expression of genes as well as for emerging CRISPR-based genome engineering approaches. Comparative genomic analyses suggest that anaerobic fungi may contain cellular machinery to aid in sexual reproduction, yet a complete mating pathway was not identified. Predicted proteomes of the anaerobic fungi also contain an unusually large fraction of proteins with homopolymeric amino acid runs consisting of five or more identical consecutive amino acids. In particular, threonine runs are especially enriched in anaerobic fungal carbohydrate active enzymes (CAZymes) and this, together with a high abundance of predicted N-glycosylation motifs, suggests that gut fungal CAZymes are heavily glycosylated, which may impact heterologous production of these biotechnologically useful enzymes. Finally, we present a codon optimization strategy to aid in the development of genetic engineering tools tailored to these early-branching anaerobic fungi.

Plasmodium falciparum translational machinery condones polyadenosine repeats

Plasmodium falciparum is a causative agent of human malaria. Sixty percent of mRNAs from its extremely AT-rich (81%) genome harbor long polyadenosine (polyA) runs within their ORFs, distinguishing the parasite from its hosts and other sequenced organisms. Recent studies indicate polyA runs cause ribosome stalling and frameshifting, triggering mRNA surveillance pathways and attenuating protein synthesis. Here, we show that P. falciparum is an exception to this rule. We demonstrate that both endogenous genes and reporter sequences containing long polyA runs are efficiently and accurately translated in P. falciparum cells. We show that polyA runs do not elicit any response from No Go Decay (NGD) or result in the production of frameshifted proteins. This is in stark contrast to what we observe in human cells or T. thermophila, an organism with similar AT-content. Finally, using stalling reporters we show that Plasmodium cells evolved not to have a fully functional NGD pathway.

Differential gene expression of Plasmodium homocircumflexum (lineage pCOLL4) across two experimentally infected passerine bird species

Plasmodium parasites are present in a wide range of host species, some of which tend to be more susceptible than others, potentially as an outcome of evolved tolerance or resistance. Common starlings seem to cope with malaria infection while common crossbills are more susceptible to the same infections. That raises the question if the parasites rely on the same molecular mechanisms regardless of host species or do Plasmodium parasites change gene-expressions in accordance to the environment different hosts might provide? We used RNA-sequencing from starlings and crossbills, experimentally infected with Plasmodium homocircumflexum (lineage pCOLL4). The assembled transcriptome contained a total of 26,733 contigs. Parasite expression patterns differed between bird species. Parasites had higher expression of cell-invasion genes when infecting crossbills compared to starlings whereas in starlings genes related to apoptosis or/and oxidative stress showed higher expression levels. This article reveals how a Plasmodium parasite might adjust its expression and gene function depending on the host species infected.

Conserved associations between G-quadruplex-forming DNA motifs and virulence gene families in malaria parasites

Background

The Plasmodium genus of malaria parasites encodes several families of antigen-encoding genes. These genes tend to be hyper-variable, highly recombinogenic and variantly expressed. The best-characterized family is the var genes, exclusively found in the Laveranian subgenus of malaria parasites infecting humans and great apes. Var genes encode major virulence factors involved in immune evasion and the maintenance of chronic infections. In the human parasite P. falciparum, var gene recombination and diversification appear to be promoted by G-quadruplex (G4) DNA motifs, which are strongly associated with var genes in P. falciparum. Here, we investigated how this association might have evolved across Plasmodium species – both Laverania and also more distantly related species which lack vars but encode other, more ancient variant gene families.

Results

The association between var genes and G4-forming motifs was conserved across Laverania, spanning ~ 1 million years of evolutionary time, with suggestive evidence for evolution of the association occurring within this subgenus. In rodent malaria species, G4-forming motifs were somewhat associated with pir genes, but this was not conserved in the Laverania, nor did we find a strong association of these motifs with any gene family in a second outgroup of avian malaria parasites. Secondly, we compared two different G4 prediction algorithms in their performance on extremely A/T-rich Plasmodium genomes, and also compared these predictions with experimental data from G4-seq, a DNA sequencing method for identifying G4-forming motifs. We found a surprising lack of concordance between the two algorithms and also between the algorithms and G4-seq data.

Conclusions

G4-forming motifs are uniquely strongly associated with Plasmodium var genes, suggesting a particular role for G4s in recombination and diversification of these genes. Secondly, in the A/T-rich genomes of Plasmodium species, the choice of prediction algorithm may be particularly influential when studying G4s in these important protozoan pathogens.

Towards deorphanizing G protein-coupled receptors of Schistosoma mansoni using the MALAR yeast two-hybrid system

Schistosomiasis is an acute and chronic disease caused by parasitic worms of the genus Schistosoma. Treatment is solely dependent on praziquantel. In the face of the worldwide dimension, projects have been initiated to develop new chemotherapies. Due to their proven druggability, G protein-coupled receptors (GPCRs) are promising targets for anthelmintics. However, to identify candidate receptors, a deeper understanding of GPCR signalling in schistosome biology is essential. Comparative transcriptomics of paired and unpaired worms and their gonads revealed 59 differentially regulated GPCR-coding genes putatively involved in neuronal processes. In general, the diversity among GPCRs and their integral membrane topology make it difficult to characterize and deorphanize these receptors. To overcome existing limitations, we performed a pilot approach and utilized the innovative Membrane-Anchored Ligand And Receptor yeast two-hybrid system (MALAR-Y2H) to associate potential neuropeptide ligands with their cognate receptors. Here, we demonstrated the ability to express full-length GPCRs of Schistosoma mansoni in a heterologous yeast-based system. Additionally, we localized GPCRs and chimeras of neuropeptides fused to the WBP1 transmembrane domain of yeast to the plasma membrane of yeast cells. Reporter gene assays indicated ligand-receptor binding, which allowed us to identify certain neuropeptides as potential ligands for two GPCRs, which had been found before to be differentially expressed in schistosomes in a pairing-dependent manner. Thus, the MALAR-Y2H system appears suitable to unravel schistosome GPCR–ligand interactions. Besides its relevance for understanding schistosome biology, identifying and characterizing GPCR–ligand interaction will also contribute to applied research aspects.

Adaptation of Translational Machinery in Malaria Parasites to Accommodate Translation of Poly-Adenosine Stretches Throughout Its Life Cycle

Malaria is caused by unicellular apicomplexan parasites of the genus Plasmodium, which includes the major human parasite Plasmodium falciparum. The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively. There is tight control of gene expression in each developmental stage, and at every level of gene synthesis: from RNA transcription, to its subsequent translation, and finally post-translational modifications of the resulting protein. Whole-genome sequencing of P. falciparum has laid the foundation for significant biological advances by revealing surprising genomic information. The P. falciparum genome is extremely AT-rich (∼80%), with a substantial portion of genes encoding intragenic polyadenosine (polyA) tracks being expressed throughout the entire parasite life cycle. In most eukaryotes, intragenic polyA runs act as negative regulators of gene expression. Recent studies have shown that translation of mRNAs containing 12 or more consecutive adenosines results in ribosomal stalling and frameshifting; activating mRNA surveillance mechanisms. In contrast, P. falciparum translational machinery can efficiently and accurately translate polyA tracks without activating mRNA surveillance pathways. This unique feature of P. falciparum raises interesting questions: (1) How is P. falciparum able to efficiently and correctly translate polyA track transcripts, and (2) What are the specifics of the translational machinery and mRNA surveillance mechanisms that separate P. falciparum from other organisms? In this review, we analyze possible evolutionary shifts in P. falciparum protein synthesis machinery that allow efficient translation of an AU rich-transcriptome. We focus on physiological and structural differences of P. falciparum stage specific ribosomes, ribosome-associated proteins, and changes in mRNA surveillance mechanisms throughout the complete parasite life cycle, with an emphasis on the mosquito and liver stages.

De novo transcriptome assembly and preliminary analyses of two avian malaria parasites, Plasmodium delichoni and Plasmodium homocircumflexum

Parasites of the genus Plasmodium infect a wide array of hosts, causing malaria in all major groups of terrestrial vertebrates including primates, reptiles, and birds. Molecular mechanisms explaining why some Plasmodium species are virulent, while other closely related malaria pathogens are relatively benign in the same hosts, remain unclear. Here, we present the RNA sequencing and subsequent transcriptome assembly of two avian Plasmodium parasites which can eventually be used to better understand the genetic mechanisms underlying Plasmodium species pathogenicity in an avian host. Plasmodium homocircumflexum, a cryptic, pathogenic species that often causes mortality and Plasmodium delichoni, a newly described, relatively benign malaria parasite that does not kill its hosts, were used to experimentally infect two Eurasian siskins (Carduelis spinus). RNA extractions were performed and RNA sequencing was completed using high throughput Illumina sequencing. Using established bioinformatics pipelines, the sequencing data from both species were used to generate transcriptomes using published Plasmodium species genomes as a scaffold. The finalized transcriptome of P. homocircumflexum contained 21,612 total contigs while that of P. delichoni contained 12,048 contigs. We were able to identify many genes implicated in erythrocyte invasion actively expressed in both P. homocircumflexum and P. delichoni, including the well described vaccine candidates Apical Membrane Antigen-1 (AMA-1) and Merozoite Surface Protein 1 (MSP1). This work acts as a stepping stone to increase available data on avian Plasmodium parasites, thus enabling future research into the evolution of pathogenicity in malaria.

Genomic Advances in Avian Malaria Research

Haemosporidian parasites causing malaria-like diseases in birds are globally distributed and have been associated with reduced host fitness and mortality in susceptible bird species. This group of parasites has not only enabled a greater understanding of host specificity, virulence, and parasite dispersal, but has also been crucial in restructuring the evolutionary history of apicomplexans. Despite their importance, genomic resources of avian haemosporidians have proved difficult to obtain, and they have, as a result, been lagging behind the congeneric Plasmodium species infecting mammals. In this review, I discuss recent genomic advances in the field of avian malaria research, and outline outstanding questions that will become possible to investigate with the continued successful efforts to generate avian haemosporidian genomic data.

Delineation of the Genera Haemoproteus and Plasmodium Using RNA-Seq and Multi-gene Phylogenetics

Members of the order Haemosporida are protist parasites that infect mammals, reptiles and birds. This group includes the causal agents of malaria, Plasmodium parasites, the genera Leucocytozoon and Fallisia, as well as the species rich genus Haemoproteus with its two subgenera Haemoproteus and Parahaemoproteus. Some species of Haemoproteus cause severe disease in avian hosts, and these parasites display high levels of diversity worldwide. This diversity emphasizes the need for accurate evolutionary information. Most molecular studies of wildlife haemosporidians use a bar coding approach by sequencing a fragment of the mitochondrial cytochrome b gene. This method is efficient at differentiating parasite lineages but insufficient for accurate phylogenetic inferences in highly diverse taxa such as haemosporidians. Recent studies have utilized multiple mitochondrial genes (cyt b, cox1 and cox3), sometimes combined with a few apicoplast and nuclear genes. These studies have been highly successful with one notable exception: the evolutionary relationships of the genus Haemoproteus remain unresolved. Here we describe the transcriptome of Haemoproteus columbae and investigate its phylogenetic position recovered from a multi-gene dataset (600 genes). This genomic approach restricts the taxon sampling to 18 species of apicomplexan parasites. We employed Bayesian inference and maximum likelihood methods of phylogenetic analyses and found H. columbae and a representative from the subgenus Parahaemoproteus to be sister taxa. This result strengthens the hypothesis of genus Haemoproteus being monophyletic; however, resolving this question will require sequences of orthologs from, in particular, representatives of Leucocytozoon species.