Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

The rate and pattern of nucleotide substitution in Drosophila mitochondrial DNA.

The nucleotide sequences of a segment of mitochondrial DNA (mtDNA) have been determined for nine species or subspecies of the subgenus Drosophila of the genus Drosophila. This segment contains two complete protein-coding genes (i.e., NADH dehydrogenase subunit 1 and cytochrome b) and a transfer RNA gene (tRNA(ser)). The G+C content at third-codon positions for the two protein-coding genes was 1.5 times higher than that in the D. melanogaster species group, which belongs to the subgenus Sophophora. However, there was a substantial difference between the nucleotide frequencies of G and C. The number of nucleotide substitutions per silent site was more than three times higher than that for nuclear DNA, although it was only 60% of that for mammalian mtDNA. Both parametric and nonparametric analyses revealed a strong transition-transversion bias in nucleotide substitution, as was observed in mammalian mtDNA. Moreover, the rate of substitution of A and T for G and C is higher than that for the opposite direction. This bias seems to be responsible for the extremely A+T-rich base composition of Drosophila mtDNA. It is also noted that the rate of transitional change between A and G is higher than that between T and C.     

Phylogenetic utility of mitochondrial COI and nuclear Gpdh genes in Drosophila

Phylogenetic utility of the mitochondrial COI (cytochrome oxidase subunit I) and nuclear Gpdh (glycerol-3-phosphate dehydrogenase) genes was studied in the Drosophila melanogaster species group. The rate of substitution was higher in the COI gene than in the Gpdh gene. In addition, multiple substitutions, not only for transitional but also for transversional substitutions, occurred faster in the COI gene. None of the trees obtained using the COI gene supported the well-established monophyly of the ananassae subgroup. In addition, the incongruence length difference test, Templeton test, and partitioned Bremer support revealed that the trees based on the COI data are considerably different from those based on the Gpdh and the combined data set. Thus, the COI gene did not show good phylogenetic performance in the melanogaster group. The present analyses based on the Gpdh gene and the combined data set revealed that the ananassae subgroup branched off first in the melanogaster group followed by the montium subgroup and further by the melanogaster subgroup in contrast to the most recent phylogenetic hypothesis based on Amy multigenes.     

Drosophila Melanogaster Mitochondrial DNA: Gene Organization and Evolutionary Considerations

The sequence of a 8351-nucleotide mitochondrial DNA (mtDNA) fragment has been obtained extending the knowledge of the Drosophila melanogaster mitochondrial genome to 90% of its coding region. The sequence encodes seven polypeptides, 12 tRNAs and the 3' end of the 16S rRNA and CO III genes. The gene organization is strictly conserved with respect to the Drosophila yakuba mitochondrial genome, and different from that found in mammals and Xenopus. The high A + T content of D. melanogaster mitochondrial DNA is reflected in a reiterative codon usage, with more than 90% of the codons ending in T or A, G + C rich codons being practically absent. The average level of homology between the D. melanogaster and D. yakuba sequences is very high (roughly 94%), although insertion and deletions have been detected in protein, tRNA and large ribosomal genes. The analysis of nucleotide changes reveals a similar frequency for transitions and transversions, and reflects a strong bias against G + C on both strands. The predominant type of transition is strand specific.     

The complete nucleotide sequence and gene organization of the mitochondrial genome of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae)

The complete nucleotide sequences of the mitochondrial genome (mitogenome) of the oriental mole cricket, Gryllotalpa orientalis (Orthoptera: Gryllotalpidae), were determined. The 15,521-bp-long G. orientalis mitogenome contains typical gene complement, base composition, and codon usage found in metazoan mitogenomes. The G. orientalis mitogenome contains the third lowest A+T content (70.5%) among the complete insects mt genome sequences. The initiation codon for the G. orientalis COI gene appears to be ATG, instead of the tetranucleotides, which have been postulated to act as initiation codon for Locusta migratoria and some lepidopteran COI genes. The initiation codon for ND2 appears to be GTG, which is rare, but has been designated as an initiator of Tricholepidion gertschi ND2. All anticodons of G. orientalis tRNAs were identical to Drosophila yakuba and L. migratoria. The tRNA(Ser)(AGN) could not form a stable stem loop structure in the DHU arm as shown in many other insect tRNA(Ser)(AGN). Phylogenetic analysis of nucleotide sequence information from all mt genes supported a monophyletic Diptera, a monophyletic Lepidoptera, a monophyletic Coleoptera, a monophyletic Mecopterida (Diptera+Lepidoptera), and a monophyletic Endopterygota (Diptera+Lepidoptera+Coleoptera), suggesting that the complete insect mitogenome sequence has a resolving power to the diversification events within Endopterygota. However, the relationships of ancient insect orders were unstable, indicating the limited use of mitogenome information at deeper phylogenetic depth.     

Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

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.      

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.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.     

