Discrepancy in divergence of the mitochondrial and nuclear genomes ofDrosophila teissieri andDrosophila yakuba

Summary Restriction sites were compared in the mitochondrial DNA (mtDNA) molecules from representatives of two closely related species of fruit flies: nine strains ofDrosophila teissieri and eight strains ofDrosophila yakuba. Nucleotide diversities amongD. teissieri strains and amongD. yakuba strains were 0.07% and 0.03%, respectively, and the nucleotide distance between the species was 0.22%. Also determined was the nucleotide sequence of a 2305-nucleotide pari (ntp) segment of the mtDNA molecule ofD. teissieri that contains the noncoding adenine+thymine (A+T)-rich region (1091 ntp) as well as the genes for the mitochondrial small-subunit rRNA, tRNA^f-met, tRNA^gln, and tRNA^ile, and portions of the ND2 and tRNA^val genes. This sequence differs from the corresponding segment of theD. yakuba mtDNA by base substitutions at 0.1% and 0.8% of the positions in the coding and noncoding regions, respectively. The higher divergence due to base substitutions in the A+T-rich region is accompanied by a greater number of insertions/deletions than in the coding regions. From alignment of theD. teissieri A+T-rich sequence with those ofD. yakuba andDrosophila virilis, it appears that the 40% of this sequence that lies adjacent to the tRNA^ile gene has been highly conserved. Divergence between the entireD. teissieri andD. yakuba mtDNA molecules, estimated from the sequences, was 0.3%; this value is close to the value (0.22%) obtained from the restriction analysis, but 10 times lower than the value estimated from published DNA hybridization results. From consideration of the relationships of mitochondrial nucleotide distance and allozyme genetic distance found among seven species of theDrosophila melanogaster subgroup, the mitochondrial nucleotide distance observed forD. teissieri andD. yakuba is anomalously low in relation to the nuclear genetic distance.     

Evolution of thedec-1 eggshell locus inDrosophila. I. Restriction site mapping and limited sequence comparison in themelanogaster species subgroup

Summary We have analyzed 18 kb of DNA in and upstream of thedefective chorion-1 (dec-1) locus of the eight known species of themelanogaster species subgroup ofDrosophila. The restriction maps ofD. simulans, D. mauritiana, D. sechellia, D. erecta, andD. orena are shown to have basically the restriction map ofD. melanogaster, whereas the maps ofD. teissieri andD. yakuba were more difficult to align. However, the basic amount of DNA and sequence arrangement appear to have been conserved in these species. A small deletion of varying length (65–200 bp) is found in a repeated sequence of the central transcribed region ofD. melanogaster, D. simulans, andD. erecta. Restriction site mapping indicated that thedec-1 gene is highly conserved in themelanogaster species subgroup. However, sequence comparison revealed that the amount of nucleotide and amino acid substitution in the repeated region is much larger than in the 5 translated region. The 5 flanking region showed noticeable restriction site polymorphisms between species. Based on calculations from the restriction maps a dendrogram was derived that supports earlier published phylogenetic relationships within themelanogaster species subgroup except that theerecta-orena pair is placed closer to themelanogaster complex than toD. teissieri andD. yakuba.     

Drosophila virilis histone gene clusters lacking H1 coding segments

Summary Approximately 30–40% ofDrosophila virilis DNA complementary to clonedDrosophila histone genes is reduced to 3.4-kilobase-pair (kbp) segments by Bgl I or Bgl II digestion. The core histone genes of a 3.4-kbp Bgl II segment cloned in the plasmid pDv3/3.4 have the same order as theD. melanogaster core histone genes in the plasmid cDm500: . Nonetheless, pDv3/3.4 and cDm500 have different histone gene configurations: In pDv3/3.4, the region between the H2B and H3 genes contains 0.35 kbp and cannot encode histone H1; in cDm500, the region contains 2.0 kbp and encodes histone H1. The lack of an H1 gene between the H2B and H3 genes in 30–40% ofD. virilis histone gene clusters suggests that changes in histone gene arrays have occurred during the evolution ofDrosophila. The ancestors of modernDrosophila may have possessed multiple varieties of histone gene clusters, which were subsequently lost differentially in thevirilis andmelanogaster lineages. Alternatively, they may have possessed a single variety, which was rearranged during evolution. The H1 genes ofD. virilis andD. melanogaster did not cross-hybridize in vitro under conditions that maintain stable duplexes between DNAs that are 75% homologous. Consequently,D. virilis H1 genes could not be visualized by hybridization to an H1-specific probe and thus remain unidentified. Our observations suggest that the coding segments in the H1 genes ofD. virilis andD. melanogaster are >25% divergent. Our estimate of sequence divergence in the H1 genes ofD. virilis andD. melanogaster seems high until one considers that the coding sequences of cloned H1 genes from the closely related speciesD. melanogaster andD. simulans are 5% divergent.     

