Cryptic Species in the Anopheles (Nyssorhynchus) albitarsis (Diptera: Culicidae) Complex: Incongruence Between Random Amplified Polymorphic DNA-Polymerase Chain Reaction Identification and Analysis of Mitochondrial DNA COI Gene Sequences

Random amplified polymorphic DNA (RAPD) diagnostic bands are one tool used to differentiate cryptic mosquito species in the Anopheles albitarsis Complex. Monophyly of four species (A. albitarsis Lynch-Arribálzaga, A. albitarsis B, A. deaneorum Rosa-Freitas, and A. marajoara Galv?o & Damasceno) currently identified with the RAPD technique was assessed using sequences of the cytochrome oxidase I (COI) mitochondrial DNA (mtDNA) gene. Maximum parsimony, maximum likelihood, and Bayesian analyses support monophyly for A. albitarsis s.s., A. albitarsis B, and A. deaneorum. Anopheles marajoara, as identified by RAPD banding patterns, was either polyphyletic or paraphyletic in all phylogenetic analyses. The phylogenetic pattern and within-species genetic distances observed in A. marajoara suggest the existence of a previously unidentified species (species E) in northern Brazil and Venezuela. Diagnostic RAPD bands were unable to distinguish between A. marajoara and species E, probably because of the low number of correlated bands used to identify species and weaknesses of the RAPD technique, in particular, violations of the untested assumption of homology of comigrating bands. A. marajoara (even without species E) is paraphyletic with respect to A. deaneorum; if A. deaneorum is a separate species from A. marajoara, then A. marajoara may consist of two or more species in Amazonian Brazil. Based on mtDNA COI sequences, there are at least four phylogenetic species within the Albitarsis Complex: A. albitarsis s.s., A. albitarsis B, A. marajoara, and species E; the species status of A. deaneorum is ambiguous.     

Intragenomic rDNA ITS2 variation in the neotropical Anopheles (Nyssorhynchus) albitarsis complex (Diptera: Culicidae)

We cloned and sequenced the rDNA internal transcribed spacer 2 (ITS2) of 4 species belonging to the neotropical Anopheles (Nyssorhynchus) albitarsis complex, that is, A. albitarsis; A. albitarsis B; Anopheles marajoara, a proven malaria vector; and Anopheles deaneorum, a suspected vector. Even though the ITS2 sequences of these species were very similar (< or =1.17% divergence), we found differences suitable for species identification and intragenomic variation of possible consequence in phylogenetic reconstruction. Variation came from 2 microsatellite regions and a number of indels and base substitutions. The existence of partially correlated subsets of clones in A. albitarsis is hypothesized either to be separate rDNA loci or to be semi-independently evolving portions of a single rDNA locus. No differences were found between males and females, suggesting that similar rDNA arrays exist on both the X and Y chromosomes. In addition, highly variant clones, possibly pseudogenes, were found in A. marajoara from Venezuela.     

Identification of Anopheles (Nyssorhynchus) albitarsis complex species (Diptera: Culicidae) using rDNA internal transcribed spacer 2-based polymerase chain reaction primes

Anopheles (Nyssorhynchus) marajoara is a proven primary vector of malaria parasites in Northeast Brazil, and An. deaneorum is a suspected vector in Western Brazil. Both are members of the morphologically similar Albitarsis Complex, which also includes An. albitarsis and an undescribed species, An. albitarsis "B". These four species were recognized and can be identified using random amplified polymorphic DNA (RAPD) markers, but various other methodologies also point to multiple species under the name An. albitarsis. We describe here a technique for identification of these species employing polymerase chain reaction (PCR) primers based on ribosomal DNA internal transcribed spacer 2 (rDNA ITS2) sequence. Since this method is based on known sequence it is simpler than the sometimes problematical RAPD-PCR. Primers were tested on samples previously identified using RAPD markers with complete correlation.     

The taxonomic status of Digramma (Pseudophyllidea: Ligulidae) inferred from DNA sequences

