Highly conserved gene arrangement of the mitochondrial genomes of 23 Plasmodium species

Hepatozoon (Eucoccidiorida: Hepatozoidae) in wild mammals of the Americas: a systematic review

BACKGROUND

The study of parasites provides insight into intricate ecological relationships in ecosystem dynamics, food web structures, and evolution on multiple scales. Hepatozoon Eucoccidiorida: Hepatozoidae) is a genus of protozoan hemoparasites with heteroxenous life cycles that switch infections between vertebrates and blood-feeding invertebrates. The most comprehensive review of the genus was published 26 years ago, and currently there are no harmonized data on the epizootiology, diagnostics, genotyping methods, evolutionary relationships, and genetic diversity of Hepatozoon in the Americas.

METHODS

Here, we provide a comprehensive review based on the PRISMA method regarding Hepatozoon in wild mammals within the American continent, in order to generate a framework for future research.

RESULTS

11 out of the 35 countries of the Americas (31.4%) had data on Hepatozoon, with Carnivora and Rodentia orders having the most characterizations. Bats, ungulates, and shrews were the least affected groups. While Hepatozoon americanum, H. americanum-like, H. canis, H. didelphydis, H. felis, H. milleri, H. griseisciuri, and H. procyonis correspond to the identified species, a plethora of genospecies is pending for a formal description combining morphology and genetics. Most of the vectors of Hepatozoon in the Americas are unknown, but some flea, mite, and tick species have been confirmed. The detection of Hepatozoon has relied mostly on conventional polymerase chain reaction (PCR), and the implementation of specific real time PCR for the genus needs to be employed to improve its diagnosis in wild animals in the future. From a genetic perspective, the V4 region of the 18S rRNA gene has been widely sequenced for the identification of Hepatozoon in wild animals. However, mitochondrial and apicoplast markers should also be targeted to truly determine different species in the genus. A phylogenetic analysis of herein retrieved 18S ribosomal DNA (rDNA) sequences showed two main clades of Hepatozoon: Clade I associated with small mammals, birds, and herpetozoa, and Clade II associated with Carnivora. The topology of the tree is also reflected in the haplotype network.

CONCLUSIONS

Finally, our review emphasizes Hepatozoon as a potential disease agent in threatened wild mammals and the role of wild canids as spreaders of Hepatozoon infections in the Americas.

Development of Mitochondria Targeting AIE-Active Cyclometalated Iridium Complexes as Potent Antimalarial Agents

The emergence of resistance to conventional antimalarial treatments remains a major cause for concern. New drugs that target the distinct development stages of Plasmodium parasites are required to address this risk. Herein, water-soluble aggregation-induced emission active cyclometalated iridium(III) polypyridyl complexes (Ir1-Ir12) are developed for the elimination of malaria parasites. Remarkably, these complexes show potent antimalarial activity in low nanomolar range against 3D7 (chloroquine and artemisinin sensitive strain), RKL9 (chloroquine resistant strain), and R539T (artemisinin resistant strains) strains of Plasmodium falciparum with faster killing rate of malaria parasites. Concomitantly, these complexes exhibit efficient in vivo antimalarial activity against both the asexual and gametocyte stages of Plasmodium berghei malaria parasite, suggesting promising transmission-blocking potential. The complexes tend to localize into mitochondria of P. falciparum determined by image and cell-based assay. The mechanistic studies reveal that these complexes exert their antimalarial activity by increasing reactive oxygen species levels and disrupting its mitochondrial membrane potential. Furthermore, the mitochondrial-dependent antimalarial activity of these complexes is confirmed in yeast model. Thus, this study for the first time highlights the potential role of targeting P. falciparum mitochondria by iridium complexes in discovering and developing the next-generation antimalarial agents for treating multidrug resistant malaria parasites.

Inheritance through the cytoplasm

Most heritable information in eukaryotic cells is encoded in the nuclear genome, with inheritance patterns following classic Mendelian segregation. Genomes residing in the cytoplasm, however, prove to be a peculiar exception to this rule. Cytoplasmic genetic elements are generally maternally inherited, although there are several exceptions where these are paternally, biparentally or doubly-uniparentally inherited. In this review, we examine the diversity and peculiarities of cytoplasmically inherited genomes, and the broad evolutionary consequences that non-Mendelian inheritance brings. We first explore the origins of vertical transmission and uniparental inheritance, before detailing the vast diversity of cytoplasmic inheritance systems across Eukaryota. We then describe the evolution of genomic organisation across lineages, how this process has been shaped by interactions with the nuclear genome and population genetics dynamics. Finally, we discuss how both nuclear and cytoplasmic genomes have evolved to co-inhabit the same host cell via one of the longest symbiotic processes, and all the opportunities for intergenomic conflict that arise due to divergence in inheritance patterns. In sum, we cannot understand the evolution of eukaryotes without understanding hereditary symbiosis.

Functional genomics of RAP proteins and their role in mitoribosome regulation in Plasmodium falciparum

