Concatenated mitochondrial DNA of the coccidian parasite Eimeria tenella

Unique Properties of Apicomplexan Mitochondria

Apicomplexan parasites constitute more than 6,000 species infecting a wide range of hosts. These include important pathogens such as those causing malaria and toxoplasmosis. Their evolutionary emergence coincided with the dawn of animals. Mitochondrial genomes of apicomplexan parasites have undergone dramatic reduction in their coding capacity, with genes for only three proteins and ribosomal RNA genes present in scrambled fragments originating from both strands. Different branches of the apicomplexans have undergone rearrangements of these genes, with Toxoplasma having massive variations in gene arrangements spread over multiple copies. The vast evolutionary distance between the parasite and the host mitochondria has been exploited for the development of antiparasitic drugs, especially those used to treat malaria, wherein inhibition of the parasite mitochondrial respiratory chain is selectively targeted with little toxicity to the host mitochondria. We describe additional unique characteristics of the parasite mitochondria that are being investigated and provide greater insights into these deep-branching eukaryotic pathogens.

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.

The complete mitochondrial genome sequence of Eimeria leuckarti (Eimeriidae, Coccidia, Apicomplexa) infecting domestic horses (Equus ferus caballus)

The complete mitochondrial genome of Eimeria leuckarti (Eimeriidae, Coccidia, Apicomplexa) was obtained. This morphologically distinctive coccidium is considered to be the only valid Eimeria species of equids and it infects a range of both domestic and wild horses and their relatives. Despite the distinctive appearance of the oocysts of E. leuckarti, the mitochondrial genome organization and gene contents were comparable to other Eimeria spp. and related eimeriid coccidia infecting a range of mammals and birds. The greatly reduced 6242 bp genome is circular-mapping and contains three protein-coding genes (COI, COIII, CytB), 18 fragments encoding the large subunit rRNA (LSU), and 13 fragments encoding the small subunit (SSU) rRNA. No tRNA was encoded similar to other Apicomplexa. A Bayesian inference tree based on aligned CDS and rDNA fragments from Eimeria leuckarti and 34 other coccidia demonstrated that this mt genome has close phylogenetic affinities to Eimeria and Isospora species, and related eimeriid coccidia.

Exploiting digital droplet PCR and Next Generation Sequencing technologies to determine the relative abundance of individual Eimeria species in a DNA sample

DNA-based diagnostic assays for detecting infections with Eimeria species have been limited to providing identification and presence/absence data for samples containing oocysts. Modern technologies that generate quantitative data, such as droplet digital PCR (ddPCR) and Next Generation Sequencing (NGS), utilize a relatively short amplicon size containing sufficient species-specific variation for reliable species level identification. Targeting the cytochrome c oxidase subunit III gene in the mitochondrial genome, we established protocols using these technologies to determine the relative abundance of the number of copies/μL of Eimeria species in a sample. Samples from chickens of known and unknown Eimeria species composition were analyzed to determine the suitability of these technologies as diagnostic assays. All technologies demonstrated robust capability of identifying and quantifying the Eimeria species in samples. The new quantitative assays described herein will produce invaluable detail of Eimeria species infections for an array of situations in commercial chicken production systems, enabling further characterization of the disease profile and allowing for the development or enhancement of new intervention strategies.

A novel fragmented mitochondrial genome in the protist pathogen Toxoplasma gondii and related tissue coccidia

Mitochondrial genome content and structure vary widely across the eukaryotic tree of life, with protists displaying extreme examples. Apicomplexan and dinoflagellate protists have evolved highly reduced mitochondrial genome sequences, mtDNA, consisting of only three cytochrome genes and fragmented rRNA genes. Here, we report the independent evolution of fragmented cytochrome genes in Toxoplasma and related tissue coccidia and evolution of a novel genome architecture consisting minimally of 21 sequence blocks (SBs) totaling 5.9 kb that exist as nonrandom concatemers. Single-molecule Nanopore reads consisting entirely of SBs ranging from 0.1 to 23.6 kb reveal both whole and fragmented cytochrome genes. Full-length cytochrome transcripts including a divergent coxIII are detected. The topology of the mitochondrial genome remains an enigma. Analysis of a cob point mutation reveals that homoplasmy of SBs is maintained. Tissue coccidia are important pathogens of man and animals, and the mitochondrion represents an important therapeutic target. The mtDNA sequence has been elucidated, but a definitive genome architecture remains elusive.