The traditional classification of the ligulid tapeworms into 2 genera, Ligula Bloch, 1782 and Digramma Cholodkovsky, 1914, remains controversial. Molecular data of sequences for the 5' end of the nuclear 28S ribosomal ribonucleic acid (rRNA) gene, the mitochondrial cytochrome c oxidase subunit I (COI) gene, and the nicotinamide adenine dinucleotide dehydrogenase subunit 1 (ND1) gene, as well as the first internal transcribed spacer (ITS1) of the nuclear ribosomal deoxyribonucleic acid (DNA), were used to characterize Digramma and to investigate its relationship with Ligula. Digramma spp. exhibited identical sequences with Ligula intestinalis both in the 28S rRNA and the COI gene and differed from L. intestinalis by 0.7% in the ITS1 region and 7.4% in the ND1 gene, respectively. A high degree of genetic conservation within 28S ribosomal DNA, COI, ITS1, and even ND1 genes, was found in Ligula and Digramma. The low genetic divergence in the 4 genes between Ligula and Digramma indicates that Digramma is probably not an independent genus. Therefore, it is proposed that Ligula and Digramma should be considered as 2 species within the genus Ligula and the tapeworms of Digramma collected from diverse localities in China belong to the same species. The present study also suggests that ITS1 and ND1 sequences can act as useful genetic markers to distinguish Ligula and Digramma.     

Combined molecular and morphological phylogenetic analyses of the New Zealand wolf spider genus Anoteropsis (Araneae: Lycosidae)

Datasets from the mitochondrial gene regions NADH dehydrogenase subunit I (ND1) and cytochrome c oxidase subunit I (COI) of the 20 species in the New Zealand wolf spider (Lycosidae) genus Anoteropsis were generated. Sequence data were phylogenetically analysed using parsimony and maximum likelihood analyses. The phylogenies generated from the ND1 and COI sequence data and a previously generated morphological dataset were significantly congruent (p<0.001). Sequence data were combined with morphological data and phylogenetically analysed using parsimony. The ND1 region sequenced included part of tRNA(Leu(CUN)), which appears to have an unstable amino-acyl arm and no TpsiC arm in lycosids. Analyses supported the existence of five species groups within Anoteropsis and the monophyly of species represented by multiple samples. A radiation of Anoteropsis species within the last five million years is inferred from the ND1 and COI likelihood phylograms, habitat and geological data, which also indicates that Anoteropsis arrived in New Zealand some time after it separated from Gondwana.     

Ixodes philipi (Acari: Ixodidae): phylogenetic status inferred from mitochondrial cytochrome oxidase subunit I gene sequence comparison

Ixodes philipi ticks were collected from the nest burrows of streaked shearwaters, Calonectris luecomelas, on 3 different islands of Japan (Awashima: 38 degrees 45'N, 139 degrees 24'E; Mikurajima: 33 degrees 52'N, 139 degrees 36'E; and Omorijima: 36 degrees 8'N, 133 degrees 10'E). The mitochondrial cytochrome oxidase subunit I (COI) gene sequence was determined for each tick. The COI sequences of 9 other ixodid tick species also were determined, and they were used for taxonomic positioning of I. philipi. A metastriata tick, Amblyomma triguttatum, was used as an outgroup reference for the analysis. Phylogenetic examination indicated that the I. philipi ticks are on the branch with Ixodes turdus and Ixodes acutitarsus weakly, and the bootstrap value of this branching was low. Three different analyses, maximum parsimony, genetic distance, and maximum likelihood, support this conclusion. To further refine this analysis, 2761 base pairs (bp) of sequence, which included the genes for tRNA(Met), NADH dehydrogenase subunit 2 (ND2), tRNA(Trp), tRNA(Cys), tRNA(Tyr), and COI, were determined and compared for 6 I. philipi ticks from the 3 different collection sites. Although a base substitution (T to C in the ND2 gene for an Awashima tick) and 2 transitions (G to A in the COI gene for 1 Omorijima tick) have occurred, the overall sequences were highly conserved. Preserved mitochondrial sequences in the ticks from 3 widely separated locations suggest the possibility of gene flow, which was probably accomplished by migratory seabirds.     

Molecular systematics of the Philippine malaria vector Anopheles flavirostris

Allozyme and molecular sequence data from the malaria vector Anopheles flavirostris (Ludlow) (Diptera: Culicidae) were analysed from 34 sites throughout the Philippines, including the type locality, to test the hypothesis that this taxon is a single panmictic species. A finer-scaled allozyme study, of mainly Luzon samples, revealed no fixed genetic differences in sympatric sites and only low levels of variation. We obtained data from partial sequences for the internal transcribed spacer 2 (ITS2) (483 bp), the third domain (D3) (330 bp) of the 28S ribosomal DNA subunit and cytochrome c oxidase subunit I (COI) of mitochondrial DNA (261 bp). No sequence variation was observed for ITS2, only a one base pair difference was observed between Philippine and Indonesian D3 sequences and An. flavirostris sequences were unique, confirming their diagnostic value for this taxon. Sixteen COI haplotypes were identified, giving 25 parsimony informative sites. Neighbour-Joining, Maximum Parsimony, Maximum Likelihood and Bayesian phylogenetic analysis of COI sequences for An. flavirostris and outgroup taxa revealed strong branch support for the monophyly of An. flavirostris, thus confirming that Philippine populations of this taxon comprise a single separate species within the Minimus Subgroup of the Funestus Group. Variation in the behaviour of An. flavirostris is likely to be intraspecific rather than interspecific in origin.     

