Domestication of olive fly through a multi-regional host shift to cultivated olives: comparative dating using complete mitochondrial genomes.

The evolutionary history of the olive fly, Bactrocera oleae, was reconstructed in a phylogenetic and coalescent framework using full mitochondrial genome data from 21 individuals covering the entire worldwide distribution of the species. Special attention was given to reconstructing the timing of the processes under study. The early subdivision of the olive fly reflects the Quaternary differentiation between Olea europea subsp. europea in the Mediterranean area and the two lineages of Olea europea subsp. cuspidata in Africa and Asia, pointing to an early and close association between the olive fly and its host. The geographic structure and timing of olive fly differentiation in the Mediterranean indicates a clear connection with the post-glacial recolonization of wild olives in the area, and is irreconcilable with the early historical process of domestication and spread of the cultivated olive from its Levantine origin. Therefore, we suggest an early co-history of the olive fly with its wild host during the Quaternary and post-glacial periods and a multi-regional shift of olive flies to cultivated olives as these cultivars gradually replaced wild olives in historical times.

The mitochondrial genomes of Campodea fragilis and Campodea lubbocki (Hexapoda: Diplura): High genetic divergence in a morphologically uniform taxon.

Complete mitochondrial genome sequences are presented from two dipluran hexapods (i.e., a group of "primarily wingless insects") of the genus Campodea and compared to those of other arthropods. Their gene order is the same as in most other hexapods and crustaceans. Structural changes have occurred in tRNA-C, tRNA-R, tRNA-S1 and tRNA-S2 as well as in both ribosomal RNAs. These mtDNAs have striking biases in nucleotide and amino acid composition. Although the two Campodea species are morphologically highly similar, their genetic divergence is larger than expected, suggesting a long evolutionary history, perhaps under stable ecological conditions.

Genes without frontiers?

For bacteria, the primary genetic barrier against the genetic exchange of DNA that is not self-transmissible is dissimilarity in the bacterial DNA sequences concerned. Genetic exchange by homologous recombination is frequent among close bacterial relatives and recent experiments have shown that it can enable the uptake of closely linked nonhomologous foreign DNA. Artificial vectors are mosaics of mobile DNA elements from free-living bacterial isolates and so bear a residual similarity to their ubiquitous natural progenitors. This homology is tightly linked to the multitude of different DNA sequences that are inserted into synthetic vectors. Can homology between vector and bacterial DNA enable the uptake of these foreign DNA inserts? In this review we investigate pUC18 as an example of an artificial vector and consider whether its homology to broad host-range antibiotic resistance transposons and plasmid origins of replication could enable the uptake of insert DNA in the light of studies of homology-facilitated foreign DNA uptake. We also discuss the disposal of recombinant DNA, its persistence in the environment and whether homologies to pUC18 may exist in naturally competent bacteria. Most DNA that is inserted into the cloning site of artificial vectors would be of little use to a bacterium, but perhaps not all.

The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa.

We sequenced the complete mitochondrial DNA (mtDNA) of the articulate brachiopod Terebratalia transversa. The circular genome is 14,291 bp in size, relatively small compared with other published metazoan mtDNAs. The 37 genes commonly found in animal mtDNA are present; the size decrease is due to the truncation of several tRNA, rRNA, and protein genes, to some nucleotide overlaps, and to a paucity of noncoding nucleotides. Although the gene arrangement differs radically from those reported for other metazoans, some gene junctions are shared with two other articulate brachiopods, Laqueus rubellus and Terebratulina retusa. All genes in the T. transversa mtDNA, unlike those in most metazoan mtDNAs reported, are encoded by the same strand. The A+T content (59.1%) is low for a metazoan mtDNA, and there is a high propensity for homopolymer runs and a strong base-compositional strand bias. The coding strand is quite G+T-rich, a skew that is shared by the confamilial (laqueid) species L. rubellus but is the opposite of that found in T. retusa, a cancellothyridid. These compositional skews are strongly reflected in the codon usage patterns and the amino acid compositions of the mitochondrial proteins, with markedly different usages being observed between T. retusa and the two laqueids. This observation, plus the similarity of the laqueid noncoding regions to the reverse complement of the noncoding region of the cancellothyridid, suggests that an inversion that resulted in a reversal in the direction of first-strand replication has occurred in one of the two lineages. In addition to the presence of one noncoding region in T. transversa that is comparable with those in the other brachiopod mtDNAs, there are two others with the potential to form secondary structures; one or both of these may be involved in the process of transcript cleavage.