The RAP (RNA-binding domain abundant in Apicomplexans) protein family has been identified in various organisms. Despite expansion of this protein family in apicomplexan parasites, their main biological functions remain unknown. In this study, we use inducible knockdown studies in the human malaria parasite, Plasmodium falciparum, to show that two RAP proteins, PF3D7_0105200 (PfRAP01) and PF3D7_1470600 (PfRAP21), are essential for parasite survival and localize to the mitochondrion. Using transcriptomics, metabolomics, and proteomics profiling experiments, we further demonstrate that these RAP proteins are involved in mitochondrial RNA metabolism. Using high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (eCLIP-seq), we validate that PfRAP01 and PfRAP21 are true RNA-binding proteins and interact specifically with mitochondrial rRNAs. Finally, mitochondrial enrichment experiments followed by deep sequencing of small RNAs demonstrate that PfRAP21 controls mitochondrial rRNA expression. Collectively, our results establish the role of these RAP proteins in mitoribosome activity and contribute to further understanding this protein family in malaria parasites.

A Mitochondria-Penetrating Peptide Exerts Potent Anti-Plasmodium Activity and Localizes at Parasites' Mitochondria

Mitochondria are considered a novel drug target as they play a key role in energy production and programmed cell death of eukaryotic cells. The mitochondria of malaria parasites differ from those of their vertebrate hosts, contributing to the drug selectivity and the development of antimalarial drugs. (F_xr)₃, a mitochondria-penetrating peptide or MPP, entered malaria-infected red cells without disrupting the membrane and subsequently killed the blood stage of P. falciparum parasites. The effects were more potent on the late stages than on the younger stages. Confocal microscopy showed that the (F_xr)₃ intensely localized at the parasite mitochondria. (F_xr)₃ highly affected both the lab-strain, chloroquine-resistant K1, and freshly isolated malaria parasites. (F_xr)₃ (1 ng/mL to 10 μg/mL) was rarely toxic towards various mammalian cells, i.e., mouse fibroblasts (L929), human leukocytes and erythrocytes. At a thousand times higher concentration (100 μg/mL) than that of the antimalarial activity, cytotoxicity and hemolytic activity of (F_xr)₃ were observed. Compared with the known antimalarial drug, atovaquone, (F_xr)₃ exhibited more rapid killing activity. This is the first report on antimalarial activity of (F_xr)₃, showing localization at parasites’ mitochondria.

The Elusive Mitochondrial Genomes of Apicomplexa: Where Are We Now?

Mitochondria are vital organelles of eukaryotic cells, participating in key metabolic pathways such as cellular respiration, thermogenesis, maintenance of cellular redox potential, calcium homeostasis, cell signaling, and cell death. The phylum Apicomplexa is entirely composed of obligate intracellular parasites, causing a plethora of severe diseases in humans, wild and domestic animals. These pathogens include the causative agents of malaria, cryptosporidiosis, neosporosis, East Coast fever and toxoplasmosis, among others. The mitochondria in Apicomplexa has been put forward as a promising source of undiscovered drug targets, and it has been validated as the target of atovaquone, a drug currently used in the clinic to counter malaria. Apicomplexans present a single tubular mitochondria that varies widely both in structure and in genomic content across the phylum. The organelle is characterized by massive gene migrations to the nucleus, sequence rearrangements and drastic functional reductions in some species. Recent third generation sequencing studies have reignited an interest for elucidating the extensive diversity displayed by the mitochondrial genomes of apicomplexans and their intriguing genomic features. The underlying mechanisms of gene transcription and translation are also ill-understood. In this review, we present the state of the art on mitochondrial genome structure, composition and organization in the apicomplexan phylum revisiting topological and biochemical information gathered through classical techniques. We contextualize this in light of the genomic insight gained by second and, more recently, third generation sequencing technologies. We discuss the mitochondrial genomic and mechanistic features found in evolutionarily related alveolates, and discuss the common and distinct origins of the apicomplexan mitochondria peculiarities.

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.

Identification and phylogenetic analysis of RNA binding domain abundant in apicomplexans or RAP proteins

The RNA binding domain abundant in apicomplexans (RAP) is a protein domain identified in a diverse group of proteins, called RAP proteins, many of which have been shown to be involved in RNA binding. To understand the expansion and potential function of the RAP proteins, we conducted a hidden Markov model based screen among the proteomes of 54 eukaryotes, 17 bacteria and 12 archaea. We demonstrated that the domain is present in closely and distantly related organisms with particular expansions in Alveolata and Chlorophyta, and are not unique to Apicomplexa as previously believed. All RAP proteins identified can be decomposed into two parts. In the N-terminal region, the presence of variable helical repeats seems to participate in the specific targeting of diverse RNAs, while the RAP domain is mostly identified in the C-terminal region and is highly conserved across the different phylogenetic groups studied. Several conserved residues defining the signature motif could be crucial to ensure the function(s) of the RAP proteins. Modelling of RAP domains in apicomplexan parasites confirmed an ⍺/β structure of a restriction endonuclease-like fold. The phylogenetic trees generated from multiple alignment of RAP domains and full-length proteins from various distantly related eukaryotes indicated a complex evolutionary history of this family. We further discuss these results to assess the potential function of this protein family in apicomplexan parasites.