Reevaluation of the Toxoplasma gondii and Neospora caninum genomes reveals misassembly, karyotype differences, and chromosomal rearrangements

Neosporacaninum primarily infects cattle, causing abortions, with an estimated impact of a billion dollars on the worldwide economy annually. However, the study of its biology has been unheeded by the established paradigm that it is virtually identical to its close relative, the widely studied human pathogen Toxoplasma gondii. By revisiting the genome sequence, assembly, and annotation using third-generation sequencing technologies, here we show that the N. caninum genome was originally incorrectly assembled under the presumption of synteny with T. gondii. We show that major chromosomal rearrangements have occurred between these species. Importantly, we show that chromosomes originally named Chr VIIb and VIII are indeed fused, reducing the karyotype of both N. caninum and T. gondii to 13 chromosomes. We reannotate the N. caninum genome, revealing more than 500 new genes. We sequence and annotate the nonphotosynthetic plastid and mitochondrial genomes and show that although apicoplast genomes are virtually identical, high levels of gene fragmentation and reshuffling exist between species and strains. Our results correct assembly artifacts that are currently widely distributed in the genome database of N. caninum and T. gondii and, more importantly, highlight the mitochondria as a previously oversighted source of variability and pave the way for a change in the paradigm of synteny, encouraging rethinking the genome as basis of the comparative unique biology of these pathogens.

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.

Mitochondrial Junction Region as Genotyping Marker for Cyclospora cayetanensis

Cyclosporiasis is an infection caused by Cyclospora cayetanensis, which is acquired by consumption of contaminated fresh food or water. In the United States, cases of cyclosporiasis are often associated with foodborne outbreaks linked to imported fresh produce or travel to disease-endemic countries. Epidemiologic investigation has been the primary method for linking outbreak cases. A molecular typing marker that can identify genetically related samples would be helpful in tracking outbreaks. We evaluated the mitochondrial junction region as a potential genotyping marker. We tested stool samples from 134 laboratory-confirmed cases in the United States by using PCR and Sanger sequencing. All but 2 samples were successfully typed and divided into 14 sequence types. Typing results were identical among samples within each epidemiologically defined case cluster for 7 of 10 clusters. These findings suggest that this marker can distinguish between distinct case clusters and might be helpful during cyclosporiasis outbreak investigations.

Insights into the evolution and drug susceptibility of Babesia duncani from the sequence of its mitochondrial and apicoplast genomes

Babesia microti and Babesia duncani are the main causative agents of human babesiosis in the United States. While significant knowledge about B. microti has been gained over the past few years, nothing is known about B. duncani biology, pathogenesis, mode of transmission or sensitivity to currently recommended therapies. Studies in immunocompetent wild type mice and hamsters have shown that unlike B. microti, infection with B. duncani results in severe pathology and ultimately death. The parasite factors involved in B. duncani virulence remain unknown. Here we report the first known completed sequence and annotation of the apicoplast and mitochondrial genomes of B. duncani. We found that the apicoplast genome of this parasite consists of a 34 kb monocistronic circular molecule encoding functions that are important for apicoplast gene transcription as well as translation and maturation of the organelle's proteins. The mitochondrial genome of B. duncani consists of a 5.9 kb monocistronic linear molecule with two inverted repeats of 48 bp at both ends. Using the conserved cytochrome b (Cytb) and cytochrome c oxidase subunit I (coxI) proteins encoded by the mitochondrial genome, phylogenetic analysis revealed that B. duncani defines a new lineage among apicomplexan parasites distinct from B. microti, Babesia bovis, Theileria spp. and Plasmodium spp. Annotation of the apicoplast and mitochondrial genomes of B. duncani identified targets for development of effective therapies. Our studies set the stage for evaluation of the efficacy of these drugs alone or in combination against B. duncani in culture as well as in animal models.

