- Open Access
The complete mitochondrial genome of Anoplocephala perfoliata, the first representative for the family Anoplocephalidae
© Guo. 2015
- Received: 8 August 2015
- Accepted: 15 October 2015
- Published: 21 October 2015
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.
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.
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.
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.
- Anoplocephala perfoliata
- Mitochondrial genome
- Phylogenetic analysis
The Cyclophyllidea order with 15 families is the biggest taxonomic group in Cestoda. The phylogenies of cyclophyllidean tapeworms have been widely studied based on morphological characters [1, 2] which are often convergent and thus of little significance [3, 4]. For example, the apparent similarities are observed between the intermicapsiferine anoplocephalids and the davaineids with fibrous polyovular egg capsules. The Anoplocephalidae is diagnosed in part by the absence of a rostellum when compared with the davaineids. However, the rostellum has been secondarily lost in numerous species or genera within the Davaineidae .
Molecular tools, using specific genes and genomic regions (12S rDNA, 28S rDNA, cox1, etc.) [3, 5–10], have been used to study phylogenetic relationships among groups of parasitic flatworms. However, some specific genes and genomic regions such as rDNA genes, may not completely provide the useful phylogenetic information at the generic level [6, 11], or may not be the most suitable marker for all taxa [12, 13]. Nevertheless, complete mitochondrial DNA sequences have become effective analytical tools in ecology, evolutionary biology and population genetics [14–18]. The complete mt genomes provide not only individual and combined mtDNA molecular markers for within- and between-species studies but also the information of gene order changes . Platyhelminthe mt genomes comprise 12 protein-coding genes, 2 ribosomal and 22 transfer RNA genes , but lack the atp8 gene. In general, the mitochondrial gene order is rather stable and gene rearrangements are rare evolutionary events. When any common change is discovered within a group, it is considered to be a phylogenetic signal [15, 20, 21]. Although mt genome data provide useful molecular markers for population genetics and phylogenetic studies, there are some families in the cestode order Cyclophyllidea for which mtDNAs are still unavailable to date. Such, few studies on the phylogeny of cyclophyllidean cestodes are conducted at the family level based on complete mitochondrial DNA sequences [11, 14].
Complete mitochondrial genome is unavailable for the Anoplocephalidae family as yet and its systemic relationships with other cestodes in the same order of Cyclophyllidea are currently poorly understood. Several species that belong to the family Anoplocephalidae can cause tapeworm infections in horses and donkeys . Among them, Anoplocephala perfoliata is a notable one with a worldwide distribution . In the present study, the complete mitochondrial genome of A. perfoliata (Anoplocephalidae) has been sequenced and annotated, which is the first representative of the family Anoplocephalidae. The A. perfoliata mt DNA sequence was used along with those of nine published mt genomes of species of cestode families in the Cyclophyllidea, including Hymenolepis diminuta (the only complete mt genome available in Hymenolepididae), Dipylidium caninum (the only complete mt genome available in Dipylidiidae), Versteria mustelae, Echinococcus granulosus, Echinococcus multilocularis, Taenia solium, Taenia asiatica, Hydatigera krepkogorski and Hydatigera parva (seven representative mt genome sequences from the four genera of Taeniidae) to explore family-level phylogenetic relationships of the four families Anoplocephalidae, Taeniidae, Dipylidiidae and Hymenolepididae.
Sampling and DNA extraction
Single A. perfoliata adult was collected from the digestive tract of a donkey slaughtered at a commercial abattoir in China for DNA extraction and sequencing. The tapeworm was washed with 0.75 % NaCl. The genomic DNA was extracted using Multisource Genomic Miniprep DNA extraction kit (AXYGEN, USA) following the manufacturer’s instructions.
Amplification of mtDNA fragments
The complete mitochondrial genomic DNA was amplified by PCR in three overlapping fragments (nad1-rrnS, rrnS-nad5, nad5-nad1) using three pairs of primers (Additional file 1: Table S4), which were designed according to the published primers  used for amplifying short fragments of mtDNAs for all Taenia species. PCR reactions were carried out in a 50 μl reaction mixture consisting of 25 μl 2 × buffer (Takara Biotechnology Co. Ltd., Dalian, China), 22.5 μl sterile deionized water, 0.5 μl each primer (50 pmol/μl), and 1.5 μl DNA template (40 ng/μl) in the following conditions: 94 °C for 5 min, then 94 °C for 30 s, 50 °C for 30 s, 68 °C for 8 min for 30 cycles, followed by 68 °C for 10 min. PCR products were purified and underwent sequencing. Each amplicon was sequenced in both directions by primer walking  with ABI PRISM 3730 (Applied Biosystems Inc.) at Takara Biotechnology Co. and Shanghai Sangon Co.
The complete mt genomic sequence of A. perfoliata was assembled using CAP3 Sequence Assembly Program  (http://doua.prabi.fr/software/cap3) and annotated by sequence comparison with those of other cestodes [14, 16] available in the GenBank database. Protein-coding genes were identified using BLAST searches and ORF Finder software at NCBI with the genetic code set for Flatworm (Translation table 9). Transfer RNA (tRNA) genes were annotated using tRNAscan-SE, combining with recognizing anticodon sequences and potential secondary structures by visual inspection. The ribosomal RNA genes (rRNA) were identified by sequence comparisons with the published cestode mt rRNA genes. Palindromes and inverted sequences in the non-coding regions were determined using Einverted and Palindrome in EMBOSS software . Secondary structures were predicted using Mfold software . Tandem repetitive elements were detected with Tandem Repeats Finder program . Nucleotide composition (%) of each gene, non-coding region and the complete mt genome was calculated using DNASTAR’s Lasergene sequence analysis software .
A total of 16 mt genome sequences were used for analysis, including seven representative mt genome sequences from the four genera of Taeniidae, and two cyclophyllidean cestode mt genomes available in GenBank (Hymenolepis diminuta (Hymenolepididae), and Dipylidium caninum (Dipylidiidae)), the mt genome sequence for A. perfoliata identified in this study, five mt genomes from Pseudophyiidea (Spirometra erinaceieuropaei, Diphyllobothrium nihonkaiense, Kiphyllobothrium latum, Diplogonoporus balaenopterae and Diplogonoporus grandis), and mt sequences from the trematode Schistosoma japonicum (outgroup) (Additional file 1: Table S1). Phylogenetic analyses were carried out using amino acid sequences, which were predicted from the nucleotide sequences of each protein-coding gene using the flatworm mt genetic code set (Translation table 9). The amino acid sequences for each protein-coding gene of 16 species were aligned individually using MAFFT 7.122 . Poorly aligned positions were removed by Gblocks server (http://molevol.cmima.csic.es/castresana/Gblocks_server.html)  using the option for a less stringent selection. A single alignment for phylogenetic analysis was created by concatenating all amino acid alignments of the 12 protein-coding genes. The optimal model of protein evolution was selected by ProtTest . The phylogenetic trees were constructed using Bayesian analysis and Maximum-likelihood (ML). Bayesian analysis was performed by MrBayes 3.2  with mtZoa model (rates = gamma, ngammacat = 5), as suggested by Rota-Stabelli et al. . Two chains (temp = 0.2) were run for 5,000,000 generations and sampled every 1,000 generations. The first 25 % of the trees were treated as burn-in. ML analysis were conducted with PhyML 3.0  using the MtArt + I + G model, and 100 bootstrap replicates were selected to calculated bootstrap support for ML trees.
The complete mt genome organization of A. perfoliata
List of annotated mitochondrial genome of Anoplocephala perfoliata
Initiation and termination codons
Number of aa
Intergenic nucleotides (bp)a
Non-coding region (NC1)
Non-coding region (NC2)
For the 12 protein-coding genes identified in A. perfoliata, the cytb (1,101 bp) is relatively large in size, whereas cox1 (1,593 bp) is among the smallest and the remaining ten genes are comparable with the normative sizes, as with the other cestodes mt genomes thus far characterized (Additional file 1: Table S3). Initiation and termination codons inferred from each protein-coding gene are shown in Table 1. Of 12 protein-coding genes, eleven genes were inferred to initiate with the ATG codon except that cytb is probably initiated with GTG. Seven were inferred to use TAG as stop codons (cox2, nad6, nad5, cytb, nad4L, atp6 and nad2), four were inferred to terminate with TAA (cox1, nad4, nad1 and nad3) and one gene (cox3) was terminated with an incomplete stop codon (TA). The truncated incomplete stop codon has also been found in cox3 gene for D. caninum. The initiation and complete or incomplete termination codons for A. perfoliata are as normal for flatworm mitochondrial protein-coding genes (Additional file 1: Table S3). Gene overlaps were found between tRNA(F) and tRNA(Q), between tRNA(Q) and tRNA(M), and between tRNA(V) and tRNA(A). The nad4 gene is inferred to overlap with the upstream nad4L gene by 34 nucleotides, which is common for these two genes among metazoan mtDNAs .
Ribosomal RNA and transfer RNA genes
The two ribosomal RNA genes, rrnS and rrnL, were predicted from the mt genome of A. perfoliata. rrnS is located between tRNA(C) and cox2, rrnL is found between tRNA(T) and tRNA(C), the same position as reported in other cestodes [14, 16]. The length of rrnL (992 bp) is relatively larger and rrnS (724 bp) is normal, compared to those found in some other cestode species (Additional file 1: Table S3). The A + T contents of rrnS and rrnL are 69.34 % and 70.67 %, respectively.
The mt genome of A. perfoliata contained two major non-coding regions (NC1 or NC2) (Table 1). NC1 (875 bp) is located between tRNA(Y) and tRNA(S2). NC2 (279 bp) is found between nad5 and tRNA(G). These are at the same relative positions as in the mt genome of H. diminuta , but NC1 is larger than NC2 in the A. perfoliata mt genome, whereas NC1 is smaller than NC2 in the H. diminuta mt genome (Additional file 1: Table S3). These differences in length between the two non-coding regions are less pronounced in cestodes and may not be used as phylogenetic signals . For the other reported taeniid mtDNAs , however, the two non-coding regions are placed between tRNA(Y) and tRNA(L1), and between nad5 and tRNA(G), respectively.
The second non-coding region (NC2) of the mt genome of A. perfoliata comprises six identical repeat elements of a 32 bp nucleotide plus one shuffling repeat element (Fig. 3b). Except for these repeat sequences, the other non-coding sequence with 52 bp in NC2 can form a stem-loop structure with 15 standard base pairs and a loop of 16 bp nucleotides (Fig. 3b). This structure, including identical repeats and a stem-loop in its non-coding region, is only observed in the published mt genome sequences of H. diminuta and Meloidogyne javanica, which have been predicted for 13 repeat elements with a stem-loop and 11 repeat elements with a stem-loop, respectively [16, 38].
The lengths of non-coding regions in other cestode mt genomes are shown in Additional file 1: Table S3. Similar and stable stem-loop structures can be also found in the non-coding regions of taeniid mtDNAs , though the functions of these non-coding regions remain unclear. In mammals, however, similar stable hairpin structures have been shown to initiate replication and transcription [39, 40]. Whether these non-coding regions of cestode mtDNAs serve a similar function remains to be investigated. Other small intergenic regions were discovered between tRNA(S2) and tRNA(L1) (30 bp) and between nad1 and tRNA(N) (26 bp). Whether any of these non-coding regions serve any function awaits investigation.
The present study is the first one to investigate the relationship among these four families using concatenated mt amino acid sequences of 12 protein-coding genes. The relationship among these four families obtained here shows that the Dipylidiidae is more closely related to the Hymenolepididae and the Anoplocephalidae than to the Taeniidae. This can also be supported by the analysis based on similarities of their mt gene arrangement, where the gene arrangement of A. perfoliata (Anoplocephalidae) mt genome is identical to those of the species of H. diminuta (Hymenolepididae) and D. caninum (Dipylidiidae), but different from those of other taeniids. The identical mt gene arrangement among the species of the families Anoplocephalidae, Hymenolepididae and Dipylidiidae supported the result from the Bayesian phylogeny based on concatenated amino acids that these three families might have a closer relationship. Additionally, sequence organization for one of the two non-coding regions in A. perfoliata is very similar to that of H. diminuta, where the front part of the non-coding region comprises 6 identical repeat elements and the remaining non-coding region sequence can form a stem-loop structure in A. perfoliata (Fig. 3b), and 13 repeat elements with a stem-loop structure in H. diminuta [16, 38]. Whether the similar feature of the non-coding region between A. perfoliata and H. diminuta can be used as phylogenetic marker to support the closer relationship between the Anoplocephalidae and the Hymenolepididae is not known yet, a wider sampling will be required to determine. In comparison with other cyclophyllidean tapeworms described earlier [14, 16, 18, 43], mitochondrial genome organization and structure of A. perfoliata are more similar with them, supporting the result from the Bayesian phylogeny that the family Anoplocephalidae belongs to the order Cyclophyllidea. Considering the limited mt genome data available for these cestode parasites, additional taxon sampling of the families in the order Cyclophyllidea including the Anoplocephalidae, the Hymenolepididae and the Dipylidiidae is required to fully assess of the inferred phylogenetic relationship.
In conclusion, this study details the complete sequencing and annotating of the mt genome of A. perfoliata. This represents the first complete mtDNA sequence from the family Anoplocephalidae. The complete mtDNA genome contains 36 genes, including 12 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes. Phylogenetic analysis inferred from amino acid sequences of 12 protein-coding genes from ten representatives belonging to four families (Hymenolepididae, Dipylidiidae, Taeniidae and Anoplocephalidae) showed that the Dipylidiidae is more closely related to the Hymenolepididae and the Anoplocephalidae than to the Taeniidae. This can also be inferred from the similarities of their mt gene arrangement. This finding provides novel insights into the phylogenetic hypothesis of the four major families of the order Cyclophyllidea. The complete A. perfoliata mt genome determined in the present study will provide novel mt genetic markers for studies of the equine tapeworm in epidemiology, molecular identification and diagnosis, and also help to reassess the systematic relationships of cestodes in the order Cyclophyllidea.
This study was supported by the “National Key Basic Research Program (973 Program) of China” (Grant No. 2015CB150300). The author would like to express her appreciation to Dr. Shuai Wang, Ms Blinda Brice, and Dr. Guo-Hua Liu for critical reading and advice on this manuscript.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hoberg EP, Jones A, Bray RA. Phylogenetic analysis among the families of the Cyclophyllidea (Eucestoda) based on comparative morphology, with new hypotheses for co-evolution in vertebrates. Syst Parasitol. 1999;42(1):51–73.View ArticlePubMedGoogle Scholar
- Hoberg EP, Gardner SL, Campbell RA. Systematics of the Eucestoda: advances toward a new phylogenetic paradigm, and observations on the early diversification of tapeworms and vertebrates. Syst Parasitol. 1999;42(1):1–12.View ArticlePubMedGoogle Scholar
- von Nickisch-Rosenegk M, Lucius R, Loos-Frank B. Contributions to the phylogeny of the Cyclophyllidea (Cestoda) inferred from mitochondrial 12S rDNA. J Mol Evol. 1999;48(5):586–96.View ArticleGoogle Scholar
- Mariaux J. Cestode systematics: any progress? Int J Parasitol. 1996;26(3):231–43.View ArticlePubMedGoogle Scholar
- Haukisalmi V, Hardman LM, Hardman M, Laakkonen J, Niemimaa J, Henttonen H. Morphological and molecular characterisation of Paranoplocephala buryatiensis n. sp. and P. longivaginata Chechulin & Gulyaev, 1998 (Cestoda: Anoplocephalidae) in voles of the genus Clethrionomys. Syst Parasitol. 2007;66(1):55–71.View ArticlePubMedGoogle Scholar
- Taleb-Hossenkhan N, Bhagwant S. Molecular characterization of the parasitic tapeworm Bertiella studeri from the island of Mauritius. Parasitol Res. 2012;110(2):759–68.View ArticlePubMedGoogle Scholar
- Mariaux J. A molecular phylogeny of the Eucestoda. J Parasitol. 1998;84(1):114–24.View ArticlePubMedGoogle Scholar
- Olson PD, Littlewood DT, Bray RA, Mariaux J. Interrelationships and evolution of the tapeworms (Platyhelminthes: Cestoda). Mol Phylogenet Evol. 2001;19(3):443–67.View ArticlePubMedGoogle Scholar
- Olson PD, Tkach VV. Advances and trends in the molecular systematics of the parasitic Platyhelminthes. Adv Parasitol. 2005;60:165–243.View ArticlePubMedGoogle Scholar
- Lee SU, Chun HC, Huh S. Molecular phylogeny of parasitic Platyhelminthes based on sequences of partial 28S rDNA D1 and mitochondrial cytochrome c oxidase subunit I. Korean J Parasitol. 2007;45(3):181–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakao M, Lavikainen A, Iwaki T, Haukisalmi V, Konyaev S, Oku Y, et al. Molecular phylogeny of the genus Taenia (Cestoda: Taeniidae): proposals for the resurrection of Hydatigera Lamarck, 1816 and the creation of a new genus Versteria. Int J Parasitol. 2013;43(6):427–37.View ArticlePubMedGoogle Scholar
- Will KW, Mishler BD, Wheeler QD. The perils of DNA barcoding and the need for integrative taxonomy. Syst Biol. 2005;54(5):844–51.View ArticlePubMedGoogle Scholar
- Zarowiecki MZ, Huyse T, Littlewood DT. Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea). Int J Parasitol. 2007;37(12):1401–18.View ArticlePubMedGoogle Scholar
- Jia WZ, Yan HB, Guo AJ, Zhu XQ, Wang YC, Shi WG, et al. Complete mitochondrial genomes of Taenia multiceps, T. hydatigena and T. pisiformis: additional molecular markers for a tapeworm genus of human and animal health significance. BMC Genomics. 2010;11:447.PubMed CentralView ArticlePubMedGoogle Scholar
- Littlewood DT, Lockyer AE, Webster BL, Johnston DA, Le TH. The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms. Mol Phylogenet Evol. 2006;39(2):452–67.View ArticlePubMedGoogle Scholar
- von Nickisch-Rosenegk M, Brown WM, Boore JL. Complete sequence of the mitochondrial genome of the tapeworm Hymenolepis diminuta: gene arrangements indicate that Platyhelminths are Eutrochozoans. Mol Biol Evol. 2001;18(5):721–30.View ArticleGoogle Scholar
- Jeon HK, Lee KH, Kim KH, Hwang UW, Eom KS. Complete sequence and structure of the mitochondrial genome of the human tapeworm, Taenia asiatica (Platyhelminthes; Cestoda). Parasitology. 2005;130(Pt 6):717–26.View ArticlePubMedGoogle Scholar
- Jeon HK, Kim KH, Eom KS. Complete sequence of the mitochondrial genome of Taenia saginata: comparison with T. solium and T. asiatica. Parasitol Int. 2007;56(3):243–6.View ArticlePubMedGoogle Scholar
- Le TH, Blair D, McManus DP. Mitochondrial genomes of parasitic flatworms. Trends Parasitol. 2002;18(5):206–13.View ArticlePubMedGoogle Scholar
- Boore JL, Brown WM. Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998;8(6):668–74.View ArticlePubMedGoogle Scholar
- Rokas A, Holland PW. Rare genomic changes as a tool for phylogenetics. Trends Ecol Evol. 2000;15(11):454–9.View ArticlePubMedGoogle Scholar
- Bohorquez GA, Luzon M, Martin-Hernandez R, Meana A. New multiplex PCR method for the simultaneous diagnosis of the three known species of equine tapeworm. Vet Parasitol. 2015;207(1–2):56–63.View ArticlePubMedGoogle Scholar
- Bohorquez A, Meana A, Pato NF, Luzon M. Coprologically diagnosing Anoplocephala perfoliata in the presence of A. magna. Vet Parasitol. 2014;204(3–4):396–401.View ArticlePubMedGoogle Scholar
- Hu M, Jex AR, Campbell BE, Gasser RB. Long PCR amplification of the entire mitochondrial genome from individual helminths for direct sequencing. Nat Protoc. 2007;2(10):2339–44.View ArticlePubMedGoogle Scholar
- Huang X, Madan A. CAP3: A DNA sequence assembly program. Genome Res. 1999;9(9):868–77.PubMed CentralView ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16(6):276–7.View ArticlePubMedGoogle Scholar
- Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31(13):3406–15.PubMed CentralView ArticlePubMedGoogle Scholar
- Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Burland TG. DNASTAR's Lasergene sequence analysis software. Methods Mol Biol. 2000;132:71–91.PubMedGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56(4):564–77.View ArticlePubMedGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27(8):1164–5.View ArticlePubMedGoogle Scholar
- Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42.PubMed CentralView ArticlePubMedGoogle Scholar
- Rota-Stabelli O, Yang Z, Telford MJ. MtZoa: a general mitochondrial amino acid substitutions model for animal evolutionary studies. Mol Phylogenet Evol. 2009;52(1):268–72.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.View ArticlePubMedGoogle Scholar
- Lupi R, de Meo PD, Picardi E, D'Antonio M, Paoletti D, Castrignano T, et al. MitoZoa: a curated mitochondrial genome database of metazoans for comparative genomics studies. Mitochondrion. 2010;10(2):192–9.View ArticlePubMedGoogle Scholar
- Boore JL, Brown WM. Complete sequence of the mitochondrial DNA of the annelid worm Lumbricus terrestris. Genetics. 1995;141(1):305–19.PubMed CentralPubMedGoogle Scholar
- Okimoto R, Chamberlin HM, Macfarlane JL, Wolstenholme DR. Repeated sequence sets in mitochondrial DNA molecules of root knot nematodes (Meloidogyne): nucleotide sequences, genome location and potential for host-race identification. Nucleic Acids Res. 1991;19(7):1619–26.PubMed CentralView ArticlePubMedGoogle Scholar
- Monnerot M, Solignac M, Wolstenholme DR. Discrepancy in divergence of the mitochondrial and nuclear genomes of Drosophila teissieri and Drosophila yakuba. J Mol Evol. 1990;30(6):500–8.View ArticlePubMedGoogle Scholar
- Wolstenholme DR. Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992;141:173–216.View ArticlePubMedGoogle Scholar
- Foronda P, Casanova JC, Valladares B, Martinez E, Feliu C. Molecular systematics of several cyclophyllid families (Cestoda) based on the analysis of 18S ribosomal DNA gene sequences. Parasitol Res. 2004;93(4):279–82.View ArticlePubMedGoogle Scholar
- Yamasaki H, Ohmae H, Kuramochi T. Complete mitochondrial genomes of Diplogonoporus balaenopterae and Diplogonoporus grandis (Cestoda: Diphyllobothriidae) and clarification of their taxonomic relationships. Parasitol Int. 2012;61(2):260–6.View ArticlePubMedGoogle Scholar
- Nakao M, Yokoyama N, Sako Y, Fukunaga M, Ito A. The complete mitochondrial DNA sequence of the cestode Echinococcus multilocularis (Cyclophyllidea: Taeniidae). Mitochondrion. 2002;1(6):497–509.View ArticlePubMedGoogle Scholar