Multilocus sequencing of Hepatozoon cf. griseisciuri infections in Ontario eastern gray squirrels (Sciurus carolinensis) uncovers two genotypically distinct sympatric parasite species

Intra-leukocytic gamonts consistent with the description of Hepatozoon griseisciuri Clark, 1958 are reported for the first time in Canadian eastern gray squirrels (Sciurus carolinensis Gmelin, 1788). Polymerase chain reaction (PCR) amplification and direct Sanger sequencing identified a pair of distinct genotypes at both a nuclear and mitochondrial locus; two 18S ribosomal RNA gene sequences (rDNA; genotype A and genotype B: 1816 base pairs (bp); 98.8% pairwise identity) and 2 distinct complete mitochondrial genome sequences (genotype A: 6311 bp; genotype B: 6114 bp; 89.1% pairwise identity) were obtained from 3 H. griseisciuri-infected squirrels sampled in Guelph, Ontario. The genetic content of both circular-mapping mitochondrial genomes was conventional for apicomplexan protists; each encoded for 3 protein-coding genes (cytochrome c oxidase subunit I (COI); cytochrome c oxidase subunit III (COIII); and cytochrome B (CytB)), 14 fragmented large subunit rDNA, 10 fragmented small subunit rDNA, and 8 unassigned rDNA. These genotypes, based on sequences obtained from a pair of loci from two parasite genomes, confirm the presence of at least two Hepatozoon species infecting Ontario eastern gray squirrels, one of which is likely to be conspecific with H. griseisciuri.

The Identification of Scientific Communities and Their Approach to Worldwide Malaria Research

It is essential to establish a pattern to detect the strengths and weaknesses of working groups publishing on malaria, to promote coordination to facilitate the eradication of the disease. Given the complexity of the scientific network of groups and institutions studying malaria, it is necessary to use a mathematical algorithm that allows us to know the real structure of research on the disease in the world. In this work, articles with the word “malaria” in the title or author keywords gathered from Elsevier Scopus database were analyzed. By means of specific software, graphs were created. The analysis of the data allowed established different scientific communities, among which two were very diverse: one formed by those groups concerned about the vector transmission and control, and another one focused on the drug resistance of the parasite. Basic, applied, and operational research to eradicate malaria is an ambitious goal of the international institutions and the scientific community. The combination of effort and the establishment of a worldwide-scientific network that allows an effective interconnection (exchange) of knowledge, infrastructure technology, collaborators, financial resources, and datasets will contribute more effectively to end the disease.

Babesia gibsoni endemic to Wuhan, China: mitochondrial genome sequencing, annotation, and comparison with apicomplexan parasites

Babesia gibsoni (B. gibsoni), an intracellular apicomplexan protozoan, poses great threat to canine health. Currently, little information is available about the B. gibsoni (WH58) endemic to Wuhan, China. Here, the mitochondrial (mt) genome of B. gibsoni (WH58) was amplified by five pairs of primers and sequenced and annotated by alignment with the reported mt genome sequences of Babesia canis (B. canis, KC207822), Babesia orientalis (KF218819), Babesia bovis (AB499088), and Theileria equi (AB499091). The evolutionary relationships were analyzed with the amino acid sequences of cytochrome c oxidase I (cox1) and cytochrome b (cob) genes in apicomplexan parasite species. Additionally, the mt genomes of Babesia, Theileria, and Plasmodium spp. were compared in size, host infection, form, distribution, and direction of the protein-coding genes. The full size of the mt genome of B. gibsoni (WH58) was 5865 bp with a linear form, containing terminal-inverted repeats on both ends, six large subunit ribosomal RNA fragments, and three protein-coding genes: cox1, cob, and cytochrome c oxidase III (cox3). Babesia, Theileria, and Plasmodium spp. had a similar mt genome size of about 6000 bp. The mt genomes of parasites that cause canine babesiosis showed a slightly smaller size than the other species. Moreover, Babesia microti (R1 strain) was about 11,100 bp in size, which was twice larger than that of the other species. The mt form was linear for Babesia and Theileria spp. but circular for Plasmodium falciparum and Plasmodium knowlesi. Additionally, all the species contained the three protein-coding genes of cox1, cox3, and cob except Toxoplasma gondii (RH strain) which only contained the cox1 and cob genes. The phylogenetic analysis indicated that B. gibsoni (WH58) was more identical to B. gibsoni (AB499087), B. canis (KC207822), and Babesia rossi (KC207823) and most divergent from Babesia conradae in Babesia spp. Despite the highest similarity to B. gibsoni (AB499087) reported in Japan, B. gibsoni (WH58) showed notable differences in the sequence of nucleotides and amino acids and the property in virulence to host and in vitro cultivation. This study compared the mt genomes of the two B. gibsoni isolates and other parasites in the phylum Apicomplexa and provided new insights into their differences and evolutionary relationships.

Plasmodium genomics: an approach for learning about and ending human malaria

