Bat-associated microbes: Opportunities and perils, an overview

The potential biotechnological uses of bat-associated bacteria are discussed briefly, indicating avenues for biotechnological applications of bat-associated microbes. The uniqueness of bats in terms of their lifestyle, genomes and molecular immunology may predispose bats to act as disease reservoirs. Molecular phylogenetic analysis has shown several instances of bats harbouring the ancestral lineages of bacterial (Bartonella), protozoal (Plasmodium, Trypanosoma cruzi) and viral (SARS-CoV2) pathogens infecting humans. Along with the transmission of viruses from bats, we also discuss the potential roles of bat-associated bacteria, fungi, and protozoan parasites in emerging diseases. Current evidence suggests that environmental changes and interactions between wildlife, livestock, and humans contribute to the spill-over of infectious agents from bats to other hosts. Domestic animals including livestock may act as intermediate amplifying hosts for bat-origin pathogens to transmit to humans. An increasing number of studies investigating bat pathogen diversity and infection dynamics have been published. However, whether or how these infectious agents are transmitted both within bat populations and to other hosts, including humans, often remains unknown. Metagenomic approaches are uncovering the dynamics and distribution of potential pathogens in bat microbiomes, which might improve the understanding of disease emergence and transmission. Here, we summarize the current knowledge on bat zoonoses of public health concern and flag the gaps in the knowledge to enable further research and allocation of resources for tackling future outbreaks.

Ribosome Specialization in Protozoa Parasites

Ribosomes, in general, are viewed as constitutive macromolecular machines where protein synthesis takes place; however, this view has been recently challenged, supporting the hypothesis of ribosome specialization and opening a completely new field of research. Recent studies have demonstrated that ribosomes are heterogenous in their nature and can provide another layer of gene expression control by regulating translation. Heterogeneities in ribosomal RNA and ribosomal proteins that compose them favor the selective translation of different sub-pools of mRNAs and functional specialization. In recent years, the heterogeneity and specialization of ribosomes have been widely reported in different eukaryotic study models; however, few reports on this topic have been made on protozoa and even less on protozoa parasites of medical importance. This review analyzes heterogeneities of ribosomes in protozoa parasites highlighting the specialization in their functions and their importance in parasitism, in the transition between stages in their life cycle, in the change of host and in response to environmental conditions.

Plasmodium malariae: the persisting mysteries of a persistent parasite

Plasmodium malariae is a 'neglected malaria parasite' in as much as the amount of research conducted on it pales into insignificance when compared to that pertaining to Plasmodium falciparum and Plasmodium vivax, its more notorious and pathogenic cousins. There has, however, been an increase in interest in this parasite over the past decade. Principally, this is because of the increasing use of sensitive molecular detection techniques that have revealed a wider than previously recorded prevalence in some regions (particularly in Africa), and high numbers of chronic, asymptomatic infections.

Mobilome of Apicomplexa Parasites

Transposable elements (TEs) are mobile genetic elements found in the majority of eukaryotic genomes. Genomic studies of protozoan parasites from the phylum Apicomplexa have only reported a handful of TEs in some species and a complete absence in others. Here, we studied sixty-four Apicomplexa genomes available in public databases, using a ‘de novo’ approach to build candidate TE models and multiple strategies from known TE sequence databases, pattern recognition of TEs, and protein domain databases, to identify possible TEs. We offer an insight into the distribution and the type of TEs that are present in these genomes, aiming to shed some light on the process of gains and losses of TEs in this phylum. We found that TEs comprise a very small portion in these genomes compared to other organisms, and in many cases, there are no apparent traces of TEs. We were able to build and classify 151 models from the TE consensus sequences obtained with RepeatModeler, 96 LTR TEs with LTRpred, and 44 LINE TEs with MGEScan. We found LTR Gypsy-like TEs in Eimeria, Gregarines, Haemoproteus, and Plasmodium genera. Additionally, we described LINE-like TEs in some species from the genera Babesia and Theileria. Finally, we confirmed the absence of TEs in the genus Cryptosporidium. Interestingly, Apicomplexa seem to be devoid of Class II transposons.

The primate malaria parasites Plasmodium malariae, Plasmodium brasilianum and Plasmodium ovale spp.: genomic insights into distribution, dispersal and host transitions

During the twentieth century, there was an explosion in understanding of the malaria parasites infecting humans and wild primates. This was built on three main data sources: from detailed descriptive morphology, from observational histories of induced infections in captive primates, syphilis patients, prison inmates and volunteers, and from clinical and epidemiological studies in the field. All three were wholly dependent on parasitological information from blood-film microscopy, and The Primate Malarias” by Coatney and colleagues (1971) provides an overview of this knowledge available at that time. Here, 50 years on, a perspective from the third decade of the twenty-first century is presented on two pairs of primate malaria parasite species. Included is a near-exhaustive summary of the recent and current geographical distribution for each of these four species, and of the underlying molecular and genomic evidence for each. The important role of host transitions in the radiation of Plasmodium spp. is discussed, as are any implications for the desired elimination of all malaria species in human populations. Two important questions are posed, requiring further work on these often ignored taxa. Is Plasmodium brasilianum, circulating among wild simian hosts in the Americas, a distinct species from Plasmodium malariae? Can new insights into the genomic differences between Plasmodium ovale curtisi and Plasmodium ovale wallikeri be linked to any important differences in parasite morphology, cell biology or clinical and epidemiological features?