Platyhelminth mitochondrial DNA: Evidence for early evolutionary origin of a tRNAserAGN that contains a dihydrouridine arm replacement loop,and of serine-specifying AGA and AGG codons

Summary The nucleotide sequence of a segment of the mitochondrial DNA (mtDNA) molecule of the liver flukeFasciola hepatica (phylum Platyhelminthes, class Trematoda) has been determined, within which have been identified the genes for tRNA^ala, tRNA^asp, respiratory chain NADH dehydrogenase subunit I (ND1), tRNA^asn, tRNA^pro, tRNA^ile, tRNA^lys, ND3, tRNA^serAGN, tRNA^trp, and cytochromec oxidase subunit I (COI). The 11 genes are arranged in the order given and are all transcribed from the same strand of the molecule. The overall order of theF. hepatica mitochondrial genes differs from what is found in other metazoan mtDNAs. All of the sequenced tRNA genes except the one for tRNA^serAGN can be folded into a secondary structure with four arms resembling most other metazoan mitochondrial tRNAs, rather than the tRNAs that contain a TψC arm replacement loop, found in nematode mtDNAs. TheF. hepatica mitochondrial tRNA^serAGN gene contains a dihydrouridine arm replacement loop, as is the case in all other metazoan mtDNAs examined to date. AGA and AGG are found in theF. hepatica mitochondrial protein genes and both codons appear to specify serine. These findings concerningF. hepatica mtDNA indicate that both a dihydrouridine arm replacement loop-containing tRNA^serAGN gene and the use of AGA and AGG codons to specify serine must first have occurred very early in, or before, the evolution of metazoa.     

The sequence,organization, and evolution of the Locusta migratoria mitochondrial genome

The sequencing of the cloned Locusta migratoria mitochondrial genome has been completed. The sequence is 15,722 by in length and contains 75.3% A+T, the lowest value in any of the five insect mitochondrial sequences so far determined. The protein coding genes have a similar A+T content (74.1%) but are distinguished by a high cytosine content at the third codon position. The gene content and organization are the same as in Drosophila yakuba except for a rearrangement of the two tRNA genes tRNA^lys and tRNA^asp. The A+T-rich region has a lower A+T nucleotide content than in other insects, and this is largely due to the presence of two G+C-rich 155-bp repetitive sequences at the 5 end of this section and the beginning of the adjacent small rRNA gene. The sizes of the large and small rRNA genes are 1,314 and 827 bp, respectively, and both sequences can be folded to form secondary structures similar to those previously predicted for Drosophila. The tRNA genes have also been modeled and these show a strong resemblance to the dipteran tRNAs, all anticodons apparently being conserved between the two species. A comparison of the protein coding nucleotide sequences of the locust DNA with the homologous sequences of five other arthropods (Drosophila yakuba, Anopheles quadrimaculatus, Anopheles gambiae, Apis mellifera, and Artemia franciscana) was performed. The amino acid composition of the encoded proteins in Locusta is similar to that of Drosophila, with a Dayhoff distance twice that of the distance between the fruit fly and the mosquitoes. A phylogenetic analysis revealed the locust genes to be more similar to those of the Dipterans than to those of the honeybee at both the nucleotide and amino acid levels. A comparative analysis of tRNA orders, using crustacean mtDNAs as outgroups, supported this. This high level of divergence in the Apis genome has been noted elsewhere and is possibly an effect of directional mutation pressure having resulted in an accelerated pattern of sequence evolution. If the general assumption that the Holometabola are monophyletic holds, then these results emphasize the difficulties of reconstructing phylogenies that include lineages with variable substitution rates and base composition biases. The need to exercise caution in using information about tRNA gene orders in phylogenetic analysis is also illustrated. However, if the honeybee sequence is excluded, the correspondence between the other five arthropod sequences supports the findings of previous studies which have endorsed the use of mtDNA sequences for studies of phylogeny at deep levels of taxonomy when mutation rates are equivalent. Correspondence to: P.K. Flook     