Malaria causes high levels of morbidity and mortality in human beings worldwide. According to the World Health Organization (WHO), about half a million people die of this disease each year. Malaria is caused by six species of parasites belonging to the Plasmodium genus: P. falciparum, P. knowlesi, P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. Currently, malaria is being kept under control with varying levels of elimination success in different countries. The development of new molecular tools as well as the use of next-generation sequencing (NGS) technologies and novel bioinformatic approaches has improved our knowledge of malarial epidemiology, diagnosis, treatment, vaccine development, and surveillance strategies. In this work, the genetics and genomics of human malarias have been analyzed. Since the first P. falciparum genome was sequenced in 2002, various population-level genetic and genomic surveys, together with transcriptomic and proteomic studies, have shown the importance of molecular approaches in supporting malaria elimination.

PCR-based detection of Plasmodium falciparum in saliva using mitochondrial cox3 and varATS primers

Background

Sampling of saliva for diagnosing Plasmodium falciparum infections is a safe, non-invasive alternative to sampling of blood. However, the use of saliva presents a challenge because lower concentrations of parasite DNA are present in saliva compared to peripheral blood. Therefore, a sensitive method is needed for detection of parasite DNA in saliva. This study utilized two recently reported “ultra-sensitive” PCR assays based on detection of the P. falciparum mitochondrial cox3 gene and the multi-copy nuclear varATS gene. The ultra-sensitive assays have an advantage over standard 18S rRNA gene-based PCR assay as they target genes with higher copy numbers per parasite genome. Stored saliva DNA samples from 60 Cameroonian individuals with infections previously confirmed by 18S rRNA gene PCR in peripheral blood were tested with assays targeting the cox3 and varATS genes.

Results

Overall, the standard 18S rRNA gene-based PCR assay detected P. falciparum DNA in 62% of the stored saliva DNA samples, whereas 77 and 68% of the samples were positive with assays that target the cox3 and varATS genes, respectively. Interestingly, the ultra-sensitive assays detected more P. falciparum infections in stored saliva samples than were originally detected by thick-film microscopy (41/60 = 68%). When stratified by number of parasites in the blood, the cox3 assay successfully detected more than 90% of infections using saliva when individuals had > 1000 parasites/μl of peripheral blood, but sensitivity was reduced at submicroscopic parasitemia levels. Bands on electrophoresis gels were distinct for the cox3 assay, whereas faint or non-specific bands were sometimes observed for varATS and 18S rRNA that made interpretation of results difficult. Assays could be completed in 3.5 and 3 h for the cox3 and varATS assays, respectively, whereas the 18S rRNA gene assays required at least 7 h.

Conclusions

This study demonstrates that a PCR assay targeting the cox3 gene detected P. falciparum DNA in more saliva samples than primers for the 18S rRNA gene. Non-invasive collection of saliva in combination with the proposed cox3 primer-based PCR assay could potentially enhance routine testing of P. falciparum during disease surveillance, monitoring, and evaluation of interventions for malaria elimination.

Electronic supplementary material

The online version of this article (10.1186/s41182-018-0100-2) contains supplementary material, which is available to authorized users.

Rapid and sensitive multiplex single-tube nested PCR for the identification of five human Plasmodium species

Malaria is caused by five species of Plasmodium in humans. Microscopy is currently used for pathogen detection, requiring considerable training and technical expertise as the parasites are often difficult to differentiate morphologically. Rapid diagnostic tests are as reliable as microscopy and offer faster diagnoses but possess lower detection limits and are incapable of distinguishing among the parasitic species. To improve global health efforts towards malaria control, a rapid, sensitive, species-specific, and economically viable diagnostic method is needed. In this study, we designed a malaria diagnostic method involving a multiplex single-tube nested PCR targeting Plasmodium mitochondrial cytochrome c oxidase III and single-stranded tag hybridization chromatographic printed-array strip. The detection sensitivity was found to be at least 40 times higher than that of agarose gel electrophoresis with ethidium bromide. This system also enables the identification of both single- and mixed-species malaria infections. The assay was validated with 152 Kenyan samples; using nested PCR as the standard, the assay's sensitivity and specificity were 88.7% and 100.0%, respectively. The turnaround time required, from PCR preparation to signal detection, is 90min. Our method should improve the diagnostic speed, treatment efficacy, and control of malaria, in addition to facilitating surveillance within global malaria eradication programs.

Utility of the cytochrome c oxidase subunit I gene for the diagnosis of toxoplasmosis using PCR

Toxoplasmosis is caused by the protozoan parasite Toxoplasma gondii, which belongs to the phylum Apicomplexa. Since this parasite causes severe clinical symptoms in immunocompromised patients, early diagnosis of toxoplasmosis is essential. PCR is currently used for early diagnosis, but there is no consensus regarding the most effective method for amplifying Toxoplasma DNA. In this study, we considered the utility of the cytochrome c subunit I (cox1) gene, which is encoded in the mitochondrial DNA of this parasite, as a novel target of PCR for the diagnosis of toxoplasmosis. To do this, we compared its copy number per haploid nuclear genome and the detection sensitivity of cox1-PCR with the previously reported target genes B1 and 18S rRNA and the AF146527 repeat element. We found that the copy number of cox1 was high and that the PCR using cox1 primers was more efficient at amplifying Toxoplasma DNA than the other PCR targets examined. In addition, PCR using clinical samples indicated that the cox1 gene would be useful for the diagnosis of toxoplasmosis. These findings suggest that use of cox1-PCR would facilitate the diagnosis of toxoplasmosis in clinical laboratories.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Mitochondrial genome sequence diversity of Indian Plasmodium falciparum isolates

