Long PCR amplification of the entire mitochondrial genome from individual helminths for direct sequencing

Exploring the evidence for mitochondrial dysfunction and genetic abnormalities in the etiopathogenesis of tropical ataxic neuropathy

Tropical ataxic neuropathy (TAN) is characterised by ataxic polyneuropathy, degeneration of the posterior columns and pyramidal tracts, optic atrophy, and sensorineural hearing loss. It has been attributed to nutritional/toxic etiologies, but evidence for the same has been equivocal. TAN shares common clinical features with inherited neuropathies and mitochondrial disorders, it may be hypothesised that genetic abnormalities may underlie the pathophysiology of TAN. This study aimed to establish evidence for mitochondrial dysfunction by adopting an integrated biochemical and multipronged genetic analysis. Patients (n = 65) with chronic progressive ataxic neuropathy with involvement of visual and/or auditory pathways underwent deep phenotyping, genetic studies including mitochondrial DNA (mtDNA) deletion analysis, mtDNA and clinical exome sequencing (CES), and respiratory chain complex (RCC) assay. The phenotypic characteristics included dysfunction of visual (n = 14), auditory (n = 12) and visual + auditory pathways (n = 29). Reduced RCC activity was present in 13 patients. Mitochondrial DNA deletions were noted in five patients. Sequencing of mtDNA (n = 45) identified a homoplasmic variant (MT-ND6) and a heteroplasmic variant (MT-COI) in one patient each. CES (n = 45) revealed 55 variants in nuclear genes that are associated with neuropathy (n = 27), deafness (n = 7), ataxia (n = 4), and mitochondrial phenotypes (n = 5) in 36 patients. This study provides preliminary evidence that TAN is associated with a spectrum of genetic abnormalities, including those associated with mitochondrial dysfunction, which is in contradistinction from the prevailing hypothesis that TAN is related to dietary toxins. Analysing the functional relevance of these genetic variants may improve the understanding of the pathogenesis of TAN.

Comparative Mitogenomics of Jumping Spiders with First Complete Mitochondrial Genomes of Euophryini (Araneae: Salticidae)

Salticidae is the most species-rich family of spiders with diverse morphology, ecology and behavior. However, the characteristics of the mitogenomes within this group are poorly understood with relatively few well-characterized complete mitochondrial genomes. In this study, we provide completely annotated mitogenomes for Corythalia opima and Parabathippus shelfordi, which represent the first complete mitogenomes of the tribe Euophryini of Salticidae. The features and characteristics of the mitochondrial genomes are elucidated for Salticidae by thoroughly comparing the known well-characterized mitogenomes. The gene rearrangement between trnL2 and trnN was found in two jumping spider species, Corythalia opima and Heliophanus lineiventris Simon, 1868. Additionally, the rearrangement of nad1 to between trnE and trnF found in Asemonea sichuanensis Song & Chai, 1992 is the first protein-coding gene rearrangement in Salticidae, which may have an important phylogenetic implication for the family. Tandem repeats of various copy numbers and lengths were discovered in three jumping spider species. The codon usage analyses showed that the evolution of codon usage bias in salticid mitogenomes was affected by both selection and mutational pressure, but selection may have played a more important role. The phylogenetic analyses provided insight into the taxonomy of Colopsus longipalpis (Żabka, 1985). The data presented in this study will improve our understanding of the evolution of mitochondrial genomes within Salticidae.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae

BACKGROUND

Arcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2-3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear.

RESULTS

We assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes.

CONCLUSIONS

The presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

Comprehensive Molecular Characterization of the Mitochondrial Genome of the Takin Lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae)

The takin lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae) causes lethal bronchopneumonia and represents severe threats to captive and wild populations. However, until now there has been very limited information available concerning the molecular epidemiology and evolutionary biology of V. eleguneniensis. Mitochondrial genomes (mtDNAs) can provide resources for investigations in these areas and, therefore, can assist with the surveillance and control of this lungworm. Herein, the complete mtDNA of V. eleguneniensis was sequenced and characterized with Illumina pipeline analyses. This circular genome (13,625 bp) encoded twelve protein-coding genes (PCGs), two rRNAs, and twenty-two tRNAs, with notable levels of AT and GC skews. Comparative genomics revealed a purifying selection among PCGs, with cox1 and nad6 having the lowest and the highest evolutionary rate, respectively. Genome-wide phylogenies showed a close relationship between V. eleguneniensis and Protostrongylus rufescens in Strongylida. Single gene (PCGs or rRNAs)-based phylogenies indicated that cox1 and nad5 genes shared the same family-level topology with that inferred from genomic datasets, suggesting that both genes could be suitable genetic markers for evolutionary and phylogenetic studies of Strongylida species. This was the first mtDNA of any member of the genus Varestrongylus, and its comprehensive molecular characterization represents a new resource for systematic, population genetic and evolutionary biological studies of Varestrongylus lungworms in wildlife.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

Sequencing and Reconstructing Helminth Mitochondrial Genomes Directly from Genomic Next-Generation Sequencing Data

We present a detailed method for extraction of high-molecular weight genomic DNA suitable for numerous DNA sequencing applications, and a straightforward in silico approach for reconstructing novel mitochondrial (mt) genomes directly from total genomic DNA extracts derived from next-generation sequencing (NGS) data sets. The in silico post-sequencing pipeline described is fast, accurate, and highly efficient, with modest memory requirements that can be performed using a standard desktop computer. The approach is particularly effective for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information currently available and overcomes many of the limitations of traditional strategies. The described methodologies are also applicable for metagenomics sequencing from mixed or pooled samples containing multiple species and subsequent specific assembly of specific mitochondrial genomes.

Mitochondrial Genome of Nonmodel Marine Metazoans by Next-Generation Sequencing (NGS)

Mitochondrial genomes (mtgenome) represent an important source of information for addressing fundamental evolutionary, phylogeographic, systematic, and ecological questions in marine organisms. In the last two decades the advent of high-throughput next-generation sequencing (NGS) has provided an unprecedented possibility to access large amount of genomic data and, as such, there has been a rapid growth in mtgenome resources and studies. In particular, NGS strategies represent a great advantage for investigating nonmodel marine organisms for which no or limited genomic resources are available. Here, we describe a routinely used standardized protocol to obtain mtgenome of nonmodel marine organisms by NGS. The protocol is composed of five main steps, including DNA extraction, DNA fragmentation, library preparation, high-throughput sequencing, and bioinformatic analyses. Each of the first three steps is followed by size/quality and concentration validations. The advantages of the described protocol rely on the assumption that no a priori information on mtgenome of the studied organism is needed and on its versatility as researchers may choose several kits for DNA extraction and library preparation and adopt different methods for DNA fragmentation depending on their needs, experience, and suppliers.