The vectors of Plasmodium knowlesi and other simian malarias Southeast Asia: challenges in malaria elimination

Plasmodium knowlesi, a simian malaria parasite of great public health concern has been reported from most countries in Southeast Asia and exported to various countries around the world. Currently P. knowlesi is the predominant species infecting humans in Malaysia. Besides this species, other simian malaria parasites such as P. cynomolgi and P. inui are also infecting humans in the region. The vectors of P. knowlesi and other Asian simian malarias belong to the Leucosphyrus Group of Anopheles mosquitoes which are generally forest dwelling species. Continual deforestation has resulted in these species moving into forest fringes, farms, plantations and human settlements along with their macaque hosts. Limited studies have shown that mosquito vectors are attracted to both humans and macaque hosts, preferring to bite outdoors and in the early part of the night. We here review the current status of simian malaria vectors and their parasites, knowledge of vector competence from experimental infections and discuss possible vector control measures. The challenges encountered in simian malaria elimination are also discussed. We highlight key knowledge gaps on vector distribution and ecology that may impede effective control strategies.

The apicoplast of Haemoproteus columbae: A comparative study of this organelle genome in Haemosporida

Apicomplexa is a phylum of parasitic protozoa; among them are the order Haemosporida, vector-borne parasites that include those that cause malaria (genus Plasmodium). Most Apicomplexa species have a non-photosynthetic plastid or apicoplast. Given its unique metabolic pathways, this organelle is considered a target for malaria therapeutics. Regardless of its importance, there is a paucity of complete apicoplast genome data hindering comparative studies. Here, the Haemoproteus (Haemoproteus) columbae apicoplast genome (lineage HAECOL1) was obtained using next-generation sequencing. This genome was included in a comparative analysis with other plastids. This 29.8 kb circular genome shares the same structure found in Plasmodium parasites. It is A + T rich (87.7%), comparable but at the higher end of A + T content observed in Plasmodium species (85.5-87.2%). As expected, considering its high A + T content, the synonymous codon usage (RSCU) and the effective number of codons (ENc) showed a moderate codon bias. Several apicoplast genes have a phylogenetic signal. However, unlike mitochondrial genes, single-gene phylogenies have low support in haemosporidian clades that diverged recently. The H. columbae apicoplast genome suggests that the apicoplast function may be conserved across Haemosporida. This parasite could be a model to study this organelle in a non-mammalian system.

A PCR-Based Technique to Track the Geographic Origin of Plasmodium falciparum With 23-SNP Barcode Analysis

Increased population movement has increased the risk of reintroducing parasites to elimination areas and also dispersing drug-resistant parasites to new regions. Therefore, reliable and repeatable methods to trace back to the source of imported infections are essential. The recently developed 23-single-nucleotide polymorphism (SNP) barcode from organellar genomes of mitochondrion (mt) and apicoplast (apico) provides a valuable tool to locate the geographic origin of Plasmodium falciparum. This study aims to explore the feasibility of using the 23-SNP barcode for tracking P. falciparum by polymerase chain reaction and sequencing, while providing geographical haplotypes of isolates that originated from Central Africa. Based on 23-SNP barcode analysis, SNPs were found at seven loci; 27 isolates were confirmed to have originated in West Africa, and this study also showed four isolates from Central Africa (Equatorial Guinea, 3; Republic of Congo, 1) that originated in East Africa. This study provides the sequence data from Central Africa and fills 23-SNP barcode data gaps of sample origins.

Comparative Plastid Genomics of Cryptomonas Species Reveals Fine-Scale Genomic Responses to Loss of Photosynthesis

Loss of photosynthesis is a recurring theme in eukaryotic evolution. In organisms that have lost the ability to photosynthesize, nonphotosynthetic plastids are retained because they play essential roles in processes other than photosynthesis. The unicellular algal genus Cryptomonas contains both photosynthetic and nonphotosynthetic members, the latter having lost the ability to photosynthesize on at least three separate occasions. To elucidate the evolutionary processes underlying the loss of photosynthesis, we sequenced the plastid genomes of two nonphotosynthetic strains, Cryptomonas sp. CCAC1634B and SAG977-2f, as well as the genome of the phototroph Cryptomonas curvata CCAP979/52. These three genome sequences were compared with the previously sequenced plastid genome of the nonphotosynthetic species Cryptomonas paramecium CCAP977/2a as well as photosynthetic members of the Cryptomonadales, including C. curvata FBCC300012D. Intraspecies comparison between the two C. curvata strains showed that although their genome structures are stable, the substitution rates of their genes are relatively high. Although most photosynthesis-related genes, such as the psa and psb gene families, were found to have disappeared from the nonphotosynthetic strains, at least ten pseudogenes are retained in SAG977-2f. Although gene order is roughly shared among the plastid genomes of photosynthetic Cryptomonadales, genome rearrangements are seen more frequently in the smaller genomes of the nonphotosynthetic strains. Intriguingly, the light-independent protochlorophyllide reductase comprising chlB, L, and N is retained in nonphotosynthetic SAG977-2f and CCAC1634B. On the other hand, whereas CCAP977/2a retains ribulose-1,5-bisphosphate carboxylase/oxygenase-related genes, including rbcL, rbcS, and cbbX, the plastid genomes of the other two nonphotosynthetic strains have lost the ribulose-1,5-bisphosphate carboxylase/oxygenase protein-coding genes.