Host species and pathogenicity effects in the evolution of the mitochondrial genomes of Eimeria species (Apicomplexa; Coccidia; Eimeriidae)

Background

Mitochondria are fundamental organelles responsible for cellular metabolism and energy production in eukaryotes via the oxidative phosphorylation pathway. Mitochondrial DNA is often used in population and species studies with the assumption of neutral evolution. However, evidence of positive selection in mitochondrial coding genes of various animal species has accumulated suggesting that amino acid changes in mtDNA might be adaptive. The functional and physiological implications of the inferred positively selected sites are usually unknown and are only evaluated based on available structural and functional models. Such studies are absent in unicellular organisms that show several crucial differences to the electron transport chain of animal mitochondria. In the present study, we explored Eimeria mitogenomes for positive selection. We also tested for association between mtDNA polymorphism and environmental variation (i.e. host species), parasite life cycle (i.e. sporulation period), and efficient host cell invasion (i.e. pathogenicity, prepatent period).

Findings

We used site- and branch-site tests to estimate the extent of purifying and positive selection at each site and each lineage of several Eimeria parasite mitogenomes retrieved from GenBank. We founded sixteen codons in the three mtDNA-encoded proteins to be under positive selection compared to a strong purifying selection. Variation in the ratios of non-synonymous to synonymous changes of the studied parasites was associated with their different host species (F = 13.748; p < 0.001), whereas pathogenicity levels were associated with both synonymous and non-synonymous changes. This association was also confirmed by the multiple regression analysis.

Conclusions

Our results suggest that host species and pathogenicity are important factors that might shape mitochondrial variation in Eimeria parasites. This supports the important role of mtDNA variations in the evolution and adaptation of these parasites.

Mitochondrial genomes of Australian chicken Eimeria support the presence of ten species with low genetic diversity among strains

Modern molecular approaches have vastly improved diagnostic capabilities for differentiating among species of chicken infecting Eimeria. Consolidating information from multiple genetic markers, adding additional poultry Eimeria species and increasing the size of available data-sets is improving the resolving power of the DNA, and consequently our understanding of the genus. This study adds information from 25 complete mitochondrial DNA genomes from Australian chicken Eimeria isolates representing all 10 species known to occur in Australia, including OTU-X, -Y and -Z. The resulting phylogeny provides a comprehensive view of species relatedness highlighting where the OTUs align with respect to others members of the genus. All three OTUs fall within the Eimeria clade that contains only chicken-infecting species with close affinities to E. maxima, E. brunetti and E. mitis. Mitochondrial genetic diversity was low among Australian isolates likely reflecting their recent introduction to the country post-European settlement. The lack of observed genetic diversity is a promising outcome as it suggests that the currently used live vaccines should continue to offer widespread protection against Eimeria outbreaks in all states and territories. Flocks were frequently found to host multiple strains of the same species, a factor that should be considered when studying disease epidemiology in the field.

Are Eimeria Genetically Diverse, and Does It Matter?

Eimeria pose a risk to all livestock species as a cause of coccidiosis, reducing productivity and compromising animal welfare. Pressure to reduce drug use in the food chain makes the development of cost-effective vaccines against Eimeria essential. For novel vaccines to be successful, understanding genetic and antigenic diversity in field populations is key. Eimeria species that infect chickens are most significant, with Eimeria tenella among the best studied and most economically important. Genome-wide single nucleotide polymorphism (SNP)-based haplotyping has been used to determine population structure, genotype distribution, and potential for cross-fertilization between E. tenella strains. Here, we discuss recent developments in our understanding of diversity for Eimeria in relation to its specialized life cycle, distribution across the globe, and the challenges posed to vaccine development.

The complete mitochondrial genomes of five Eimeria species infecting domestic rabbits

