Sequence variability in two mitochondrial DNA regions and internal transcribed spacer among three cestodes infecting animals and humans from China

Sequence variability in two mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 4 (nad4), and internal transcribed spacer (ITS) of rDNA among and within three cestodes, Spirometra erinaceieuropaei, Taenia multiceps and Taenia hydatigena, from different geographical origins in China was examined. A portion of the cox1 (pcox1), nad4 genes (pnad4) and the ITS (ITS1+5.8S rDNA+ITS2) were amplified separately from individual cestodes by polymerase chain reaction (PCR). Representative amplicons were subjected to sequencing in order to estimate sequence variability. While the intra-specific sequence variations within each of the tapeworm species were 0-0.7% for pcox1, 0-1.7% for pnad4 and 0.1-3.6% for ITS, the inter-specific sequence differences were significantly higher, being 12.1-17.6%, 18.7-26.2% and 31-75.5% for pcox1, pnad4 and ITS, respectively. Phylogenetic analyses based on the pcox1 sequence data revealed that T. multiceps and T. hydatigena were more closely related to the other members of the Taenia genus, and S. erinaceieuropaei was more closely related to the other members of the Spirometra genus. These findings demonstrated clearly the usefulness of mtDNA and rDNA sequences for population genetic studies of these cestodes of socio-economic importance.     

Sequence variability in three mitochondrial genes between the two pig nodule worms Oesophagostomum dentatum and O. quadrispinulatum

In this study, sequence variation in three mitochondrial DNA regions, namely cytochrome c oxidase subunit (cox1) and NADH dehydrogenase subunits 1 and 4 (nad1 and nad4), between Oesophagostomum dentatum and O. quadrispinulatum isolated from pigs in different geographical origins in Mainland China was examined, and their phylogenetic relationships were reconstructed. A partial of the cox1 (pcox1), nad1, and nad4 genes (pnad1 and pnad4) were amplified separately from individual nodule worms by PCR and were subjected to direct sequencing in order to define sequence variations. While the intraspecific sequence variations within each of the two species were 0.3-5.2% for pcox1, 0-4.9% for pnad1, and 0-7.1% for pnad4, the interspecific sequence differences were significantly higher, being 10.7-13.4% for pcox1, 11-14.6% for pnad1, and 14.9-18% for pnad4, respectively. There were a number of nucleotide positions in the pcox1, pnad1, and pnad4 sequences with no apparent intraspecific variation but distinct interspecific differences among those samples of Oesophagostomum spp. examined, which may be used as genetic makers for the identification and differentiation of the Oesophagostomum spp. Phylogenetic analyses using three inference methods, namely Bayesian inference, maximum likelihood, and maximum parsimony based on the combined sequences of pcox1, pnad1, and pnad4 revealed that the O. dentatum and O. quadrispinulatum form monophyletic groups, respectively. These findings demonstrated clearly the usefulness of the three mitochondrial sequences for studying systematics, population genetic structures, and the molecular ecology of Oesophagostomum spp.     

猪鞭虫和斯氏鞭虫线粒体pnad4和pcox1基因的扩增与序列分析

以采自广东不同地区的猪鞭虫(Trichuris suis)和斯氏鞭虫(Trichuris skrjabini)为研究对象,扩增线粒体基因组的pnad4和pcox1基因,将扩增产物纯化后克隆至pMD18-T载体,对阳性菌落进行测序和序列分析。结果表明,来源于不同地区的鞭虫上述两个基因种内相对保守,差异性不大,但种间差异较大。猪鞭虫和斯氏鞭虫的pnad4长度均为468 bp,相互间有162个种间遗传标记位点;猪鞭虫和斯氏鞭虫的pcox1长度分别为600 bp与438 bp,相互间有156个种间遗传标记位点。两种鞭虫之间的pnad4基因差异率为45.1%-45.9%,pcox1基因差异率为44.1%-45.4%。上述线粒体两个基因均可作为区分两种线虫的种间遗传标记。     