Combined molecular and morphological phylogenetic analyses of the New Zealand wolf spider genus Anoteropsis (Araneae: Lycosidae)

Datasets from the mitochondrial gene regions NADH dehydrogenase subunit I (ND1) and cytochrome c oxidase subunit I (COI) of the 20 species in the New Zealand wolf spider (Lycosidae) genus Anoteropsis were generated. Sequence data were phylogenetically analysed using parsimony and maximum likelihood analyses. The phylogenies generated from the ND1 and COI sequence data and a previously generated morphological dataset were significantly congruent (p<0.001). Sequence data were combined with morphological data and phylogenetically analysed using parsimony. The ND1 region sequenced included part of tRNA^Leu(CUN), which appears to have an unstable amino-acyl arm and no TψC arm in lycosids. Analyses supported the existence of five species groups within Anoteropsis and the monophyly of species represented by multiple samples. A radiation of Anoteropsis species within the last five million years is inferred from the ND1 and COI likelihood phylograms, habitat and geological data, which also indicates that Anoteropsis arrived in New Zealand some time after it separated from Gondwana.     

Amazonian malaria vector anopheline relationships interpreted from ITS2 rDNA sequences

Species identification of anopheline mosquitoes (Diptera: Culicidae) can be problematic because many of them belong to complexes of morphologically similar species, often with contrasted ecology, behaviour and vectorial importance. The application of DNA-based diagnostics has proved to be useful for distinguishing between such species. We determined ribosomal DNA sequences of the second internal transcribed spacer (ITS2) from samples of 16 species of Anopheles captured in the Amazon Basin, Brazil. Length of the ITS2 varied from 323 to 410 base pairs, with GC content ranging from 50.7% to 66.5% and sequence identity from 25% to 99% between species. Maximum-likelihood paup analysis separated two distinct groups of species conforming with the recognized subgenera Anopheles (represented by eiseni, mattogrossensis, mediopunctatus and peryassui) and Nyssorhynchus (represented by 12 spp.). For the latter group, the neighbour-joining tree generated from rDNA sequence ITS2 relationships is compatible with the morphological taxonomic key established for these Amazonian species: albitarsis, aquasalis, benarrochi, braziliensis, darlingi, deaneorum, dunhami, evansae, nuneztovari, oswaldoi, rangeli and triannulatus. These ITS2 sequence data proved to be a useful tool for species identification and, potentially, to solve taxonomic problems.     

Molecular comparison of topotypic specimens confirms Anopheles (Nyssorhynchus) dunhami Causey (Diptera: Culicidae) in the Colombian Amazon

The presence of Anopheles (Nyssorhynchus) dunhami Causey in Colombia (Department of Amazonas) is confirmed for the first time through direct comparison of mtDNA cytochrome c oxidase I (COI) barcodes and nuclear rDNA second internal transcribed spacer (ITS2) sequences with topotypic specimens of An. dunhami from Tefé, Brazil. An. dunhami was identified through retrospective correlation of DNA sequences following misidentification as Anopheles nuneztovari s.l. using available morphological keys for Colombian mosquitoes. That An. dunhami occurs in Colombia and also possibly throughout the Amazon Basin, is of importance to vector control programs, as this non-vector species is morphologically similar to known malaria vectors including An. nuneztovari, Anopheles oswaldoi and Anopheles trinkae. Species identification of An. dunhami and differentiation from these closely related species are highly robust using either DNA ITS2 sequences or COI DNA barcode. DNA methods are advocated for future differentiation of these often sympatric taxa in South America.     

Molecular evaluation of nominal species in the Ceratitis fasciventris, C. anonae, C. rosa complex (Diptera: Tephritidae)