Rabbit coccidiosis caused by members of the genus Eimeria can cause enormous economic impact worldwide, but the genetics, epidemiology and biology of these parasites remain poorly understood. In the present study, we sequenced and annotated the complete mitochondrial (mt) genomes of five Eimeria species that commonly infect the domestic rabbits. The complete mt genomes of Eimeria intestinalis, Eimeria flavescens, Eimeria media, Eimeria vejdovskyi and Eimeria irresidua were 6261bp, 6258bp, 6168bp, 6254bp, 6259bp in length, respectively. All of the mt genomes consist of 3 genes for proteins (cytb, cox1, and cox3), 14 gene fragments for the large subunit (LSU) rRNA and 11 gene fragments for the small subunit (SSU) rRNA, but no transfer RNA (tRNA) genes. The gene order of the mt genomes is similar to that of Plasmodium, but distinct from Haemosporida and Theileria. Phylogenetic analyses based on full nucleotide sequences using Bayesian analysis revealed that the monophyly of the Eimeria of rabbits was strongly statistically supported with a Bayesian posterior probabilities. These data provide novel mtDNA markers for studying the population genetics and molecular epidemiology of the Eimeria species, and should have implications for the molecular diagnosis, prevention and control of coccidiosis in rabbits.

Genetic similarities between Cyclospora cayetanensis and cecum-infecting avian Eimeria spp. in apicoplast and mitochondrial genomes

Background

Cyclospora cayetanensis is an important cause for diarrhea in children in developing countries and foodborne outbreaks of cyclosporiasis in industrialized nations. To improve understanding of the basic biology of Cyclospora spp. and development of molecular diagnostic tools and therapeutics, we sequenced the complete apicoplast and mitochondrial genomes of C. cayetanensis.

Methods

The genome of one Chinese C. cayetanensis isolate was sequenced using Roche 454 and Illumina technologies. The assembled genomes of the apicoplast and mitochondrion were retrieved, annotated, and compared with reference genomes for other apicomplexans to infer genome organizations and phylogenetic relationships. Sequence variations in the mitochondrial genome were identified by comparison of two C. cayetanensis nucleotide sequences from this study and a recent publication.

Results

The apicoplast and mitochondrial genomes of C. cayetanensis are 34,155 and 6,229 bp in size and code for 65 and 5 genes, respectively. Comparative genomic analysis showed high similarities between C. cayetanensis and Eimeria tenella in both genomes; they have 85.6 % and 90.4 % nucleotide sequence similarities, respectively, and complete synteny in gene organization. Phylogenetic analysis of the genomic sequences confirmed the genetic similarities between cecum-infecting avian Eimeria spp. and C. cayetanensis. Like in other coccidia, both genomes of C. cayetanensis are transcribed bi-directionally. The apicoplast genome is circular, codes for the complete machinery for protein biosynthesis, and contains two inverted repeats that differ slightly in LSU rRNA gene sequences. In contrast, the mitochondrial genome has a linear concatemer or circular mapping topology. Eight single-nucleotide and one 7-bp multiple-nucleotide variants were detected between the mitochondrial genomes of C. cayetanensis from this and recent studies.

Conclusions

The apicoplast and mitochondrial genomes of C. cayetanensis are highly similar to those of cecum-infecting avian Eimeria spp. in both genome organization and sequences. The availability of sequence data beyond rRNA and heat shock protein genes could facilitate studies of C. cayetanensis biology and development of genotyping tools for investigations of cyclosporiasis outbreaks.

The complete mitochondrial genome sequence of Eimeria innocua (Eimeriidae, Coccidia, Apicomplexa)

The complete mitochondrial genome of Eimeria innocua KR strain (Eimeriidae, Coccidia, Apicomplexa) was sequenced. This coccidium infects turkeys (Meleagris gallopavo), Bobwhite quails (Colinus virginianus), and Grey partridges (Perdix perdix). Genome organization and gene contents were comparable with other Eimeria spp. infecting galliform birds. The circular-mapping mt genome of E. innocua is 6247 bp in length with three protein-coding genes (cox1, cox3, and cytb), 19 gene fragments encoding large subunit (LSU) rRNA and 14 gene fragments encoding small subunit (SSU) rRNA. Like other Apicomplexa, no tRNA was encoded. The mitochondrial genome of E. innocua confirms its close phylogenetic affinities to Eimeria dispersa.