Genetic differentiation in eastern European and western Asian populations of the liver fluke, Fasciola hepatica, as revealed by mitochondrial nad1 and cox1 genes

Partial sequences of mitochondrial genes nad1 (316 bp) and cox1 (429 bp) were analyzed to estimate the variability of the liver fluke samples collected in 20 localities in Russia, Belarus, Ukraine, Bulgaria, Armenia, Azerbaijan, Georgia, Turkey, Turkmenistan, and China. The sequences had 4.1% (nad1) and 2.3% (cox1) of variable sites, and 13 and 10 haplotypes were identified among nad1 and cox1 genes, respectively. Spatial analysis of genetic and nucleotide diversity indicated little or no structuring of genetic variation between hosts or regions. The analysis of distribution of both separate and combined (nad1 + cox1) haplotypes revealed the existence of 2 well-defined lineages with 2 main haplotypes and a number of shared divergent haplotypes. Our study showed that the first lineage included the main N1-C1 haplotype, which was found in Australia, China, Georgia, Turkey, Armenia, Azerbaijan, and in all European populations (from Russia, Belarus, Ukraine, Bulgaria). The second lineage was found in all European populations and in populations from Armenia and Azerbaijan. It was suggested that one of the lineages (I) has an Asian origin. The possible source of mtDNA variability and associations between lineage divergence of parasite and its definitive hosts (cattle and sheep) are discussed.     

Sequence analysis of ribosomal and mitochondrial genes of the giant liver fluke Fascioloides magna (Trematoda: Fasciolidae): intraspecific variation and differentiation from Fasciola hepatica

Complete sequences of ribosomal and mitochondrial genes of the giant liver fluke Fascioloides magna are presented. In particular, small subunit (18S) and internal transcribed spacers (ITS1 and ITS2) of the ribosomal gene (rDNA), as well as cytochrome c oxidase subunit I (cox1) and nicotinamide dehydrogenase subunit I (nad1) of the mitochondrial DNA (mtDNA), were analyzed. The 18S and ITS sequences were compared with previously published sequences of the liver fluke Fasciola hepatica. Fixed interspecific genetic differences were determined that allow molecular differentiation of F. magna and F. hepatica using either the PCR-RFLP method or PCR amplification of species-specific DNA regions. Additionally, intraspecific sequence polymorphism of the complete cox1 and nad1 mitochondrial genes in geographically distinct F. magna populations was determined. Based on the sequence divergences, short (< 500 bp) variable regions suitable for broader biogeographical studies of giant liver fluke were designed.     

Characterization of the complete mitochondrial genome sequence of Spirometra erinaceieuropaei (Cestoda: Diphyllobothriidae) from China

Sparganosis, caused by the plerocercoid larvae of members of the genus Spirometra, can cause significant public health problem and considerable economic losses. In the present study, the complete mitochondrial DNA (mtDNA) sequence of Spirometra erinaceieuropaei from China was determined, characterized and compared with that of S. erinaceieuropaei from Japan. The gene arrangement in the mt genome sequences of S. erinaceieuropaei from China and Japan is identical. The identity of the mt genomes was 99.1% between S. erinaceieuropaei from China and Japan, and the complete mtDNA sequence of S. erinaceieuropaei from China is slightly shorter (2 bp) than that from Japan. Phylogenetic analysis of S. erinaceieuropaei with other representative cestodes using two different computational algorithms [Bayesian inference (BI) and maximum likelihood (ML)] based on concatenated amino acid sequences of 12 protein-coding genes, revealed that S. erinaceieuropaei is closely related to Diphyllobothrium spp., supporting classification based on morphological features. The present study determined the complete mtDNA sequences of S. erinaceieuropaei from China that provides novel genetic markers for studying the population genetics and molecular epidemiology of S. erinaceieuropaei in humans and animals.     

The Mitochondrial Genome of Chlorogonium elongatum Inferred from the Complete Sequence

