Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.     

Search for characteristic structural features of mammalian mitochondrial tRNAs

A number of mitochondrial (mt) tRNAs have strong structural deviations from the classical tRNA cloverleaf secondary structure and from the conventional L-shaped tertiary structure. As a consequence, there is a general trend to consider all mitochondrial tRNAs as "bizarre" tRNAs. Here, a large sequence comparison of the 22 tRNA genes within 31 fully sequenced mammalian mt genomes has been performed to define the structural characteristics of this specific group of tRNAs. Vertical alignments define the degree of conservation/variability of primary sequences and secondary structures and search for potential tertiary interactions within each of the 22 families. Further horizontal alignments ascertain that, with the exception of serine-specific tRNAs, mammalian mt tRNAs do fold into cloverleaf structures with mostly classical features. However, deviations exist and concern large variations in size of the D- and T-loops. The predominant absence of the conserved nucleotides G18G19 and T54T55C56, respectively in these loops, suggests that classical tertiary interactions between both domains do not take place. Classification of the tRNA sequences according to their genomic origin (G-rich or G-poor DNA strand) highlight specific features such as richness/poorness in mismatches or G-T pairs in stems and extremely low G-content or C-content in the D- and T-loops. The resulting 22 "typical" mammalian mitochondrial sequences built up a phylogenetic basis for experimental structural and functional investigations. Moreover, they are expected to help in the evaluation of the possible impacts of those point mutations detected in human mitochondrial tRNA genes and correlated with pathologies.     

Mamit-tRNA, a database of mammalian mitochondrial tRNA primary and secondary structures