The Complete Mitochondrial Genome of the Foodborne Parasitic Pathogen Cyclospora cayetanensis

Cyclospora cayetanensis is a human-specific coccidian parasite responsible for several food and water-related outbreaks around the world, including the most recent ones involving over 900 persons in 2013 and 2014 outbreaks in the USA. Multicopy organellar DNA such as mitochondrion genomes have been particularly informative for detection and genetic traceback analysis in other parasites. We sequenced the C. cayetanensis genomic DNA obtained from stool samples from patients infected with Cyclospora in Nepal using the Illumina MiSeq platform. By bioinformatically filtering out the metagenomic reads of non-coccidian origin sequences and concentrating the reads by targeted alignment, we were able to obtain contigs containing Eimeria-like mitochondrial, apicoplastic and some chromosomal genomic fragments. A mitochondrial genomic sequence was assembled and confirmed by cloning and sequencing targeted PCR products amplified from Cyclospora DNA using primers based on our draft assembly sequence. The results show that the C. cayetanensis mitochondrion genome is 6274 bp in length, with 33% GC content, and likely exists in concatemeric arrays as in Eimeria mitochondrial genomes. Phylogenetic analysis of the C. cayetanensis mitochondrial genome places this organism in a tight cluster with Eimeria species. The mitochondrial genome of C. cayetanensis contains three protein coding genes, cytochrome (cytb), cytochrome C oxidase subunit 1 (cox1), and cytochrome C oxidase subunit 3 (cox3), in addition to 14 large subunit (LSU) and nine small subunit (SSU) fragmented rRNA genes.

A linear mitochondrial genome of Cyclospora cayetanensis (Eimeriidae, Eucoccidiorida, Coccidiasina, Apicomplexa) suggests the ancestral start position within mitochondrial genomes of eimeriid coccidia

The near complete mitochondrial genome for Cyclospora cayetanensis is 6184 bp in length with three protein-coding genes (Cox1, Cox3, CytB) and numerous lsrDNA and ssrDNA fragments. Gene arrangements were conserved with other coccidia in the Eimeriidae, but the C. cayetanensis mitochondrial genome is not circular-mapping. Terminal transferase tailing and nested PCR completed the 5'-terminus of the genome starting with a 21 bp A/T-only region that forms a potential stem-loop. Regions homologous to the C. cayetanensis mitochondrial genome 5'-terminus are found in all eimeriid mitochondrial genomes available and suggest this may be the ancestral start of eimeriid mitochondrial genomes.

Simultaneous identification and DNA barcoding of six Eimeria species infecting turkeys using PCR primers targeting the mitochondrial cytochrome c oxidase subunit I (mtCOI) locus

Species-specific PCR primers targeting the mitochondrial cytochrome c oxidase subunit I (mtCOI) locus were generated that allow for the specific identification of the most common Eimeria species infecting turkeys (i.e., Eimeria adenoeides, Eimeria meleagrimitis, Eimeria gallopavonis, Eimeria meleagridis, Eimeria dispersa, and Eimeria innocua). PCR reaction chemistries were optimized with respect to divalent cation (MgCl2) and dNTP concentrations, as well as PCR cycling conditions (particularly anneal temperature for primers). Genomic DNA samples from single oocyst-derived lines of six Eimeria species were tested to establish specificity and sensitivity of these newly designed primer pairs. A mixed 60-ng total DNA sample containing 10 ng of each of the six Eimeria species was used as DNA template to demonstrate specific amplification of the correct product using each of the species-specific primer pairs. Ten nanograms of each of the five non-target Eimeria species was pooled to provide a non-target, control DNA sample suitable to test the specificity of each primer pair. The amplifications of the COI region with species-specific primer pairs from pooled samples yielded products of expected sizes (209 to 1,012 bp) and no amplification of non-target Eimeria sp. DNA was detected using the non-target, control DNA samples. These primer pairs specific for Eimeria spp. of turkeys did not amplify any of the seven Eimeria species infecting chickens. The newly developed PCR primers can be used as a diagnostic tool capable of specifically identifying six turkey Eimeria species; additionally, sequencing of the PCR amplification products yields sequence-based genotyping data suitable for identification and molecular phylogenetics.