Complete sequence of the mitochondrial genome of the tapeworm Hymenolepis diminuta: gene arrangements indicate that Platyhelminths are Eutrochozoans.

Using "long-PCR," we amplified in overlapping fragments the complete mitochondrial genome of the tapeworm Hymenolepis diminuta (Platyhelminthes: Cestoda) and determined its 13,900-nt sequence. The gene content is the same as that typically found for animal mitochondrial DNA (mtDNA) except that atp8 appears to be lacking, a condition found previously for several other animals. Despite the small size of this mtDNA, there are two large noncoding regions, one of which contains 13 repeats of a 31-nt sequence and a potential stem-loop structure of 25 bp with an 11-member loop. Large potential secondary structures were identified also for the noncoding regions of two other cestode mtDNAS: Comparison of the mitochondrial gene arrangement of H. diminuta with those previously published supports a phylogenetic position of flatworms as members of the Eutrochozoa, rather than placing them basal to either a clade of protostomes or a clade of coelomates.

A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus.

We determined the complete mtDNA sequence of the centipede Lithobius forficatus and found that only one of the 22 inferred tRNA genes encodes a fully paired aminoacyl acceptor stem. The other 21 genes encode tRNAs with up to five mismatches in these stems, and some of these overlap extensively with the downstream genes. Because a well-paired acceptor stem is required for proper tRNA functioning, RNA editing in the products of these genes was suspected. We investigated this hypothesis by studying cDNA sequences from eight tRNAs and found the editing of up to 5 nt at their 3' ends. This editing appears to occur by a novel mechanism with the 5' end of the acceptor stem being used as a template for the de novo synthesis of the 3' end, presumably by an RNA-dependent RNA polymerase. In addition, unusual secondary structures for several tRNAs were found, including those lacking a TPsiC (T) or a dihydrouridine (D) arm, and having an unusual number of base pairs in the acceptor or anticodon stems.

The complete mitochondrial DNA sequence of the horseshoe crab Limulus polyphemus.

We determined the complete 14,985-nt sequence of the mitochondrial DNA of the horseshoe crab Limulus polyphemus (Arthropoda: Xiphosura). This mtDNA encodes the 13 protein, 2 rRNA, and 22 tRNA genes typical for metazoans. The arrangement of these genes and about half of the sequence was reported previously; however, the sequence contained a large number of errors, which are corrected here. The two strands of Limulus mtDNA have significantly different nucleotide compositions. The strand encoding most mitochondrial proteins has 1. 25 times as many A's as T's and 2.33 times as many C's as G's. This nucleotide bias correlates with the biases in amino acid content and synonymous codon usage in proteins encoded by different strands and with the number of non-Watson-Crick base pairs in the stem regions of encoded tRNAs. The sizes of most mitochondrial protein genes in Limulus are either identical to or slightly smaller than those of their Drosophila counterparts. The usage of the initiation and termination codons in these genes seems to follow patterns that are conserved among most arthropod and some other metazoan mitochondrial genomes. The noncoding region of Limulus mtDNA contains a potential stem-loop structure, and we found a similar structure in the noncoding region of the published mtDNA of the prostriate tick Ixodes hexagonus. A simulation study was designed to evaluate the significance of these secondary structures; it revealed that they are statistically significant. No significant, comparable structure can be identified for the metastriate ticks Rhipicephalus sanguineus and Boophilus microplus. The latter two animals also share a mitochondrial gene rearrangement and an unusual structure of mt-tRNA(C) that is exactly the same association of changes as previously reported for a group of lizards. This suggests that the changes observed are not independent and that the stem-loop structure found in the noncoding regions of Limulus and Ixodes mtDNA may play the same role as that between trnN and trnC in vertebrates, i.e., the role of lagging strand origin of replication.

Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa.