A vital function for mitochondrial DNA in the petite-negative yeastKluyveromyces lactis

Petite-negative yeasts do not form viable respiratory-deficient mutants on treatment with DNA-targeting drugs that readily eliminate the mitochondial DNA (mtDNA) from petite-positive yeasts. However, in the petite-negative yeastKluyveromyces lactis, specific mutations in the nuclear genesMGI2 andMGI5 encoding the- and-subunits of the mitochondrial F₁-ATPase, allow mtDNA to be lost. In this study we show that wild-typeK. lactis does not survive in the absence of its mitochondrial genome and that the function ofmgi mutations is to suppress lethality caused by loss of mtDNA. Firstly, we find that loss of a multicopy plasmid bearing amgi allele readily occurs from a wild-type strain with functional mtDNA but is not tolerated in the absence of mtDNA. Secondly, we cloned theK. lactis homologue of theSaccharomyces cerevisiae mitochondrial genome maintenance geneMGM101, and disrupted one of the two copies in a diploid. Following sporulation, we find that segregants containing the disrupted gene form minicolonies containing 6-8000 inviable cells. By contrast, disruption ofMGM101 is not lethal in a haploidmgi strain with a specific mutation in a subunit of the mitochondrial F₁-ATPase. These observations suggest that mtDNA inK. lactis encodes a vital function which may reside in one of the three mitochondrially encoded subunits of F₀.     

Mitochondrial DNA of the sea anemone,Metridium senile (Cnidaria): Prokaryote-like genes for tRNAf-Met and small-subunit ribosomal RNA,and standard genetic code specificities for AGR and ATA codons

The nucleotide sequence of a segment of the mitochondrial DNA (mtDNA) molecule of the sea anemone Metridium senile (phylum Cnidaria, class Anthozoa, order Actiniaria) has been determined, within which have been identified the genes for respiratory chain NADH dehydrogenase subunit 2 (ND2), the small-subunit rRNA (s-rRNA), cytochrome c oxidase subunit II(COII), ND4, ND6, cytochrome b (Cyt b), tRNA^f-Met, and the large-subunit rRNA (1-rRNA). The eight genes are arranged in the order given and are all transcribed from the same strand of the molecule. The overall order of the M. senile mt-genes differs from that of other metazoan mtDNAs. In M. senile mt-protein genes, AGA and AGG codons appear to have the standard genetic code specification of arginine, rather than serine as found for other invertebrate mt-genetic codes. Also, ATA has the standard genetic code specification of isoleucine. TGA occurs in three M. senile mt-protein genes and may specify tryptophan as in other metazoan, protozoan, and some fungal mt-genetic codes. The M. senile mt-rRNA^f-Met gene has primary and secondary structure features closely resembling those of the Escherichia coli initiator tRNA, including standard dihydrouridine and TC loop sequences and a mismatch pair at the top of the aminoacyl stem. Determinations of the 5 and 3 end nucleotides of the M. senile mt-srRNAs indicated that these molecules have a homogenous size of 1,081 ntp, larger than any other known metazoan mt-s-rRNAs. Consistent with its larger size, the M. senile mt-s-rRNA can be folded into a secondary structure that more closely resembles that of the E. coli 16S rRNA than can any other metazoan mt-s-rRNA. These findings concerning M. senile mtDNA indicate that most of the unusual features regarding metazoan mt-genetic codes, rRNAs, and probably tRNAs developed after divergence of the Cnidarian line from the ancestral line common to other metazoa.Correspondence to: D.R. Wolstenholme     

Sequence evolution of theGpdh gene in theDrosophila virilis species group