Mamit-tRNA (http://mamit-tRNA.u-strasbg.fr), a database for mammalian mitochondrial genomes, has been developed for deciphering structural features of mammalian mitochondrial tRNAs and as a helpful tool in the frame of human diseases linked to point mutations in mitochondrial tRNA genes. To accommodate the rapid growing availability of fully sequenced mammalian mitochondrial genomes, Mamit-tRNA has implemented a relational database, and all annotated tRNA genes have been curated and aligned manually. System administrative tools have been integrated to improve efficiency and to allow real-time update (from GenBank Database at NCBI) of available mammalian mitochondrial genomes. More than 3000 tRNA gene sequences from 150 organisms are classified into 22 families according to the amino acid specificity as defined by the anticodon triplets and organized according to phylogeny. Each sequence is displayed linearly with color codes indicating secondary structural domains and can be converted into a printable two-dimensional (2D) cloverleaf structure. Consensus and typical 2D structures can be extracted for any combination of primary sequences within a given tRNA specificity on the basis of phylogenetic relationships or on the basis of structural peculiarities. Mamit-tRNA further displays static individual 2D structures of human mitochondrial tRNA genes with location of polymorphisms and pathology-related point mutations. The site offers also a table allowing for an easy conversion of human mitochondrial genome nucleotide numbering into conventional tRNA numbering. The database is expected to facilitate exploration of structure/function relationships of mitochondrial tRNAs and to assist clinicians in the frame of pathology-related mutation assignments.     

Undetected antisense tRNAs in mitochondrial genomes?

Background  

The hypothesis that both mitochondrial (mt) complementary DNA strands of tRNA genes code for tRNAs (sense-antisense coding) is explored. This could explain why mt tRNA mutations are 6.5 times more frequently pathogenic than in other mt sequences. Antisense tRNA expression is plausible because tRNA punctuation signals mt sense RNA maturation: both sense and antisense tRNAs form secondary structures potentially signalling processing. Sense RNA maturation processes by default 11 antisense tRNAs neighbouring sense genes. If antisense tRNAs are expressed, processed antisense tRNAs should have adapted more for translational activity than unprocessed ones. Four tRNA properties are examined: antisense tRNA 5′ and 3′ end processing by sense RNA maturation and its accuracy, cloverleaf stability and misacylation potential.     

The Complete Mitochondrial Genomes of Six Species of Tetranychus Provide Insights into the Phylogeny and Evolution of Spider Mites

Many spider mites belonging to the genus Tetranychus are of agronomical importance. With limited morphological characters, Tetranychus mites are usually identified by a combination of morphological characteristics and molecular diagnostics. To clarify their molecular evolution and phylogeny, the mitochondrial genomes of the green and red forms of Tetranychus urticae as well as T. kanzawai, T. ludeni, T. malaysiensis, T. phaselus, T. pueraricola were sequenced and compared. The seven mitochondrial genomes are typical circular molecules of about 13,000 bp encoding and they are composed of the complete set of 37 genes that are usually found in metazoans. The order of the mitochondrial (mt) genes is the same as that in the mt genomes of Panonychus citri and P. ulmi, but very different from that in other Acari. The J-strands of the mitochondrial genomes have high (∼84%) A+T contents, negative GC-skews and positive AT-skews. The nucleotide sequence of the cox1 gene, which is commonly used as a taxon barcode and molecular marker, is more highly conserved than the nucleotide sequences of other mitochondrial genes in these seven species. Most tRNA genes in the seven genomes lose the D-arm and/or the T-arm. The functions of these tRNAs need to be evaluated. The mitochondrial genome of T. malaysiensis differs from the other six genomes in having a slightly smaller genome size, a slight difference in codon usage, and a variable loop in place of the T-arm of some tRNAs by a variable loop. A phylogenic analysis shows that T. malaysiensis first split from other Tetranychus species and that the clade of the family Tetranychoidea occupies a basal position in the Trombidiformes. The mt genomes of the green and red forms of T. urticae have limited divergence and short evolutionary distance.     

Evolution and phylogenetic information content of mitochondrial genomic structural features illustrated with acrodont lizards

DNA sequences from 195 squamate reptiles indicate that mitochondrial gene order is the most reliable phylogenetic character establishing monophyly of acrodont lizards and of the snake families Boidae, Colubridae, and Viperidae. Gene order shows no evidence of evolutionary parallelisms or reversals in these taxa. Derived secondary structures of mitochondrial tRNAs also prove to be useful phylogenetic characters showing no reversals. Parallelisms for secondary structures of tRNAs are restricted to deep lineages that are separated by at least 200 million years of independent evolution. Presence of a stem-and-loop structure between the genes encoding tRNA(Asn) and tRNA(Cys), where the replication origin for light-strand synthesis is typically located in vertebrate mitochondrial genomes, is found to undergo at least three and possibly as many as seven evolutionary shifts, most likely parallel losses. This character is therefore a less desirable phylogenetic marker than the other structural changes examined. Sequencing regions that contain multiple genes, including tRNA genes, may be preferable to the common practice of obtaining single-gene fragments for phylogenetic inference because it permits observation of major structural changes in the mitochondrial genome. Such characters may occasionally provide phylogenetic information on relatively short internal branches for which base substitutional changes are expected to be relatively uninformative.     

Two genetic codes, one genome: frameshifted primate mitochondrial genes code for additional proteins in presence of antisense antitermination tRNAs

Genomic amino acid usages coevolve with cloverleaf formation capacities of corresponding primate mitochondrial tRNAs, also for antisense tRNAs, suggesting translational function for sense and antisense tRNAs. Some antisense tRNAs are antitermination tRNAs (anticodons match stops (UAR: UAA, UAG; AGR: AGA, AGG)). Genomes possessing antitermination tRNAs avoid corresponding stops in frames 0 and +1, preventing translational antitermination. In frame +2, AGR stop frequencies and corresponding antisense antitermination tRNAs coevolve positively. This suggests expression of frameshifted overlapping genes, potentially shortening genomes, increasing metabolic efficiency. Blast analyses of hypothetical proteins translated from one and seven +1, respectively, +2 frameshifted human mitochondrial protein coding genes align with eleven GenBank sequences (31% of the mitochondrial coding regions). These putative overlap genes contain few UARs, AGRs align with arginine. Overlap gene numbers increase in presence of, and with time since evolution of antitermination tRNA AGR in 57 primate mitochondrial genomes. Numbers of putative proteins translated from antisense protein coding sequences and detected by blast also coevolve positively with antitermination tRNAs; expression of two of these ‘antisense’ mRNAs increases under low resource availability. Although more direct evidence is still lacking for the existence of proteins translated from overlapping mitochondrial genes and for antisense tRNAs activity, coevolutions between predicted overlap genes and the antitermination tRNAs required to translate them suggest expression of overlapping genes by an overlapping genetic code. Functions of overlapping genes remain unknown, perhaps originating from dual lifestyles of ancestral free living-parasitic mitochondria. Their amino acid composition suggests expression under anaerobic conditions.     

Genes for tRNAAsp,tRNAPro, tRNATyr and two tRNAsSer in wheat mitochondrial DNA

We have begun a systematic search for potential tRNA genes in wheat mtDNA, and present here the sequences of regions of the wheat mitochondrial genome that encode genes for tRNA^Asp (anticodon GUC), tRNA^Pro (UGG), tRNA^Tyr (GUA), and two tRNAs^Ser (UGA and GCU). These genes are all solitary, not immediately adjacent to other tRNA or known protein coding genes. Each of the encoded tRNAs can assume a secondary structure that conforms to the standard cloverleaf model, and that displays none of the structural aberrations peculiar to some of the corresponding mitochondrial tRNAs from other eukaryotes. The wheat mitochondrial tRNA sequences are, on average, substantially more similar to their eubacterial and chloroplast counterparts than to their homologues in fungal and animal mitochondria. However, an analysis of regions 150 nucleotides upstream and 100 nucleotides downstream of the tRNA coding regions has revealed no obvious conserved sequences that resemble the promoter and terminator motifs that regulate the expression of eubacterial and some chloroplast tRNA genes. When restriction digests of wheat mtDNA are probed with ³²P-labelled wheat mitochondrial tRNAs, <20 hybridizing bands are detected, whether enzymes with 4 bp or 6 bp recognition sites are used. This suggests that the wheat mitochondrial genome, despite its large size, may carry a relatively small number of tRNA genes.     

Evolution of mitochondrial gene orders in echinoderms

A comprehensive analysis of the mitochondrial gene orders of all previously published and two novel Antedon mediterranea (Crinoidea) and Ophiura albida (Ophiuroidea) complete echinoderm mitochondrial genomes shows that all major types of rearrangement operations are necessary to explain the evolution of mitochondrial genomes. In addition to protein coding genes we include all tRNA genes as well as the control region in our analysis. Surprisingly, 7 of the 16 genomes published in the GenBank database contain misannotations, mostly unannotated tRNAs and/or mistakes in the orientation of tRNAs, which we have corrected here. Although the gene orders of mt genomes appear very different, only 8 events are necessary to explain the evolutionary history of echinoderms with the exception of the ophiuroids. Only two of these rearrangements are inversions, while we identify three tandem-duplication-random-loss events and three transpositions.     

A correlation between N2-dimethylguanosine presence and alternate tRNA conformers.

Even though the evolutionary conservation of the cloverleaf model is strongly suggestive of powerful constraints on the secondary structure of functional tRNAs, some mitochondrial tRNAs cannot be folded into this form. From the optimal base pairing pattern of these recalcitrant tRNAs, structural correlations between the length of the anticodon stem and the lengths of connector regions between the two helical domains, formed by the coaxial stacking of the anticodon and D-stems and the acceptor and T-stems, have been derived and used to scan the tRNA and tRNA gene database. We show here that some cytosolic tRNA gene sequences that are compatible with the cloverleaf model can also be folded into patterns proposed for the unusual mitochondrial tRNAs. Furthermore, the ability to be folded into these atypical structures correlates in the mature RNA sequences with the presence of dimethylguanosine, whose role may be to prevent the unusual mitochondrial tRNA pattern folding.     

Pathogenic mutations in antisense mitochondrial tRNAs

Pathogenic mutations in mitochondrial tRNAs are 6.5 times more frequent than in other mitochondrial genes. This suggests that tRNA mutations perturb more than one function. A potential additional tRNA gene function is that of templating for antisense tRNAs. Pathogenic mutations weaken cloverleaf secondary structures of sense tRNAs. Analyses here show similar effects for most antisense tRNAs, especially after adjusting for associations between sense and antisense cloverleaf stabilities. These results imply translational activity by antisense tRNAs. For sense tRNAs Ala and Ser UCN, pathogenicity associates as much with sense as with antisense cloverleaf formation. For tRNA Pro, pathogenicity seems associated only with antisense, not sense tRNA cloverleaf formation. Translational activity by antisense tRNAs is expected for the 11 antisense tRNAs processed by regular sense RNA maturation, those recognized by their cognate amino acid’s tRNA synthetase, and those forming relatively stable cloverleaves as compared to their sense counterpart. Most antisense tRNAs probably function routinely in translation and extend the tRNA pool (extension hypothesis); others do not (avoidance hypothesis). The greater the expected translational activity of an antisense tRNA, the more pathogenic mutations weaken its cloverleaf secondary structure. Some evidence for RNA interference, a more classical role for antisense tRNAs, exists only for tRNA Ser UCN. Mutation pathogenicity probably frequently results from a mixture of effects due to sense and antisense tRNA translational activity for many mitochondrial tRNAs. Genomic studies should routinely explore for translational activity by antisense tRNAs.     

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.     

Three LIF-dependent signatures and gene clusters with atypical expression profiles, identified by transcriptome studies in mouse ES cells and early derivatives

Background

The apparent scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). This subclass encompasses over 48,000 species and forms the largest group within the Arachnida. Although mitochondrial genomes are widely utilised for phylogenetic and population genetic studies, only 20 mitochondrial genomes of Acari have been determined, of which only one belongs to the diverse order of the Sarcoptiformes. In this study, we describe the mitochondrial genome of the European house dust mite Dermatophagoides pteronyssinus, the most important member of this largely neglected group.