The initiator tRNA genes of Drosophila melanogaster: evidence for a tRNA pseudogene

We have isolated four segments of Drosophila melanogaster DNA that hybridize to homologous initiator tRNAMet. Three of the cloned fragments contain initiator tRNA genes, each of which can be transcribed in vitro. The fourth clone, pPW568, contains an initiator tRNA pseudogene which is not transcribed in vitro by RNA polymerase III. The pseudogene is contained in a 1.15 kb DNA fragment. This fragment has the characteristics of dispersed repetitive DNA and hybridizes in situ to at least 30 sites in the Drosophila genome. The arrangement of the initiator tRNA genes we have isolated, is different to that of other Drosophila tRNA gene families. The initiator tRNA genes are not clustered nor intermingled with other tRNA genes. They occur as single copies within an approximately 415-bp repeat segment, which is separated from other initiator tRNA genes by a mean distance of 17 kb. In situ hybridization to polytene chromosomes localizes these genes to the 61D region of the Drosophila genome. Hybridization analysis of genomic DNA indicates the presence of 8-9 non-allelic initiator tRNA genes in Drosophila melanogaster.     

"Silent" sites in Drosophila genes are not neutral: evidence of selection among synonymous codons

The patterns of synonymous codon usage in 91 Drosophila melanogaster genes have been examined. Codon usage varies strikingly among genes. This variation is associated with differences in G+C content at silent sites, but (unlike the situation in mammalian genes) these differences are not correlated with variation in intron base composition and so are not easily explicable in terms of mutational biases. Instead, those genes with high G+C content at silent sites, resulting from a strong "preference" for a particular subset of the codons that are mostly C- ending, appear to be the more highly expressed genes. This suggests that G+C content is reduced in sequences where selective constraints are weaker, as indeed seen in a pseudogene. These and other data discussed are consistent with the effects of translational selection among synonymous codons, as seen in unicellular organisms. The existence of selective constraints on silent substitutions, which may vary in strength among genes, has implications for the use of silent molecular clocks.      

Ixodes philipi (Acari: Ixodidae): phylogenetic status inferred from mitochondrial cytochrome oxidase subunit I gene sequence comparison

Ixodes philipi ticks were collected from the nest burrows of streaked shearwaters, Calonectris luecomelas, on 3 different islands of Japan (Awashima: 38 degrees 45'N, 139 degrees 24'E; Mikurajima: 33 degrees 52'N, 139 degrees 36'E; and Omorijima: 36 degrees 8'N, 133 degrees 10'E). The mitochondrial cytochrome oxidase subunit I (COI) gene sequence was determined for each tick. The COI sequences of 9 other ixodid tick species also were determined, and they were used for taxonomic positioning of I. philipi. A metastriata tick, Amblyomma triguttatum, was used as an outgroup reference for the analysis. Phylogenetic examination indicated that the I. philipi ticks are on the branch with Ixodes turdus and Ixodes acutitarsus weakly, and the bootstrap value of this branching was low. Three different analyses, maximum parsimony, genetic distance, and maximum likelihood, support this conclusion. To further refine this analysis, 2761 base pairs (bp) of sequence, which included the genes for tRNA(Met), NADH dehydrogenase subunit 2 (ND2), tRNA(Trp), tRNA(Cys), tRNA(Tyr), and COI, were determined and compared for 6 I. philipi ticks from the 3 different collection sites. Although a base substitution (T to C in the ND2 gene for an Awashima tick) and 2 transitions (G to A in the COI gene for 1 Omorijima tick) have occurred, the overall sequences were highly conserved. Preserved mitochondrial sequences in the ticks from 3 widely separated locations suggest the possibility of gene flow, which was probably accomplished by migratory seabirds.     

Variable strength of translational selection among 12 Drosophila species

Codon usage bias in Drosophila melanogaster genes has been attributed to negative selection of those codons whose cellular tRNA abundance restricts rates of mRNA translation. Previous studies, which involved limited numbers of genes, can now be compared against analyses of the entire gene complements of 12 Drosophila species whose genome sequences have become available. Using large numbers (6138) of orthologs represented in all 12 species, we establish that the codon preferences of more closely related species are better correlated. Differences between codon usage biases are attributed, in part, to changes in mutational biases. These biases are apparent from the strong correlation (r = 0.92, P < 0.001) among these genomes' intronic G + C contents and exonic G + C contents at degenerate third codon positions. To perform a cross-species comparison of selection on codon usage, while accounting for changes in mutational biases, we calibrated each genome in turn using the codon usage bias indices of highly expressed ribosomal protein genes. The strength of translational selection was predicted to have varied between species largely according to their phylogeny, with the D. melanogaster group species exhibiting the strongest degree of selection.     