Diagnosis of primary mitochondrial disorders -Emphasis on myopathological aspects

Mitochondrial disorders are one of the most common neurometabolic disorders affecting all age groups. The phenotype-genotype heterogeneity in these disorders can be attributed to the dual genetic control on mitochondrial functions, posing a challenge for diagnosis. Though the advancement in the high-throughput sequencing and other omics platforms resulted in a "genetics-first" approach, the muscle biopsy remains the benchmark in most of the mitochondrial disorders. This review focuses on the myopathological aspects of primary mitochondrial disorders. The utility of muscle biopsy is not limited to analyse the structural abnormalities; rather it also proves to be a potential tool to understand the deranged sub-cellular functions.

A Novel L1 Linker Mutation in DES Resulted in Total Absence of Protein

Desminopathies (MIM*601419) are clinically heterogeneous, manifesting with myopathy and/or cardiomyopathy and with intra-sarcoplasmic desmin-positive deposits. They have either an autosomal dominant (AD) or recessive (AR) pattern of inheritance. Desmin is a crucial intermediate filament protein regulating various cellular functions in muscle cells. Here, we report a 13-year-old girl, born of second-degree consanguineous parents, with normal developmental milestones, who presented with dilated cardiomyopathy, respiratory insufficiency and predominant distal upper limb weakness. A striking feature on muscle biopsy was the presence of a peripheral chain of nuclei in addition to myopathic features. Immunostaining showed complete lack of desmin expression, further confirmed by western blot analysis. Ultrastructurally, subsarcolemmal granular material, expanded Z-band aggregation, distortion of myofilaments, focal Z-band streaming, lobed and clustered myonuclei were observed. Next-generation sequencing revealed a novel homozygous nonsense mutation c.448C>T, p.R150X in the patient, while the parents were heterozygous carriers. Single mitochondrial DNA deletion and isolated complex IV deficiency were noted. Our findings add to the ever-expanding phenotype and molecular spectrum of desminopathies.

Nanopore Sequencing Resolves Elusive Long Tandem-Repeat Regions in Mitochondrial Genomes

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

Clinico-pathological and Molecular Spectrum of Mitochondrial Polymerase γ Mutations in a Cohort from India

Polymerase γ catalytic subunit (POLG), a nuclear gene, encodes the enzyme responsible for mitochondrial DNA (mtDNA) replication. POLG mutations are a major cause of inherited mitochondrial diseases. They present with varied phenotypes, age of onset, and severity. Reports on POLG mutations from India are limited. Hence, this study aimed to describe the clinico-pathological and molecular observations of POLG mutations. A total of 446 patients with clinical diagnosis of mitochondrial disorders were sequenced for all exons and intron-exon boundaries of POLG. Of these, 19 (4.26%) patients (M:F: 10:9) had POLG mutations. The age of onset ranged from 5 to 55 years with an overlapping phenotypic spectrum. Ptosis, peripheral neuropathy, seizures, and ataxia were the common neurological features observed. The most common clinical phenotype was chronic progressive external ophthalmoplegia (CPEO) and CPEO plus (n = 14). Muscle biopsy showed characteristic features of mitochondrial myopathy in fourteen patients (14/19) and respiratory chain enzyme deficiency in eleven patients (11/19). Multiple mtDNA deletions were seen in 47.36% (9/19) patients. Eight pathogenic POLG variations including two novel variations (p.G132R and p.V1106A) were identified. The common pathogenic mutation identified was p.L304R, being present in eight patients (42.1%) predominantly in the younger age group followed by p.W748S in four patients (21%). To the best of our knowledge, this is the first extensive study from India, highlights the clinico-pathological and molecular spectrum of POLG mutations.

Characterization of the Complete Mitochondrial Genomes of Two Sibling Species of Parasitic Roundworms, Haemonchus contortus and Teladorsagia circumcincta

Haemonchus contortus and Teladorsagia circumcincta are among the two most pathogenic internal parasitic nematodes infecting small ruminants, such as sheep and goats, and are a global animal health issue. Accurate identification and delineation of Haemonchidae species is essential for development of diagnostic and control strategies with high resolution for Trichostrongyloidea infection in ruminants. Here, we describe in detail and compare the complete mitochondrial (mt) genomes of the New Zealand H. contortus and T. circumcincta field strains to improve our understanding of species- and strain-level evolution in these closely related roundworms. In the present study, we performed extensive comparative bioinformatics analyses on the recently sequenced complete mt genomes of the New Zealand H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP field strains. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP). The AT-rich mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP are 14,001 bp (A+T content of 77.4%) and 14,081 bp (A+T content of 77.3%) in size, respectively. All 36 of the typical nematode mt genes are transcribed in the forward direction in both species and comprise of 12 protein-encoding genes (PCGs), 2 ribosomal RNA (rrn) genes, and 22 transfer RNA (trn) genes. The secondary structures for the 22 trn genes and two rrn genes differ between H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP, however the gene arrangements of both are consistent with other Trichostrongylidea sequenced to date. Comparative analyses of the complete mitochondrial nucleotide sequences, PCGs, A+T rich and non-coding repeat regions of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP further reinforces the high levels of diversity and gene flow observed among Trichostrongylidea, and supports their potential as ideal markers for strain-level identification from different hosts and geographical regions with high resolution for future studies. The complete mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP presented here provide useful novel markers for further studies of the meta-population connectivity and the genetic mechanisms driving evolution in nematode species.

Vortex fluidics-mediated DNA rescue from formalin-fixed museum specimens

DNA from formalin-preserved tissue could unlock a vast repository of genetic information stored in museums worldwide. However, formaldehyde crosslinks proteins and DNA, and prevents ready amplification and DNA sequencing. Formaldehyde acylation also fragments the DNA. Treatment with proteinase K proteolyzes crosslinked proteins to rescue the DNA, though the process is quite slow. To reduce processing time and improve rescue efficiency, we applied the mechanical energy of a vortex fluidic device (VFD) to drive the catalytic activity of proteinase K and recover DNA from American lobster tissue (Homarus americanus) fixed in 3.7% formalin for >1-year. A scan of VFD rotational speeds identified the optimal rotational speed for recovery of PCR-amplifiable DNA and while 500+ base pairs were sequenced, shorter read lengths were more consistently obtained. This VFD-based method also effectively recovered DNA from formalin-preserved samples. The results provide a roadmap for exploring DNA from millions of historical and even extinct species.