The nucleotide sequence of theGpdh gene from six taxa,D. virilis, D. lummei, D. novamexicana, D. a. americana, D. a. texana andD. ezoana, belonging to thevirilis species group was determined to examine details of evolutionary change in the structure of theGpdh gene. TheGpdh gene is comprised of one 5 non-translated region, eight exons, seven introns and three 3 non-translated regions. Exon/intron organization was identical in all the species examined, but different from that of mammals. Interspecific nucleotide divergence in the entireGpdh gene followed the common pattern: it was low in the exon, high in the intron and intermediate in the non-translated regions. The degree of nucleotide divergence differed within these regions, suggesting that selection exerts constraints differentially on nucleotide change of theGpdh gene. A phylogenetic tree of thevirilis phylad constructed from nucleotide variation of total sequence was consistent with those obtained from other data.Nucleotide sequences for theGpdh gene ofD. lummei, D. novamexicana, D. a. americana, D. a. texana andD. ezoana have been submitted to GenBank with accession numbers D50087, D50088, D50089, D50090 and D50091.     

The mitochondrial ribosomal RNA genes of the nematodes Caenorhabditis elegans and Ascaris suum: Consensus secondary-structure models and conserved nucleotide sets for phylogenetic analysis

The small- and large-subunit mitochondrial ribosomal RNA genes (mt-s-rRNA and mt-1-rRNA) of the nematode worms Caenorhabditis elegans and Ascaris suum encode the smallest rRNAs so far reported for metazoa. These size reductions correlate with the previously described, smaller, structurally anomalous mt-tRNAs of C. elegans and A. suum. Using primer extension analysis, the 5 end nucleotides of the mt-s-rRNA and mt-1-rRNA genes were determined to be adjacent to the 3 end nucleotides of the tRNA^Glu and tRNA^His genes, respectively. Detailed, consensus secondary-structure models were constructed for the mt-s-rRNA genes and the 3 64% of mt-1-rRNA genes of the two nematodes. The mt-s-rRNA secondary-structure model bears a remarkable resemblance to the previously defined universal core structure of E. coli 16S rRNA: most of the nucleotides that have been classified as variable or semiconserved in the E. coli model appear to have been eliminated from the C. elegans and A. suum sequences. Also, the secondary structure model constructed for the 3 64% of the mt-1-rRNA is similar to the corresponding portion of the previously defined E. coli 23S rRNA core secondary structure. The proposed C. elegans/A. suum mt-s-rRNA and mt-1-rRNA models include all of the secondary-structure element-forming sequences that in E. coli rRNAs contain nucleotides important for A-site and P-site (but not E-site) interactions with tRNAs. Sets of apparently homologous sequences within the mt-s-rRNA and mt-1-rRNA core structures, derived by alignment of the C. elegans and A. suum mt-rRNAs to the corresponding mt-rRNAs of other eukaryotes, and E. coli rRNAs were used in maximum-likelihood analyses. The patterns of divergence of metazoan phyla obtained show considerable agreement with the most prevalent metazoan divergence patterns derived from more classical, morphological, and developmental data.Correspondence to: D.R. Wolstenholme     

Functional conservation of a glucose-repressible amylase gene promoter from Drosophila virilis in Drosophila melanogaster

Summary Previous studies have demonstrated that the expression of the -amylase gene is repressed by dietary glucose in Drosophila melanogaster. Here, we show that the -amylase gene of a distantly related species, D. virilis, is also subject to glucose repression. Moreover, the cloned amylase gene of D. virilis is shown to be glucose repressible when it is transiently expressed in D. melanogaster larvae. This cross-species, functional conservation is mediated by a 330-bp promoter region of the D. virilis amylase gene. These results indicate that the promoter elements required for glucose repression are conserved between distantly related Drosophila species. A sequence comparison between the amylase genes of D. virilis and D. melanogaster shows that the promoter sequences diverge to a much greater degree than the coding sequences. The amylase promoters of the two species do, however, share small clusters of sequence similarity, suggesting that these conserved cis-acting elements are sufficient to control the glucose-regulated expression of the amylase gene in the genus Drosophila.Offprint requests to: D.A. Hickey     

Five transfer RNA genes lacking CCA termini are clustered in the chromosome of Streptomyces rimosus

Summary The nucleotide sequence of a 1105 by Streptomyces rimosus DNA fragment containing five transfer RNA genes was determined. Two tRNA^Gln (CUG) genes, differing by 1 by in the aminoacyl stem, and three identical tRNA^Glu (CUC) genes were identified. The five tRNA genes, arranged in the order: Gln1-Glul-Glu2-Gln2-Glu3, were separated by short, nonhomologous intergenic regions. Surprisingly, none of these tRNA genes encoded the CCA 3 terminus of mature tRNAs. All five encoded tRNAs for the translation of GC rich codons, which are preferentially used in Streptomyces genes (CAG and GAG, respectively). We recently reported nucleotide sequences of two initiator tRNA genes from S. rimosus, which also do not encode the CCA end of mature tRNAs. It is therefore very likely that S. rimosus represents an example of those eubacteria in which the majority of tRNA genes do not encode the 3 terminal CCA end of mature tRNAs. Evolutionary implications of this finding remain to be elucidated.     

