- Open Access
The complete mitochondrial genome of Anoplocephala perfoliata, the first representative for the family Anoplocephalidae
Parasites & Vectors volume 8, Article number: 549 (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.
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.
Results and discussion
The complete mt genome organization of A. perfoliata
The complete mt genome of A. perfoliata (GenBank: KR054960) is a circular, double-stranded DNA molecule with 14,459 bp in size (Fig. 1), which is slightly larger than those of other cestodes (13,482 bp–14,296 bp) [14, 36] (Additional file 1: Table S1). High A + T content (70.99 %) and low G + C content (29.01 %) were observed in the mt genome sequence of A. perfoliata (Additional file 1: Table S2). The overall A + T content is very similar to those of other cestode species reported so far (H. diminuta [71.04 %], D. caninum [72.89 %], T. solium [72.1 %], T. asiatica [71.4 %], E. multilocularis [69.04 %] and E. granulosus [67.02 %]). The 12 protein-coding genes (cox1-3, atp6, nad1-6, nad4L and cytb), 22 tRNA genes and two rRNA genes (rrnS and rrnL) were identified in A. perfoliata mt genome by comparing their sequence identities and/or secondary structures with those of other flatworms. The atp8 gene that appears in most animal mtDNAs is missing in this A. perfoliata isolate, as in other known flatworm mt genomes. All the genes are encoded by one strand in the same direction (Fig. 1). A list of gene arrangement, gene length and intergenic spacer regions of mt genome for A. perfoliata are shown in Table 1. It is noteworthy that the gene order of A. perfoliata mt genome is identical to those of H. diminuta and D. caninum, but slightly different from those of other taeniids due to an order switch between tRNA(S2) and tRNA(L1).
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 22 tRNA genes were identified in the mtDNA of A. perfoliata. The sizes of the tRNA genes are 57–73 bp (Table 1). Predicted secondary structures of 18 tRNAs have the standard cloverleaf structure in A. perfoliata, but that of four tRNA genes (tRNA(C), tRNA(R) and two tRNA(S)) lack the DHU arm (Fig. 2), as reported in some Taenia species mtDNAs . However, in the tRNAs of the H. diminuta mt genome, only three tRNAs (the tRNA(R) and the two tRNA(S) genes) lack the DHU arm. Each of the four tRNAs has the potential for an extra base-pairing, probably allowing a longer anticodon stem in A. perfoliata (Fig. 2). Similar tRNA structures with a longer anticodon stem are also found in Taenia species  and H. diminuta mtDNA . Additionally, it is noteworthy that an additional base-pairing for a longer acceptor stem may also be possible for the tRNA(C) in A. perfoliata (Fig. 2), yet no similar structure has been reported in other cestode mtDNAs.
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.
In the A. perfoliata mt genome, NC1 is mainly composed of four repeat elements of 169 bp, plus a part of the repeat element of 85 bp. Only two bases vary among these four repeat elements. The presence of long repeats is a common feature of non-coding regions for cestodes . A number of potential secondary structures in the NC1 of A. perfoliata were also predicted. For example, a stable potential secondary structure with two arms and two loops could be folded in each of these repeat regions (Fig. 3a). This secondary structure contains 11 and five canonical base pairs for the stems, with 3 bp and 5 bp loops, respectively (Fig. 3a). For NC1, three two-stem structures are feasible, but a large potential secondary structure presented as a pentamer is also possible with 24 canonical base pairs plus long loops (Fig. 3a). The large potential secondary structures predicted in triplex have been observed for non-coding regions between nad5 and tRNA(G) or between tRNA(Y) and tRNA(L1) in other cestode mtDNAs . Competitive and mutually exclusive hairpins could be present in this case, as reported in other mtDNAs [16, 37]. Besides the aforementioned major repeat regions, there were two sets of short inverted sequences at the NC1 terminals (TAATACTATTA and ATAAATTTAAAAATTTAT), and each of them could be folded into a hairpin structure.
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.
In order to evaluate the phylogenetic relationship among the families of the order Cyclophyllidea based on mt genome data, phylogeny was reconstructed using inferred amino acid sequences of the 12 protein-coding genes from 16 representatives including ten species from the order Cyclophyllidea, five species from the order Pseudophyllidea and one trematode outgroup. An identical phylogenetic tree topology was obtained through the Maximum-likelihood (ML) and Bayesian analyses (Fig. 4). The tree inferred using the ML method was the same as that of the Bayesian analysis and was omitted. The trees showed that Pseudophyllidea and Cyclophyllidea individuals formed a monophyletic group which is consistent with previous molecular analyses [11, 41, 42]. Within the Taeniidae group, four valid genera were robustly supported with high nodal support. This result is in agreement with the study by Nakao et al. , but different from those of other studies suggesting that the cestode family Taeniidae has only two valid genera . The tree revealed sister-group relationship between the Anoplocephalidae and the Hymenolepididae, which agrees with what has been reported using the 18S gene based on distance-based Neighbour-Joining method, but is not consistent with the parsimony-based phylogeny . This inferred relationship between the Anoplocephalidae and the Hymenolepididae is also consistent with some earlier morphological and molecular studies (12S rDNA and cox1) suggesting both families as very close parents [3, 5]. Furthermore, the trees indicated that the Dipylidiidae was phylogenetically closer to the Hymenolepididae and the Anoplocephalidae than to the Taeniidae. The evolutionary relationship between these four families revealed in the present study is incongruent with the hypothesis based on the phylogenetic tree derived from 12S rDNA  which suggests that the Dipylidiidae is more closely related to the Taeniidae than to either the Hymenolepididae or the Anoplocephalidae, and also disagrees with Hoberg’s proposal  that the Hymenolepididae is closer to the Dipylidiidae and the Taeniidae than to the Anoplocephalidae based on comparative morphology analysis.
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.
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.
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.
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.
Mariaux J. Cestode systematics: any progress? Int J Parasitol. 1996;26(3):231–43.
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.
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.
Mariaux J. A molecular phylogeny of the Eucestoda. J Parasitol. 1998;84(1):114–24.
Olson PD, Littlewood DT, Bray RA, Mariaux J. Interrelationships and evolution of the tapeworms (Platyhelminthes: Cestoda). Mol Phylogenet Evol. 2001;19(3):443–67.
Olson PD, Tkach VV. Advances and trends in the molecular systematics of the parasitic Platyhelminthes. Adv Parasitol. 2005;60:165–243.
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.
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.
Will KW, Mishler BD, Wheeler QD. The perils of DNA barcoding and the need for integrative taxonomy. Syst Biol. 2005;54(5):844–51.
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.
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.
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.
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.
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.
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.
Le TH, Blair D, McManus DP. Mitochondrial genomes of parasitic flatworms. Trends Parasitol. 2002;18(5):206–13.
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.
Rokas A, Holland PW. Rare genomic changes as a tool for phylogenetics. Trends Ecol Evol. 2000;15(11):454–9.
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.
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.
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.
Huang X, Madan A. CAP3: A DNA sequence assembly program. Genome Res. 1999;9(9):868–77.
Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16(6):276–7.
Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31(13):3406–15.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.
Burland TG. DNASTAR's Lasergene sequence analysis software. Methods Mol Biol. 2000;132:71–91.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
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.
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.
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.
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.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.
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.
Boore JL, Brown WM. Complete sequence of the mitochondrial DNA of the annelid worm Lumbricus terrestris. Genetics. 1995;141(1):305–19.
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.
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.
Wolstenholme DR. Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992;141:173–216.
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.
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.
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.
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.
The author declares that she has no competing interests.
Cestode species used for comparative analysis with Anoplocephala. perfoliata. Table S2. A + T content (%) of the protein-coding, tRNA, rRNA genes and non-coding regions of mitochondrial genome of Anoplocephala perfoliata. Table S3. Properties of mtDNA protein-coding genes, rRNA genes and non-coding regions of Anoplocephala perfoliata and other cestode species. Table S4. Primers used to amplify and sequence mitochondrial genome from Anoplocephala perfoliata. (DOC 209 kb)