Results

The mitochondrial genome of D. pteronyssinus is a circular DNA molecule of 14,203 bp. It contains the complete set of 37 genes (13 protein coding genes, 2 rRNA genes and 22 tRNA genes), usually present in metazoan mitochondrial genomes. The mitochondrial gene order differs considerably from that of other Acari mitochondrial genomes. Compared to the mitochondrial genome of Limulus polyphemus, considered as the ancestral arthropod pattern, only 11 of the 38 gene boundaries are conserved. The majority strand has a 72.6% AT-content but a GC-skew of 0.194. This skew is the reverse of that normally observed for typical animal mitochondrial genomes. A microsatellite was detected in a large non-coding region (286 bp), which probably functions as the control region. Almost all tRNA genes lack a T-arm, provoking the formation of canonical cloverleaf tRNA-structures, and both rRNA genes are considerably reduced in size. Finally, the genomic sequence was used to perform a phylogenetic study. Both maximum likelihood and Bayesian inference analysis clustered D. pteronyssinus with Steganacarus magnus, forming a sistergroup of the Trombidiformes.

Conclusion

Although the mitochondrial genome of D. pteronyssinus shares different features with previously characterised Acari mitochondrial genomes, it is unique in many ways. Gene order is extremely rearranged and represents a new pattern within the Acari. Both tRNAs and rRNAs are truncated, corroborating the theory of the functional co-evolution of these molecules. Furthermore, the strong and reversed GC- and AT-skews suggest the inversion of the control region as an evolutionary event. Finally, phylogenetic analysis using concatenated mt gene sequences succeeded in recovering Acari relationships concordant with traditional views of phylogeny of Acari.     