Ceratitis fasciventris, C. anonae and C. rosa form a complex of economically important fruit fly pests infesting a variety of crops in African countries. Hitherto only adult males of these species can be distinguished easily by morphological characters. Other stages cannot, and for some taxa the taxonomic interpretation and species boundaries remain unclear. In order to clarify phylogenetic relationships and taxonomic status of these species, sequences of mitochondrial (16S, COI, ND6) and nuclear markers (period, ITS1) were analysed in specimens of the three morphospecies throughout the distribution of the complex. Maximum likelihood trees did not recover monophyletic groups corresponding to the morphospecies. Conversely, ND6 and COI divided West African C. fasciventris specimens in two consistent and bootstrap supported clades, involving specimens from Benin and from Mali/Ivory Coast, while the nuclear gene fragments per and ITS1 recovered a well-supported clade corresponding to C. fasciventris from Kenya/Uganda. Hence, the phylogenetic relationships and taxonomic interpretation of the complex appear more intricate than previously hypothesised. The current molecular data do not allow to identify C. fasciventris, C. anonae and C. rosa as distinct phylogenetic species but rather suggest that the morphospecies C. fasciventris is itself a complex of cryptic taxa.     

Phylogenetic reconstruction and DNA barcoding for closely related pine moth species (Dendrolimus) in China with multiple gene markers

Unlike distinct species, closely related species offer a great challenge for phylogeny reconstruction and species identification with DNA barcoding due to their often overlapping genetic variation. We tested a sibling species group of pine moth pests in China with a standard cytochrome c oxidase subunit I (COI) gene and two alternative internal transcribed spacer (ITS) genes (ITS1 and ITS2). Five different phylogenetic/DNA barcoding analysis methods (Maximum likelihood (ML)/Neighbor-joining (NJ), "best close match" (BCM), Minimum distance (MD), and BP-based method (BP)), representing commonly used methodology (tree-based and non-tree based) in the field, were applied to both single-gene and multiple-gene analyses. Our results demonstrated clear reciprocal species monophyly for three relatively distant related species, Dendrolimus superans, D. houi, D. kikuchii, as recovered by both single and multiple genes while the phylogenetic relationship of three closely related species, D. punctatus, D. tabulaeformis, D. spectabilis, could not be resolved with the traditional tree-building methods. Additionally, we find the standard COI barcode outperforms two nuclear ITS genes, whatever the methods used. On average, the COI barcode achieved a success rate of 94.10-97.40%, while ITS1 and ITS2 obtained a success rate of 64.70-81.60%, indicating ITS genes are less suitable for species identification in this case. We propose the use of an overall success rate of species identification that takes both sequencing success and assignation success into account, since species identification success rates with multiple-gene barcoding system were generally overestimated, especially by tree-based methods, where only successfully sequenced DNA sequences were used to construct a phylogenetic tree. Non-tree based methods, such as MD, BCM, and BP approaches, presented advantages over tree-based methods by reporting the overall success rates with statistical significance. In addition, our results indicate that the most closely related species D. punctatus, D. tabulaeformis, and D. spectabilis, may be still in the process of incomplete lineage sorting, with occasional hybridizations occurring among them.     

Genetic relationships between Oeneis urda and O. mongolica (Nymphalidae: Lepidoptera)

The species status of Oeneis urda (Eversmann) and O. mongolica (Oberthür) has been argued based on morphological characters. Reexamination of their major morphological characters has shown a slight differentiation in the two species. Sequences of three mitochondrial genes (COI, ND6, and ND1) and one nuclear region (internal transcribed spacer 2, ITS2) from two O. urda populations (Yangyang and Mt. Hanla) and one O. mongolica population (Uljin) were performed for phylogenetic and population genetic inferences. Sharing of identical sequences in the ND6 gene and ITS2, minimal sequence divergence in the COI and ND1 genes, and phylogenetically undividable sequence types in all mitochondrial genes and ITS2 suggest genetic continuity between the two species. Nevertheless, significant F_ST estimates (P < 0.05) were found for the COI gene in comparisons between Yangyang (O. urda) and Uljin (O. mongolica), between Yangyang (O. urda) and Mt. Hanla (O. urda), and between Uljin (O. mongolica) and Mt. Hanla (O. urda) populations. These F_ST estimates, along with other gene‐based analyses collectively suggest isolation of the two species at some point in the past, but incomplete separation between the two species on the mainland (Yangyang and Uljin) and biogeographically forced isolation of the O. urda population on Mt. Hanla collectively appear to complicate species status of these two species that were once further clearly separated.     

ITS2 sequence variations among members of Diachasmimorpha longicaudata complex (Hymenoptera: Braconidae) in Thailand