MitoZ: a toolkit for animal mitochondrial genome assembly, annotation and visualization

Mitochondrial genome (mitogenome) plays important roles in evolutionary and ecological studies. It becomes routine to utilize multiple genes on mitogenome or the entire mitogenomes to investigate phylogeny and biodiversity of focal groups with the onset of High Throughput Sequencing (HTS) technologies. We developed a mitogenome toolkit MitoZ, consisting of independent modules of de novo assembly, findMitoScaf (find Mitochondrial Scaffolds), annotation and visualization, that can generate mitogenome assembly together with annotation and visualization results from HTS raw reads. We evaluated its performance using a total of 50 samples of which mitogenomes are publicly available. The results showed that MitoZ can recover more full-length mitogenomes with higher accuracy compared to the other available mitogenome assemblers. Overall, MitoZ provides a one-click solution to construct the annotated mitogenome from HTS raw data and will facilitate large scale ecological and evolutionary studies. MitoZ is free open source software distributed under GPLv3 license and available at https://github.com/linzhi2013/MitoZ.

Mitochondrial Genomes of Two Thaparocleidus Species (Platyhelminthes: Monogenea) Reveal the First rRNA Gene Rearrangement among the Neodermata

Phylogenetic framework for the closely related Ancylodiscoidinae and Ancyrocephalinae subfamilies remains contentious. As this issue was never studied using a large molecular marker, we sequenced the first two Ancylodiscoidinae mitogenomes: Thaparocleidus asoti and Thaparocleidus varicus. Both mitogenomes had two non-coding regions (NCRs) that contained a number of repetitive hairpin-forming elements (RHE). Due to these, the mitogenome of T. asoti (16,074 bp) is the longest among the Monogenea; especially large is its major NCR, with 3500 bp, approximately 1500 bp of which could not be sequenced (thus, the total mitogenome size is ≈ 17,600 bp). Although RHEs have been identified in other monopisthocotyleans, they appear to be independently derived in different taxa. The presence of RHEs may have contributed to the high gene order rearrangement rate observed in the two mitogenomes, including the first report of a transposition of rRNA genes within the Neodermata. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrinae embedded within the Ancyrocephalinae (paraphyly), whereas Ancylodiscoidinae formed a sister-group with them. This was also supported by the gene order analysis. 28S rDNA dataset produced polyphyletic Dactylogyridae and Ancyrocephalinae. The phylogeny of the two subfamilies shall have to be further evaluated with more data.

Characterization of the Complete Mitochondrial Genome of Ostertagia trifurcata of Small Ruminants and its Phylogenetic Associations for the Trichostrongyloidea Superfamily

The complete mitochondrial (mt) genome of Ostertagia trifurcata, a parasitic nematode of small ruminants, has been sequenced and its phylogenetic relationship with selected members from the superfamily Trichostrongyloidea was investigated on the basis of deduced datasets of mt amino acid sequences. The entire mt genome of Ostertagia trifurcata is circular and 14,151 bp in length. It consists of a total of 36 genes comprising 12 genes coding for proteins (PCGs), 2 genes for ribosomal RNA (rRNA), 22 transfer RNA (tRNA) genes and 2 non-coding regions, since all genes are transcribed in the same direction. The phylogenetic analysis based on the concatenated datasets of predicted amino acid sequences of the 12 protein coding genes supported monophylies of the Haemonchidae, Dictyocaulidae and Molineidae families, but rejected monophylies of the Trichostrongylidae family. The complete characterization and provision of the mtDNA sequence of Ostertagia trifurcata provides novel genetic markers for molecular epidemiological investigations, systematics, diagnostics and population genetics of Ostertagia trifurcata and its correspondents.

Three new Diplozoidae mitogenomes expose unusual compositional biases within the Monogenea class: implications for phylogenetic studies

BACKGROUND

As the topologies produced by previous molecular and morphological studies were contradictory and unstable (polytomy), evolutionary relationships within the Diplozoidae family and the Monogenea class (controversial relationships among the Discocotylinea, Microcotylinea and Gastrocotylinea suborders) remain unresolved. Complete mitogenomes carry a relatively large amount of information, sufficient to provide a much higher phylogenetic resolution than traditionally used morphological traits and/or single molecular markers. However, their implementation is hampered by the scarcity of available monogenean mitogenomes. Therefore, we sequenced and characterized mitogenomes belonging to three Diplozoidae family species, and conducted comparative genomic and phylogenomic analyses for the entire Monogenea class.

RESULTS

Taxonomic identification was inconclusive, so two of the species were identified merely to the genus level. The complete mitogenomes of Sindiplozoon sp. and Eudiplozoon sp. are 14,334 bp and 15,239 bp in size, respectively. Paradiplozoon opsariichthydis (15,385 bp) is incomplete: an approximately 2000 bp-long gap within a non-coding region could not be sequenced. Each genome contains the standard 36 genes (atp8 is missing). G + T content and the degree of GC- and AT-skews of these three mitogenome (and their individual elements) were higher than in other monogeneans. nad2, atp6 and nad6 were the most variable PCGs, whereas cox1, nad1 and cytb were the most conserved. Mitochondrial phylogenomics analysis, conducted using concatenated amino acid sequences of all PCGs, indicates that evolutionary relationships of the three genera are: (Eudiplozoon, (Paradiplozoon, Sindiplozoon)); and of the three suborders: (Discocotylinea, (Microcotylinea, Gastrocotylinea)). These intergeneric relationships were also supported by the skewness and principal component analyses.

CONCLUSIONS

Our results show that nad2, atp6 and nad6 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in Diplozoidae. Nucleotide bias and codon and amino acid usage patterns of the three diplozoid mitogenomes are more similar to cestodes and trematodes than to other monogenean flatworms. This unusual mutational bias was reflected in disproportionately long branches in the phylogram. Our study offsets the scarcity of molecular data for the subclass Polyopisthocotylea to some extent, and might provide important new insights into the evolutionary history of the three genera and three suborders.

Extensive gene rearrangements in the mitochondrial genomes of two egg parasitoids, Trichogramma japonicum and Trichogramma ostriniae (Hymenoptera: Chalcidoidea: Trichogrammatidae)