Two ancient bacterial endosymbionts have coevolved with the planthoppers (Insecta: Hemiptera: Fulgoroidea)

Background

Pseudoscorpions are chelicerates and have historically been viewed as being most closely related to solifuges, harvestmen, and scorpions. No mitochondrial genomes of pseudoscorpions have been published, but the mitochondrial genomes of some lineages of Chelicerata possess unusual features, including short rRNA genes and tRNA genes that lack sequence to encode arms of the canonical cloverleaf-shaped tRNA. Additionally, some chelicerates possess an atypical guanine-thymine nucleotide bias on the major coding strand of their mitochondrial genomes.

Results

We sequenced the mitochondrial genomes of two divergent taxa from the chelicerate order Pseudoscorpiones. We find that these genomes possess unusually short tRNA genes that do not encode cloverleaf-shaped tRNA structures. Indeed, in one genome, all 22 tRNA genes lack sequence to encode canonical cloverleaf structures. We also find that the large ribosomal RNA genes are substantially shorter than those of most arthropods. We inferred secondary structures of the LSU rRNAs from both pseudoscorpions, and find that they have lost multiple helices. Based on comparisons with the crystal structure of the bacterial ribosome, two of these helices were likely contact points with tRNA T-arms or D-arms as they pass through the ribosome during protein synthesis. The mitochondrial gene arrangements of both pseudoscorpions differ from the ancestral chelicerate gene arrangement. One genome is rearranged with respect to the location of protein-coding genes, the small rRNA gene, and at least 8 tRNA genes. The other genome contains 6 tRNA genes in novel locations. Most chelicerates with rearranged mitochondrial genes show a genome-wide reversal of the CA nucleotide bias typical for arthropods on their major coding strand, and instead possess a GT bias. Yet despite their extensive rearrangement, these pseudoscorpion mitochondrial genomes possess a CA bias on the major coding strand. Phylogenetic analyses of all 13 mitochondrial protein-coding gene sequences consistently yield trees that place pseudoscorpions as sister to acariform mites.

Conclusion

The well-supported phylogenetic placement of pseudoscorpions as sister to Acariformes differs from some previous analyses based on morphology. However, these two lineages share multiple molecular evolutionary traits, including substantial mitochondrial genome rearrangements, extensive nucleotide substitution, and loss of helices in their inferred tRNA and rRNA structures.     