The proofreading activity of Pfprex from Plasmodium falciparum can prevent mutagenesis of the apicoplast genome by oxidized nucleotides

The DNA polymerase module of the Pfprex enzyme (PfpPol) is responsible for duplication of the genome of the apicoplast organelle in the malaria parasite. We show that PfpPol can misincorporate oxidized nucleotides such as 8oxodGTP opposite dA. This event gives rise to transversion mutations that are known to lead to adverse physiological outcomes. The apicoplast genome is particularly vulnerable to the harmful effects of 8oxodGTP due to very high AT content (~ 87%). We show that the proofreading activity of PfpPol has the unique ability to remove the oxidized nucleotide from the primer terminus. Due to this property, the proofreading domain of PfpPol is able to prevent mutagenesis of the AT-rich apicoplast genome and neutralize the deleterious genotoxic effects of ROS generated in the apicoplast due to normal metabolic processes. The proofreading activity of the Pfprex enzyme may, therefore, represent an attractive target for therapeutic intervention. Also, a survey of DNA repair pathways shows that the observed property of Pfprex constitutes a novel form of dynamic error correction wherein the repair of promutagenic damaged nucleotides is concomitant with DNA replication.

Annotation and characterization of Babesia gibsoni apicoplast genome

Background

Babesia gibsoni is an apicomplexan parasite transmitted by ticks, which can infect canine species and cause babesiosis. The apicoplast is an organelle associated with isoprenoids metabolism, is widely present in apicomplexan parasites, except for Cryptosporidium. Available data indicate that the apicoplast is essential for the survival of apicomplexan parasites.

Methods

Here, the apicoplast genome of B. gibsoni was investigated by high-throughput genome sequencing, bioinformatics analysis, and conventional PCR.

Results

The apicoplast genome of B. gibsoni-Wuhan strain (B. gibsoni-WH) consists of a 28.4 kb circular molecule, with A + T content of 86.33%, similar to that of B. microti. Specifically, this genome encodes genes involved in maintenance of the apicoplast DNA, transcription, translation and maturation of organellar proteins, which contains 2 subunits of ribosomal RNAs, 17 ribosomal proteins, 1 EF-Tu elongation factor (tufA), 5 DNA-dependent RNA polymerase beta subunits, 2 Clp protease chaperones, 23 tRNA genes and 5 unknown open reading frames (hypothetical proteins). Phylogenetic analysis revealed high similarity of B. gibsoni apicoplast genome to that of B. orientalis and B. bovis.

Conclusions

To our knowledge, this is the first report of annotation and characterization of B. gibsoni-WH apicoplast genome. The results will facilitate the development of new anti-Babesia drug targets.

There Is Treasure Everywhere: Reductive Plastid Evolution in Apicomplexa in Light of Their Close Relatives

The phylum Apicomplexa (Alveolates) comprises a group of host-associated protists, predominately intracellular parasites, including devastating parasites like Plasmodium falciparum, the causative agent of malaria. One of the more fascinating characteristics of Apicomplexa is their highly reduced (and occasionally lost) remnant plastid, termed the apicoplast. Four core metabolic pathways are retained in the apicoplast: heme synthesis, iron–sulfur cluster synthesis, isoprenoid synthesis, and fatty acid synthesis. It has been suggested that one or more of these pathways are essential for plastid and plastid genome retention. The past decade has witnessed the discovery of several apicomplexan relatives, and next-generation sequencing efforts are revealing that they retain variable plastid metabolic capacities. These data are providing clues about the core genes and pathways of reduced plastids, while at the same time further confounding our view on the evolutionary history of the apicoplast. Here, we examine the evolutionary history of the apicoplast, explore plastid metabolism in Apicomplexa and their close relatives, and propose that the differences among reduced plastids result from a game of endosymbiotic roulette. Continued exploration of the Apicomplexa and their relatives is sure to provide new insights into the evolution of the apicoplast and apicomplexans as a whole.

Novel Molecular Synapomorphies Demarcate Different Main Groups/Subgroups of Plasmodium and Piroplasmida Species Clarifying Their Evolutionary Relationships

