Mitochondrial genome sequences of Nematocera (lower Diptera): evidence of rearrangement following a complete genome duplication in a winter crane fly.

The complete mitochondrial DNA sequences of eight representatives of lower Diptera, suborder Nematocera, along with nearly complete sequences from two other species, are presented. These taxa represent eight families not previously represented by complete mitochondrial DNA sequences. Most of the sequences retain the ancestral dipteran mitochondrial gene arrangement, while one sequence, that of the midge Arachnocampa flava (family Keroplatidae), has an inversion of the trnE gene. The most unusual result is the extensive rearrangement of the mitochondrial genome of a winter crane fly, Paracladura trichoptera (family Trichocera). The pattern of rearrangement indicates that the mechanism of rearrangement involved a tandem duplication of the entire mitochondrial genome, followed by random and nonrandom loss of one copy of each gene. Another winter crane fly retains the ancestral diperan gene arrangement. A preliminary mitochondrial phylogeny of the Diptera is also presented.

Mitochondrial genome sequences of representatives of three families of scorpionflies (Order Mecoptera) and evolution in a major duplication of coding sequence.

The complete mitochondrial DNA sequences of a hangingfly, Bittacus pilicornis (Mecoptera: Bittacidae), a snow scorpion fly, Boreus elegans (Mecoptera: Boreidae), and a nearly complete sequence from another scorpionfly species, Microchorista philpotti (Mecoptera: Nannochoristidae) were determined. The coding sequence of all three genomes includes the 37 genes normally found in insect mtDNAs, in the same gene order as first described in Drosophila. In addition to the standard set of genes, the Microchorista sequence includes a large duplication of the coding region. The duplication is at least 4 kb (and may be much larger) and includes the remnants of three protein-coding genes and seven tRNA genes. The duplication evidently arose as a single event, and the duplicated region can be aligned in its entirety with the corresponding region of the functional genome. Although most of the genes contain defects that render them nonfunctional, analysis shows that the protein-coding genes in the duplicated region evolved for a considerable period under constraints expected of functional protein-coding genes. It is evident, therefore, that for a period two copies of some of the mitochondrial genes were functional in this species, including genes coding for proteins.

Episodic radiations in the fly tree of life.

Flies are one of four superradiations of insects (along with beetles, wasps, and moths) that account for the majority of animal life on Earth. Diptera includes species known for their ubiquity (Musca domestica house fly), their role as pests (Anopheles gambiae malaria mosquito), and their value as model organisms across the biological sciences (Drosophila melanogaster). A resolved phylogeny for flies provides a framework for genomic, developmental, and evolutionary studies by facilitating comparisons across model organisms, yet recent research has suggested that fly relationships have been obscured by multiple episodes of rapid diversification. We provide a phylogenomic estimate of fly relationships based on molecules and morphology from 149 of 157 families, including 30 kb from 14 nuclear loci and complete mitochondrial genomes combined with 371 morphological characters. Multiple analyses show support for traditional groups (Brachycera, Cyclorrhapha, and Schizophora) and corroborate contentious findings, such as the anomalous Deuterophlebiidae as the sister group to all remaining Diptera. Our findings reveal that the closest relatives of the Drosophilidae are highly modified parasites (including the wingless Braulidae) of bees and other insects. Furthermore, we use micro-RNAs to resolve a node with implications for the evolution of embryonic development in Diptera. We demonstrate that flies experienced three episodes of rapid radiation--lower Diptera (220 Ma), lower Brachycera (180 Ma), and Schizophora (65 Ma)--and a number of life history transitions to hematophagy, phytophagy, and parasitism in the history of fly evolution over 260 million y.

Characterization of mature mitochondrial transcripts in Drosophila, and the implications for the tRNA punctuation model in arthropods.

The mitochondrial genome of Drosophila melanogaster is thought to be transcribed in about five polycistronic primary transcripts, which are processed into 11 mRNA, 22 tRNA and two rRNA species required for the genome's function. The tRNA punctuation model has been proposed to predict the cleavage sites used in this process. In this model, tRNAs are removed from the primary transcripts and the fragments that remain become mRNA and rRNA transcripts. Thus the 5' and 3' ends of the major gene transcripts are defined by the endpoints of the intervening tRNA sequences. We used 5' and 3' RACE, and circularization and RT-PCR methods to determine the sequences of both ends of all major gene transcripts of the D. melanogaster mitochondrial genome. In general, the tRNA punctuation model accurately predicts the 3' ends of most mRNA and rRNA molecules, even where there are non-coding residues present. Non-coding residues at the 5' end are evidently removed during RNA processing. The mRNAs begin precisely at the start codon for each gene. In particular, the 5' end of the cox1 gene is the first in-frame sense codon, UCG, implying that this codon serves as the start. In-frame TAA stop codons immediately preceding both the cox1 and cox2 genes may serve to prevent in-frame translation of these genes prior to the completion of processing, and are removed from the mature transcripts. Where multiple tRNA genes are present, as between nad3 and nad5, they are removed sequentially in a 3' to 5' direction.

