Analysis of the complete Fischoederius elongatus (Paramphistomidae, Trematoda) mitochondrial genome
- Xin Yang†1,
- Yunyang Zhao†1,
- Lixia Wang2,
- Hanli Feng3,
- Li Tan1,
- Weiqiang Lei1,
- Pengfei Zhao1,
- Min Hu1 and
- Rui Fang1Email author
© Yang et al.; licensee BioMed Central. 2015
Received: 17 March 2015
Accepted: 9 May 2015
Published: 20 May 2015
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.
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.
The complete mt genome sequence of F. elongates should provide information for phylogenetic and epidemiological studies for F. elongates and the family Paramphistomidae.
KeywordsFischoederius elongates Mitochondrial genome
Paramphistomes are distributed worldwide and have been reported in many countries, such as Bulgaria, France, Poland, Hungary, Italy, India, Russia, Sardinia and Yugoslavia . The paramphistome can infect fishes, reptiles, birds and mammals, some of which can lead to huge economic losses related to seriously gastrointestinal diseases, low producitivity or death in ruminants . In Arumeru District, the prevalence rate of paramphistomes is as high as 56.7 % in cattle .
Fischoederius elongates is an important member of paramphistomes, the parasite usually inhabits the rumen of cattle, buffaloes, sheep and goats. Ruminants are usually infected by ingesting snails, such as Lymnaea acuminata, Lymnaea succinea or Gyraulus euphraticus . Ruminants infected with F. elongates show weakness, mental fatigue and eventually death. More seriously, F. elongates maybe a zoonotic trematode, a Chinese woman from Guangdong Province was reported to be the first human infection case , but it is still unknown how she was infected.
Untill now, the most common diagnostic method for F. elongates is the microscopical examination, but it’s time-consuming, and hard to distinguish with other paramphistomes. As a useful marker, mt genome has been widely used for species identification [6–10]. The complete mt genome of F. elongates can provide alternative molecular markers for the species identification, epidemiology and genetic diversity of paramphistomes.
In the present study, we got the full sequence and gene arrangement of mt genome of F. elongates and compared it with selected trematodes. We found that F. elongates had the closest relationship with P. cervi.
The study was performed under the instructions and approval of Laboratory Animals Research Centre of Hubei province in P. R. China and the ethics committee of Huazhong Agricultural University (Permit number: 4200695757).
Parasite collection and DNA isolation
F. elongates adults were collected from the rumen and reticulum of naturally infected cattle in Zhanggang, Tianmen, Hubei province, PR China, according to the Animal Ethics Guidelines of Huazhong Agricultural University. Then, the adult worms were washed extensively in 0.9 % sodium chloride solution, and identified through morphological examinations . Subsequently, one worm was stained for identification , and the rest were fixed in 75 % alcohol (V/V) and stored at −20 °C until use . Total genomic DNA was isolated from one worm . The ITS-2 region of F. elongates was amplified and sequenced as reported previously , it was 100 % similar to that of F. elongates (GenBank accession no. JQ688410.1).
Amplification and sequencing of F. elongates mt genome
Primers used in the present study
Then, 12 additional primers (Table 1) were designed based on the obtained sequencing results to amplify six regions from genomic DNA (~40-80 ng) by long-PCR. PCRs (50 μl) were performed in reactions containing 0.4 mM each of dNTPs, 5 μl 10× LA Taq buffer II(Mg2+ Plus), 2.5 μM of each primer, 2.5 U LA Taq polymerase (Takara) and 2.5 μl genomic DNA. And the reactions were run under the following program: 94 °C for 5 min, followed by 35 cycles of 94 °C/30 s, 50 °C/30 s and 72 °C/1-5 min (depending on the size of F. hepatica). Amplicons were cloned into pGEM-T-Easy vector (Promega, USA) and then sequenced using a primer-walking strategy .
F. elongates mt genome sequences were assembled manually and then aligned with the mt genome sequences of F. hepatica, C. sinensis and P. cervi using the program Clustal X 1.83 . Open reading frames were identified by ORF Finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html) using the echinoderm and flatworm mitochondrial code. Initiation and termination codons of the 12 protein-coding genes were identified as reported . The 22 tRNA genes were predicted using tRNAscan-SE or manual adjustments [20,21]. The two rRNA genes were predicted by comparison with those of F. hepatica , C. sinensis  and P. cervi . Amino acid sequences of 12 protein-coding genes were inferred using ExPASy Translate tool (http://web.expasy.org/translate/) using the echinoderm and flatworm mitochondrial codes, and aligned using MEGA 5.0 with default settings .
Nucleotide variation analysis
The nucleotide variation between F. elongates and P. cervi was analysed by sliding window analysis as reported .
Amino acid sequences translated from individual genes of the mt genome of F. elongates were aligned with those predicted from mt genomes of selected trematodes, including C. sinensis (NC_012147) , Dicrocoelium dendriticum (NC_025280.1) , F. hepatica (NC_002546) , Haplorchis taichui (NC_022433.1) , Metagonimus yokogawai (KC330755.1), Opisthorchis viverrini (JF739555.1) , P. cervi (NC_023095.1) , Schistosoma haematobium (NC_008074) , Schistosoma japonicum (AF215860) , Schistosoma mekongi (NC_002529) , Schistosoma spindale (NC_008067) , and the cestode Taenia solium (outgroup) (NC_004022.1) . The amino acid sequences of selected trematodes were aligned using MEGA 5.0 , and phylogenetic analysis of the aligned amino acid sequences was conducted in MEGA 5.0 using the Maximum Likelihood (ML) method.
Results and discussion
Features of the mt genome of F. elongates
The complete mitochondrial genome of F. elongates (GenBank accession no. KM_397348) is 14,120 bp in length. The length of the F. elongates mt genome is larger than the mtDNA genomes of C. sinensis (13,875 bp) and S. japonicum (14,085 bp), but smaller than D. dendriticum (14,884 bp), F. hepatica (14,462 bp), H. taichui (15,130 bp), M. yokogawai (15,258 bp), S. haematobium (15,003 bp), S. mekongi (14,072 bp) and S. spindale (16,901 bp).
The organization of the mitochondrial genome of Fischoederius elongatus
Number of aa1
Nucleotide contents of genes and the non-coding region within the mitochondrial genome of Fischoederius elongatus
A + T(%)
The present F. elongates mt genome can provide useful information for the studies of epidemiology, species identification and genetic diversity of Fischoederius spp. At the same, it will also make contribution to the taxonomy study of Fischoederius spp. With the full mt genome of F. elongates, we can undertake a study within F. elongates from different regions or among Fischoederius spp. by combining the morphological features with genetic analyses (with molecular markers from mitochondria or ribosome, such as cox1, nad4, 18S, ITS-1 and ITS-2). Meanwhile, the mt genome of F. elongates may also provide information for the prevention and diagnosis of Fischoederius spp. and perhaps, this mt genome information may assist in the new drug, since mitochondria is the target of some drugs, such as decoquinate.
Codon usage for 12 protein-coding genes in the mitochondrial genome of Fischoederius elongatus
Transfer RNA and ribosomal RNA genes
The F. elongates mt genome encodes 22 tRNAs, and the length of 22 tRNA genes ranged from 60 bp to 71 bp (Table 2). There are two non-coding regions in F. elongates mt genome, rrnS (751 bp) and rrnL (995 bp) (Table 2). The location of rrnS is between tRNA-Cys and cox2 and the rrnL is between tRNA-Thr and tRNA-Cys, which is the same as other trematodes, such as F. hepatica , C. sinensis  and P. cervi .
Many flatworms have non-coding regions, it’s common to find two non-coding regions in trematodes: one long non-coding region (LNR) and one short non-coding region (SNR). In F. elongates, there is a short non-coding region (SNR: 62 nucleotides), which is located between cytb and nad4L. In addition, there is also a long non-coding region (LNR: 468 nucleotides) between tRNA-Phe and cox3 (Table 2), the LNR has two obvious features, one is microsatellite-like sequences, such as (TA)n (n <5); the other is homopolymer sequences, such as (T)n (n <7). People still don’t understand clearly why the non-coding regions exist, and the function of them, people just knew the non-coding regions may participate in the replication of mitochondria .
Nucleotide variability between F. elongates and P. cervi
Sincere thanks to Professor Bang Shen for comments on the manuscript. This work was supported in part by the “National Key Basic Research Program (973 Program) of China” (Grant No. 2015CB150300), the “Special Fund for Agro-scientific Research in the Public Interest” (Grant No. 201303037) and “Huazhong Agricultural University Students Research Fund” (Grant No. 2015054).
- Horak IG. Paramphistomiasis of domestic ruminants. Adv Parasitol. 1971;9:33–72.View ArticlePubMedGoogle Scholar
- Li XR. Color atlas of animal parasitosis (Second Edition). Beijing: China Agriculture Press; 2011.Google Scholar
- Nzalawahe J, Kassuku AA, Stothard JR, Coles GC, Eisler MC. Trematode infections in cattle in Arumeru District, Tanzania are associated with irrigation. Parasit Vectors. 2014;7:107.View ArticlePubMed CentralPubMedGoogle Scholar
- Yamaguti S. A synoptical review of life histories of digenetic trematodes of vertebrates. Tokyo: Keigaku Publishing Co.; 1975.Google Scholar
- Li D. A case of Fischoederius elongatus infection in China. Annual Bull Soc Parasitol Guangdong Province. 1991;12(11–13):155–6.Google Scholar
- Ramesh A, Small ST, Kloos ZA, Kazura JW, Nutman TB, Serre D, et al. The complete mitochondrial genome sequence of the filarial nematode Wuchereria bancrofti from three geographic isolates provides evidence of complex demographic history. Mol Biochem Parasitol. 2012;183(1):32–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Gasser RB, Jabbar A, Mohandas N, Hoglund J, Hall RS, Littlewood DT, et al. Assessment of the genetic relationship between Dictyocaulus species from Bos taurus and Cervus elaphus using complete mitochondrial genomic datasets. Parasit Vectors. 2012;5:241.View ArticlePubMed CentralPubMedGoogle Scholar
- Choi KS, Koekemoer LL, Coetzee M. Population genetic structure of the major malaria vector Anopheles funestus s.s. and allied species in southern Africa. Parasit Vectors. 2012;5:283.View ArticlePubMed CentralPubMedGoogle Scholar
- Jane EH, Bernard F. Echinostoma and Echinostomiasis. Adv Parasitol. 1990;29:215–69.View ArticleGoogle Scholar
- Ghatani S, Shylla JA, Roy B, Tandon V. Multilocus sequence evaluation for differentiating species of the trematode Family Gastrothylacidae, with a note on the utility of mitochondrial COI motifs in species identification. Gene. 2014;548(2):277–84.View ArticlePubMedGoogle Scholar
- Wang GL. Laboratory diagnostic techniques of parasites. Feeding Livestock. 2013;3:39–43.Google Scholar
- Liu GH, Wang Y, Song HQ, Li MW, Ai L, Yu XL, et al. Characterization of the complete mitochondrial genome of Spirocerca lupi: sequence, gene organization and phylogenetic implications. Parasit Vectors. 2013;6:45.View ArticlePubMed CentralPubMedGoogle Scholar
- Gasser RB, Hu M, Chilton NB, Campbell BE, Jex AJ, Otranto D, et al. Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation. Nat Protoc. 2006;1(6):3121–8.View ArticlePubMedGoogle Scholar
- Itagaki T, Tsumagari N, Tsutsumi K, Chinone S. Discrimination of three amphistome species by PCR-RFLP based on rDNA ITS2 markers. J Vet Med Sci. 2003;65(8):931–3.View ArticlePubMedGoogle Scholar
- Le TH, Blair D, McManus DP. Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda). Parasitology. 2001;123(Pt 6):609–21.PubMedGoogle Scholar
- Shekhovtsov SV, Katokhin AV, Kolchanov NA, Mordvinov VA. The complete mitochondrial genomes of the liver flukes Opisthorchis felineus and Clonorchis sinensis (Trematoda). Parasitol Int. 2010;59(1):100–3.View ArticlePubMedGoogle Scholar
- Yan HB, Wang XY, Lou ZZ, Li L, Blair D, Yin H, et al. The mitochondrial genome of Paramphistomum cervi (Digenea), the first representative for the family Paramphistomidae. PLoS One. 2013;8(8):e71300.View ArticlePubMed CentralPubMedGoogle 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
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu M, Chilton NB, Gasser RB. The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea). Int J Parasitol. 2002;32(2):145–58.View ArticlePubMedGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu GH, Yan HB, Otranto D, Wang XY, Zhao GH, Jia WZ, et al. Dicrocoelium chinensis and Dicrocoelium dendriticum (Trematoda: Digenea) are distinct lancet fluke species based on mitochondrial and nuclear ribosomal DNA sequences. Mol Phylogenet Evol. 2014;79:325–31.View ArticlePubMedGoogle Scholar
- Lee D, Choe S, Park H, Jeon HK, Chai JY, Sohn WM, et al. Complete mitochondrial genome of Haplorchis taichui and comparative analysis with other trematodes. Korean J Parasitol. 2013;51(6):719–26.View ArticlePubMed CentralPubMedGoogle Scholar
- Cai XQ, Liu GH, Song HQ, Wu CY, Zou FC, Yan HK, et al. Sequences and gene organization of the mitochondrial genomes of the liver flukes Opisthorchis viverrini and Clonorchis sinensis (Trematoda). Parasitol Res. 2012;110(1):235–43.View 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
- Le TH, Blair D, Agatsuma T, Humair PF, Campbell NJ, Iwagami M, et al. Phylogenies inferred from mitochondrial gene orders-a cautionary tale from the parasitic flatworms. Mol Biol Evol. 2000;17(7):1123–5.View ArticlePubMedGoogle Scholar
- Nakao M, Sako Y, Ito A. The mitochondrial genome of the tapeworm Taenia solium: a finding of the abbreviated stop codon U. J Parasitol. 2003;89(3):633–5.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.