The parasitoid Diachasmimorpha longicaudata complex in Thailand contains at least 3 cryptic species informally designated as species D. longicaudata A, B and BB. DNA sequence data of nuclear ITS2 (second internal transcribed spacer) were used to characterize members of this D. longicaudata complex. The polymerase chain reaction (PCR) amplicon of ITS2 region of D. longicaudata B (≈ 650 bp) clearly differentiated this species from A and BB (amplicon of ≈ 590 bp). Sequence alignment of individual parasitoids revealed that low intraspecies differences ranged from 0.457 to 3.991%, but interspecies differences ranged from 7.566 to 12.989%. Phylogenetic trees constructed using Neighbor-Joining (NJ) and Maximum Parsimony (MP) methods, taking the parasitoid Psyttalia concolor complex as an outgroup, revealed that D. longicaudata A, B and BB formed a monophyletic group, with species A and BB being more closely related than species B. ITS2 characterization of D. longicaudata complex has revealed an interesting divergence of the three cryptic sibling species in Thailand.     

Molecular phylogenetics and mitochondrial genomic evolution in the chamaeleonidae (Reptilia,Squamata)

A phylogenetic hypothesis for the lizard family Chamaeleonidae is generated from 1503 aligned base positions (883 parsimony-informative) of mitochondrial DNA for specimens representing 59 species (57 ingroup and two outgroup). Sequences are reported for a genomic segment encoding eight transfer RNAs, NADH dehydrogenase component 2 (ND2), and portions of NADH dehydrogenase component 1 (ND1) and cytochrome c oxidase subunit 1 (COI). Newly reported genomic rearrangements and duplications support the hypothesis that mitochondrial gene order and content are destabilized by phylogenetic loss of a functional origin for light-strand replication between the genes encoding tRNA(Asn) and tRNA(Cys). A novel gene order characterizes all sampled Brookesia except B. nasus. Brookesia nasus, the apparent sister taxon of a clade formed by all other Brookesia, has the ancestral gene order but contains a large tandem duplication. An apparently noncoding 220 base pair insertion between the genes encoding ND2 and tRNA(Trp) is reported for Bradypodion tavetanum. Phylogenetic analysis identifies nine clades whose ancestral lineages diverged early in chamaeleonid evolutionary history: (1) Brookesia (possibly excluding B. nasus), (2) Chamaeleo subgenus Chamaeleo (excluding C. namaquensis), (3) Chamaeleo subgenus Trioceros, (4) viviparous Bradypodion, (5) oviparous Bradypodion, (6) genus Furcifer (except F. balteatus), and (7-9) three distinct clades of Calumma. Chamaeleo namaquensis, Brookesia nasus, Furcifer balteatus, Rhampholeon brevicaudatus, and R. spectrum represent ancient lineages dating to approximately the same time. Multiple independent losses and a possible secondary gain of horns are inferred for Trioceros. Viviparity has at least two separate origins in chameleons, one in Bradypodion and     

Molecular phylogeny in mytilids supports the wooden steps to deep-sea vents hypothesis

Molecular data were used to study the diversity of mytilids associated with sunken-woods sampled in the Solomon Islands and discuss the 'wooden steps to deep-sea vent' hypothesis proposed by Distel et al. First, COI data used in a barcoding approach confirm the presence of four distinct species. Analyses of the 18S rDNA and COI dataset then confirmed that these sunken-wood mytilids belonged to a monophyletic group including all species from deep-sea reducing environments. Finally, we analyzed the relationships within this monophyletic group that include the Bathymodiolinae using a COI dataset and a combined analysis of mitochondrial COI and ND4 genes and nuclear rDNA 18S and 28S. Our study supported the 'wooden steps to deep-sea vent' hypothesis: one of the sunken-wood species had a basal position within the Bathymodiolionae, and all described vent and seep mussels included in our analyses were derived taxa within Bathymodiolinae.     

Taxonomic status of two South American sympatric squat lobsters, Munida gregaria and Munida subrugosa (Crustacea: Decapoda: Galatheidae), challenged by DNA sequence information