Insect mitochondrial genomics 3: the complete mitochondrial genome sequences of representatives from two neuropteroid orders: a dobsonfly (order Megaloptera) and a giant lacewing and an owlfly (order Neuroptera).

We describe the complete mitochondrial genomes from representatives of two orders of the Neuropterida: a dobsonfly, Corydalus cornutus (Megaloptera: Corydalidae, GenBank Accession No. FJ171323), a giant lacewing Polystoechotes punctatus (Neuroptera: Polystoechotidae, FJ171325), and an owlfly, Ascaloptynx appendiculatus (Neuroptera: Ascalaphidae, FJ171324). The dobsonfly sequence is 15,687 base pairs with a major noncoding (A+T rich) region of approximately 967 bp. The gene content and organization of the dobsonfly is identical to that of most insects. The giant lacewing sequence is 16 036 bp with a major noncoding region of about 1123 bp, while the owlfly sequence is 15,877 bp with a major noncoding region of about 1066 bp. The two Neuroptera sequences include a transposition of two tRNA genes, tRNATrp and tRNACys. These tRNA genes are coded on opposite strands and overlap by seven residues in the standard insect mitochondrial gene arrangement. Thus, the transposition required a duplication of at least the region of overlap. It is likely that the transposition occurred by a duplication of both genes followed by deletion of one copy of each gene. Examination of this region in two other neuropteroid species, a snakefly, Agulla sp. (Raphidioptera: Raphidiidae), and an antlion, Myrmeleon immaculatus (Neuroptera: Myrmeleontidae), shows that the rearrangement is widespread in the order Neuroptera but not present in either of the other two orders of Neuropterida.

Evolution of the mitochondrial genomes of gall midges (Diptera: Cecidomyiidae): rearrangement and severe truncation of tRNA genes.

We determined the complete mitochondrial genome sequences of two species of gall midges (Diptera: Cecidomyiidae), as well as partial sequence from a third cecidomyiid and a species from a related family, the Sciaridae. The sciarid sequence has a number of rearrangements of tRNA genes, relative to other dipterans, but is otherwise unremarkable. In contrast, the cecidomyiid genomes possess a number of very unusual features. First, the two complete sequences are very small compared with other dipteran mitochondrial genomes. The genome of Mayetiola destructor is only 14,759 bp while that of Rhopalomyia pomum is only 14,503 bp, comparable with genome sizes observed in some arachnids. Second, all three cecidomyiid species have very high A + T content--more than 83% for the coding region. Third, all three cecidomyiid species possess a number of rearrangements of tRNA genes, including variations within the family. Fourth, the most extraordinary feature of cecidomyiids examined in this study is an extreme truncation of all tRNA genes, including the loss of TPsiC arms and apparent absence of the 3' part of the aminoacyl stems.The truncated tRNA genes of cecidomyiids are very similar to those previously reported for spiders and appear to represent a second, independent origin of these structural features. It is likely that they are made functional through RNA editing, perhaps using the 5' end of the aminoacyl stem as a template for the construction of the required 3' end.

Recoding of translation in turtle mitochondrial genomes: programmed frameshift mutations and evidence of a modified genetic code.

A +1 frameshift insertion has been documented in the mitochondrial gene nad3 in some birds and reptiles. By sequencing polyadenylated mRNA of the chicken (Gallus gallus), we have shown that the extra nucleotide is transcribed and is present in mature mRNA. Evidence from other animal mitochondrial genomes has led us to hypothesize that certain mitochondrial translation systems have the ability to tolerate frameshift insertions using programmed translational frameshifting. To investigate this, we sequenced the mitochondrial genome of the red-eared slider turtle (Trachemys scripta), where both the widespread nad3 frameshift insertion and a novel site in nad4l were found. Sequencing the region surrounding the insertion in nad3 in a number of other turtles and tortoises reveal general mitochondrial +1 programmed frameshift site features as well as the apparent redefinition of a stop codon in Parker's snake-neck turtle (Chelodina parkeri), the first known example of this in vertebrate mitochondria.

