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.     

Sequence variability in three mitochondrial DNA regions of Spirometra erinaceieuropaei spargana of human and animal health significance

Sequence variability in three mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 3 (cox3), NADH dehydrogenase subunits 1 and 4 (nad1 and nad4) in Spirometra erinaceieuropaei spargana from different geographical regions in China was examined. A portion of each of the cox3 (pcox3), nad1 (pnad1) and nad4 genes (pnad4) were amplified separately from individual S. erinaceieuropaei spargana by polymerase chain reaction (PCR). Representative amplicons were subjected to sequencing in order to estimate sequence variability. The sequences of pcox3, pnad1 and pnad4 were 541, 607 and 847?bp in length, respectively. The A+T contents of the sequences were 68.39-68.76% (pcox3), 63.76-64.91% (pnad1) and 67.18-67.77% (pnad4), respectively, while the intra-specific sequence variations within each of the S. erinaceieuropaei spargana were 0-1.5% for pcox3, 0-2.8% for pnad1 and 0-2.7% for pnad4. Phylogenetic analysis using neighbour joining (NJ), maximum likelihood (ML) and maximum parsimony (MP) methods, indicated that all the spargana isolates in Hunan Province represented S. erinaceieuropaei. These findings demonstrated clearly the usefulness of the three mtDNA sequences for population genetics studies of S. erinaceieuropaei spargana of human and animal health significance.     

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.     

Nonneutral Mitochondrial DNA Variation in Humans and Chimpanzees

We sequenced the NADH dehydrogenase subunit 3 (ND3) gene from a sample of 61 humans, five common chimpanzees, and one gorilla to test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution. Within humans and within chimpanzees, the ratio of replacement to silent nucleotide substitutions was higher than observed in comparisons between species, contrary to neutral expectations. To test the generality of this result, we reanalyzed published human RFLP data from the entire mitochondrial genome. Gains of restriction sites relative to a known human mtDNA sequence were used to infer unambiguous nucleotide substitutions. We also compared the complete mtDNA sequences of three humans. Both the RFLP data and the sequence data reveal a higher ratio of replacement to silent nucleotide substitutions within humans than is seen between species. This pattern is observed at most or all human mitochondrial genes and is inconsistent with a strictly neutral model. These data suggest that many mitochondrial protein polymorphisms are slightly deleterious, consistent with studies of human mitochondrial diseases.     

Affinities between Asian non-human Schistosoma species,the S. indicum group,and the African human schistosomes

Schistosoma species have traditionally been arranged in groups based on egg morphology, geographical origins, and the genus or family of snail intermediate host. One of these groups is the 'S. indicum group' comprising species from Asia that use pulmonate snails as intermediate hosts. DNA sequences were obtained from the four members of this group (S. indicum, S. spindale, S. nasale and S. incognitum) to provide information concerning their phylogenetic relationships with other Asian and African species and species groups. The sequences came from the second internal transcribed spacer (ITS2) of the ribosomal gene repeat, part of the 28S ribosomal RNA gene (28S), and part of the mitochondrial cytochrome c oxidase subunit 1 (CO1) gene. Tree analyses using both distance and parsimony methods showed the S. indicum group not to be monophyletic. Schistosoma indicum, S. spindale and S. nasale were clustered among African schistosomes, while S. incognitum was placed as sister to the African species (using ITS2 and 28S nucleotide sequences and CO1 amino acid sequences), or as sister to all other species of Schistosoma (CO1 nucleotide sequences). Based on the present molecular data, a scenario for the evolution of the S. indicum group is discussed.     

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.     

Paramecium mitochondrial genes. I. Small subunit rRNA gene sequence and microevolution