The class Hematozoa encompasses several clinically important genera, including Plasmodium, whose members cause the major life-threating disease malaria. Hence, a good understanding of the interrelationships of organisms from this class and reliable means for distinguishing them are of much importance. This study reports comprehensive phylogenetic and comparative analyses on protein sequences on the genomes of 28 hematozoa species to understand their interrelationships. In addition to phylogenetic trees based on two large datasets of protein sequences, detailed comparative analyses were carried out on the genomes of hematozoa species to identify novel molecular synapomorphies consisting of conserved signature indels (CSIs) in protein sequences. These studies have identified 79 CSIs that are exclusively present in specific groups of Hematozoa/Plasmodium species, also supported by phylogenetic analysis, providing reliable means for the identification of these species groups and understanding their interrelationships. Of these CSIs, six CSIs are specifically shared by all hematozoa species, two CSIs serve to distinguish members of the order Piroplasmida, five CSIs are uniquely found in all Piroplasmida species except B. microti and two CSIs are specific for the genus Theileria. Additionally, we also describe 23 CSIs that are exclusively present in all genome-sequenced Plasmodium species and two, nine, ten and eight CSIs which are specific for members of the Plasmodium subgenera Haemamoeba, Laverania, Vinckeia and Plasmodium (excluding P. ovale and P. malariae), respectively. Additionally, our work has identified several CSIs that support species relationships which are not evident from phylogenetic analysis. Of these CSIs, one CSI supports the ancestral nature of the avian-Plasmodium species in comparison to the mammalian-infecting groups of Plasmodium species, four CSIs strongly support a specific relationship of species between the subgenera Plasmodium and Vinckeia and three CSIs each that reliably group P. malariae with members of the subgenus Plasmodium and P. ovale within the subgenus Vinckeia, respectively. These results provide a reliable framework for understanding the evolutionary relationships among the Plasmodium/Piroplasmida species. Further, in view of the exclusivity of the described molecular markers for the indicated groups of hematozoa species, particularly large numbers of unique characteristics that are specific for all Plasmodium species, they provide important molecular tools for biochemical/genetic studies and for developing novel diagnostics and therapeutics for these organisms.

Comparative analysis of apicoplast genomes of Babesia infective to small ruminants in China

Background

Babesiosis is an economically important disease caused by tick-borne apicomplexan protists of the genus Babesia. Most apicomplexan parasites, including Babesia, have a plastid-derived organelle termed an apicoplast, which is involved in critical metabolic pathways such as fatty acid, iron-sulphur, haem and isoprenoid biosynthesis. Apicoplast genomic data can provide significant information for understanding and exploring the biological features, taxonomic and evolutionary relationships of apicomplexan parasites, and identify targets for anti-parasitic drugs. However, there are limited data on the apicoplast genomes of Babesia species infective to small ruminants.

Methods

PCR primers were designed based on the previously reported apicoplast genome sequences of Babesia motasi Lintan and Babesia sp. Xinjiang using Illumina technology. The overlapped apicoplast genomic fragments of six ovine Babesia isolates were amplified and sequenced using the Sanger dideoxy chain-termination method. The full-length sequences of the apicoplast genomes were assembled and annotated using bioinformatics software. The gene contents and order of apicoplast genomes obtained in this study were defined and compared with those of other apicomplexan parasites. Phylogenetic trees were constructed on the concatenated amino acid sequences of 13 gene products using MEGA v.6.06.

Results

The results showed that the six ovine Babesia apicoplast genomes consisted of circular DNA. The genome sizes were 29,916–30,846 bp with 78.7–81.0% A + T content, 29–31 open reading frames (ORF) and 23–24 transport RNAs. The ORFs encoded four DNA-directed RNA polymerase subunits (rpoB, rpoCl, rpoC2a and rpoC2b), 13 ribosomal proteins, one elongation factor TU (tufA), two ATP-dependent Clp proteases (ClpC) and 7–11 hypothetical proteins. Babesia sp. has three more genes than Babesia motasi (rpl5, rps8 and rpoB). Phylogenetic analysis showed that Babesia sp. is located in a separate clade. Babesia motasi Lintan/Tianzhu and B. motasi Ningxian/Hebei were divided into two subclades.

Conclusions

To our knowledge, this study is the first to elucidate the whole apicoplast genomic structural features of six Babesia isolates infective to small ruminants in China using Sanger sequencing. The data provide useful information confirming the taxonomic relationships of these parasites and identifying targets for anti-apicomplexan parasite drugs.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3581-x) contains supplementary material, which is available to authorized users.

Apicoplast phylogeny reveals the position of Plasmodium vivax basal to the Asian primate malaria parasite clade

The malaria parasite species, Plasmodium vivax infects not only humans, but also African apes. Human specific P. vivax has evolved from a single ancestor that originated from a parasite of African apes. Although previous studies have proposed phylogenetic trees positioning P. vivax (the common ancestor of human and African ape P. vivax) within the assemblages of Asian primate parasites, its position has not yet been robustly confirmed. We determined nearly complete apicoplast genome sequences from seven Asian primate parasites, Plasmodium cynomolgi (strains Ceylonensis and Berok), P. knowlesi P. fragile, P. fieldi, P. simiovale, P. hylobati, P. inui, and an African primate parasite, P. gonderi, that infects African guenon. Phylogenetic relationships of the Plasmodium species were analyzed using newly and previously determined apicoplast genome sequences. Multigene maximum likelihood analysis of 30 protein coding genes did not position P. vivax within the Asian primate parasite clade but positioned it basal to the clade, after the branching of an African guenon parasite, P. gonderi. The result does not contradict with the emerging notion that P. vivax phylogenetically originated from Africa. The result is also supported by phylogenetic analyses performed using massive nuclear genome data of seven primate Plasmodium species.