Animal mitochondrial genomes usually exhibit conserved gene arrangement across major lineages, while those in the Hymenoptera are known to possess frequent rearrangements, as are those of several other orders of insects. Here, we sequenced two complete mitochondrial genomes of Trichogramma japonicum and Trichogramma ostriniae (Hymenoptera: Chalcidoidea: Trichogrammatidae). In total, 37 mitochondrial genes were identified in both species. The same gene arrangement pattern was found in the two species, with extensive gene rearrangement compared with the ancestral insect mitochondrial genome. Most tRNA genes and all protein-coding genes were encoded on the minority strand. In total, 15 tRNA genes and seven protein-coding genes were rearranged. The rearrangements of cox1 and nad2 as well as most tRNA genes were novel. Phylogenetic analysis based on nucleotide sequences of protein-coding genes and on gene arrangement patterns produced identical topologies that support the relationship of (Agaonidae + Pteromalidae) + Trichogrammatidae in Chalcidoidea. CREx analysis revealed eight rearrangement operations occurred from presumed ancestral gene order of Chalcidoidea to form the derived gene order of Trichogramma. Our study shows that gene rearrangement information in Chalcidoidea can potentially contribute to the phylogeny of Chalcidoidea when more mitochondrial genome sequences are available.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Sequencing of the complete mitochondrial genomes of eight freshwater snail species exposes pervasive paraphyly within the Viviparidae family (Caenogastropoda)

Phylogenetic relationships among snails (Caenogastropoda) are still unresolved, and many taxonomic categories remain non-monophyletic. Paraphyly has been reported within a large family of freshwater snails, Viviparidae, where the taxonomic status of several species remains questionable. As many endemic Chinese viviparid species have become endangered during the last few decades, this presents a major obstacle for conservation efforts. Mitochondrial genomes (mitogenomes) carry a large amount of data, so they can often provide a much higher resolution for phylogenetic analyses in comparison to the traditionally used molecular markers. To help resolve their phylogenetic relationships, the complete mitogenomes of eight Chinese viviparid snails, Viviparus chui, Cipangopaludina chinensis, C. ussuriensis, C. dianchiensis (endangered), Margarya melanioides (endangered), M. monodi (critically endangered), Bellamya quadrata and B. aeruginosa, were sequenced and compared to almost all of the available caenogastropod mitogenomes. Viviparidae possess the largest mitogenomes (16 392 to 18 544 bp), exhibit the highest A+T bias (72.5% on average), and some exhibit unique gene orders (a rearrangement of the standard MYCWQGE box), among the Caenogastropoda. Apart from the Vermetidae family and Cerithioidea superfamily, which possessed unique gene orders, the remaining studied caenogastropod mitogenomes exhibited highly conserved gene order, with minimal variations. Maximum likelihood and Bayesian inference analyses, used to reconstruct the phylogenetic relationships among 49 almost complete (all 37 genes) caenogastropod mitogenomes, produced almost identical tree topologies. Viviparidae were divided into three clades: a) Margarya and Cipangopaludina (except C. ussuriensis), b) Bellamya and C. ussuriensis, c) Viviparus chui. Our results present evidence that some Cipangopaludina species (dianchiensis and cathayensis) should be renamed into the senior genus Margarya. The phylogenetic resolution obtained in this study is insufficient to fully resolve the relationships within the 'b' clade, but if C. chinensis proves to be a valid representative of the genus, C. ussuriensis may have to be reassigned a different genus (possibly Bellamya, or even a new genus). Non-monophyly also remains pervasive among the higher (above the family-level) Caenogastropod taxonomic classes. Gene order distance matrix produced a different phylogenetic signal from the nucleotide sequences, which indicates a limited usability of this approach for inferring caenogastropod phylogenies. As phenotypic homoplasy appears to be widespread among some viviparid genera, in order to effectively protect the rapidly diminishing endemic Viviparid populations in China, further detailed molecular phylogenetic studies are urgently needed to resolve the taxonomic status of several species.

C. Costa M, A. Andow D, Paula DP. Mitogenome sequence accuracy using different elucidation methods

Mitogenome sequences are highly desired because they are used in several biological disciplines. Their elucidation has been facilitated through the development of massive parallel sequencing, accelerating their deposition in public databases. However, sequencing, assembly and annotation methods might induce variability in their quality, raising concerns about the accuracy of the sequences that have been deposited in public databases. In this work we show that different sequencing methods (number of species pooled in a library, insert size and platform) and assembly and annotation methods generated variable completeness and similarity of the resulting mitogenome sequences, using three species of predaceous ladybird beetles as models. The identity of the sequences varied considerably depending on the method used and ranged from 38.19 to 90.1% for Cycloneda sanguinea, 72.85 to 91.06% for Harmonia axyridis and 41.15 to 93.60% for Hippodamia convergens. Dissimilarities were frequently found in the non-coding A+T rich region, but were also common in coding regions, and were not associated with low coverage. Mitogenome completeness and sequence identity were affected by the sequencing and assembly/annotation methods, and high within-species variation was also found for other mitogenome depositions in GenBank. This indicates a need for methods to confirm sequence accuracy, and guidelines for verifying mitogenomes should be discussed and developed by the scientific community.

Opportunities and challenges in metabarcoding approaches for helminth community identification in wild mammals

Despite metabarcoding being widely used to analyse bacterial community composition, its application in parasitological research remains limited. What interest there has been has focused on previously intractable research settings where traditional methods are inappropriate, for example, in longitudinal studies and studies involving endangered species. In settings such as these, non-invasive sampling combined with metabarcoding can provide a fast and accurate assessment of component communities. In this paper we review the use of metabarcoding in the study of helminth communities in wild mammals, outlining the necessary procedures from sample collection to statistical analysis. We highlight the limitations of the metabarcoding approach and speculate on what type of parasitological study would benefit from such methods in the future.

Complete Mitochondrial Genome of a Tongue Worm Armillifer agkistrodontis

Armillifer agkistrodontis (Ichthyostraca: Pantastomida) is a parasitic pathogen, only reported in China, which can cause a zoonotic disease, pentastomiasis. A complete mitochondrial (mt) genome was 16,521 bp comprising 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and 1 non-coding region (NCR). A phylogenetic tree drawn with the concatenated amino acid sequences of the 6 conserved PCGs (atp6, cox1-3, and nad2) showed that A. agkistrodontis and Armillifer armillatus constituted a clade Pentastomida which was a sister group of the Branchiura. The complete mt genome sequence of A. agkistrodontis provides important genetic markers for both phylogenetic and epidemiological studies of pentastomids.