线粒体ND4-ND4L基因在黑腹果蝇种组中的进化特征

本实验对黑腹果蝇种组(melanogaster species group)中8个种亚组33个样品两个线粒体基因ND4和ND4L进行了测序,并分析了ND4基因的序列差异和碱基替换特点,发现近缘物种中存在很明显的转换倾向,而在远缘物种中由于重复替换导致转换数处于饱和状态,我们的实验数据证实了线粒体基因较核基因有较快的进化速度。最后根据D．melanogaster与D．yakuba的遗传距离推算了8个种亚组的分化时间,ananassae种亚组最先分化,然后依次是montium,melanogaster,ficsphila,eugracilis,elegans,suzukii和takahashii最后分化。     

Rapid evolution and gene-specific patterns of selection for three genes of spermatogenesis in Drosophila

Hybrid males resulting from crosses between closely related species of Drosophila are sterile. The F1 hybrid sterility phenotype is mainly due to defects occurring during late stages of development that relate to sperm individualization, and so genes controlling sperm development may have been subjected to selective diversification between species. It is also possible that genes of spermatogenesis experience selective constraints given their role in a developmental pathway. We analyzed the molecular evolution of three genes playing a role during the sperm developmental pathway in Drosophila at an early (bam), a mid (aly), and a late (dj) stage. The complete coding region of these genes was sequenced in different strains of Drosophila melanogaster and Drosophila simulans. All three genes showed rapid divergence between species, with larger numbers of nonsynonymous to synonymous differences between species than polymorphisms. Although this could be interpreted as evidence for positive selection at all three genes, formal tests of selection do not support such a conclusion. Departures from neutrality were detected only for dj and bam but not aly. The role played by selection is unique and determined by gene-specific characteristics rather than site of expression. In dj, the departure was due to a high proportion of neutral synonymous polymorphisms in D. simulans, and there was evidence of purifying selection maintaining a high lysine amino acid protein content that is characteristic of other DNA-binding proteins. The earliest spermatogenesis gene surveyed, which plays a role in both male and female gametogenesis, was bam, and its significant departure from neutrality was due to an excess of nonsynonymous substitutions between species. Bam is degraded at the end of mitosis, and rapid evolutionary changes among species might be a characteristic shared with other degradable transient proteins. However, the large number of nonsynonymous changes between D. melanogaster and D. simulans and a phylogenetic comparative analysis among species confirms evidence of positive selection driving the evolution of Bam and suggests an yet unknown germ cell line developmental adaptive change between these two species.     

Drosophila mitochondrial DNA: Conserved sequences in the A+T-rich region and supporting evidence for a secondary structure model of the small ribosomal RNA

Summary The sequence of a segment of theDrosophila virilis mitochondrial DNA (mtDNA) molecule that contains the A+T-rich region, the small rRNA gene, the tRNA^f-met, tRNA^gln, and tRNA^ile genes, and portions of the ND2 and tRNA^val genes is presented and compared with the corresponding segment of theD. yakuba mtDNA molecule. The A+T-rich regions ofD. virilis andD. yakuba contain two correspondingly located sequences of 49 and 276/274 nucleotides that appear to have been conserved during evolution. In each species the replication origin of the mtDNA molecule is calculated to lie within a region that overlaps the larger conserved sequence, and within this overlap is found a potential hairpin structure. Substitutions between the larger conserved sequences of the A+T-rich regions, the small mt-rRNA genes, and the ND2 genes are biased in favor of transversions, 71–97% of which are AT changes. There is a 13.8 times higher frequency of nucleotide differences between the 5 halves than between the 3 halves of theD. virilis andD. yakuba small mt-rRNA genes. Considerations of the effects of observed substitutions and deletion/insertions on possible nucleotide pairing within the small mt-rRNA genes ofD. virilis andD. yakuba strongly support the secondary structure model for theDrosophila small mt-rRNA that we previously proposed.     

Conservatism of sites of tRNA loci among the linkage groups of severalDrosophila species

Summary The sites of seven tRNA genes (Arg-2, Lys-2, Ser-2b, Ser-7, Thr-3, Thr-4, Val-3b) were studied by in situ hybridization.¹²⁵I-labeled tRNA probes fromDrosophila melanogaster were hybridized to spreads of polytene chromosomes prepared from fourDrosophila species representing different evolutionary lineages (D. melanogaster, Drosophila hydei, Drosophila pseudoobscura, andDrosophila virilis). Most tRNA loci occurred on homologous chromosomal elements of all four species. In some cases the number of hybridization sites within an element varied and sites on nonhomologous elements were found. It was observed that both tRNA ₂ ^Arg and tRNA ₂ ^Lys hybridized to the same site on homologous elements in several species. These data suggest a limited amount of exchange among different linkage groups during the evolution ofDrosophila species.