The sequences of the small subunit mitochondrial rRNA genes from two divergent species of Paramecium (primaurelia and tetraurelia) were determined. The gene lies near the center of the linear mitochondrial genome, on the same strand as are all other currently identified genes. The sequences generally resemble their counterparts found in cytoplasmic, procaryotic, and other mitochondrial sources. The rDNA gene boundaries were located by nuclease S1 protection. Small subunit rDNA spans about 1680 nucleotides, including an extraneous 83-base pair sequence very near the 3' end which is unique to Paramecium mitochondria. This "insert" occurs at the apex of the highly variable in length penultimate helix, according to proposed models for small subunit rRNA secondary structure. A discontinuity occurs in isolated rRNA near the start of the insert, resulting in a stable 13 S RNA species and a small segment containing the remaining 3' portion of the gene. The overall rRNA gene sequence was 94% conserved between the two species, and the nucleotide differences consisted of 53% transitions, 37% transversions, and 9% insertions plus deletions. These substitutions were somewhat clustered, and the two most divergent regions coincided with the gene boundaries. The sequence was aligned with Escherichia coli 16 S rRNA for direct comparison of sequence and structure.     

Molecular differentiation of the microgastrine species commonly found in paddy fields from Southeast Asia,with additional data on their phylogeny （Hymenoptera：Braconidae）

Partial DNA sequences of three genes, that is, mitochondrial large ribosomal subunit (16S), nuclear large ribosomal subunit (28S D2) and mitochondrial NADH1 dehydrogenase (NADH1) gene, were sequenced from different microgas trine species(Braconidae: Microgastrinae) collected fresh from paddy fields. The DNA sequences were used to determine the extent of sequence variation among species in order to evaluate the specific status of each species. Cladistic analysis was also used to infer a phylogenetic relationship among these species. The results showed that sequence divergence among species of the same genus Cotesia was much lower than those among different genera, such as Cotesia, Exoryza and Apanteles; the sequence similarity of 16S rDNA and NADH 1 genes between Cotesia sp. and C. chilonis was higher than that between C. sp. and C. ruficrus.Phylogenetic analyses suggested that four species of Cotesia were always grouped in the same clade regardless of using different analysis methods; Cotesia sp. and C. chilonis are more closely related to each other than to C. ruficrus, different from previous morphological results. Additionally, sequence analyses indicated that NADH1 gene has more parsimony informative characters than 28S rDNA D2 and 16S rDNA at the species-level analysis,indicating that NADH1 gene might be a useful marker for species-level analysis.     

Molecular taxonomic position of the elephant schistosome, Bivitellobilharzia nairi, newly discovered in Sri Lanka

Bivitellobilharzia nairi (Mudaliar and Ramanujachar, 1945) Dutt and Srivastava, 1955 was first recorded in India. A number of adult worm specimens of this schistosome species were recovered from a domestic elephant, which died in 1999 in Sri Lanka. This is the first report of this schistosome from Sri Lanka. In the present study, in order to clarify the phylogenetic relationship with other species of schistosomes, sequences from the second internal transcribed spacer (ITS2) of the ribosomal gene repeat, part of the 28S ribosomal RNA gene (28S), and part of the mitochondrial cytochrome c oxidase subunit 1 (CO1) gene from B. nairi were analyzed. Two intraspecific variations were seen within 13 individuals in the ITS2 region. In the CO1 region of the mitochondrial DNA, there were four haplotypes in the nucleotide sequences and two haplotypes in the amino acid sequences. Phylogenetic analysis using the nuclear DNA showed that B. nairi was basal to all of species of the genus Schistosoma. The 28S tree also showed that the mammalian lineage was monophyletic. However, phylogenetic analysis using the mitochondrial DNA showed that B. nairi was nested within the genus Schistosoma. The taxonomical position for this species as well as the contradiction between the results from the nuclear and mitochondrial genes were discussed.     

Mitochondrial sequence evolution in spiders: intraspecific variation in tRNAs lacking the TPsiC Arm