We have analysed the whole mitochondrial (mt) genome sequences (each ~6 kilo nucleotide base pairs in length) of four field isolates of the malaria parasite Plasmodium falciparum collected from different locations in India. Comparative genomic analyses of mt genome sequences revealed three novel India-specific single nucleotide polymorphisms. In general, high mt genome diversity was found in Indian P. falciparum, at a level comparable to African isolates. A population phylogenetic tree placed the presently sequenced Indian P. falciparum with the global isolates, while a previously sequenced Indian isolate was an outlier. Although this preliminary study is limited to a few numbers of isolates, the data have provided fundamental evidence of the mt genome diversity and evolutionary relationships of Indian P. falciparum with that of global isolates.

Mitochondrial population genomic analyses reveal population structure and demography of Indian Plasmodium falciparum

Inference on the genetic diversity of Plasmodium falciparum populations could help in better management of malaria. A very recent study with mitochondrial (mt) genomes in global P. falciparum had revealed interesting evolutionary genetic patterns of Indian isolates in comparison to global ones. However, no population genetic study using the whole mt genome sequences of P. falciparum isolates collected in the entire distribution range in India has yet been performed. We herewith have analyzed 85 whole mt genomes (48 already published and 37 entirely new) sampled from eight differentially endemic Indian locations to estimate genetic diversity and infer population structure and historical demography of Indian P. falciparum. We found 19 novel Indian-specific Single Nucleotide Polymorphisms (SNPs) and 22 novel haplotypes segregating in Indian P. falciparum. Accordingly, high haplotype and nucleotide diversities were detected in Indian P. falciparum in comparison to many other global isolates. Indian P. falciparum populations were found to be moderately sub-structured with four different genetic clusters. Interestingly, group of local populations aggregate to form each cluster; while samples from Jharkhand and Odisha formed a single cluster, P. falciparum isolates from Asom formed an independent one. Similarly, Surat, Bilaspur and Betul formed a single cluster and Goa and Mangalore formed another. Interestingly, P. falciparum isolates from the two later populations were significantly genetically differentiated from isolates collected in other six Indian locations. Signature of historical population expansion was evident in five population samples, and the onset of expansion event was found to be very similar to African P. falciparum. In agreement with the previous finding, the estimated Time to Most Recent Common Ancestor (TMRCA) and the effective population size were high in Indian P. falciparum. All these genetic features of Indian P. falciparum with high mt genome diversity are somehow similar to Africa, but quite different from other Asian population samples.

Improved detection of malaria cases in island settings of Vanuatu and Kenya by PCR that targets the Plasmodium mitochondrial cytochrome c oxidase III (cox3) gene

Detection of sub-microscopic parasitemia is crucial for all malaria elimination programs. PCR-based methods have proven to be sensitive, but two rounds of amplification (nested PCR) are often needed to detect the presence of Plasmodium DNA. To simplify the detection process, we designed a nested PCR method whereby only the primary PCR is required for the detection of the four major human Plasmodium species. Primers designed for the detection of the fifth species, Plasmodium knowlesi, were not included in this study due to the absence of appropriate field samples. Compared to the standard 18S rDNA PCR method, our cytochrome c oxidase III (cox3) method detected 10-50% more cases while maintaining high sensitivities (1.00) for all four Plasmodium species in our samples from Vanuatu (n=77) and Kenya (n=76). Improvement in detection efficiency was more substantial for samples with sub-microscopic parasitemia (54%) than those with observable parasitemia (10-16%). Our method will contribute to improved malaria surveillance in low endemicity settings.

New insights into the evolutionary history of Plasmodium falciparum from mitochondrial genome sequence analyses of Indian isolates

Estimating genetic diversity and inferring the evolutionary history of Plasmodium falciparum could be helpful in understanding origin and spread of virulent and drug-resistant forms of the malaria pathogen and therefore contribute to malaria control programme. Genetic diversity of the whole mitochondrial (mt) genome of P. falciparum sampled across the major distribution ranges had been reported, but no Indian P. falciparum isolate had been analysed so far, even though India is highly endemic to P. falciparum malaria. We have sequenced the whole mt genome of 44 Indian field isolates and utilized published data set of 96 genome sequences to present global genetic diversity and to revisit the evolutionary history of P. falciparum. Indian P. falciparum presents high genetic diversity with several characteristics of ancestral populations and shares many of the genetic features with African and to some extent Papua New Guinean (PNG) isolates. Similar to African isolates, Indian P. falciparum populations have maintained high effective population size and undergone rapid expansion in the past with oldest time to the most recent common ancestor (TMRCA). Interestingly, one of the four single nucleotide polymorphisms (SNPs) that differentiates P. falciparum from P. falciparum-like isolates (infecting non-human primates in Africa) was found to be segregating in five Indian P. falciparum isolates. This SNP was in tight linkage with other two novel SNPs that were found exclusively in these five Indian isolates. The results on the mt genome sequence analyses of Indian isolates on the whole add to the current understanding on the evolutionary history of P. falciparum.