The 22,704-bp circular mitochondrial DNA (mtDNA) of the chlamydomonad alga Chlorogonium elongatum was completely cloned and sequenced. The genome encodes seven proteins of the respiratory electron transport chain, subunit 1 of the cytochrome oxidase complex (cox1), apocytochrome b (cob), five subunits of the NADH dehydrogenase complex (nad1, nad2, nad4, nad5, and nad6), a set of three tRNAs (Q, W, M), and the large (LSU)- and small (SSU)-subunit ribosomal RNAs. Six group-I introns were found, two each in the cox1, cob, and nad5 genes. In each intron an open reading frame (ORF) related to maturases or endonucleases was identified. Both the LSU and the SSU rRNA genes are split into fragments intermingled with each other and with other genes. Although the average A + T content is 62.2%, GC-rich clusters were detected in intergenic regions, in variable domains of the rRNA genes, and in introns and intron-encoded ORFs. A comparison of the genome maps reveals that C. elongatum and Chlamydomonas eugametos mtDNAs are more closely related to one another than either is to Chlamydomonas reinhardtii mtDNA. Received: 3 November 1997 / Accepted: 12 January 1998     

Mitochondrial DNA sequence variation among geographical isolates of Opisthorchis viverrini in Thailand and Lao PDR, and phylogenetic relationships with other trematodes

The present study compared the genetic variation among 14 different geographical isolates of Opisthorchis viverrini sensu lato from Thailand and Lao PDR using sequence data for 2 mitochondrial DNA genes, the subunit 1 of NADH dehydrogenase gene (nad1) and cytochrome c oxidase gene (cox1). Four different nad1 haplotypes were detected among isolates, all of which were identical at the amino acid sequence level. Nucleotide sequence variation among 14 isolates ranged from 0 to 0.3% for nad1. Two different cox1 haplotypes were detected among isolates. These two haplotypes differed at 2 nucleotide positions, one of which resulted in a change in the amino acid sequence. Nucleotide sequence variation among isolates for cox1 ranged from 0 to 0.5%. Comparison of cox1 sequences of O. viverrini to those of other trematodes revealed nucleotide differences of 13-31%. A phylogenetic analysis of the cox1 sequence data revealed strong statistical support for a clade containing O. viverrini and 2 other species of opisthorchid trematodes; O. felineus and Clonorchis sinsensis.     

The mitochondrial genome of the thermal dimorphic fungus Penicillium marneffei is more closely related to those of molds than yeasts

We report the complete sequence of the mitochondrial genome of Penicillium marneffei, the first complete mitochondrial DNA sequence of a thermal dimorphic fungus. This 35 kb mitochondrial genome contains the genes encoding ATP synthase subunits 6, 8, and 9 (atp6, atp8, and atp9), cytochrome oxidase subunits I, II, and III (cox1, cox2, and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxireductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), ribosomal protein of the small ribosomal subunit (rps), 28 tRNAs, and small and large ribosomal RNAs. Analysis of gene contents, gene orders, and gene sequences revealed that the mitochondrial genome of P. marneffei is more closely related to those of molds than yeasts.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.     

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).     

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.     

Maternal inheritance of mitochondria in Eucalyptus globulus

It is important to verify mitochondrial inheritance in plant species in which mitochondrial DNA (mtDNA) will be used as a source of molecular markers. We used a polymerase chain reaction (PCR)/restriction fragment length polymorphism (RFLP) approach to amplify mitochondrial introns from subunits 1, 4, 5, and 7 of NADH dehydrogenase (nad) and cytochrome oxidase subunit II (cox2) in Eucalyptus globulus. PCR fragments were then either sequenced or cut with restriction enzymes to reveal polymorphism. Sequencing cox2 showed that eucalypts lack the intron between exons 1 and 2. One polymorphism was found in intron 2-3 of nad7 following restriction digests with HphI. Fifty-four F1 progeny from seven families with parents distinguishable in their mitochondrial nad7 were screened to show that mitochondria were maternally inherited in E. globulus. These results constitute the first report of mitochondrial inheritance in the family Myrtaceae.     

Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.     

Mitochondrial DNA divergence in populations of the tapeworm Diphyllobothrium nihonkaiense and its phylogenetic relationship with Diphyllobothrium klebanovskii