Evolution of the noncoding regions inDrosophila mitochondrial DNA: Two intergenic regions

The sequences of the mitochondrial DNA (mtDNA) segment containing the two intergenic regions were determined for six species belonging to theDrosophila immigrans species group and compared to the corresponding segments ofDrosophila species which had been studied previously. We found remarkable differences in the evolutionary rates of the two intergenic regions. The Intergenic I region, which lies between thetRNA ^gln and thetRNA ^ile genes, was found to be highly conserved in terms of both size (30 ntp) and nucleotide sequence among the species studied. In contrast, the sequences of the Intergenic II region, which lies between thetRNA ^f-met and thetRNA ^ile genes, showed considerable variation. The size of the Intergenic II region ranged from 0 to 88 ntp, and accurate alignment was possible only among sequences from geographical strains or very closely related species in thenasuta species subgroup. The observed differences in conservation of the two mtDNA intergenic regions are discussed in light of functional constraints on mtDNA sequences.     

Mitochondrial DNA of the blowflyPhormia regina: restriction analysis and gene localization

A study of an invertebrate mitochondrial genome, that of the blowflyPhormia regina, has been initiated to compare its structural and functional relatedness to other metazoan mitochondrial genomes. A restriction map of mitochondrial DNA (mtDNA) isolated from sucrose gradient-purified mitochondria has been established using a combination of single and double restriction endonuclease digestions and hybridizations with isolated mtDNA fragments, revealing a genome size of 17.5 kilobases (kb). A number of mitochondrial genes including those encoding the 12 S and 16 S ribosomal RNA, the cytochromec oxidase I subunit (COI) and an unidentified open reading frame (URF2) have been located on thePhormia mtDNA by Southern blot analysis using as probes both isolated mtDNA fragments and oligonucleotides derived from the sequences of previously characterized genes from rat andDrosophila yakuba mtDNAs. These data indicate that for those regions examined, the mitochondrial genome organization of blowfly mtDNA is the same as that ofDrosophila yakuba, the order being COI-URF2-12 S-16 S. These data also report the presence of an A + T-rich region, located as a 2.5-kb region between the URF2 and the 12 S rRNA genes, and its amplification by the polymerase chain reaction is described.     

Fertile somatic hybrids between transgenicNicotiana tabacum and transgenicN. debneyi selected by dual-antibiotic resistance

Summary A simple, yet effective selection system was used to produce fertile somatic hybrids betweenNicotiana tabacum andN. debneyi. This approach utilized transgenic antibiotic-resistantN. tabacum andN. Debneyi as donor plants for mesophyll protoplast fusions. Thirteen somatic hybrid plants were regenerated from calli capable of growth on medium containing both antibiotics. The majority of the hybrids displayed a range of leaf and floral morphologies and growth habits that were intermediate to those of the parental species, and had chromosome numbers varying from amphidiploid (2n = 96) to hypoaneuploid (2n = 60). Isoenzyme and RFLP analysis demonstrated the presence and expression of nuclear genes from both parents in all of the hybrids. Most plants are fully fertile. Thus, these plants differ from the malesterile tobacco cybrids and alloplasmic lines produced by transferring theN. debneyi cytoplasm to tobacco. A nonrandom pattern of cytoplasmic segregation in the fusion products occurred with a bias towards the presence ofN. debneyi cp and mtDNA. Evidence for the presence of rearranged or recombinant cp and mtDNA in some of the hybrids was obtained. The somatic hybrids were successfully backcrossed to theN. tabacum parent and are now being tested for immunity to black root rot, a trait specific toN. debneyi, but not existent in theN. tabacum parental line.     

Sequence analysis of three tRNAPhe nuclear genes and a mutated gene,and one gene for tRNAAla from Arabidopsis thaliana