We investigated the taxonomic status of two sympatric morphospecies of squat lobsters from southern South America (Beagle Channel, Strait of Magellan, and Burdwood Bank), Munida gregaria and Munida subrugosa , by DNA sequence analysis of three mitochondrial (mt)DNA gene fragments [416 bp of 16S rDNA(165), 566 bp of cytochrome c oxidase subunit I(COI) and 418 bp of NADH dehydrogenase subunit 1 (ND1)]; and the nuclear rDNA internal transcribed spacer (ITS) 1 (883–952 bp). We obtained a total of 79 sequences from 32 individuals. The 16S sequences of all M. gregaria and M. subrugosa were invariant and identical, whereas COI and ND1 showed 12 and 15 variable sites, respectively. These polymorphisms were shared between morphospecies. Interspecific Tamura–Nei distances for COI and ND1 sequences were 0.0024 and 0.0032, respectively, and were not significantly different from intraspecific distances (Kruskal–Wallis tests: P  = 0.58 and P  = 0.69, for COI and ND1, respectively). Similar to the results obtained from the mtDNA sequences, no relationship was found between the ITS1 maximum parsimony tree topology and the morphologic classification of specimens in M. gregaria and M. subrugosa . We conclude that M. gregaria and M. subrugosa from southern South America may either represent a case of a dimorphic species, or a case of incomplete lineage sorting. The fact that these two morphospecies did not show fixed differences over a total of 1947 bp analysed reinforces the hypothesis of a single dimorphic species.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 421–434.     

Phylogenetic relationships among the Caribbean members of the Cliona viridis complex (Porifera, Demospongiae, Hadromerida) using nuclear and mitochondrial DNA sequences

Species complexes - groups of closely related species in which intraspecific and interspecific variability overlap - have generated considerable interest and study. Frequently, members of a species complex do not have complete reproductive isolation; therefore, the complex may go through extensive gene flow. In the Caribbean Sea, some encrusting and excavating sponges of the genus Cliona (Porifera, Hadromerida, Clionaidae) are grouped within the great "Cliona viridis" complex because of their morphological similarities. This study examined the evolutionary relationships of the Caribbean members of this complex (C. caribbaea, C. tenuis, C. aprica and C. varians) and related taxa based on nuclear (ITS1 and ITS2) and mitochondrial (3' end of ND6) DNA sequences. The intragenomic ITS variation and its secondary structures were evaluated using a mixed approach of Denaturing Gradient Gel Electrophoresis (DGGE), DNA sequencing and secondary structure prediction. Considerable intragenomic variation was found in all the species, with apparently functional ITS1 and ITS2 secondary structures. Despite the subtle but clear morphological differentiation in these excavating sponges, the intragenomic copies of C. caribbaea, C. tenuis and C. aprica had a polyphyletic placement in the ITS1 and ITS2 genealogies and very low divergence. Therefore, it is clear that these species constitute a species complex (herein called Ct-complex). Genetic distances within the Ct-complex revealed that an important part of the interspecific variation overlapped with intraspecific variation, suggesting either incomplete lineage sorting or extensive gene flow. In contrast, C. varians and an unidentified "Pione" species emerged as monophyletic clades, being the closest sister groups to the Ct-complex. Additionally, our results support that C. laticavicola and C. delitrix conform a monophyletic group, but absence of reciprocal monophyly in these species suggests they may be life stages or ecophenotypes of a single species or they have diverged recently. Our work showed that the 3' end of the ND6 mitochondrial gene was highly conserved and not suitable for phylogenetic analysis at the interspecific level.     

Nucleotide sequence of nine protein-coding genes and 22 tRNAs in the mitochondrial DNA of the sea starPisaster ochraceus

Summary We have cloned and sequenced over 9 kb of the mitochondrial genome from the sea starPisaster ochraceus. Within a continuous 8.0-kb fragment are located the genes for NADH dehydrogenase subunits 1, 2, 3, and 4L (ND1, ND2, ND3, and ND4L), cytochrome oxidase subunits I, II, and III (COI, COII, and COIII), and adenosine triphosphatase subunits 6 and 8 (ATPase 6 and ATPase 8). This large fragment also contains a cluster of 13 tRNA genes between ND1 and COI as well as the genes for isoleucine tRNA between ND1 and ND2, arginine tRNA between COI and ND4L, lysine tRNA between COII and ATPase 8, and the serine (UCN) tRNA between COIII and ND3. The genes for the other five tRNAs lie outside this fragment. The gene for phenylalanine tRNA is located between cytochrome b and the 12S ribosomal genes. The genes for tRNA^glu and tRNA^thr are 3 to the 12S ribosomal gene. The tRNAs for histidine and serine (AGN) are adjacent to each other and lie between ND4 and ND5. These data confirm the novel gene order in mitochondrial DNA (mtDNA) of sea stars and delineate additional distinctions between the sea star and other mtDNA molecules.