Novel high resolution melting (HRM) and snapback assays for simultaneous detection and differentiation of Plasmodium ovale spp

Plasmodium ovale spp. are two of the six species of apicomplexan parasites belonging to the genus Plasmodium commonly causing disease in humans. A recent phylogeny study has identified both Plasmodium ovale species (P. ovale curtisi and P. ovale wallikeri) as two sympatric occurring species. The actual prevalence and clinical relevance of P. ovale spp. are likely underestimated due to low parasitemia and mixed infections, which pose a major challenge to microscopic diagnosis and are frequently undetectable using malaria Rapid Diagnostic Tests (RDTs). The aim of this work is to develop a HRM-based assay for simultaneous detection and differentiation of P. ovale wallikeri and P. ovale curtisi. Thirty three well-documented P. ovale spp. samples from previous studies were used for this study. The newly developed High Resolution Melting (HRM) assay targeting the apicoplast genome was highly specific to both P. ovale species. Adding a snapback tail at the 5' end of the forward primer for a nested HRM PCR, increased the melting temperature (Tm) difference between the two species. To our knowledge this study reports the first direct HRM assay developed on the apicoplast genome, specific for both P. ovale species. This method provides added value to the WHO open request of developing new practical malaria diagnostic methods for the malaria elimination program and could contribute to a quick and efficient diagnosis of low-level parasitemia, symptomatic or asymptomatic, as well as mixed or single P. ovale infections.

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.

Multiple origins of endosymbionts in Chlorellaceae with no reductive effects on the plastid or mitochondrial genomes

Ancient endosymbiotic relationships have led to extreme genomic reduction in many bacterial and eukaryotic algal endosymbionts. Endosymbionts in more recent and/or facultative relationships can also experience genomic reduction to a lesser extent, but little is known about the effects of the endosymbiotic transition on the organellar genomes of eukaryotes. To understand how the endosymbiotic lifestyle has affected the organellar genomes of photosynthetic green algae, we generated the complete plastid genome (plastome) and mitochondrial genome (mitogenome) sequences from three green algal endosymbionts (Chlorella heliozoae, Chlorella variabilis and Micractinium conductrix). The mitogenomes and plastomes of the three newly sequenced endosymbionts have a standard set of genes compared with free-living trebouxiophytes, providing no evidence for functional genomic reduction. Instead, their organellar genomes are generally larger and more intron rich. Intron content is highly variable among the members of Chlorella, suggesting very high rates of gain and/or loss of introns during evolution. Phylogenetic analysis of plastid and mitochondrial genes demonstrated that the three endosymbionts do not form a monophyletic group, indicating that the endosymbiotic lifestyle has evolved multiple times in Chlorellaceae. In addition, M. conductrix is deeply nested within the Chlorella clade, suggesting that taxonomic revision is needed for one or both genera.

Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution

Elucidation of the evolutionary history and interrelatedness of Plasmodium species that infect humans has been hampered by a lack of genetic information for three human-infective species: P. malariae and two P. ovale species (P. o. curtisi and P. o. wallikeri). These species are prevalent across most regions in which malaria is endemic and are often undetectable by light microscopy, rendering their study in human populations difficult. The exact evolutionary relationship of these species to the other human-infective species has been contested. Using a new reference genome for P. malariae and a manually curated draft P. o. curtisi genome, we are now able to accurately place these species within the Plasmodium phylogeny. Sequencing of a P. malariae relative that infects chimpanzees reveals similar signatures of selection in the P. malariae lineage to another Plasmodium lineage shown to be capable of colonization of both human and chimpanzee hosts. Molecular dating suggests that these host adaptations occurred over similar evolutionary timescales. In addition to the core genome that is conserved between species, differences in gene content can be linked to their specific biology. The genome suggests that P. malariae expresses a family of heterodimeric proteins on its surface that have structural similarities to a protein crucial for invasion of red blood cells. The data presented here provide insight into the evolution of the Plasmodium genus as a whole.

Genome-scale comparison of expanded gene families in Plasmodium ovale wallikeri and Plasmodium ovale curtisi with Plasmodium malariae and with other Plasmodium species

Malaria in humans is caused by six species of Plasmodium parasites, of which the nuclear genome sequences for the two Plasmodium ovale spp., P. ovale curtisi and P. ovale wallikeri, and Plasmodium malariae have not yet been analyzed. Here we present an analysis of the nuclear genome sequences of these three parasites, and describe gene family expansions therein. Plasmodium ovale curtisi and P. ovale wallikeri are genetically distinct but morphologically indistinguishable and have sympatric ranges through the tropics of Africa, Asia and Oceania. Both P. ovale spp. show expansion of the surfin variant gene family, and an amplification of the Plasmodium interspersed repeat (pir) superfamily which results in an approximately 30% increase in genome size. For comparison, we have also analyzed the draft nuclear genome of P. malariae, a malaria parasite causing mild malaria symptoms with a quartan life cycle, long-term chronic infections, and wide geographic distribution. Plasmodium malariae shows only a moderate level of expansion of pir genes, and unique expansions of a highly diverged transmembrane protein family with over 550 members and the gamete P25/27 gene family. The observed diversity in the P. ovale wallikeri and P. ovale curtisi surface antigens, combined with their phylogenetic separation, supports consideration that the two parasites be given species status.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Molecular characterization of misidentified Plasmodium ovale imported cases in Singapore

