The complete mitochondrial genome sequence ofEimeria mitis(Apicomplexa: Coccidia)

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.

Comprehensive molecular characterization of complete mitogenome assemblies of 33 Eimeria isolates infecting domestic chickens

BACKGROUND

Coccidiosis caused by Eimeria is one of the most severe chicken diseases and poses a great economic threat to the poultry industry. Understanding the evolutionary biology of chicken Eimeria parasites underpins development of new interactions toward the improved prevention and control of this poultry disease.

METHODS

We presented an evolutionary blueprint of chicken coccidia by genetically characterizing complete mitogenome assemblies of 33 isolates representing all seven known Eimeria species infecting chickens in China. Further genome- and gene-level phylogenies were also achieved to better understand the evolutionary relationships of these chicken Eimeria at the species level.

RESULTS

33 mitogenomes of chicken eimerian parasites ranged from 6148 bp to 6480 bp in size and encoded typical mitochondrial compositions of apicomplexan parasites including three protein-coding genes (PCGs), seven fragmented small subunit (SSU) and 12/13 fragmented large subunit (LSU) rRNAs. Comparative genomics provided an evolutionary scenario for the genetic diversity of PCGs-cytochrome c oxidase subunits 1 and 3 (cox1 and cox3) and cytochrome b (cytb); all were under purifying selection with cox1 and cox3 being the lowest and highest evolutionary rates, respectively. Genome-wide phylogenies classified the 33 Eimeria isolates into seven subgroups, and furthermore Eimeria tenella and Eimeria necatrix were determined to be more closely related to each other than to the other eight congenic species. Single/concatenated mitochondrial protein gene-based phylogenies supported cox1 as the genetic marker for evolutionary and phylogenetic studies for avain coccidia.

CONCLUSIONS

To our knowledge, these are the first population-level mitogenomic data on the genus Eimeria, and its comprehensive molecular characterization provides valuable resources for systematic, population genetic and evolutionary biological studies of apicomplexan parasites in poultry.

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.

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.

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.

The full mitochondrial genome sequence of Raillietina tetragona from chicken (Cestoda: Davaineidae)

In the present study, the complete mitochondrial DNA (mtDNA) sequence of Raillietina tetragona was sequenced and its gene contents and genome organizations was compared with that of other tapeworm. The complete mt genome sequence of R. tetragona is 14,444 bp in length. It contains 12 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and two non-coding region. All genes are transcribed in the same direction and have a nucleotide composition high in A and T. The contents of A + T of the complete mt genome are 71.4% for R. tetragona. The R. tetragona mt genome sequence provides novel mtDNA marker for studying the molecular epidemiology and population genetics of Raillietina and has implications for the molecular diagnosis of chicken cestodosis caused by Raillietina.

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.

A selective review of advances in coccidiosis research

Coccidiosis is a widespread and economically significant disease of livestock caused by protozoan parasites of the genus Eimeria. This disease is worldwide in occurrence and costs the animal agricultural industry many millions of dollars to control. In recent years, the modern tools of molecular biology, biochemistry, cell biology and immunology have been used to expand greatly our knowledge of these parasites and the disease they cause. Such studies are essential if we are to develop new means for the control of coccidiosis. In this chapter, selective aspects of the biology of these organisms, with emphasis on recent research in poultry, are reviewed. Topics considered include taxonomy, systematics, genetics, genomics, transcriptomics, proteomics, transfection, oocyst biogenesis, host cell invasion, immunobiology, diagnostics and control.

Sequence variability in four mitochondrial genes among pinworm Aspicularis tetraptera isolates from laboratory mice in four provinces, China

This research aimed at exploring sequence variability in four mitochondrial (mt) genes, namely, cytochrome c oxidase subunit 1 (cox1), cytochrome b (cytb) and NADH dehydrogenase subunits 1 and 5 (nad1 and nad5), among pinworm Aspicularis tetraptera isolates from laboratory mice in four different provinces, China. A part of the cox1 (pcox1), cytb (pcytb), nad1 and nad5 genes (pnad1 and pnad5) were amplified separately from individual pinworms by polymerase chain reaction (PCR) and sequenced to determine sequence variations and examine their phylogenetic relationships. Herein, the intra-specific sequence variations within A. tetraptera were 0-0.5% for pcox1, 0-1.4% for pcytb, 0-1.8% for pnad1 and 0-1.7% for pnad5, respectively. In contrast, the inter-specific sequence differences among members of the Oxyuridae were significantly higher, being 13.7-17.0% for pcox1, 24.5-34.7% for pcytb, 26.6-29.6% for pnad1 and 24.4-25.5% for pnad5, respectively. Three methods, namely, Bayesian inference (BI), maximum likelihood (ML) and maximum parsimony (MP), were used for phylogenetic analyses based on the combined sequences of the four mt gene sequences, and the results indicated that all A. tetraptera samples form monophyletic groups, but samples from the same geographical origin did not always cluster together. This study demonstrated the existence of low-level intra-specific variation in four mtDNA sequences among A. tetraptera isolates from laboratory mice in different geographic regions in China, indicating no obvious geographical distinction among A. tetraptera isolates in China. These findings have important implications for studying systematics, molecular epidemiology and population genetics of A. tetraptera.

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.