Insect mitochondrial genomics 2: The complete mitochondrial genome sequence of a giant stonefly, Pteronarcys princeps, asymmetric directional mutation bias, and conserved plecopteran A+T-region elements.

Mitochondrial (mt) genome sequences of insects are receiving renewed attention in molecular phylogentic studies, studies of mt-genome rearrangement, and other unusual molecular phenomena, such as translational frameshifting. At present, the basal neopteran lineages are poorly represented by mt-genome sequences. Complete mt-genome sequences are available in the databases for only the Orthoptera and Blatteria; 9 orders are unrepresented. Here, we present the complete mt-genome sequence of a giant stonefly, Pteronarcys princeps (Plecoptera; Pteronarcyidae). The 16,004 bp genome is typical in its genome content, gene organisation, and nucleotide composition. The genome shows evidence of strand-specific mutational biases, correlated with the time between the initiation of leading and the initiation of lagging strand replication. Comparisons with other insects reveal that this trend is seen in other insect groups, but is not universally consistent among sampled mt-genomes. The A+T region is compared with that of 2 stoneflies in the family Peltoperlidae. Conserved stem-loop structures and sequence blocks are noted between these distantly related families.

The phylogeny of the subgroups within the melanogaster species group: likelihood tests on COI and COII sequences and a Bayesian estimate of phylogeny.

The relationships among the majority of the subgroups in the Drosophila melanogaster species group remain unresolved. We present a 2223basepair dataset for mitochondrial cytochrome oxidase I and cytochrome oxidase II for 43 species (including new data from 11 species), sampled to include the major subgroups. After a brief review of competing hypotheses for the ananassae, montium, suzukii, and takahashii subgroups, we combine the two genes based on a new use of the SH test and present KH and SH likelihood comparisons (Kishino and Hasegawa, 1989. J. Mol. Evol. 29, 170-179; Shimodaira and Hasegawa, 1999) to test the monophyly and placement of these subgroups within the larger species group. Although we find insignificant differences between the two suggested placements for the ananassae subgroup, the ananassae is sister to the rest of the subgroups in the melanogaster species group in every investigation. For the takahashii subgroup, although we cannot reject monophyly, the species are so closely related to the suzukii subgroup for these data that the two subgroups often form one clade. Finally, we present a Bayesian estimate of the phylogeny for both genes combined, utilizing a recently published method that allows for different models of evolution for different sites.

Single nucleotide +1 frameshifts in an apparently functional mitochondrial cytochrome b gene in ants of the genus Polyrhachis.

Twelve of 30 species examined in the ant genus Polyrhachis carry single nucleotide insertions at one or two positions within the mitochondrial cytochrome b (cytb) gene. Two of the sites are present in more than one species. Nucleotide substitutions in taxa carrying insertions show the strong codon position bias expected of functional protein coding genes, with substitutions concentrated in the third positions of the original reading frame. This pattern of evolution of the sequences strongly suggests that they are functional cytb sequences. This result is not the first report of +1 frameshift insertions in animal mitochondrial genes. A similar site was discovered in vertebrates, where single nucleotide frameshift insertions in many birds and a turtle were reported by Mindell et al. (Mol Biol Evol 15:1568, 1998). They hypothesized that the genes are correctly decoded by a programmed frameshift during translation. The discovery of four additional sites gives us the opportunity to look for common features that may explain how programmed frameshifts can arise. The common feature appears to be the presence of two consecutive rare codons at the insertion site. We hypothesize that the second of these codons is not efficiently translated, causing a pause in the translation process. During the stall the weak wobble pairing of the tRNA bound in the peptidyl site of the ribosome, together with an exact Watson-Crick codon-anticodon pairing in the +1 position, allows translation to continue in the +1 reading frame. The result of these events is an adequate level of translation of a full-length and fully functional protein. A model is presented for decoding of these mitochondrial genes, consistent with known features of programmed translational frameshifting in the yeast TY1 and TY3 retrotransposons.

Insect mitochondrial genomics: the complete mitochondrial genome sequence of the meadow spittlebug Philaenus spumarius (Hemiptera: Auchenorrhyncha: Cercopoidae).