BACKGROUND

Plasmodium ovale, considered the rarest of the malaria parasites of humans, consists of two morphologically identical but genetically distinct sympatric species, Plasmodium ovale curtisi and Plasmodium ovale wallikeri. These parasites resemble morphologically to Plasmodium vivax with which they also share a tertian periodicity and the ability to cause relapses, making them easily misidentified as P. vivax. Plasmodium ovale infections are rarely reported, but given the likelihood of misidentification, their prevalence might be underestimated.

METHODS

Morphological and molecular analysis of confirmed malaria cases admitted in Singapore in 2012-2014 detected nine imported P. ovale cases that had been misidentified as P. vivax. Since P. ovale had not been previously officially reported in Singapore, a retrospective analysis of available, frozen, archival blood samples was performed and returned two additional misidentified P. ovale cases in 2003 and 2006. These eleven P. ovale samples were characterized with respect to seven molecular markers (ssrRNA, Potra, Porbp2, Pog3p, dhfr-ts, cytb, cox1) used in recent studies to distinguish between the two sympatric species, and to a further three genes (tufa, clpC and asl).

RESULTS

The morphological features of P. ovale and the differential diagnosis with P. vivax were reviewed and illustrated by microphotographs. The genetic dimorphism between P. ovale curtisi and P. ovale wallikeri was assessed by ten molecular markers distributed across the three genomes of the parasite (Genbank KP050361-KP050470). The data obtained for seven of these markers were compared with those published and confirmed that both P. ovale species were present. This dimorphism was also confirmed for the first time on: (1) two genes from the apicoplast genome (tufA and clpC genes); and, (2) the asl gene that was used for phylogenetic analyses of other Plasmodium species, and that was found to harbour the highest number of dimorphic loci between the two P. ovale species.

CONCLUSION

Misidentified P. ovale infections are reported for the first time among imported malaria cases in Singapore. Genetic dimorphism between P. ovale curtisi and P. ovale wallikeri was confirmed using markers from the parasites' three genomes. The apparent increase of imported P. ovale since 2012 (with yearly detection of cases) is puzzling. Given decrease in the overall number of malaria cases recorded in Singapore since 2010 the 'resurgence' of this neglected species raises public health concerns.

Characterization and annotation of Babesia orientalis apicoplast genome

Background

Babesia orientalis is an obligate intraerythrocytic protozoan parasite of the buffalo (Bubalus bubalis, Linnaeus, 1758) transmitted by the tick Rhipicephalus heamaphysaloides. It is the causative agent of water buffalo babesiosis, one of the most important pathogens of water buffalo in central and southern China. As a member of the phylum Apicomplexa, B. orientalis possesses a relatively independent and alga originated organelle the apicoplast. Apicoplasts in other apicomplexa parasites are involved in the biosynthesis of haem, fatty acids, iron-sulphur clusters and isoprenoids. Some of these metabolic pathways were shown to be essential for parasite survival, therefore can serve as potential drug targets.

Methods

30 pairs of primers were designed based on the full genome sequence of B. orientalis (unpublished data) and by aligning reported apicoplast genomes of Babesia bovis and Theileria parva. Conventional PCRs was performed to obtain overlapped fragments to cover the whole apicoplast genome. Then the apicoplast genome of B.orientalis was sequenced, assembled and aligned with reported apicoplast genomes of B. bovis and T. parva. The obtained apicoplast genome was annotated by using Artemis and comparing with published apicomplexan apicoplast genomes. The SSU and LSU nucleotide sequences generated were used in a phylogenetic analysis using the maximum likelihood implemented in MAGE 6.0.

Results

We have obtained and analyzed the complete genome sequence of the B. orientalis apicoplast. It consisted of a 33.2 kb circular DNA (78.9 % A + T). The apicoplast genome unidirectionally encodes one large and one small subunit ribosomal RNAs, 24 tRNA genes, 4 DNA-dependent RNA polymerase beta subunits (rpoB, rpoC1, rpoC2a and rpoC2b), 17 ribosomal proteins, one EF-Tu elongation factor, 2 Clp protease chaperones, and 14 hypothetical proteins. In addition, it includes two copies of the clpC gene. The structure and organization of the B. orientalis apicoplast genome are most similar to those of the B. bovis apicoplast.

Conclusions

This is the first report of the complete sequence of the B. orientalis apicoplast genome. This information should be useful in the development of safe and efficient treatment against buffalo babesiosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1158-x) contains supplementary material, which is available to authorized users.

Plasmodium Apicoplast Gln-tRNAGln Biosynthesis Utilizes a Unique GatAB Amidotransferase Essential for Erythrocytic Stage Parasites