Mitochondrial Phylogenomics yields Strongly Supported Hypotheses for Ascaridomorph Nematodes

Ascaridomorph nematodes threaten the health of humans and other animals worldwide. Despite their medical, veterinary and economic importance, the identification of species lineages and establishing their phylogenetic relationships have proved difficult in some cases. Many working hypotheses regarding the phylogeny of ascaridomorphs have been based on single-locus data, most typically nuclear ribosomal RNA. Such single-locus hypotheses lack independent corroboration, and for nuclear rRNA typically lack resolution for deep relationships. As an alternative approach, we analyzed the mitochondrial (mt) genomes of anisakids (~14 kb) from different fish hosts in multiple countries, in combination with those of other ascaridomorphs available in the GenBank database. The circular mt genomes range from 13,948-14,019 bp in size and encode 12 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNA genes. Our analysis showed that the Pseudoterranova decipiens complex consists of at least six cryptic species. In contrast, the hypothesis that Contracaecum ogmorhini represents a complex of cryptic species is not supported by mt genome data. Our analysis recovered several fundamental and uncontroversial ascaridomorph clades, including the monophyly of superfamilies and families, except for Ascaridiidae, which was consistent with the results based on nuclear rRNA analysis. In conclusion, mt genome analysis provided new insights into the phylogeny and taxonomy of ascaridomorph nematodes.

The complete mitochondrial genome of Orthocoelium streptocoelium (Digenea: Paramphistomidae) for comparison with other digeneans

Orthocoelium streptocoelium is a common paramphistome species parasitizing the rumen and/or reticulum of small ruminants, leading to significant losses. This study first determined the complete mitochondrial (mt) genome of O. streptocoelium. The complete mt genome of O. streptocoelium was amplified, sequenced, assembled, analysed and then compared with those of other digeneans. The entire mt genome of O. streptocoelium is 13,800 bp in length, which is smaller than those of other digeneans except for Opisthorchis viverrini. This mt genome contains 12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and two non-coding regions. The arrangement of the O. streptocoelium mt genome is the same as those of other digeneans except for Schistosoma haematobium and Schistosoma spindale. Phylogenetic analyses based on concatenated amino acid sequences of the 12 protein-coding genes representing 16 digeneans were conducted to assess the relationship of O. streptocoelium with other digeneans. The result indicated that O. streptocoelium is closely related to Paramphistomum cervi and Fischoederius elongates, which is in accordance with their relationships by taxonomy. This complete mt genome of O. streptocoelium enriched the mitochondrial genome data of paramphistomes and provided important molecular markers for diagnostics and studies of population variation, epidemiology, ecology and evolution of O. streptocoelium and other digeneans.

The complete mitochondrial genome of Gastrothylax crumenifer (Gastrothylacidae, Trematoda) and comparative analyses with selected trematodes

In the present study, we sequenced and analyzed the mitochondrial (mt) genome of Gastrothylax crumenifer and compared it with other selected trematodes. The full mt genome of G. crumenifer was amplified, sequenced, assembled, analyzed and then subjected to phylogenetic analysis. The complete mt genome of G. crumenifer is 14,801 bp in length and contains two rRNA genes, two non-coding regions (LNR and SNR), 12 protein-coding genes, and 22 transfer RNA genes. The gene organization of the G. crumenifer mt genome is the same as that of other trematodes, except for Schistosoma haematobium and Schistosoma spindale. All the genes are transcribed in the same direction and rich in "A + T", which is in accordance with other trematodes, such as Fasciola hepatica, Paramphistomum cervi, and Fischoederius elongatus. Phylogenetic analysis using concatenated amino acid sequences of the 12 protein-coding genes showed that G. crumenifer is closely related to F. elongatus. The availability of mt genome sequence of G. crumenifer can provide useful DNA markers for studying the molecular epidemiology and population genetics of this parasite and other paramphistomes.

Harnessing the Toxocara Genome to Underpin Toxocariasis Research and New Interventions

Parasitic worms, such as flatworms (platyhelminths) and roundworms (nematodes), cause substantial morbidity and mortality in animals and people globally. The ascaridoid nematode Toxocara canis is a zoonotic parasite of socioeconomic significance worldwide. In humans, this worm causes toxocariasis (disease) mainly in underprivileged communities in both the developed and developing worlds. While reasonably well studied from clinical and epidemiological perspectives, little is understood about the molecular biology of T. canis, its relationship with its hosts and the disease that it causes. However, a recent report of the draft genome and transcriptomes of T. canis should underpin many fundamental and applied research areas in the future. The present article gives a background on Toxocara and toxocariasis, a brief account of diagnostic approaches for specific identification and genetic analysis, and gives a perspective on the impact that the genome of T. canis and advanced molecular technologies could have on our understanding of the parasite and the diseases that it causes as well as the design of new and improved approaches for the diagnosis, treatment and control of toxocariasis.

Mitochondrial genome of Ogmocotyle sikae and implications for phylogenetic studies of the Notocotylidae trematodes

Ogmocotyle spp. (Trematoda: Digenea: Notocotylidae) are neglected but important trematodes that can infect numerous mammal species, causing significant economic losses to livestock industries. However, there have been few studies on the molecular ecology of these trematodes. We amplified and sequenced the complete mitochondrial (mt) genome of Ogmocotyle sikae (14,307 bp). The gene content and gene arrangement of O. sikae mt genome was similar to that of Dicrocoelium chinensis, except that trnE and trnG were reversed. Phylogenetic analysis of O. sikae and selected parasites using Bayesian inference was performed based on concatenated amino acid sequence datasets conceptually translated from the 12 protein-coding genes. The results indicated that the family Notocotylidae is related to the family Paramphistomatidae. Our description of O. sikae mt genome provides a significant resource of molecular markers for future comparative studies of the Notocotylidae and other trematodes.

The complete mitochondrial genome of Anoplocephala perfoliata, the first representative for the family Anoplocephalidae

Background

Mitochondrial (mt) genome sequences are widely used to understand phylogenetic relationships among parasites. However, no complete mt genome sequence is available in the family Anoplocephalidae to date. This study sequenced and annotated the complete mt genome of Anoplocephala perfoliata (Anoplocephalidae), and investigated its phylogenetic relationships with other species from the families Hymenolepididae, Dipylidiidae and Taeniidae of the order Cyclophyllidea using the amino acid sequences of the 12 proteins in their mt genomes.

Methods