We present the complete mitochondrial genome sequence of the meadow spittlebug Philaenus spumarius (Auchenorrhyncha: Cercopoidae). This contribution represents the second mitochondrial genome from the Hemiptera and the second of the three hemipteran suborders sampled. The genome is a circular molecule of 16 324 bp with a total A+T content of 77.0% and 76.7% for coding regions only. The gene content, order, and structure are consistent with the Drosophila yakuba genome structure (Clary and Wolstenholme 1985) and the hypothesized ancestral arthropod genome arrangement (Crease 1999). Nucleotide composition and codon usage are near the means observed in other insect mitochondria sequenced to date but have a higher A+T richness compared with the other hemipteran example, the kissing bug Triatoma dimidiata (Dotson and Beard. 2001. Insect Mol. Biol. 10: 205-215). The major noncoding region (the A+T rich region or putative control region) between the small ribosomal subunit and the tRNAIle gene includes two extensive repeat regions. The first repeat region includes 19 tandem repeats of a 46-bp sequence, whereas the second contains a longer sequence (146 bp) tandemly repeated four times.

Molecular analysis of the biting midges (Diptera: Ceratopogonidae), based on mitochondrial cytochrome oxidase subunit 2.

Sequences from the mitochondrial cytochrome oxidase subunit 2 gene (cox2) were determined for 14 species from the family Ceratopogonidae, representing 12 genera and all five subfamilies, along with six representatives of other nematoceran families. The purpose was to develop a molecular phylogeny of the Ceratopogonidae, and interpret the phylogenetic position of the family within the infraorder Culicomorpha. These taxa have been analysed using cladistic methodology which, in combination with an excellent fossil record, provides a well established morphological phylogeny. Sequence analysis of cox2 revealed a high degree of sequence divergence among the species, reflecting in part the antiquity of the family, but also a significant acceleration of sequence evolution in the ceratopogonids compared to other nematoceran Diptera. Phylogenetic reconstruction by neighbor-joining and maximum parsimony gave strong support for an early separation of an ancient lineage that includes the two genera, Austroconops and Leptoconops, from the remainder of the family. The results support the existence of a clade that includes two subfamilies, Dasyheleinae and Forcipomyiinae, and this clade appears as sister to the remaining subfamily, Ceratopogoninae. The molecular phylogeny also supports monophyly of the Ceratopogonidae, and either a sister or paraphyletic relationship of this family with the Chironomidae.

Phylogenetic and genomic analysis of the complete mitochondrial DNA sequence of the spotted asparagus beetle Crioceris duodecimpunctata.

We report the complete mitochondrial DNA sequence of the spotted asparagus beetle, Crioceris duodecimpunctata. The genome complement, gene order, and nucleotide composition of this beetle's mitochondrial genome were found to be typical of those reported for other insects. Unusual features of this genome include the substitution of UCU for GCU as the anticodon for tRNA(Ser), an unusual TpsiC loop for the tRNA(Ile) gene, and the identification of a putative ATT start codon for cox1. The utility of complete mitochondrial genome data for phylogenetic inference of the insect orders was tested, and compared to that of cox1 and combined mitochondrial ribosomal DNA sequences. Even though the number of insect orders represented by complete mitochondrial genomes is still limited, several well-established relationships are evident in the phylogenetic analysis of the complete sequences. Monophyly of the orders Diptera, Lepidoptera, and Coleoptera were consistently recovered. Monophyly of the Holometabola was also observed in some (though not all) analyses. The accumulation of complete mitochondrial sequences from a broader array of insect orders holds the promise of clarifying the early diversification of insects.

SELECTION AND THE "SEX-RATIO" POLYMORPHISM IN NATURAL POPULATIONS OF DROSOPHILA PSEUDOOBSCURA.

"Sex-ratio" (SR) is a naturally occurring X-linked meiotic drive system, where the SR-X chromosome is transmitted to nearly all progeny of SR males. It occurs at frequencies of up to 25% in some populations of Drosophila pseudoobscura. Because of the twofold drive advantage, SR should rapidly fix in populations, causing the extinction of the species, unless opposed by strong selection. I examine several of the adult components of fitness, including the frequencies of all genotypic mating combinations, fertilities, and fecundities of flies from two populations in southeastern Arizona. Significant reduction of fecundity of SR/SR females was observed in the Tucson population. No evidence was found for either lower fertility or reduced mating success of SR males, relative to standard males. Most selection opposing SR appears to be operating at the larval stages in nature.