The malaria parasite Plasmodium falciparum apicoplast indirect aminoacylation pathway utilizes a non-discriminating glutamyl-tRNA synthetase to synthesize Glu-tRNA(Gln) and a glutaminyl-tRNA amidotransferase to convert Glu-tRNA(Gln) to Gln-tRNA(Gln). Here, we show that Plasmodium falciparum and other apicomplexans possess a unique heterodimeric glutamyl-tRNA amidotransferase consisting of GatA and GatB subunits (GatAB). We localized the P. falciparum GatA and GatB subunits to the apicoplast in blood stage parasites and demonstrated that recombinant GatAB converts Glu-tRNA(Gln) to Gln-tRNA(Gln) in vitro. We demonstrate that the apicoplast GatAB-catalyzed reaction is essential to the parasite blood stages because we could not delete the Plasmodium berghei gene encoding GatA in blood stage parasites in vivo. A phylogenetic analysis placed the split between Plasmodium GatB, archaeal GatE, and bacterial GatB prior to the phylogenetic divide between bacteria and archaea. Moreover, Plasmodium GatA also appears to have emerged prior to the bacterial-archaeal phylogenetic divide. Thus, although GatAB is found in Plasmodium, it emerged prior to the phylogenetic separation of archaea and bacteria.

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.

A barcode of organellar genome polymorphisms identifies the geographic origin of Plasmodium falciparum strains

Malaria is a major public health problem that is actively being addressed in a global eradication campaign. Increased population mobility through international air travel has elevated the risk of re-introducing parasites to elimination areas and dispersing drug-resistant parasites to new regions. A simple genetic marker that quickly and accurately identifies the geographic origin of infections would be a valuable public health tool for locating the source of imported outbreaks. Here we analyse the mitochondrion and apicoplast genomes of 711 Plasmodium falciparum isolates from 14 countries, and find evidence that they are non-recombining and co-inherited. The high degree of linkage produces a panel of relatively few single-nucleotide polymorphisms (SNPs) that is geographically informative. We design a 23-SNP barcode that is highly predictive (~92%) and easily adapted to aid case management in the field and survey parasite migration worldwide.

Tracing the source of malarial infections is an important step towards monitoring and controlling the disease. Here, Preston et al. analyse sequence data from 711 isolates and design a genetic barcode based on combined mitochondrial and apicoplast genomes that is able to distinguish between malaria parasites isolated from different geographical regions.

Transcriptomic analysis reveals evidence for a cryptic plastid in the colpodellid Voromonas pontica, a close relative of chromerids and apicomplexan parasites

Colpodellids are free-living, predatory flagellates, but their close relationship to photosynthetic chromerids and plastid-bearing apicomplexan parasites suggests they were ancestrally photosynthetic. Colpodellids may therefore retain a cryptic plastid, or they may have lost their plastids entirely, like the apicomplexan Cryptosporidium. To find out, we generated transcriptomic data from Voromonas pontica ATCC 50640 and searched for homologs of genes encoding proteins known to function in the apicoplast, the non-photosynthetic plastid of apicomplexans. We found candidate genes from multiple plastid-associated pathways including iron-sulfur cluster assembly, isoprenoid biosynthesis, and tetrapyrrole biosynthesis, along with a plastid-type phosphate transporter gene. Four of these sequences include the 5' end of the coding region and are predicted to encode a signal peptide and a transit peptide-like region. This is highly suggestive of targeting to a cryptic plastid. We also performed a taxon-rich phylogenetic analysis of small subunit ribosomal RNA sequences from colpodellids and their relatives, which suggests that photosynthesis was lost more than once in colpodellids, and independently in V. pontica and apicomplexans. Colpodellids therefore represent a valuable source of comparative data for understanding the process of plastid reduction in humanity's most deadly parasite.

Ycf93 (Orf105), a small apicoplast-encoded membrane protein in the relict plastid of the malaria parasite Plasmodium falciparum that is conserved in Apicomplexa

Malaria parasites retain a relict plastid (apicoplast) from a photosynthetic ancestor shared with dinoflagellate algae. The apicoplast is a useful drug target; blocking housekeeping pathways such as genome replication and translation in the organelle kills parasites and protects against malaria. The apicoplast of Plasmodium falciparum encodes 30 proteins and a suite of rRNAs and tRNAs that facilitate their expression. orf105 is a hypothetical apicoplast gene that would encode a small protein (PfOrf105) with a predicted C-terminal transmembrane domain. We produced antisera to a predicted peptide within PfOrf105. Western blot analysis confirmed expression of orf105 and immunofluorescence localised the gene product to the apicoplast. Pforf105 encodes a membrane protein that has an apparent mass of 17.5 kDa and undergoes substantial turnover during the 48-hour asexual life cycle of the parasite in blood stages. The effect of actinonin, an antimalarial with a putative impact on post-translational modification of apicoplast proteins like PfOrf105, was examined. Unlike other drugs perturbing apicoplast housekeeping that induce delayed death, actinonin kills parasites immediately and has an identical drug exposure phenotype to the isopentenyl diphosphate synthesis blocker fosmidomycin. Open reading frames of similar size to PfOrf105, which also have predicted C-terminal trans membrane domains, occur in syntenic positions in all sequenced apicoplast genomes from Phylum Apicomplexa. We therefore propose to name these genes ycf93 (hypothetical chloroplast reading frame 93) according to plastid gene nomenclature convention for conserved proteins of unknown function.

