The Mitochondrial Genome of the Honeybee Apis Mellifera: Complete Sequence and Genome Organization

The complete sequence of honeybee (Apis mellifera) mitochondrial DNA is reported being 16,343 bp long in the strain sequenced. Relative to their positions in the Drosophila map, 11 of the tRNA genes are in altered positions, but the other genes and regions are in the same relative positions. Comparisons of the predicted protein sequences indicate that the honeybee mitochondrial genetic code is the same as that for Drosophila; but the anticodons of two tRNAs differ between these two insects. The base composition shows extreme bias, being 84.9% AT (cf. 78.6% in Drosophila yakuba). In protein-encoding genes, the AT bias is strongest at the third codon positions (which in some cases lack guanines altogether), and least in second codon positions. Multiple stepwise regression analysis of the predicted products of the protein-encoding genes shows a significant association between the numbers of occurrences of amino acids and %T in codon family, but not with the number of codons per codon family or other parameters associated with codon family base composition. Differences in amino acid abundances are apparent between the predicted Apis and Drosophila proteins, with a relative abundance in the Apis proteins of lysine and a relative deficiency of alanine. Drosophila alanine residues are as often replaced by serine as conserved in Apis. The differences in abundances between Drosophila and Apis are associated with %AT in the codon families, and the degree of divergence in amino acid composition between proteins correlates with the divergence in %AT at the second codon positions. Overall, transversions are about twice as abundant as transitions when comparing Drosophila and Apis protein-encoding genes, but this ratio varies between codon positions. Marked excesses of transitions over chance expectation are seen for the third positions of protein-coding genes and for the gene for the small subunit of ribosomal RNA. For the third codon positions the excess of transitions is adequately explained as due to the restriction of observable substitutions to transitions for conserved amino acids with two-codon families; the excess of transitions over expectation for the small ribosomal subunit suggests that the conservation of nucleotide size is favored by selection.     

Nucleotide sequence and evolution of the rat mitochondrial cytochrome b gene containing the ochre termination codon

The nucleotide sequences of the genes for cytochrome b and three potential transfer RNAs (tRNAPro, tRNAThr and tRNAGlu) in cloned rat mitochondrial DNA were determined. The derived amino acid sequence of the cytochrome b protein from the light strand indicated that the C-terminal amino acid is asparagine and the ochre termination codon is encoded in the DNA, in contrast to the the lack of termination codon in the reading frame of human [Anderson et al., Nature 290 (1981) 457] or mouse [Bibb et al., Cell 26 (1981) 167] mitochondrial DNA. The first ATG codon of the cytochrome b gene was spaced five nucleotides from the 5'-end of the tRNAGlu gene on the heavy strand. There was a single nucleotide spacing between the termination codon of the cytochrome b gene and the 5' end of the tRNAThr gene in the light strand. There was also a single nucleotide spacing between the 3'-end of the tRNAThr gene and the 3'-end of the tRNAPro gene on the heavy strand. The amino acid and nucleotide sequences of the cytochrome b genes of mammals and yeast [Nobrega and Tzagoloff, J. Biol. Chem. 255 (1980) 9828] were compared to reveal structural differences in two very different species. At the same time, amino acid substitutions in particular regions of the mammalian gene corresponding to the exon-intron boundaries in the yeast gene were noted. These genetic features are discussed in relation to the extreme compression of genetic information in the mammalian mitochondrial genome as related to the evolution of the gene organization and its sequence.     

Evolution of the mitochondrial cytochrome oxidase II gene among 10 orders of insects.

We examine the complete nucleotide sequences of the mitochondrial cytochrome oxidase II gene of 13 species of insects, representing 10 orders. The genes range from 673 to 690 bp in length, encoding 226 to 229 amino acids. Several insertion or deletion events, each involving one or two codons, can be observed. The 3' end of the gene is extremely variable in both length and sequence, making alignment of the ends unreliable. Using the first 639 nucleotide positions, for which unambiguous alignments could be obtained, we examine the neighbor-joining trees based on nucleotide divergences and based on conserved subsets of that data, including transversion and amino acid and second codon position divergences. Each of these subsets produces different trees, none of which can be easily reconciled with trees constructed using morphology and the fossil record. Bootstrap analysis using second codon positions strongly supports affinities between the order Blatteria (cockroaches) and the order Isoptera (termites) and between a wasp and the published honeybee sequence (Order Hymenoptera). The divergence of insect orders is very ancient and may have occurred too rapidly for easy resolution using mitochondrial protein sequences. Unambiguous resolution of insect orders will probably require analysis of many additional taxa, using the COII gene and other conserved sequences.     