The complete mitochondrial genome sequence of Hepatozoon catesbianae (Apicomplexa: Coccidia: Adeleorina), a blood parasite of the green frog, Lithobates (formerly Rana) clamitans

A complete mitochondrial genome for the blood parasite Hepatozoon catesbianae (Alveolata; Apicomplexa; Coccidia; Adeleorina; Hepatozoidae) was obtained through PCR amplification and direct sequencing of resulting PCR products. The mitochondrial genome of H. catesbianae is 6,397 bp in length and contains 3 protein-coding genes (cytochrome c oxidase subunit I [COI]; cytochrome c oxidase subunit III [COIII]; and cytochrome B [CytB]). Sequence similarities to previously published mitochondrial genomes of other apicomplexan parasites permitted annotation of 23 putative rDNA fragments in the mitochondrial genome of H. catesbianae, 14 large subunit rDNA fragments, and 9 small subunit rDNA fragments. Sequences corresponding to rDNA fragments RNA5, RNA8, RNA11, and RNA19 of Plasmodium falciparum were not identified in the mitrochondrial genome sequence of H. catesbianae. Although the presence of 3 protein-coding regions and numerous putative rDNA fragments is a feature typical for apicomplexan mitochondrial genomes, the mitochondrial genome of H. catesbianae possesses a structure and gene organization that is distinct among the Apicomplexa. This is the first complete mitochondrial genome sequence obtained from any apicomplexan parasite in the suborder Adeleorina.

Mitochondrial genome of Babesia orientalis, apicomplexan parasite of water buffalo (Bubalus babalis, Linnaeus, 1758) endemic in China

BACKGROUND

Apicomplexan parasites of the genus Babesia, Theileria and Plasmodium are very closely related organisms. Interestingly, their mitochondrial (mt) genomes are highly divergent. Among Babesia, Babesia orientalis is a new species recently identified and specifically epidemic to the southern part of China, causing severe disease to water buffalo. However, no information on the mt genome of B. orientalis was available.

METHODS

Four pairs of primers were designed based on the full genome sequence of B. orientalis (unpublished data) and by aligning reported mt genomes of B. bovis, B. bigemina, and T. parva. The entire mt genome was amplified by four sets of PCR. The obtained mt genome was annotated by aligning with published apicomplexan mt genomes and Artemis software v11. Phylogenetic analysis was performed by using cox1 and cob amino acid sequences.

RESULTS

The complete mt genome of B. orientalis (Wuhan strain) was sequenced and characterized. The entire mt genome is 5996 bp in length with a linear form, containing three protein-coding genes including cytochrome c oxidase I (cox1), cytochrome b (cob) and cytochrome c oxidase III (cox3) and six rRNA large subunit gene fragments. The gene arrangement in B. orientalis mt genome is similar to those of B. bovis, B. gibsoni and Theileria parva, but different from those of T. orientalis, T. equi and Plasmodium falciparum. Comparative analysis indicated that cox1 and cob genes were more conserved than cox3. Phylogenetic analysis based on amino acid sequences of cox1, cob and cox1 + cob, respectively, revealed that B. orientalis fell into Babesia clade with the closest relationship to B. bovis.

CONCLUSIONS

The availability of the entire mt genome sequences of B. orientalis provides valuable information for future phylogenetic, population genetics and molecular epidemiological studies of apicomplexan parasites.

Whole mitochondrial genome sequence of an Indian Plasmodium falciparum field isolate

Mitochondrial genome sequence of malaria parasites has served as a potential marker for inferring evolutionary history of the Plasmodium genus. In Plasmodium falciparum, the mitochondrial genome sequences from around the globe have provided important evolutionary understanding, but no Indian sequence has yet been utilized. We have sequenced the whole mitochondrial genome of a single P. falciparum field isolate from India using novel primers and compared with the 3D7 reference sequence and 1 previously reported Indian sequence. While the 2 Indian sequences were highly divergent from each other, the presently sequenced isolate was highly similar to the reference 3D7 strain.

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.

Whole genome mapping and re-organization of the nuclear and mitochondrial genomes of Babesia microti isolates

Babesia microti is the primary causative agent of human babesiosis, an emerging pathogen that causes a malaria-like illness with possible fatal outcome in immunocompromised patients. The genome sequence of the B. microti R1 strain was reported in 2012 and revealed a distinct evolutionary path for this pathogen relative to that of other apicomplexa. Lacking from the first genome assembly and initial molecular analyses was information about the terminal ends of each chromosome, and both the exact number of chromosomes in the nuclear genome and the organization of the mitochondrial genome remained ambiguous. We have now performed various molecular analyses to characterize the nuclear and mitochondrial genomes of the B. microti R1 and Gray strains and generated high-resolution Whole Genome maps. These analyses show that the genome of B. microti consists of four nuclear chromosomes and a linear mitochondrial genome present in four different structural types. Furthermore, Whole Genome mapping allowed resolution of the chromosomal ends, identification of areas of misassembly in the R1 genome, and genomic differences between the R1 and Gray strains, which occur primarily in the telomeric regions. These studies set the stage for a better understanding of the evolution and diversity of this important human pathogen.