The complete mt genome of A. perfoliata was amplified by Long-range PCR, sequenced using primer walking and annotated by comparing with those of other cestodes. Its phylogenetic relationship with the species from the families Hymenolepididae, Dipylidiidae and Taeniidae was inferred using the 12 protein sequences based on Maximum likelihood and Bayesian methods.

Results

The complete circular mt genome sequence for A. perfoliata is 14,459 bp in size, and includes 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes. The mt gene arrangement of A. perfoliata is identical to those of previously reported Hymenolepis diminuta (Hymenolepididae) and Dipylidium caninum (Dipylidiidae), but slightly different from those of other taeniids due to an order switch between tRNA(S2) and tRNA(L1). The phylogenetic analyses showed that the Dipylidiidae was more closely related to Anoplocephalidae and Hymenolepididae than to Taeniidae. The relationship among the four families obtained by Maximum likelihood and Bayesian inferences based on predicted amino acid sequences of protein-coding genes is consistent with that based on their mt gene arrangement similarities.

Conclusions

This study determined the first mt genome for the family Anoplocephalidae, providing rich sources for selecting useful molecular markers for ecological and phylogenetic studies. Analyses on mt genome sequences of the four families of cestodes provide novel insights into their phylogenetic relationships. Of couse, more taxon sampling is necessary for future phylogenetic studies of these cestodes using mt genome sequences.

Electronic supplementary material

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

The complete mitochondrial genome of rabbit pinworm Passalurus ambiguus: genome characterization and phylogenetic analysis

Passalurus ambiguus (Nematda: Oxyuridae) is a common pinworm which parasitizes in the caecum and colon of rabbits. Despite its significance as a pathogen, the epidemiology, genetics, systematics, and biology of this pinworm remain poorly understood. In the present study, we sequenced the complete mitochondrial (mt) genome of P. ambiguus. The circular mt genome is 14,023 bp in size and encodes of 36 genes, including 12 protein-coding, two ribosomal RNA, and 22 transfer RNA genes. The mt gene order of P. ambiguus is the same as that of Wellcomia siamensis, but distinct from that of Enterobius vermicularis. Phylogenetic analyses based on concatenated amino acid sequences of 12 protein-coding genes by Bayesian inference (BI) showed that P. ambiguus was more closely related to W. siamensis than to E. vermicularis. This mt genome provides novel genetic markers for studying the molecular epidemiology, population genetics, systematics of pinworm of animals and humans, and should have implications for the diagnosis, prevention, and control of passaluriasis in rabbits and other animals.

The mitochondrial genome of Angiostrongylus mackerrasae as a basis for molecular, epidemiological and population genetic studies

Background

Angiostrongylus mackerrasae is a metastrongyloid nematode endemic to Australia, where it infects the native bush rat, Rattus fuscipes. This lungworm has an identical life cycle to that of Angiostrongylus cantonensis, a leading cause of eosinophilic meningitis in humans. The ability of A. mackerrasae to infect non-rodent hosts, specifically the black flying fox, raises concerns as to its zoonotic potential. To date, data on the taxonomy, epidemiology and population genetics of A. mackerrasae are unknown. Here, we describe the mitochondrial (mt) genome of A. mackerrasae with the aim of starting to address these knowledge gaps.

Methods

The complete mitochondrial (mt) genome of A. mackerrasae was amplified from a single morphologically identified adult worm, by long-PCR in two overlapping amplicons (8 kb and 10 kb). The amplicons were sequenced using the MiSeq Illumina platform and annotated using an in-house pipeline. Amino acid sequences inferred from individual protein coding genes of the mt genomes were concatenated and then subjected to phylogenetic analysis using Bayesian inference.

Results

The mt genome of A. mackerrasae is 13,640 bp in size and contains 12 protein coding genes (cox1-3, nad1-6, nad4L, atp6 and cob), and two ribosomal RNA (rRNA) and 22 transfer RNA (tRNA) genes.

Conclusions

The mt genome of A. mackerrasae has similar characteristics to those of other Angiostrongylus species. Sequence comparisons reveal that A. mackerrasae is closely related to A. cantonensis and the two sibling species may have recently diverged compared with all other species in the genus with a highly specific host selection. This mt genome will provide a source of genetic markers for explorations of the epidemiology, biology and population genetics of A. mackerrasae.

An RNA-based approach to sequence the mitogenome of Hypoptopoma incognitum (Siluriformes: Loricariidae)

Hypoptopoma incognitum is a fish of the fifth most species-rich family of vertebrates and abundant in rivers from the Brazilian Amazon. Only two species of Loricariidae fish have their complete mitogenomes sequence deposited in the Genbank. An innovative RNA-based approach was used to assemble the complete mitogenome of H. incognitum with an average coverage depth of 5292×. The typical vertebrate mitochondrial features were found; 22 tRNA genes, two rRNA genes, 13 protein-coding genes, and a non-coding control region. Moreover, the use of this approach allowed the measurement of mtRNA expression levels, the punctuation pattern of editing, and the detection of heteroplasmies.

Gnathostoma spinigerum Mitochondrial Genome Sequence: a Novel Gene Arrangement and its Phylogenetic Position within the Class Chromadorea

Human gnathostomiasis is an emerging food-borne parasitic disease caused by nematodes in the genus Gnathostoma. In spite of their significance as pathogens, these parasites remain poorly understood at the molecular level. In the present study, we sequenced the mitochondrial (mt) genome of G. spinigerum, which infects a range of definitive hosts including dogs, cats, tigers, leopards and humans. The mt genome of G. spinigerum is 14,079 bp in size and shows substantial changes in gene order compared to other nematodes studied to date. Phylogenetic analyses of mt genome sequences by Bayesian inference (BI) revealed that the infraorder Gnathostomatomorpha (represented by G. spinigerum) is closely related to the infraorder Ascaridomorpha. G. spinigerum is the first species from the infraorder Gnathostomatomorpha for which a complete mt genome has been sequenced. The new data will help understand the evolution, population genetics and systematics of this medically important group of parasites.

Mitochondrial genomes of Trichinella species and genotypes – a basis for diagnosis, and systematic and epidemiological explorations

In the present study we sequenced or re-sequenced, assembled and annotated 15 mitochondrial genomes representing the 12 currently recognised taxa of Trichinella using a deep sequencing-coupled approach. We then defined and compared the gene order in individual mitochondrial genomes (14 to 17.7 kb), evaluated genetic differences among species/genotypes and re-assessed the relationships among these taxa using the mitochondrial nucleic acid or amino acid sequence data sets. In addition, a rich source of mitochondrial genetic markers was defined that could be used in future systematic, epidemiological and population genetic studies of Trichinella. The sequencing-bioinformatic approach employed herein should be applicable to a wide range of eukaryotic parasites.