The CO-I and CO-II region of honeybee mitochondrial DNA: evidence for variation in insect mitochondrial evolutionary rates

The sequence of a region of honeybee (Apis mellifera ligustica) mitochondrial DNA, which contains the genes for cytochrome c oxidase subunits I and II (CO-I and CO-II) and inferred genes for tRNA(Asp), tRNA(Leu)UUR, tRNA(Lys), and tRNA(Trp), is presented. The region includes the segment previously identified as incurring a length increase in some other bee strains, including Africanized bees. The sequence information of this study and of that by Vlasak et al. shows that several shifts of tRNA genes have occurred between Apis and Drosophila, but shifts of other kinds of genes have yet to be demonstrated. The CO-I and CO-II gene sequences are both more A+T rich than are the corresponding Drosophila genes. Parsimony analyses using the mouse and Xenopus sequences as outgroups show significantly more amino acid substitutions on the branch to Apis (120) than on that to Drosophila (44), indicating a difference in the long-term evolutionary rates of hymenopteran and dipteran mtDNA.     

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.     

Current Intraspecific Dynamics of Sequence Evolution Differs from Long-Term Trends and Can Account for the AT-Richness of Honeybee Mitochondrial DNA

The very high AT content of hymenopteran mtDNA has warranted speculation about nucleotide substitution processes in this group. Here we investigate the pattern of honeybee, Apis mellifera, mtDNA nucleotide polymorphisms inferred from phylogeny in terms of differences between the ATPase6, COI, COII, COIII, cytochrome b, and ND2 genes and strand asymmetry in mutation rates. The observed transition/transversion ratios and the distribution of nonsynonymous substitutions between regions differed significantly. The pattern of differences between genes leading to these heterogeneities (the ATPase6 and COIII genes group apart from the rest) differed markedly from that predicted on the basis of long-term evolutionary change and may indicate differences between current and long-term dynamics of sequence evolution. Also, there is strong strand asymmetry in substitutions, which probably results in a mutability of G and C sufficiently high to account for the AT-richness of honeybee mtDNA. Received: 21 October 1998 / Accepted: 27 January 1999     

Phylogeny of the Ascaridoidea (Nematoda: Ascaridida) based on three genes and morphology: hypotheses of structural and sequence evolution

Ascaridoid nematodes parasitize the gastrointestinal tract of vertebrate definitive hosts and are represented by more than 50 described genera. We used 582 nucleotides (83% of the coding sequence) of the mitochondrial gene cytochrome oxidase subunit 2, in combination with published small- and large-subunit nuclear rDNA sequences (2,557 characters) and morphological data (20 characters), to produce a phylogenetic hypothesis for representatives of this superfamily. This combined evidence phylogeny strongly supported clades that, with 1 exception, were consistent with Fagerholm's 1991 classification. Parsimony mapping of character states on the combined evidence tree was used to develop hypotheses for the evolution of morphological, life history, and amino acid characters. This analysis of character evolution revealed that certain key features that have been used by previous workers for developing taxonomic and evolutionary hypotheses represent plesiomorphic states. Cytochrome oxidase subunit 2 nucleotides show a strong compositional bias to A+T and a substitution bias to thymine. These biases are most apparent at third positions of codons and 4-fold degenerate sites, which is consistent with the nonrandom substitution pattern of A+T pressure. Despite nucleotide bias, cytochrome oxidase amino acid sequences show conservation and retention of critical functional residues, as inferred from comparisons to other organisms.     