Plasmodium falciparum mitochondrial genetic diversity exhibits isolation-by-distance patterns supporting a sub-Saharan African origin

The geographical distribution of single nucleotide polymorphism (SNP) in the mitochondrial genome of the human malaria parasite Plasmodium falciparum was investigated. We identified 88 SNPs in 516 isolates from seven parasite populations in Africa, Southeast Asia and Oceania. Analysis of the SNPs postulated a sub-Saharan African origin and recovered a strong negative correlation between within-population SNP diversity and geographic distance from the putative African origin over Southeast Asia and Oceania. These results are consistent with those previously obtained for nuclear genome-encoded housekeeping genes, indicating that the pattern of inheritance does not substantially affect the geographical distribution of SNPs.

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.

Characterisation of full-length mitochondrial copies and partial nuclear copies (numts) of the cytochrome b and cytochrome c oxidase subunit I genes of Toxoplasma gondii, Neospora caninum, Hammondia heydorni and Hammondia triffittae (Apicomplexa: Sarcocystidae)

Genomic DNA was extracted from three oocyst isolates of Hammondia triffittae from foxes and two oocyst isolates of Hammondia heydorni from dogs, as well as from cell culture-derived tachyzoites of Toxoplasma gondii (RH strain) and Neospora caninum (NC-Liverpool strain), and examined by PCR with primers targeting the cytochrome b (cytb) and the cytochrome c oxidase subunit I (cox1) genes in order to characterise both genes and, if possible, the remainder of the mitochondrial genome of these species. Several primers were designed and used in various combinations to amplify regions within and between both genes and to determine gene order. When certain forward primers targeting cytb were used in combination with certain reverse primers targeting cox1, two overlapping sequences were obtained for each species and isolate studied, which showed that a full-length copy of cytb was followed 36-37 bp downstream by a full-length copy of cox1, and these sequences are believed to represent the true mitochondrial genes and the gene order in the mitochondrial genome of the four species examined. The cytb of T. gondii, N. caninum, H. heydorni and H. triffittae comprised a total of 1,080 bp (359 amino acids) and used ATG and TAA as start and stop codon, respectively. The cox1 of these species also used TAA as stop codon, whereas the most likely start codon was ATG, resulting in a gene comprising 1,491 bp (496 amino acids). Pair-wise sequence comparisons based on either cytb or cox1 clearly separated T. gondii from N. caninum and both of these species from the two Hammondia species, whereas the latter two species were 100 % identical at cytb and shared 99.3 % identity at cox1. Phylogenetic analyses using the maximum-likelihood method confirmed these findings and placed T. gondii in a clade separate from the three other species and all four Toxoplasmatinae in a sister clade to Eimeria spp. PCR with other primers and/or primer pairs than those used to obtain the full-length mitochondrial genes yielded several types of about 1-1.5 kb long sequences, which comprised stretches of the primer-targeted genes at both ends and an intervening non-coding sequence of various length and composition. Thus, portions of cytb could be found both upstream and downstream from portions of cox1 and portions of the same gene could be found adjacent to each other (cytb→cox1; cox1→cytb; cytb→cytb; cox1→cox1). Sequence comparisons revealed that some of these gene fragments were truncated genes, whereas others included the putative start or stop codon of the full-length mitochondrial genes. From the nature of the gene fragments and/or their flanking sequences, they are assumed to be located on the chromosomes of the nuclear genome and to represent nuclear mitochondrial DNA segments (numts) or pseudogenes. In the four species examined, there were no nucleotide differences between the full-length mitochondrial copies of cytb and cox1 and their various incomplete nuclear counterparts. With a few exceptions, identical numt types and closely similar flanking sequences were obtained for all four species, which would indicate that the original transfer of these mitochondrial genes to the nuclear genome and/or the majority of any subsequent rearrangements of these gene fragments within the nuclear genome happened before the four species diverged. Yet, there were species-specific differences in the nucleotide composition of the nuclear gene fragments, identical to the differences in the mitochondrial genes, which would indicate that the incomplete nuclear copies of cytb and cox1 have been continuously updated during evolution to conform to their mitochondrial parent genes. The PCR-based findings of numts were further supported by Basic Local Alignment Search Tool (BLAST) searches against genome sequences of T. gondii and N. caninum using the concatenated mitochondrial cytb/cox1 sequences as queries. These searches revealed the presence of numerous numts of eighth distinct types in both species, with each one having a fixed starting and end point with respect to the nucleotide positions in the full-length mitochondrial genes. Four numt types were completely homologous between both species, whereas four other types differed with respect to their end point and/or the absence/presence of a 96-bp deletion. Each starting and end point was associated with a unique 100-200-bp long flanking sequence, which further revealed the presence of numts. For both species, the numt types and their various arrangements with respect to each other were identical or similar to those obtained by PCR in all four species examined. None of the identified numts covered a full-length gene, but together, the various numts covered the entire mitochondrial cytb and cox1 genes in an overlapping manner. In addition, they were fairly closely spaced on the chromosomes, and these features may explain why the nuclear copies were preferentially amplified to the exclusion of the true mitochondrial genes with most primers and primer pairs used in the present study. The possibility of a similar high prevalence of numts occurring in the nuclear genome of dinoflagellates is discussed.