Analyses of mitochondrial DNA sequences from three species of Habronattus jumping spiders (Chelicerata: Arachnida: Araneae) reveal unusual inferred tRNA secondary structures and gene arrangements, providing new information on tRNA evolution within chelicerate arthropods. Sequences from the protein-coding genes NADH dehydrogenase subunit 1 (ND1), cytochrome oxidase subunit I (COI), and subunit II (COII) were obtained, along with tRNA, tRNA, and large-subunit ribosomal RNA (16S) sequences; these revealed several peculiar features. First, inferred secondary structures of tRNA and, likely, tRNA, lack the TPsiC arm and the variable arm and therefore do not form standard cloverleaf structures. In place of these arms is a 5-6-nt T arm-variable loop (TV) replacement loop such as that originally described from nematode mitochondrial tRNAs. Intraspecific variation occurs in the acceptor stem sequences in both tRNAs. Second, while the proposed secondary structure of the 3' end of 16S is similar to that reported for insects, the sequence at the 5' end is extremely divergent, and the entire gene is truncated about 300 nt with respect to Drosophila yakuba. Third, initiation codons appear to consist of ATY (ATT and ATC) and TTG for ND1 and COII, respectively. Finally, Habronattus shares the same ND1-tRNA-16S gene arrangement as insects and crustaceans, thus illustrating variation in a tRNA gene arrangement previously proposed as a character distinguishing chelicerates from insects and crustaceans.     

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.     

Mitochondrial DNA divergence in yoshinobori gobies (Rhinogobius species complex) between the Bonin Islands and the Japan–Ryukyu Archipelago

To investigate the genetic differentiation between the Bonin (Ogasawara) Islands' freshwater goby Rhinogobius sp. Bonin Island (BI) form (Ogasawara-yoshinobori) and the Japan–Ryukyu Archipelago relatives, the mitochondrial DNA (mtDNA) phylogeny of Japanese Rhinogobius species was inferred from partial nucleotide sequences of the mitochondrial NADH dehydrogenase 5 subunit (ND5) gene (945 bp). The resultant tree showed that the Bonin Islands group separated first from the other Japanese lineage, and a test calculation indicated the divergence date to be approximately 3 million years BP. Although it is necessary to use a more reliable estimate to confirm the divergence date, Rhinogobius sp. BI has retained its mtDNA lineage in the islands for millions of years.     

Molecular phylogeny of parasitic Platyhelminthes based on sequences of partial 28S rDNA D1 and mitochondrial cytochrome c oxidase subunit I

The phylogenic relationships existing among 14 parasitic Platyhelminthes in the Republic of Korea were investigated via the use of the partial 28S ribosomal DNA (rDNA) D1 region and the partial mitochondrial cytochrome c oxidase subunit 1 (mCOI) DNA sequences. The nucleotide sequences were analyzed by length, G + C %, nucleotide differences and gaps in order to determine the analyzed phylogenic relationships. The phylogenic patterns of the 28S rDNA D1 and mCOI regions were closely related within the same class and order as analyzed by the PAUP 4.0 program, with the exception of a few species. These findings indicate that the 28S rDNA gene sequence is more highly conserved than are the mCOI gene sequences. The 28S rDNA gene may prove useful in studies of the systematics and population genetic structures of parasitic Platyhelminthes.     

A novel mitochondrial gene order in the crinoid echinoderm Florometra serratissima

The complete nucleotide sequence of the mitochondrial genome of the crinoid Florometra serratissima has been determined. It is a circular DNA molecule, 16,005 bp in length, containing the genes for 13 proteins, small and large ribosomal RNAs, and 22 transfer RNAs (tRNAs). Three regions of unassigned sequence (UAS) greater than 73 bp have been located. The largest, UAS I, is 432 bp long and exhibits sequence similarity to the putative mitochondrial control regions seen in other animals. UAS II (77 bp) and UAS III (73 bp) are located between the 5' ends of coding sequences and may play roles as bidirectional promoters. Analyses of nucleotide composition revealed that the major peptide-encoding strand is high in T and low in C. This bias is reflected in a specific pattern of codon usage. Molecular phylogenetic analyses based on cytochrome c oxidase (COI, COII, and COIII) amino acid and nucleotide sequences did not resolve all the relationships between echinoderm classes. The overall animal mitochondrial gene content has been maintained in the crinoid, but there is extensive rearrangement with respect to both the echinoid and the asteroid mtDNA gene maps. Florometra serratissima has a novel genome organization in a segment containing most of the tRNA genes, large and small rRNA genes, and the NADH dehydrogenase subunit 1 and 2 genes. Potential pathways and mechanisms for gene rearrangements between mitochondrial gene maps of echinoderm classes and vertebrates are discussed as indicators of early deuterostome phylogeny.     

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.     