The axial ligands of heme in cytochromes: a near-infrared magnetic circular dichroism study of yeast cytochromes c, c1, and b and spinach cytochrome f

Room temperature near-infrared magnetic circular dichroism and low-temperature electron paramagnetic resonance measurements have been used to characterize the ligands of the heme iron in mitochondrial cytochromes c, c1, and b and in cytochrome f of the photosynthetic electron transport chain. The MCD data show that methionine is the sixth ligand of the heme of oxidized yeast cytochrome c1; the identify of this residue is inferred to be the single conserved methionine identified from a partial alignment of the available cytochrome c1 amino acid sequences. A different residue, which is most likely lysine, is the sixth heme ligand in oxidized spinach cytochrome f. The data for oxidized yeast cytochrome b are consistent with bis-histidine coordination of both hemes although the possibility that one of the hemes is ligated by histidine and lysine cannot be rigorously excluded. The neutral and alkaline forms of oxidized yeast cytochrome c have spectroscopic properties very similar to those of the horse heart proteins, and thus, by analogy, the sixth ligands are methionine and lysine, respectively.     

Increased genetic diversity in mitochondrial genes is correlated with the evolution of parasitism in the Hymenoptera

A higher AT content and rate of mtDNA sequence divergence was found in parasitic wasps (Apocrita) compared with nonparasitic wasps (Symphyta). The compositional bias was reflected in extreme codon bias for a cytochrome oxidase I protein coding gene fragment as well as in the types of amino acid substitutions that have occurred during the evolution of this gene fragment. In some instances, compositional bias influenced the definition of a conservative amino acid change. The increased rate of mtDNA sequence evolution probably arose during the early Jurassic, coincident with the first appearance of parasitic wasps in the fossil record. Our results suggest a causal link between the rate of sequence divergence and the parasitic lifestyle.Abbreviations AT adenosine-thymine - CO-1 cytochrome oxidase 1 - mtDNA mitochondrial DNA - Myr million years Correspondence to: M. Dowton     

Evolution of the mitochondrial DNA cytochrome b gene in Tetraonidae birds

Mitochondrial fragments containing the cytochrome b gene (1020 bp in size) of four bird species belonging to four genera of the family Tetraonidae (Tetrao parvirostris, Bonasa umbellus, Lagopus lagopus scoticus, and Falcipennis falcipennis) were directly sequenced. Of the 1020 nucleotide positions, 186 were variable and uniformly distributed over the gene and only 46 were parsimony informative. Most substitutions were synonymous. Replacement substitutions were detected for 15 out of 340 amino acid sites; only four replacements were parsimony informative. The greatest codon bias was found for leucine and serine. The C-T transitions and the G-C transversions were, respectively, the most common (60.7%) and the most rare (5.9%). The mutation frequencies were high at the third codon position (85.2%) and relatively low at the first and the second position. At the third codon position of the species examined, the guanine content was the lowest (3.3%) and the cytosine content was the highest (44.5%). Based on the cytochrome b gene sequences, phylogenetic relationships in the order Galliformes are inferred.     

Evolution of the WANCY region in amniote mitochondrial DNA

In most vertebrate mitochondrial genomes, the site for initiation of light-strand replication, OL, is found within a cluster of five transfer RNA (tRNA) genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), and tRNA(Tyr)). This region and part of the adjacent cytochrome c oxydase subunit I (COI) gene were sequenced for two crocodilian, two turtle, and one snake species and for Sphenodon punctatus; part of the adjacent nicotinamide adenine dinucleotide dehydrogenase subunit 2 (ND2) gene was also sequenced for the crocodilian and turtle species. All had the typical vertebrate gene order. The turtles and the snake have a lengthy noncoding sequence between the tRNA(Asn) and tRNA(Cys) genes that we assumed to be homologous to the mammalian OL. The crocodilians and Sphenodon lack such a sequence, a condition they share with birds. Most proposed phylogenies for the amniotes require that OL at this position was lost at least twice during their diversification or was evolved independently more than once. Within the five tRNA genes, frequencies of substitutions are much higher in loops than in stems. Many loops vary dramatically in size among the species; in the most extreme case, the D-arm of the Sphenodon tRNA(Cys) is a "D-arm replacement" loop of seven nucleotides. Frequency of transitions in stems is relatively uniform across tRNAs, but frequency of transversions varies greatly. Mismatches in stems are infrequent, and their relative frequency in a specific tRNA is unrelated to the frequency of substitution in the corresponding gene. Several features of mammalian mitochondrial tRNAs are conserved in WANCY tRNAs throughout amniotes. The inferred initiation codon for COI is GTG in crocodilians, turtles, and the snake, a condition they share with fishes, certain amphibians, and birds. TTG appears to be the initiation codon for COI in Sphenodon; if correct, this would be a novel initiation codon for vertebrate mitochondrial DNA. Phylogenetic analyses of the inferred amino acid sequences of ND2 and COI support the sister-group relationship of birds and crocodilians and suggest that mammals are an early derived lineage within the amniotes.      