Characterization of active ribosomal RNA harboring MITEs insertion in microsporidian Nosema bombycis genome

Microsporidia are a group of obligate intracellular parasites of medical and agricultural importance, which can infect almost all animals, including human beings. Using the genome data of Nosema bombycis, four families of miniature inverted-repeat transposable elements (MITEs) in ribosomal DNA (rDNA) were characterized in the microsporidian N. bombycis and were named LSUME1, ITSME1, SSUME1, and SSUME2, respectively. The genome-wide investigation of these MITEs shows that these MITEs families distribute randomly in N. bombycis genome. All insertion sequences have conserved characteristics of MITEs, the direct repeat sequence and terminal inverted-repeat sequence at both ends of each MITEs sequence. Additionally, using the CLC RNA Workbench Software, secondary structures of rRNA containing MITEs sequence have been predicted and were located in variable region or expansion segment. Furthermore, using two different probes, one is prepared by MITE sequence only (short probe) and the other is prepared by MITE sequence flanking partial rDNA sequence (long probe); northern blotting and dot blotting have been performed to detect the transcriptional and functional activity of the rDNA containing MITEs insertion. Fortunately, we found that the rDNA, which harbors the MITE, not only can be transcripted but also can form a complete ribosome. This is an interesting thing that one gene can keep active even when it has been inserted with another sequence. But the biological and structural significance of this observation is not readily apparent.

Novel type of linear mitochondrial genomes with dual flip-flop inversion system in apicomplexan parasites, Babesia microti and Babesia rodhaini

Background

Mitochondrial (mt) genomes vary considerably in size, structure and gene content. The mt genomes of the phylum Apicomplexa, which includes important human pathogens such as the malaria parasite Plasmodium, also show marked diversity of structure. Plasmodium has a concatenated linear mt genome of the smallest size (6-kb); Babesia and Theileria have a linear monomeric mt genome (6.5-kb to 8.2-kb) with terminal inverted repeats; Eimeria, which is distantly related to Plasmodium and Babesia/Theileria, possesses a mt genome (6.2-kb) with a concatemeric form similar to that of Plasmodium; Cryptosporidium, the earliest branching lineage within the phylum Apicomplexa, has no mt genome. We are interested in the evolutionary origin of linear mt genomes of Babesia/Theileria, and have investigated mt genome structures in members of archaeopiroplasmid, a lineage branched off earlier from Babesia/Theileria.

Results

The complete mt genomes of archaeopiroplasmid parasites, Babesia microti and Babesia rodhaini, were sequenced. The mt genomes of B. microti (11.1-kb) and B. rodhaini (6.9-kb) possess two pairs of unique inverted repeats, IR-A and IR-B. Flip-flop inversions between two IR-As and between two IR-Bs appear to generate four distinct genome structures that are present at an equi-molar ratio. An individual parasite contained multiple mt genome structures, with 20 copies and 2 – 3 copies per haploid nuclear genome in B. microti and B. rodhaini, respectively.

Conclusion

We found a novel linear monomeric mt genome structure of B. microti and B. rhodhaini equipped with dual flip-flop inversion system, by which four distinct genome structures are readily generated. To our knowledge, this study is the first to report the presence of two pairs of distinct IR sequences within a monomeric linear mt genome. The present finding provides insight into further understanding of evolution of mt genome structure.

The fragmented mitochondrial ribosomal RNAs of Plasmodium falciparum

Background

The mitochondrial genome in the human malaria parasite Plasmodium falciparum is most unusual. Over half the genome is composed of the genes for three classic mitochondrial proteins: cytochrome oxidase subunits I and III and apocytochrome b. The remainder encodes numerous small RNAs, ranging in size from 23 to 190 nt. Previous analysis revealed that some of these transcripts have significant sequence identity with highly conserved regions of large and small subunit rRNAs, and can form the expected secondary structures. However, these rRNA fragments are not encoded in linear order; instead, they are intermixed with one another and the protein coding genes, and are coded on both strands of the genome. This unorthodox arrangement hindered the identification of transcripts corresponding to other regions of rRNA that are highly conserved and/or are known to participate directly in protein synthesis.

Principal Findings

The identification of 14 additional small mitochondrial transcripts from P. falcipaurm and the assignment of 27 small RNAs (12 SSU RNAs totaling 804 nt, 15 LSU RNAs totaling 1233 nt) to specific regions of rRNA are supported by multiple lines of evidence. The regions now represented are highly similar to those of the small but contiguous mitochondrial rRNAs of Caenorhabditis elegans. The P. falciparum rRNA fragments cluster on the interfaces of the two ribosomal subunits in the three-dimensional structure of the ribosome.

Significance

All of the rRNA fragments are now presumed to have been identified with experimental methods, and nearly all of these have been mapped onto the SSU and LSU rRNAs. Conversely, all regions of the rRNAs that are known to be directly associated with protein synthesis have been identified in the P. falciparum mitochondrial genome and RNA transcripts. The fragmentation of the rRNA in the P. falciparum mitochondrion is the most extreme example of any rRNA fragmentation discovered.