The apicoplast genome of Leucocytozoon caulleryi, a pathogenic apicomplexan parasite of the chicken

Leucocytozoon caulleryi, a haemosporidian parasite of the chicken (Gallus gallus domesticus), can be highly pathogenic and often fatal. Although this parasite is extremely relevant to veterinary science, knowledge of its genomic features is limited. To gain information applicable to developing novel control methods for the parasite, we analyzed the apicoplast genome of L. caulleryi. This extranuclear organellar DNA of 85.1 % A + T and a unit of 34,779 bp was found to encode almost the same set of genes as the plastid genome of Plasmodium falciparum, including 16 tRNA and 30 protein coding genes, and except for one open reading frame, ORF91 absent in L. caulleryi. As in P. falciparum, the L. caulleryi apicoplast DNA contains two sets of a unique inverted repeat (IR), each one 5,253 bp and encoding genes specifying one large and one small rRNA subunit and nine tRNAs but no protein, and separated by a unique 13 bp sequence. Studies of several haemosporidian apicoplast DNA sequences have identified a corresponding IR region; however, none of these studies has looked at the complete sequence, even for well-studied species such as P. falciparum. Phylogenetic studies using a concatenated amino acid sequence based on the open reading frames confirmed the close relationship between L. caulleryi and Plasmodium spp. In this study, we determined the nucleotide sequence of the entire L. caulleryi apicoplast genome, including the region connecting the two IR units. This is the first report of the complete nucleotide sequence of a haemosporidian apicoplast DNA with a canonical IR.

The unique structure of the apicoplast genome of the rodent malaria parasite Plasmodium chabaudi chabaudi

The apicoplast, a non-photosynthetic plastid of apicomplexan species, has an extremely reduced but highly conserved genome. Here, the apicoplast genome of the rodent malaria parasite Plasmodium chabaudi chabaudi (Pcc) isolate CB was characterized. Although the set of genes in the genome is identical, the copy number of some tRNA genes differs between Pcc and other Plasmodium species because the Pcc DNA has only one rRNA/tRNA gene cluster, which is normally duplicated in other species. The location of the duplicated trnR(ACG) and trnM implies that one of the duplicated clusters in the ancestral molecule has been lost due to an intramolecular recombination event. The Pcc DNA occurs in two isoforms with an internal inversion between them. The presence of a unique variant in the duplicated trnT gene suggests that the two isoforms are interconvertible. This is the first report of the complete nucleotide sequence of a Plasmodium apicoplast DNA.

Diversity of mitochondrial genome structure in the phylum Apicomplexa

Mitochondria are ubiquitous organelles in all eukaryotes that are essential for a range of cellular processes and cellular signaling. Nearly all mitochondria have their own DNA or mitochondrial (mt) genome, which varies considerably in size, structure and organization. The phylum Apicomplexa includes a variety of unicellular eukaryotes, some of which are parasites of clinical or economic importance. Recent studies have demonstrated that apicomplexan mt genomes, which include the smallest 6 kb genome of the malaria parasites, exhibit remarkably diverse structures. Apicomplexan parasites are interesting model organisms in order to understand the evolution of mt genomes. This review summarizes the structure of apicomplexan mt genomes and highlights the unique features and the evolution of the mt genome.

Effective preparation of Plasmodium vivax field isolates for high-throughput whole genome sequencing

Whole genome sequencing (WGS) of Plasmodium vivax is problematic due to the reliance on clinical isolates which are generally low in parasitaemia and sample volume. Furthermore, clinical isolates contain a significant contaminating background of host DNA which confounds efforts to map short read sequence of the target P. vivax DNA. Here, we discuss a methodology to significantly improve the success of P. vivax WGS on natural (non-adapted) patient isolates. Using 37 patient isolates from Indonesia, Thailand, and travellers, we assessed the application of CF11-based white blood cell filtration alone and in combination with short term ex vivo schizont maturation. Although CF11 filtration reduced human DNA contamination in 8 Indonesian isolates tested, additional short-term culture increased the P. vivax DNA yield from a median of 0.15 to 6.2 ng µl(-1) packed red blood cells (pRBCs) (p = 0.001) and reduced the human DNA percentage from a median of 33.9% to 6.22% (p = 0.008). Furthermore, post-CF11 and culture samples from Thailand gave a median P. vivax DNA yield of 2.34 ng µl(-1) pRBCs, and 2.65% human DNA. In 22 P. vivax patient isolates prepared with the 2-step method, we demonstrate high depth (median 654X coverage) and breadth (≥89%) of coverage on the Illumina GAII and HiSeq platforms. In contrast to the A+T-rich P. falciparum genome, negligible bias was observed in coverage depth between coding and non-coding regions of the P. vivax genome. This uniform coverage will greatly facilitate the detection of SNPs and copy number variants across the genome, enabling unbiased exploration of the natural diversity in P. vivax populations.