Analysis of the complete Fischoederius elongatus (Paramphistomidae, Trematoda) mitochondrial genome

Background

Fischoederius elongates is an important trematode of Paramphistomes in ruminants. Animals infected with F. elongates often don’t show obvious symptoms, so it is easy to be ignored. However it can cause severe economic losses to the breeding industry. Knowledge of the mitochondrial genome of F. elongates can be used for phylogenetic and epidemiological studies.

Findings

The complete mt genome sequence of F. elongates is 14,120 bp in length and contains 12 protein-coding genes, 22 tRNA genes, two rRNA genes and two non-coding regions (LNR and SNR). The gene arrangement of F. elongates is the same as other trematodes, such as Fasciola hepatica and Paramphistomum cervi. Phylogenetic analyses using concatenated amino acid sequences of the 12 protein-coding genes by Maximum-likelihood and Neighbor-joining analysis method showed that F. elongates was closely related to P. cervi.

Conclusion

The complete mt genome sequence of F. elongates should provide information for phylogenetic and epidemiological studies for F. elongates and the family Paramphistomidae.

The complete mitochondrial genome of Pseudoterranova azarasi and comparative analysis with other anisakid nematodes

Anisakiasis/anisakidosis caused by anisakid nematodes is an emerging infectious disease that can cause a wide range of clinical syndromes and are difficult to diagnose and treat in humans. In spite of their significance as pathogens, the systematics, genetics, epidemiology and biology of these parasites remain poorly understood. In the present study, we sequenced the complete mitochondrial (mt) genome of Pseudoterranova azarasi, which is one of the most important zoonotic anisakid parasites. The circular mt genome is 13,954 bp in size and encodes of 36 genes, including 12 protein-coding, 2 ribosomal RNA and 22 transfer RNA genes. The mt gene order of P. azarasi is the same as those of Ascaris spp. (Ascarididae), Toxocara spp. (Toxocaridae) and Anisakis simplex (Anisakidae), but distinct from those of Ascaridia spp. (Ascaridiidae) and Cucullanus robustus (Cucullanidae). Phylogenetic analyses based on concatenated amino acid sequences of 12 protein-coding genes by Bayesian inference (BI) showed that Pseudoterranova were more closely related to Anisakis than they were to Contracaecum with strong a posterior probability support. This mt genome provides a novel genetic markers for exploring cryptic/sibling species and host affiliations, and should have implications for the diagnosis, prevention and control of anisakidosis in humans.

A specific indel marker for the Philippines Schistosoma japonicum revealed by analysis of mitochondrial genome sequences

In the present study, near-complete mitochondrial (mt) genome sequences for Schistosoma japonicum from different regions in the Philippines and Japan were amplified and sequenced. Comparisons among S. japonicum from the Philippines, Japan, and China revealed a geographically based length difference in mt genomes, but the mt genomic organization and gene arrangement were the same. Sequence differences among samples from the Philippines and all samples from the three endemic areas were 0.57-2.12 and 0.76-3.85 %, respectively. The most variable part of the mt genome was the non-coding region. In the coding portion of the genome, protein-coding genes varied more than rRNA genes and tRNAs. The near-complete mt genome sequences for Philippine specimens were identical in length (14,091 bp) which was 4 bp longer than those of S. japonicum samples from Japan and China. This indel provides a unique genetic marker for S. japonicum samples from the Philippines. Phylogenetic analyses based on the concatenated amino acids of 12 protein-coding genes showed that samples of S. japonicum clustered according to their geographical origins. The identified mitochondrial indel marker will be useful for tracing the source of S. japonicum infection in humans and animals in Southeast Asia.

Mitochondrial and nuclear ribosomal DNA dataset supports that Paramphistomum leydeni (Trematoda: Digenea) is a distinct rumen fluke species

BACKGROUND

Rumen flukes parasitize the rumen and reticulum of ruminants, causing paramphistomiasis. Over the years, there has been considerable debate as to whether Paramphistomum leydeni and Paramphistomum cervi are the same or distant species.

METHODS

In the present study, the complete mitochondrial (mt) genome of P. leydeni was amplified using PCR-based sequencing and compared with that of P. cervi. The second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA (rDNA) of P. leydeni specimens (n = 6) and P. cervi specimens (n = 8) was amplified and then sequenced. Phylogenetic relationship of the concatenated amino acid sequence data for 12 protein-coding genes of the two rumen flukes and selected members of Trematoda was evaluated using Bayesian inference (BI).

RESULTS

The complete mt genome of P. leydeni was 14,050 bp in size. Significant nucleotide difference between the P. leydeni mt genome and that of P. cervi (14.7%) was observed. For genetic divergence in ITS-2, sequence difference between P. leydeni and P. cervi was 3.1%, while no sequence variation was detected within each of them. Phylogenetic analysis indicated that P. leydeni and P. cervi are closely-related but distinct rumen flukes.

CONCLUSIONS

Results of the present study support the proposal that P. leydeni and P. cervi represent two distinct valid species. The mt genome sequences of P. leydeni provide plentiful resources of mitochondrial markers, which can be combined with nuclear markers, for further comparative studies of the biology of P. leydeni and its congeners from China and other countries.

The complete mitochondrial genome of the gullet worm Gongylonema pulchrum: gene content, arrangement, composition and phylogenetic implications

Background

Gongylonema pulchrum (Nematoda: Gongylonematidae), a thread-like spirurid gullet worm, infects a range of mammalian definitive hosts, including cattle, pigs, equines, goats, primates and humans, and can cause gongylonemiasis.

Methods

In the present study, the complete mitochondrial (mt) genome of G. pulchrum was obtained using Long-range PCR and subsequent primer walking. The phylogenetic position of G. pulchrum within the Spiruromorpha was established using Bayesian analyses of the protein-coding genes at the amino acid level.

Results

The length of this AT-rich (75.94%) mt genome is 13,798 bp. It contains 12 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and one non-coding region. The gene arrangement is the same as those of Thelazia callipaeda (Thelaziidae) and Setaria digitata (Onchocercidae), but distinct from that of Heliconema longissimum (Physalopteridae). Phylogenetic analyses, based on the concatenated amino acid sequence data for all 12 protein-coding genes using Bayesian inference (BI) method, showed that G. pulchrum (Gongylonematidae) was more closely related to Spirocerca lupi (Spiruroidea) than other members of the infraorder Spiruromorpha.