Effects of artesunate treatment on Plasmodium gallinaceum transmission in the vectors Aedes aegypti and Culex quinquefasciatus

In the absence of vaccines, chemotherapy is an effective and economical way for controlling malaria. Development of anti-malarial drugs that target pathogenic blood stage parasites and gametocytes is preferable for the treatment as it can alleviate the host's morbidity and mortality and block transmission of the Plasmodium parasite. Recently, our laboratory has developed an in vivo transmission blocking assay that involves administration of 7 consecutive daily doses of a test compound into domestic chickens (Gallus gallus domesticus) infected with the avian malaria parasite Plasmodium gallinaceum with 10% parasitaemia and 1% gametocytaemia. To compromise the cost and time for artesunate (ATN) treatment, this study aimed to investigate effects of a 5-day consecutive administration of 10 milligrams per kilogram (mg/kg) ATN on P. gallinaceum infection in chickens and transmission to two natural vectors, Aedes aegypti and Culex quinquefasciatus. Our study showed that the treatment with 10 mg/kg ATN for 7 days, but not 5 days, completely eliminated blood stage infections, prevented recrudescence and blocked gametocyte production and transmission of P. gallinaceum to its vectors, thereby confirming the potent schizontocidal and gametocytocidal activities of ATN. This regimen should be further evaluated in field trials.

Complete mitochondrial genome sequences from five Eimeria species (Apicomplexa; Coccidia; Eimeriidae) infecting domestic turkeys

BACKGROUND

Clinical and subclinical coccidiosis is cosmopolitan and inflicts significant losses to the poultry industry globally. Seven named Eimeria species are responsible for coccidiosis in turkeys: Eimeria dispersa; Eimeria meleagrimitis; Eimeria gallopavonis; Eimeria meleagridis; Eimeria adenoeides; Eimeria innocua; and, Eimeria subrotunda. Although attempts have been made to characterize these parasites molecularly at the nuclear 18S rDNA and ITS loci, the maternally-derived and mitotically replicating mitochondrial genome may be more suited for species level molecular work; however, only limited sequence data are available for Eimeria spp. infecting turkeys. The purpose of this study was to sequence and annotate the complete mitochondrial genomes from 5 Eimeria species that commonly infect the domestic turkey (Meleagris gallopavo).

METHODS

Six single-oocyst derived cultures of five Eimeria species infecting turkeys were PCR-amplified and sequenced completely prior to detailed annotation. Resulting sequences were aligned and used in phylogenetic analyses (BI, ML, and MP) that included complete mitochondrial genomes from 16 Eimeria species or concatenated CDS sequences from each genome.

RESULTS

Complete mitochondrial genome sequences were obtained for Eimeria adenoeides Guelph, 6211 bp; Eimeria dispersa Briston, 6238 bp; Eimeria meleagridis USAR97-01, 6212 bp; Eimeria meleagrimitis USMN08-01, 6165 bp; Eimeria gallopavonis Weybridge, 6215 bp; and Eimeria gallopavonis USKS06-01, 6215 bp). The order, orientation and CDS lengths of the three protein coding genes (COI, COIII and CytB) as well as rDNA fragments encoding ribosomal large and small subunit rRNA were conserved among all sequences. Pairwise sequence identities between species ranged from 88.1% to 98.2%; sequence variability was concentrated within CDS or between rDNA fragments (where indels were common). No phylogenetic reconstruction supported monophyly of Eimeria species infecting turkeys; Eimeria dispersa may have arisen via host switching from another avian host. Phylogenetic analyses suggest E. necatrix and E. tenella are related distantly to other Eimeria of chickens.