Heterogeneity in regional GC content and differential usage of codons and amino acids in GC-poor and GC-rich regions of the genome of Apis mellifera

The honeybee (Apis mellifera) has a genome with a wide variation in GC content showing 2 clear modal GC values, in some ways reminiscent of an isochore-like structure. To gain insight into causes and consequences of this pattern, we used a comparative approach to study the genome-wide alignment of primarily coding sequence of A. mellifera with Drosophila melanogaster and Anopheles gambiae. The latter 2 species show a higher average GC content than A. mellifera and no indications of bimodality, suggesting that the GC-poor mode is a derived condition in honeybee. In A. mellifera, synonymous sites of genes generally adopt the GC content of the region in which they reside. A large proportion of genes in GC-poor regions have not been assigned to the honeybee assembly because of the low sequence complexity of their genome neighborhood. The synonymous substitution rate between A. mellifera and the other species is very close to saturation, but analyses of nonsynonymous substitutions as well as amino acid substitutions indicate that the GC-poor regions are not evolving faster than the GC-rich regions. We describe the codon usage and amino acid usage and show that they are remarkably heterogeneous within the honeybee genome between the 2 different GC regions. Specifically, the genes located in GC-poor regions show a much larger deviation in both codon usage bias and amino acid usage from the Dipterans than the genes located in the GC-rich regions.     

Numerous transposed sequences of mitochondrial cytochrome oxidase I-II in aphids of the genus Sitobion (Hemiptera: Aphididae)

Polymerase chain reaction (PCR) products corresponding to 803 bp of the cytochrome oxidase subunits I and II region of mitochondrial DNA (mtDNA COI-II) were deduced to consist of multiple haplotypes in three Sitobion species. We investigated the molecular basis of these observations. PCR products were cloned, and six clones from one individual per species were sequenced. In each individual, one sequence was found commonly, but also two or three divergent sequences were seen. The divergent sequences were shown to be nonmitochondrial by sequencing from purified mtDNA and Southern blotting experiments. All seven nonmitochondrial clones sequenced to completion were unique. Nonmitochondrial sequences have a high proportion of unique sites, and very few characters are shared between nonmitochondrial clones to the exclusion of mtDNA. From these data, we infer that fragments of mtDNA have been transposed separately (probably into aphid chromosomes), at a frequency only known to be equalled in humans. The transposition phenomenon appears to occur infrequently or not at all in closely related genera and other aphids investigated. Patterns of nucleotide substitution in mtDNA inferred over a parsimony tree are very different from those in transposed sequences. Compared with mtDNA, nonmitochondrial sequences have less codon position bias, more even exchanges between A, G, C and T, and a higher proportion of nonsynonymous replacements. Although these data are consistent with the transposed sequences being under less constraint than mtDNA, changes in the nonmitochondrial sequences are not random: there remains significant position bias, and probable excesses of synonymous replacements and of conservative inferred amino acid replacements. We conclude that a proportion of the inferred change in the nonmitochondrial sequences occurred before transposition. We believe that Sitobion aphids (and other species exhibiting mtDNA transposition) may be important for studying the molecular evolution of mtDNA and pseudogenes. However, our data highlight the need to establish the true evolutionary relationships between sequences in comparative investigations.      