Complete nucleotide sequence of the mitochondrial genome of a Malagasy poison frog Mantella madagascariensis: evolutionary implications on mitochondrial genomes of higher anuran groups

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of a Malagasy poison frog, Mantella madagascariensis (family Mantellidae), and partial sequences of two Mantella (M. baroni and M. bernhardi) and two additional mantellid species (Boophis madagascariensis and Mantidactylus cf. ulcerosus). The M. madagascariensis genome was shown to be the largest (23kbp) of all vertebrate mtDNAs investigated so far. Furthermore, the following unique features were revealed: (1) the positions of some genes and gene regions were rearranged compared to mitochondrial genomes typical for vertebrates and other anuran groups, (2) two distinct genes and a pseudogene corresponding to transfer RNA gene for methionine (tRNA-Met) were encoded, and (3) two control regions with very high sequence homology were present. These features were shared by the two other Mantella species but not the other mantellid species, indicating dynamic genome reorganization in a common ancestor linage before divergence of the Mantella genus. The reorganization pathway could be explained by a model of gene duplication and deletion. Duplication and deletion events also seem to have been responsible for concerted sequence evolution of the control regions in Mantella mt genomes. It is also suggested that the pseudo tRNA-Met gene sustained for a long time in Mantella mt genomes possibly functions as a punctuation marker for NADH dehydrogenase subunit (ND) 2 mRNA processing. Phylogenetic analyses employing a large sequence data set of mt genes supported the monophyly of Mantellidae and Rhacophoridae and other recent phylogenetic views for ranoid frogs. The resultant phylogenetic relationship also suggested parallel occurrence of two tRNA-Met genes, duplicated control regions, and ND5 gene translocation in independent ranoid lineages.     

Nucleotide sequence and evolution of the 18S ribosomal RNA gene in maize mitochondria.

The nucleotide sequence of the gene coding for the 18S ribosomal RNA of maize mitochondria has been determined and a model for the secondary structure is proposed. Dot matrix analysis has been used to compare the extent and distribution of sequence similarities of the entire maize mitochondrial 18S rRNA sequence with that of 15 other small subunit rRNA sequences. The mitochondrial gene shows great similarity to the eubacterial sequences and to the maize chloroplast, and less similarity to mitochondrial rRNA genes in animals and fungi. We propose that this similarity is due to a slow rate of nucleotide divergence in plant mtDNA compared to the mtDNA of animals. Sequence comparisons indicate that the evolution of the maize mitochondrial 18S, chloroplast 16S and nuclear 17S ribosomal genes have been essentially independent, in spite of evidence for DNA transfer between organelles and the nucleus.     

Marmot phylogeny revisited: molecular evidence for a diphyletic origin of sociality

We established the phylogeny of 11 species of the genus Marmota based on the entire sequence of the mitochondrial cytochrome b ( cyt-b ) gene (1.1 kb) and a partial sequence of the NADH dehydrogenase subunit 4 ( ND4 ) gene (1.2 kb). In three species ( Marmota caligata , Marmota olympus , and Marmota bobac ) full-sized nuclear pseudogenes of the mitochondrial cyt-b were identified. The mitochondrial cyt-b genes and the three pseudogenes form separate clusters in the maximum parsimony dendrogram. This finding suggests that the pseudogenes originated from a single transfer to the nucleus that may have occurred prior to the radiation of the genus Marmota . Notably, compared with their functional mitochondrial equivalents the pseudogenes show a much lower substitution rate. In the dendrograms deduced from the mitochondrial sequences two distinct clusters become apparent: one cluster consists of the North-west American species, the other contains the Eurasian species together with the North American species Marmota monax . The position of M. monax as a member of the Eurasian clade is in accordance with the evolution of chromosome numbers. The results are of special interest with respect to the evolution of social systems in the genus that vary from solitary species ( M. monax ) to highly social species living in family groups (e.g. Marmota marmota ). The molecular phylogeny suggests a diphyletic origin of high sociality in the genus Marmota .