Conclusions

The present study represents the first mt genome sequence for the family Gongylonematidae, which provides the opportunity to develop novel genetic markers for studies of epidemiology, population genetics and systematics of this nematode of human and animal health significance.

Electronic supplementary material

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

Mitochondrial genome of Hypoderaeum conoideum - comparison with selected trematodes

Background

Hypoderaeum conoideum is a neglected but important trematode. The life cycle of this parasite is complex: snails serve as the first intermediate hosts: bivalves, fishes or tadpoles serve as the second intermediate hosts, and poultry (such as chickens and ducks) act as definitive hosts. In recent years, H. conoideum has caused significant economic losses to the poultry industry in some Asian countries. Despite its importance, little is known about the molecular ecology and population genetics of this parasite. Knowledge of mitochondrial (mt) genome of H. conoideum can provide a foundation for phylogenetic studies as well as epidemiological investigations.

Methods

The entire mt genome of H. conoideum was amplified in five overlapping fragments by PCR and sequenced, annotated and compared with mt genomes of selected trematodes. A phylogenetic analysis of concatenated mt amino acid sequence data for H. conoideum, eight other digeneans (Clonorchis sinensis, Fasciola gigantica, F. hepatica, Opisthorchis felineus, Schistosoma haematobium, S. japonicum, S. mekongi and S. spindale) and one tapeworm (Taenia solium; outgroup) was conducted to assess their relationships.

Results

The complete mt genome of H. conoideum is 14,180 bp in length, and contains 12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and one non-coding region (NCR). The gene arrangement is the same as in Fasciola spp, with all genes being transcribed in the same direction. The phylogenetic analysis showed that H. conoideum had a relatively close relationship with F. hepatica and other members of the Fasciolidae, followed by the Opisthorchiidae, and then the Schistosomatidae.

Conclusions

The mt genome of H. conoideum should be useful as a resource for comparative mt genomic studies of trematodes and for DNA markers for systematic, population genetic and epidemiological studies of H. conoideum and congeners.

Sequencing and annotation of mitochondrial genomes from individual parasitic helminths

Mitochondrial (mt) genomics has significant implications in a range of fundamental areas of parasitology, including evolution, systematics, and population genetics as well as explorations of mt biochemistry, physiology, and function. Mt genomes also provide a rich source of markers to aid molecular epidemiological and ecological studies of key parasites. However, there is still a paucity of information on mt genomes for many metazoan organisms, particularly parasitic helminths, which has often related to challenges linked to sequencing from tiny amounts of material. The advent of next-generation sequencing (NGS) technologies has paved the way for low cost, high-throughput mt genomic research, but there have been obstacles, particularly in relation to post-sequencing assembly and analyses of large datasets. In this chapter, we describe protocols for the efficient amplification and sequencing of mt genomes from small portions of individual helminths, and highlight the utility of NGS platforms to expedite mt genomics. In addition, we recommend approaches for manual or semi-automated bioinformatic annotation and analyses to overcome the bioinformatic "bottleneck" to research in this area. Taken together, these approaches have demonstrated applicability to a range of parasites and provide prospects for using complete mt genomic sequence datasets for large-scale molecular systematic and epidemiological studies. In addition, these methods have broader utility and might be readily adapted to a range of other medium-sized molecular regions (i.e., 10-100 kb), including large genomic operons, and other organellar (e.g., plastid) and viral genomes.

The mitochondrial genome of Parascaris univalens--implications for a "forgotten" parasite

Background

Parascaris univalens is an ascaridoid nematode of equids. Little is known about its epidemiology and population genetics in domestic and wild horse populations. PCR-based methods are suited to support studies in these areas, provided that reliable genetic markers are used. Recent studies have shown that mitochondrial (mt) genomic markers are applicable in such methods, but no such markers have been defined for P. univalens.

Methods

Mt genome regions were amplified from total genomic DNA isolated from P. univalens eggs by long-PCR and sequenced using Illumina technology. The mt genome was assembled and annotated using an established bioinformatic pipeline. Amino acid sequences inferred from all protein-encoding genes of the mt genomes were compared with those from other ascaridoid nematodes, and concatenated sequences were subjected to phylogenetic analysis by Bayesian inference.

Results

The circular mt genome was 13,920 bp in length and contained two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with those of other ascaridoids. Phylogenetic analysis of the concatenated amino acid sequence data for the 12 mt proteins showed that P. univalens was most closely related to Ascaris lumbricoides and A. suum, to the exclusion of other ascaridoids.

Conclusions

This mt genome representing P. univalens now provides a rich source of genetic markers for future studies of the genetics and epidemiology of this parasite and its congener, P. equorum. This focus is significant, given that there is no published information on the specific prevalence and distribution of P. univalens infection in domestic and wild horse populations.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-3305-7-428) contains supplementary material, which is available to authorized users.

Characterization of the complete mitochondrial genomes of Nematodirus oiratianus and Nematodirus spathiger of small ruminants

BACKGROUND

Nematodirus spp. are among the most common nematodes of ruminants worldwide. N. oiratianus and N. spathiger are distributed worldwide as highly prevalent gastrointestinal nematodes, which cause emerging health problems and economic losses. Accurate identification of Nematodirus species is essential to develop effective control strategies for Nematodirus infection in ruminants. Mitochondrial DNA (mtDNA) could provide powerful genetic markers for identifying these closely related species and resolving phylogenetic relationships at different taxonomic levels.

METHODS

In the present study, the complete mitochondrial (mt) genomes of N. oiratianus and N. spathiger from small ruminants in China were obtained using Long-range PCR and sequencing.

RESULTS

The complete mt genomes of N. oiratianus and N. spathiger were 13,765 bp and 13,519 bp in length, respectively. Both mt genomes were circular and consisted of 36 genes, including 12 genes encoding proteins, 2 genes encoding rRNA, and 22 genes encoding tRNA. Phylogenetic analyses based on the concatenated amino acid sequence data of all 12 protein-coding genes by Bayesian inference (BI), Maximum likelihood (ML) and Maximum parsimony (MP) showed that the two Nematodirus species (Molineidae) were closely related to Dictyocaulidae.

CONCLUSIONS

The availability of the complete mtDNA sequences of N. oiratianus and N. spathiger not only provides new mtDNA sources for a better understanding of nematode mt genomics and phylogeny, but also provides novel and useful genetic markers for studying diagnosis, population genetics and molecular epidemiology of Nematodirus spp. in small ruminants.