We report a contiguous region of more than half (> 7,500 nt) of the mitochondrial genomes for Platynereis dumerii (Annelida: Polychaeta), Helobdella robusta (Annelida: Hirudinida), and Galathealinum brachiosum (Pogonophora: Perviata). The relative arrangements of all 22 genes identified for Helobdella and Galathealinum are identical to one another and to their arrangements in the mtDNA of the previously studied oligochaete annelid Lumbricus. In contrast, Platynereis differs from these taxa in the positions of several tRNA genes and in having two additional tRNA genes (trnC and trnM) and a large noncoding sequence in this region. Comparisons of relative gene arrangements and of the nucleotide and inferred amino acid sequences among these and other published taxa provide strong support for an annelid-mollusk clade that excludes arthropods, and for the inclusion of pogonophorans within Annelida, rather than giving them separate phylum status. Gene arrangement comparisons include the first use of a recently described method on previously unpublished data. Although a variety of alternative initiation codons are typically used by mitochondrial protein-encoding genes, ATG appears to be the initiator for all but one reported here. The large noncoding region (1,091 nt) identified in Platynereis has no significant sequence similarity to the noncoding region of Lumbricus, although each contains runs of TA dinucleotides and of homopolymers, which could potentially serve as signaling elements. There is strong bias for synonymous codon usage in Helobdella and especially in Galathealinum. In this latter taxon, 5 codons are completely unused, 13 are used three or fewer times, and G appears at third codon positions in only 26 of the 2,236 codons. Nucleotide composition bias appears to influence amino acid composition of the proteins.

Complete sequence, gene arrangement, and genetic code of mitochondrial DNA of the cephalochordate Branchiostoma floridae (Amphioxus)

We have determined the 15,083-nucleotide (nt) sequence of the mitochondrial DNA (mtDNA) of the lancelet Branchiostoma floridae (Chordata: Cephalochordata). As is typical in metazoans, the mtDNA encodes 13 protein, 2 rRNA, and 22 tRNA genes. The gene arrangement differs from the common vertebrate arrangement by only four tRNA gene positions. Three of these are unique to Branchiostoma, but the fourth is in a position that is primitive for chordates. It shares the genetic code variations found in vertebrate mtDNAs except that AGA = serine, a code variation found in many invertebrate phyla but not in vertebrates (the related codon AGG was not found). Branchiostoma mtDNA lacks a vertebrate-like control region; its largest noncoding region (129 nt) is unremarkable in sequence or base composition, and its location between ND5 and tRNAG differs from that usually found in vertebrates. It also lacks a potential hairpin DNA structure like those found in many (though not in all) vertebrates to serve as the second-strand (i.e., L-strand) origin of replication. Perhaps related to this, the sequence corresponding to the DHU arm of tRNAC cannot form a helical stem, a condition found in a few other vertebrate mtDNAs that also lack a canonical L-strand origin of replication. ATG and GTG codons appear to initiate translation in 11 and 2 of the protein-encoding genes, respectively. Protein genes end with complete (TAA or TAG) or incomplete (T or TA) stop codons; the latter are presumably converted to TAA by post-transcriptional polyadenylation.

Animal mitochondrial genomes.

Animal mitochondrial DNA is a small, extrachromosomal genome, typically approximately 16 kb in size. With few exceptions, all animal mitochondrial genomes contain the same 37 genes: two for rRNAs, 13 for proteins and 22 for tRNAs. The products of these genes, along with RNAs and proteins imported from the cytoplasm, endow mitochondria with their own systems for DNA replication, transcription, mRNA processing and translation of proteins. The study of these genomes as they function in mitochondrial systems-'mitochondrial genomics'-serves as a model for genome evolution. Furthermore, the comparison of animal mitochondrial gene arrangements has become a very powerful means for inferring ancient evolutionary relationships, since rearrangements appear to be unique, generally rare events that are unlikely to arise independently in separate evolutionary lineages. Complete mitochondrial gene arrangements have been published for 58 chordate species and 29 non-chordate species, and partial arrangements for hundreds of other taxa. This review compares and summarizes these gene arrangements and points out some of the questions that may be addressed by comparing mitochondrial systems.

Big trees from little genomes: mitochondrial gene order as a phylogenetic tool.

Gene arrangement comparisons are a powerful tool for phylogenetic studies, especially those focused on ancient relationships. Recent reports using metazoan mitochondrial genomes address evolutionary relationships as well as rates and mechanisms of rearrangement. Mitochondrial systems serve as a model for larger-scale comparisons of whole organismal genomes and a stimulus for developing methods for reconstructing the patterns of rearrangement.

Transmission of mitochondrial DNA--playing favorites?

Mitochondria are essential subcellular organelles containing an extranuclear genome (mtDNA). Mutations in mtDNA have recently been identified as causing a variety of human hereditary disease. In most of these cases, the tissues of the affected individual contain a mixture of mutant and normal mtDNA, with this ratio determining the severity of symptoms. Stochastic factors alone have generally been believed to determine this ratio. Jenuth et al.(1), however, examining mice that contain a mixture of mtDNA types, show evidence of strong selective forces at work in favoring one mtDNA variant over another in some tissues.