The largemouth bass cytochrome b gene

The cytochrome b gene is one of the protein-coding genes of the mitochondrial genome that has gained importance because of the ease with which molecular techniques can be applied to the analysis of its structure. The nucleotide sequence of the largemouth bass ( Microplerus salmoides ) cytochrome b gene was determined and the inferred amino acid structure is presented in the form of a structural model derived originally from rat cytochrome b . The inferred amino acid sequences from divergent animal species are aligned and compared in the context of this model. The data suggest that regions of the gene may be evolving at different rates due to different selection pressures associated with functional constraints. Conserved and variable regions of cytochrome b have been identified and can be targeted for species identification, the examination of intraspeciflc variation, and phylogenetic reconstructions in future research.     

Interaction of silent and replacement changes in eukaryotic coding sequences

We examined the codon usages in well-conserved and less-well-conserved regions of vertebrate protein genes and found them to be similar. Despite this similarity, there is a statistically significant decrease in codon bias in the less-well-conserved regions. Our analysis suggests that although those codon changes initially fixed under amino acid replacements tend to follow the overall codon usage pattern, they also reduce the bias in codon usage. This decrease in codon bias leads one to predict that the rate of change of synonymous codons should be greater in those regions that are less well conserved at the amino acid level than in the better-conserved regions. Our analysis supports this prediction. Furthermore, we demonstrate a significantly elevated rate of change of synonymous codons among the adjacent codons 5' to amino acid replacement positions. This provides further support for the idea that there are contextual constraints on the choice of synonymous codons in eukaryotes.     

The actin genes of Drosophila: protein coding regions are highly conserved but intron positions are not

The entire set of six closely related Drosophila actin genes was isolated using recombinant DNA methodology, and the structures of the respective coding regions were characterized by gene mapping techniques and by nucleotide sequencing of selected portions. Structural comparisons of these genes have resulted in several unexpected findings. Most striking is the nonconservation of the positions of intervening sequences within the protein-encoding regions of these genes. One of the Drosophila actin genes, DmA4, is split within a glycine codon at position 13; none of the remaining five genes is interrupted in the analogous position. Another gene, DmA6, is split within a glycine codon at position 307; at least two of the Drosophila actin genes are not split in the analogous position. Additionally, none of the Drosophila actin genes is split within codon four, where the yeast actin gene is interrupted. The six Drosophila actin genes encode several different proteins, but the amino acid sequence of each is similar to that of vertebrate cytoplasmic actins. None of the genes encodes a protein comparable in primary sequence to vertebrate skeletal muscle actin. Surprisingly, in each of these derived actin amino acid sequences in the initiator methionine is directly followed by a cysteine residue, which in turn precedes the string of three acidic amino acids characteristic of the amino termini of mature vertebrate cytoplasmic actins. We discuss these findings in the context of actin gene evolution and function.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Multisequence comparisons in protein coding genes. Search for functional constraints

A very powerful method for detecting functional constraints operative in biological macromolecules is presented. This method entails performing a base permanence analysis of protein coding genes at each codon position simultaneously in different species. It calculates the degree of permanence of subregions of the gene by dividing it into segments, c codons long, counting how many sites remain unchanged in each segment among all species compared. By comparing the base permanence among several sequences with the expectations based on a stochastic evolutionary process, gene regions showing different degrees of conservation can be selected. This means that wherever the permanence deviates significantly from the expected value generated by the simulation, the corresponding regions are considered "constrained" or "hypervariable". The constrained regions are of two types: alpha and beta. The alpha regions result from constraints at the amino acid level, whereas the beta regions are those probably involved in "control" processing. The method has been applied to mitochondrial genes coding for subunit 6 of the ATPase and subunit 1 of the cytochrome oxidase in four mammalian species: human, rat, mouse, and cow. In the two mitochondrial genes a few regions that are highly conserved in all codon positions have been identified. Among these regions a sequence, common to both genes, that is complementary to a strongly conserved region of 12S rRNA has been found. This method can also be of great help in studying molecular evolution mechanisms.