Diphyllobothrium nihonkaiense [Y. Yamane, H. Kamo, G. Bylund, J.P. Wilkgren. Diphyllobothrium nihonkaiense sp. nov (Cestoda: Diphyllobothriidae)- revised identification of Japanese broad tapeworm. Shimane J Med Sci 1986;10:29-48.] and Diphyllobothrium klebanovskii [I.V. Muratov, P.S. Posokhov. Causative agent of human diphyllobothriasis - Diphyllobothrium klebanovskii sp. n. Parazitologiia. 1988;22:165-170.] are two major species of human diphyllobothriasis in Japan and Far East Russia, respectively, but their taxonomical relationship remains unclear. In this study, we analysed the DNA sequences of 16 clinical isolates of D. nihonkaiense from Japanese people, 3 isolates of D. klebanovskii from a bear in Kamchatka, and 4 clinical isolates of D. klebanovskii from native Udygeyci people in Russia, as well as 4 plerocercoids from Oncorhynchus spp. 18S rDNA and internal transcribed spacer 1 (ITS1) sequences from D. nihonkaiense and D. klebanovskii showed a high level of similarity, indicating synonymy of the two species. Analyses of mitochondrial DNA (mtDNA) sequence polymorphisms in the cox1 and nad3 genes of D. nihonkaiense (D. klebanovskii) revealed two deeply divergent lineages, A and B, with genetic distances (Kimura-2 parameter) of 0.018-0.022. Furthermore, the distinct monophyletic groupings of cox1 haplotypes corresponded to the distinct monophyletic groupings of nad3 haplotypes. The two lineages were neither distinguished by morphological features nor defined by the localities of the samples. These results suggest that the two morphologically cryptic lineages have diverged and coexisted over a long period of time.     

Intraspecific variation of Spirometra erinaceieuropaei and phylogenetic relationship between Spirometra and Diphyllobothrium inferred from mitochondrial CO1 gene sequences

Spirometra erinaceieuropaei is a diphyllobothriid cestode whose adult stage occurs mainly in cat-like carnivores, but occasionally in canids and humans. Although it is generally accepted that the distribution of S. erinaceieuropaei is cosmopolitan, it is controversial as to whether all of S. erinaceieuropaei reported are the same species. This study determined partial sequences of the CO1 gene from several isolates in Asian countries and compared them to sequence data from the GenBank/EMBL/DDBJ nucleotide sequence database. Then intraspecific variation of S. erinaceieuropaei and its phylogenetic relationship with Diphyllobothrium were evaluated. The level of nucleotide variation in the CO1 gene sequences within S. erinaceieuropaei was less than 2.6%. Although it was a little larger than that within each species of Diphyllobothrium (0.1-1.0%), it was much smaller than the interspecific variation within the genus Diphyllobothrium (6.2-14%). These facts indicate that all isolates of S. erinaceieuropaei used in this study, which were collected from Asia, Australia and New Zealand, belong to the same species. Based on CO1 gene sequences, genus Spirometra is clearly separate from the genus Diphyllobothrium. It seems that the genus Spirometra is not a synonym of the genus Diphyllobothrium. The phylogenetic relationship between S. erinaceieuropaei and Sparganum proliferum inferred from the CO1 gene clearly confirm the previous opinion that S. proliferum is a distinct species from S. erinaceieuropaei.     

Exploitation of mitochondrial nad6 as a complementary marker for studying population variability in Lepidoptera

The applicability of mitochondrial nad6 sequences to studies of DNA and population variability in Lepidoptera was tested in four species of economically important moths and one of wild butterflies. The genetic information so obtained was compared to that of cox1 sequences for two species of Lepidoptera. nad6 primers appropriately amplified all the tested DNA targets, the generated data proving to be as informative and suitable in recovering population structures as that of cox1. The proposal is that, to obtain more robust results, this mitochondrial region can be complementarily used with other molecular sequences in studies of low level phylogeny and population genetics in Lepidoptera.