Stable tRNA-based phylogenies using only 76 nucleotides

tRNAs are among the most ancient, highly conserved sequences on earth, but are often thought to be poor phylogenetic markers because they are short, often subject to horizontal gene transfer, and easily change specificity. Here we use an algorithm now commonly used in microbial ecology, UniFrac, to cluster 175 genomes spanning all three domains of life based on the phylogenetic relationships among their complete tRNA pools. We find that the overall pattern of similarities and differences in the tRNA pools recaptures universal phylogeny to a remarkable extent, and that the resulting tree is similar to the distribution of bootstrapped rRNA trees from the same genomes. In contrast, the trees derived from tRNAs of identical specificity or of individual isoacceptors generally produced trees of lower quality. However, some tRNA isoacceptors were very good predictors of the overall pattern of organismal evolution. These results show that UniFrac can extract meaningful biological patterns from even phylogenies with high level of statistical inaccuracy and horizontal gene transfer, and that, overall, the pattern of tRNA evolution tracks universal phylogeny and provides a background against which we can test hypotheses about the evolution of individual isoacceptors.     

The effect of model choice on phylogenetic inference using mitochondrial sequence data: lessons from the scorpions

Chelicerates are a diverse group of arthropods, with around 65,000 described species occupying a wide range of habitats. Many phylogenies describing the relationships between the various chelicerate orders have been proposed. While some relationships are widely accepted, others remain contentious. To increase the taxonomic sampling of species available for phylogenetic study based on mitochondrial genomes we produced the nearly complete sequence of the mitochondrial genome of the scorpion Mesobuthus gibbosus. Mitochondrial gene order in M. gibbosus largely mirrors that in Limulus polyphemus but tRNA secondary structures are truncated. A recent analysis argued that independent reversal of mitochondrial genome strand-bias in several groups of arthropods, including spiders and scorpions, could compromise phylogenetic reconstruction and proposed an evolutionary model that excludes mutational events caused by strand-bias (Neutral Transitions Excluded, NTE). An arthropod dataset of six mitochondrial genes, when analyzed under NTE, yields strong support for scorpions as sister taxon to the rest of Chelicerata. We investigated the robustness of this result by exploring the effect of adding additional chelicerate genes and taxa and comparing the phylogenies obtained under different models. We find evidence that (1) placement of scorpions arising at the base of the Chelicerata is an artifact of model mis-specification and scorpions are strongly supported as basal arachnids and (2) an expanded chelicerate dataset finds support for several proposed interordinal relationships (ticks plus mites [Acari] and spiders plus whip spiders plus whip scorpions [Araneae+Pedipalpi]). Mitochondrial sequence data are subject to systematic bias that is positively misleading for evolutionary inference and thus extreme methodological care must be taken when using them to infer phylogenies.     

TRANSITION TO A LAND FLORA: A MOLECULAR PHYLOGENETIC PERSPECTIVE

Phylogenetic analyses of morphoultrastructural characters and DNA sequences of chloroplast, mitochondrial, and nuclear genes and investigations of genomic structural features in both organellar genomes have clarified several difficult issues in charophyte and bryophyte phylogenies. Mesostigma seems to represent the first lineage of streptophytes, and Charales are a derived member of charophytes. Liverworts and hornworts both are candidates for the basal-most lineage in land plants according to analyses of morphoultrastructural and sequence data, but mitochondrial intron distribution favors the former. Takakia is likely a basal moss. More studies using both multigene sequence analyses and genomic structural characters will significantly improve our understanding of the origin and early evolution of land plants.     

Comparing two evolutionary mechanisms of modern tRNAs

All modern tRNA gene families have a high similarity in their primary structure, and share the same cloverleaf secondary structure and an inverted L tertiary structure, which provide the clues for the study of their origin and evolution. There are two important mechanisms of the tRNA sequences evolution. One is point mutation, another is complementary duplication method. Both of them are supported by some specific examples. To find out the superior one of the two mechanisms or find out the most suitable mechanism for modern tRNAs evolution, we constructed two types of networks, parallel and antiparallel networks, based on the two mechanisms respectively, and then compared the degree distribution and clustering coefficient of networks constructed by the tRNA sequences of the single anticodon group, single isoaccepting group, and the whole tRNAs group of the two types of networks. The result of the comparison seems consistent with the idea that modern tRNA sequences evolved primarily by the mechanism of complementary method, and point mutation is an important and indispensable auxiliary mechanism during the evolutionary event.     

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