Complete sequence of the mitochondrial DNA of the annelid worm Lumbricus terrestris.

We have determined the complete nucleotide (nt) sequence of the mitochondrial genome of an oligochaete annelid, the earthworm Lumbricus terrestris. This genome contains the 37 genes typical of metazoan mitochondrial DNA (mtDNA), including ATPase8, which is missing from some invertebrate mtDNAs. ATPase8 is not immediately upstream of ATPase6, a condition found previously only in the mtDNA of snails. All genes are transcribed from the same DNA strand. The largest noncoding region is 384 nt and is characterized by several homopolymer runs, a tract of alternating TA pairs, and potential secondary structures. All protein-encoding genes either overlap the adjacent downstream gene or end at an abbreviated stop codon. In Lumbricus mitochondria, the variation of the genetic code that is typical of most invertebrate mitochondrial genomes is used. Only the codon ATG is used for translation initiation. Lumbricus mtDNA is A + T rich, which appears to affect the codon usage pattern. The DHU arm appears to be unpaired not only in tRNAser(AGN), as is typical for metazoans, but perhaps also in tRNAser(UCN), a condition found previously only in a chiton and among nematodes. Relating the Lumbricus gene organization to those of other major protostome groups requires numerous rearrangements.

Deducing the pattern of arthropod phylogeny from mitochondrial DNA rearrangements.

The origins of arthropods and the phylogenetic relationships among their three major living groups (atelocerates, crustaceans and chelicerates) are vigorously contended. To help resolve this, we determined mitochondrial gene arrangements for a chelicerate, a myriapod, two crustaceans, an onychophoran, a mollusc and an annelid, and compared them with published gene orders of other species. The result strongly supports the monophyly of Arthropoda and of Mandibulata (atelocerates plus crustaceans) and refutes the Uniramia (atelocerates plus onychophorans). Gene arrangement comparisons are emerging as a powerful new tool for resolving ancient phylogenetic relationships.

Relocation of mitochondria to the prospective dorsal marginal zone during Xenopus embryogenesis.

Dorsal-ventral axis formation in Xenopus laevis begins with a cytoplasmic rotation during the first cell cycle and culminates in a series of cell interactions and movements during gastrulation and neurulation that lead to the formation of dorsal-anterior structures. Evidence reported here indicates that mitochondria are differentially redistributed along the prospective dorsal-ventral axis as a consequence of the cortical-cytoplasmic rotation during the first cell cycle. This finding reinvigorates a possibility that has been considered for many years: asymmetries in cytoplasmic components and metabolic activities contribute to the development of morphological asymmetries.

Complete DNA sequence of the mitochondrial genome of the black chiton, Katharina tunicata.

The DNA sequence of the 15,532-base pair (bp) mitochondrial DNA (mtDNA) of the chiton Katharina tunicata has been determined. The 37 genes typical of metazoan mtDNA are present: 13 for protein subunits involved in oxidative phosphorylation, 2 for rRNAs and 22 for tRNAs. The gene arrangement resembles those of arthropods much more than that of another mollusc, the bivalve Mytilus edulis. Most genes abut directly or overlap, and abbreviated stop codons are inferred for four genes. Four junctions between adjacent pairs of protein genes lack intervening tRNA genes; however, at each of these junctions there is a sequence immediately adjacent to the start codon of the downstream gene that is capable of forming a stem-and-loop structure. Analysis of the tRNA gene sequences suggests that the D arm is unpaired in tRNA(ser)(AGN), which is typical of metazoan mtDNAs, and also in tRNA(ser)(UCN), a condition found previously only in nematode mtDNAs. There are two additional sequences in Katharina mtDNA that can be folded into structures resembling tRNAs; whether these are functional genes is unknown. All possible codons except the stop codons TAA and TAG are used in the protein-encoding genes, and Katharina mtDNA appears to use the same variation of the mitochondrial genetic code that is used in Drosophila and Mytilus. Translation initiates at the codons ATG, ATA and GTG. A + T richness appears to have affected codon usage patterns and, perhaps, the amino acid composition of the encoded proteins. A 142-bp non-coding region between tRNA(glu) and CO3 contains a 72-bp tract of alternating A and T.