Three genes and one mutant gene for tRNA^Phe (GAA) and one gene for tRNA^Ala (UGC) were isolated from a whole-cell DNA library of Arabidopsis thaliana. All three tRNA^Phe genes are identical in their nucleotide sequence, but differ in their 5 and 3 flanking regions. The mutant tRNA^Phe (GAA) gene differs from the other three genes by one nucleotide change from highly conserved G to C at the 57th nucleotide position. The primary structure of the first tRNA^Ala gene was also determined in this experiment.     

Genetics of esterases in Drosophila. V. Characteristics of the “pupal” esterase in D. virilis

From analysis of the properties of the pupal esterase (p-esterase) in Drosophila virilis, it is concluded that it is heat stable, its electrophoretic detection depends on culture density, its expression is stage specific, and it is not a variant of esterase 2. It was also demonstrated that p-esterase, like esterase 6, is activated by injections of the juvenile hormone into larvae. Heat treatment of heat-resistant D. virilis stocks led to decreased activities of the juvenile hormone dependent esterases but did not affect those of the heat-sensitive stocks. It is suggested that heat resistance in D. virilis is related to some functional features of the system of modifier genes controlling the phenotypic expression of esterases.     

The potato mitochondrial initiator methionine tRNA gene and its flanking regions: an illustration of the diversity of mitochondrial genome rearrangements among plant species

The initiator methionine transfer RNA (tRNAf ^Met) gene was identified on a 347 bpEco RI-Hind III DNA fragment of the potato mitochondrial (mt) genome. The sequence of this gene shows 1 to 7 nucleotide differences with the other plant mt tRNAsf ^Met or tRNAf ^Met genes studied so far. Whereas the tRNAf ^Met gene is present as a single copy in the potato mt genome, a tRNA pseudogene corresponding to 60% of a complete tRNA (from the 5 end to the variable region) and located at 105 nucleotides upstream of the tRNAf ^Met gene on the opposite strand was shown to be repeated at least three times. Furthermore, the physical environment of the tRNAf ^Met gene in the mt genome is very different among plants, which suggests that the tRNAf ^Met gene region has often been implicated in recombination events of plant mt genomes leading to important rearrangements in gene order.     

Indirect evidence of alteration in the expression of the rDNA genes in interspecific hybrids betweenDrosophila melanogaster andDrosophila simulans

Crosses betweenDrosophila melanogaster females andD. simulans males produce viable hybrid females, while males are lethal. These males are rescued if they carry theD. simulans Lhr gene. This paper reports that females of the wild-typeD. melanogaster population Staket do not produce viable hybrid males when crossed withD. simulans Lhr males, a phenomenon which we designate as the Staket phenotype. The agent responsible for this phenomenon was found to be the StaketX chromosome (X ^mel ,Stk). Analysis of the Staket phenotype showed that it is suppressed by extra copies ofD. melanogaster rDNA genes and that theX ^mel ,Stk chromosome manifests a weak bobbed phenotype inD. melanogaster X ^mel ,Stk/0 males. The numbers of functional rDNA genes inX ^mel ,Stk andX ^mel ,y w (control) chromosomes were found not to differ significantly. Thus a reduction in rDNA gene number cannot account for the weak bobbedX ^mel ,Stk phenotype let alone the Staket phenotype. The rRNA precursor molecules transcribed from theX ^mel ,Stk rDNA genes seem to be correctly processed in both intraspecific (melanogaster) and interspecific (melanogaster-simulans) conditions. It is therefore suggested that theX ^mel ,Stk rDNA genes are inefficiently transcribed in themelanogaster-simulans hybrids.     

Structural evolution of the Drosophila 5S ribosomal genes

We compare the 5S gene structure from nine Drosophila species. New sequence data (5S genes of D. melanogaster, D. mauritiana, D. sechellia, D. yakuba, D. erecta, D. orena, and D. takahashii) and already-published data (5S genes of D. melanogaster, D. simulans, and D. teissieri) are used in these comparisons. We show that four regions within the Drosophila 5S genes display distinct rates of evolution: the coding region (120 bp), the 5-flanking region (54–55 bp), the 3-flanking region (21–22 bp), and the internal spacer (149–206 bp). Intra- and interspecific heterogeneity is due mainly to insertions and deletions of 6–17-bp oligomers. These small rearrangements could be generated by fork slippages during replication and could produce rapid sequence divergence in a limited number of steps. Correspondence to: M. Wegnez     

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.