CONCLUSIONS

Mitochondrial genomes of Eimeria species sequenced to date are highly conserved with regard to gene content and structure. Nonetheless, complete mitochondrial genome sequences and, particularly the three CDS, possess sufficient sequence variability for differentiating Eimeria species of poultry. The mitochondrial genome sequences are highly suited for molecular diagnostics and phylogenetics of coccidia and, potentially, genetic markers for molecular epidemiology.

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.

Gene conversion shapes linear mitochondrial genome architecture

Recently, it was shown that gene conversion between the ends of linear mitochondrial chromosomes can cause telomere expansion and the duplication of subtelomeric loci. However, it is not yet known how widespread this phenomenon is and how significantly it has impacted organelle genome architecture. Using linear mitochondrial DNAs and mitochondrial plasmids from diverse eukaryotes, we argue that telomeric recombination has played a major role in fashioning linear organelle chromosomes. We find that mitochondrial telomeres frequently expand into subtelomeric regions, resulting in gene duplications, homogenizations, and/or fragmentations. We suggest that these features are a product of subtelomeric gene conversion, provide a hypothetical model for this process, and employ genetic diversity data to support the idea that the greater the effective population size the greater the potential for gene conversion between subtelomeric loci.

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.

In vivo transmission blocking activities of artesunate on the avian malaria parasite Plasmodium gallinaceum

Infection and transmission of the avian malaria parasite Plasmodium gallinaceum in domestic chickens is associated with high economic burden and presents a major challenge to poultry industry in South East Asia. Development of drugs targeting both asexual blood stage parasites and sexual stages of the avian malarias will be beneficial for malaria treatment and eradication. However, current drugs recommended for treatment of the avian malaria parasites target specifically the asexual blood stage parasites, but have little or no impact to the gametocytes, the major target for development of transmission-blocking strategies. In the present work, we established a simple procedure to evaluate gametocytocidal and transmission blocking activities in a P. gallinaceum-avian model. The assays involved administration of seven consecutive daily doses of test compounds into P. gallinaceum-infected chickens with 10% parasitaemia and 1% gametocytaemia. Our studies indicated that intramuscular injection with seven daily low doses (the minimum effective dose of 10mg/kg) of artesunate blocked the gametocyte production and transmission to the mosquito vector Aedes aegypti. This assay can be further applicable for testing new compounds against P. gallinaceum and for other parasitic protozoa infecting birds.

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.

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.

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

Mitochondrial (mt) genomes from diverse phylogenetic groups vary considerably in size, structure and organization. The genus Plasmodium, the causative agent of malaria, has the smallest mt genome in the form of a tandemly repeated, linear element of 6 kb. The Plasmodium mt genome encodes only three protein genes (cox1, cox3 and cob) and large- and small-subunit ribosomal RNA (rRNA) genes, which are highly fragmented with 19 identified rRNA pieces. The complete mt genome sequences of 21 Plasmodium species have been published but a thorough investigation of the arrangement of rRNA gene fragments has been undertaken for only Plasmodium falciparum, the human malaria parasite. In this study, we determined the arrangement of mt rRNA gene fragments in 23 Plasmodium species, including two newly determined mt genome sequences from P. gallinaceum and P. vinckei vinckei, as well as Leucocytozoon caulleryi, an outgroup of Plasmodium. Comparative analysis reveals complete conservation of the arrangement of rRNA gene fragments in the mt genomes of all the 23 Plasmodium species and L. caulleryi. Surveys for a new rRNA gene fragment using hidden Markov models enriched with recent mt genome sequences led us to suggest the mtR-26 sequence as a novel candidate LSU rRNA fragment in the mt genomes of the 24 species. Additionally, we found 22-25 bp-inverted repeat sequences, which may be involved in the generation of lineage-specific mt genome arrangements after divergence from a common ancestor of the genera Eimeria and Plasmodium/Leucocytozoon.