桔小实蝇线粒体基因组全序列及其分析

桔小实蝇Bactrocera dorsalis线粒体基因组全序列对研究实蝇分子系统进化具有重要意义。本研究通过DNA测序和克隆技术,对桔小实蝇mtDNA全序列进行了测定和分析。结果表明：桔小实蝇线粒体基因组全长15 915 bp（GenBank序列号: DQ845759）。基因组碱基组成为39.3%A,16.2%C,10.2%G,34.3%T,由13个蛋白编码基因、22个tRNA基因、2个rRNA基因以及一个非编码的控制区域（A+T-rich区）组成。7个蛋白编码基因和13个tRNA基因从J链编码,其余6个蛋白编码基因和9个tRNA基因从N链编码。位于J链上的蛋白编码基因具有近似的A、T含量,而位于N链上的蛋白编码基因的A的含量明显高于T的含量。以mtDNA COⅠ基因为例,比较了桔小实蝇与其他14种实蝇的亲缘关系,结果显示其与同亚属（果实蝇亚属Bactrocera）内的其他近缘种相互间的同源性很高。     

Seventeen new complete mtDNA sequences reveal extensive mitochondrial genome evolution within the Demospongiae

Two major transitions in animal evolution--the origins of multicellularity and bilaterality--correlate with major changes in mitochondrial DNA (mtDNA) organization. Demosponges, the largest class in the phylum Porifera, underwent only the first of these transitions and their mitochondrial genomes display a peculiar combination of ancestral and animal-specific features. To get an insight into the evolution of mitochondrial genomes within the Demospongiae, we determined 17 new mtDNA sequences from this group and analyzing them with five previously published sequences. Our analysis revealed that all demosponge mtDNAs are 16- to 25-kbp circular molecules, containing 13-15 protein genes, 2 rRNA genes, and 2-27 tRNA genes. All but four pairs of sampled genomes had unique gene orders, with the number of shared gene boundaries ranging from 1 to 41. Although most demosponge species displayed low rates of mitochondrial sequence evolution, a significant acceleration in evolutionary rates occurred in the G1 group (orders Dendroceratida, Dictyoceratida, and Verticillitida). Large variation in mtDNA organization was also observed within the G0 group (order Homosclerophorida) including gene rearrangements, loss of tRNA genes, and the presence of two introns in Plakortis angulospiculatus. While introns are rare in modern-day demosponge mtDNA, we inferred that at least one intron was present in cox1 of the common ancestor of all demosponges. Our study uncovered an extensive mitochondrial genomic diversity within the Demospongiae. Although all sampled mitochondrial genomes retained some ancestral features, including a minimally modified genetic code, conserved structures of tRNA genes, and presence of multiple non-coding regions, they vary considerably in their size, gene content, gene order, and the rates of sequence evolution. Some of the changes in demosponge mtDNA, such as the loss of tRNA genes and the appearance of hairpin-containing repetitive elements, occurred in parallel in several lineages and suggest general trends in demosponge mtDNA evolution.     

Mitochondrial tRNA cleavage by tRNA-targeting ribonuclease causes mitochondrial dysfunction observed in mitochondrial disease

Mitochondrial DNA (mtDNA) is a genome possessed by mitochondria. Since reactive oxygen species (ROS) are generated during aerobic respiration in mitochondria, mtDNA is commonly exposed to the risk of DNA damage. Mitochondrial disease is caused by mitochondrial dysfunction, and mutations or deletions on mitochondrial tRNA (mt tRNA) genes are often observed in mtDNA of patients with the disease. Hence, the correlation between mt tRNA activity and mitochondrial dysfunction has been assessed. Then, cybrid cells, which are constructed by the fusion of an enucleated cell harboring altered mtDNA with a ρ⁰ cell, have long been used for the analysis due to difficulty in mtDNA manipulation. Here, we propose a new method that involves mt tRNA cleavage by a bacterial tRNA-specific ribonuclease. The ribonuclease tagged with a mitochondrial-targeting sequence (MTS) was successfully translocated to the mitochondrial matrix. Additionally, mt tRNA cleavage, which resulted in the decrease of cytochrome c oxidase (COX) activity, was observed.     

Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer

The mitochondrial genome of grape (Vitis vinifera), the largestorganelle genome sequenced so far, is presented. The genomeis 773,279 nt long and has the highest coding capacity amongknown angiosperm mitochondrial DNAs (mtDNAs). The proportionof promiscuous DNA of plastid origin in the genome is also thelargest ever reported for an angiosperm mtDNA, both in absoluteand relative terms. In all, 42.4% of chloroplast genome of Vitishas been incorporated into its mitochondrial genome. In orderto test if horizontal gene transfer (HGT) has also contributedto the gene content of the grape mtDNA, we built phylogenetictrees with the coding sequences of mitochondrial genes of grapeand their homologs from plant mitochondrial genomes. Many incongruentgene tree topologies were obtained. However, the extent of incongruencebetween these gene trees is not significantly greater than thatobserved among optimal trees for chloroplast genes, the commonancestry of which has never been in doubt. In both cases, weattribute this incongruence to artifacts of tree reconstruction,insufficient numbers of characters, and gene paralogy. Thisfinding leads us to question the recent phylogenetic interpretationof Bergthorsson et al. (2003, 2004) and Richardson and Palmer(2007) that rampant HGT into the mtDNA of Amborella best explainsphylogenetic incongruence between mitochondrial gene trees forangiosperms. The only evidence for HGT into the Vitis mtDNAfound involves fragments of two coding sequences stemming fromtwo closteroviruses that cause the leaf roll disease of thisplant. We also report that analysis of sequences shared by bothchloroplast and mitochondrial genomes provides evidence fora previously unknown gene transfer route from the mitochondrionto the chloroplast.     

Mitochondrial DNA sequences of the Afro-Arabian spiny-tailed lizards (genus Uromastyx; family Agamidae): phylogenetic analyses and evolution of gene arrangements

Approximately 1.7 kbp of mitochondrial DNA were sequenced from 29 individuals assignable to 11 Uromastyx species or subspecies and two other agamids. U. ocellata and U. ornata had an insertion between the glutamine and isoleucine tRNA genes, which could be folded into a stable stem-and-loop structure, and the insertion for U. ornata additionally retained a sequence similar to the glutamine tRNA gene. This corroborates the role of tandem duplication in reshaping mitochondrial gene arrangements, and supports the idea that the origin of light-strand replication could be relocated within mitochondrial genomes. Molecular phylogeny from different tree-building methods consistently placed African and Arabian taxa in mutually monophyletic groups, excluding U. hardwickii inhabiting India and Pakistan. Unlike previous studies based on morphology , U. macfadyeni did not cluster with morphologically similar Arabian taxa, suggesting convergent evolution to be responsible for the morphological similarities. Divergence times estimated among the Uromastyx taxa, together with geological and palaeontological evidence, suggest that the Uromastyx agamids originated from Central Asia during the Eocene and colonized Africa after its connection with Eurasia in the early Miocene. Their radiation may have been facilitated by repeated aridification of North Africa since the middle Miocene, and geological events such as the expansion of the Red Sea and the East African Rift Valley.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 247–260.     

Novel Repetitive Structures, Deviant Protein-Encoding Sequences andUnidentified ORFs in the Mitochondrial Genome of the BrachiopodLingula anatina

Complete sequence determination of the brachiopod Lingula anatina mtDNA (28,818 bp) revealed an organization that is remarkably atypical for an animal mt-genome. In addition to the usual set of 37 animal mitochondrial genes, which make up only 57% (16,555 bp) of the entire sequence, the genome contains lengthy unassigned sequences. All the genes are encoded in the same DNA strand, generally in a compact way, whereas the overall gene order is highly divergent in comparison with known animal mtDNA. Individual genes are generally longer and deviate considerably in sequence from their homologues in other animals. The genome contains two major repeat regions, in which 11 units of unassigned sequences and six genes (atp8, trnM, trnQ, trnV, and part of cox2 and nad2) are found in repetition, in the form of nested direct repeats of unparalleled complexity. One of the repeat regions contains unassigned repeat units dispersed among several unique sequences, novel repetitive structure for animal mtDNAs. Each of those unique sequences contains an open reading frame for a polypeptide between 80 and 357 amino acids long, potentially encoding a functional molecule, but none of them has been identified with known proteins. In both repeat regions, tRNA genes or tRNA gene-like sequences flank major repeated units, supporting the view that those structures play a role in the mitochondrial gene rearrangements. Although the intricate repeated organization of this genome can be explained by recurrent tandem duplications and subsequent deletions mediated by replication errors, other mechanisms, such as nonhomologous recombinations, appear to explain certain structures more easily.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

The complete mitochondrial genome of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae: gene order and trn gene clusters reveal a common evolutionary course for all Sordariomycetes, while intergenic regions show variation

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae, with a total size of 24,673 bp, was one of the smallest known mtDNAs of Pezizomycotina. It contained the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, a single intron that harbored an intronic ORF coding for a putative ribosomal protein (rps) within the large rRNA gene (rnl), and a set of 24 tRNA genes which recognized codons for all amino acids, except proline and valine. Gene order comparison with all known mtDNAs of Sordariomycetes illustrated a highly conserved genome organization for all the protein- and rRNA-coding genes, as well as three clusters of tRNA genes. By considering all mitochondrial essential protein-coding genes as one unit a phylogenetic study of these small genomes strongly supported the common evolutionary course of Sordariomycetes (100% bootstrap support) and highlighted the advantages of analyzing small genomes (mtDNA) over single genes. In addition, comparative analysis of three intergenic regions demonstrated sequence variability that can be exploited for intra- and inter-specific identification of Metarhizium. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Nucleotide Sequence and Gene Organization of the Starfish Asterina Pectinifera Mitochondrial Genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCU)(Ser), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.     

Novel features of metazoan mtDNA revealed from sequence analysis of three mitochondrial DNA segments of the land snail Albinaria turrita (Gastropoda: Clausiliidae)

The mitochondrial DNA (mtDNA) size of the terrestrial gastropod Albinaria turrita was determined by restriction enzyme mapping and found to be approximately 14.5 kb. Its partial gene content and organization were examined by sequencing three cloned segments representing about one-fourth of the mtDNA molecule. Complete sequences of cytochrome c oxidase subunit II (COII), and ATPase subunit 8 (ATPase8), as well as partial sequences of cytochrome c oxidase subunit I (COI), NADH dehydrogenase subunit 6 (ND6), and the large ribosomal RNA (IrRNA) genes were determined. Nine putative tRNA genes were also identified by their ability to conform to typical mitochondrial tRNA secondary structures. An 82-nt sequence resembles a noncoding region of the bivalve Mytilus edulis, even though it might contain a tenth tRNA gene with an unusual 5-nt overlap with another tRNA gene. The genetic code of Albinaria turrita appears to be the same as that of Drosophila and Mytilus edulis. The structures of COI and COII are conservative, but those of ATPase8 and ND6 are diversified. The sequenced portion of thelrRNA gene (1,079 nt) is characterized by conspicuous deletions in the 5 and 3 ends; this gene represents the smallest coelomate IrRNA gene so far known. Sequence comparisons of the identified genes indicate that there is greater difference between Albinaria and Mytilus than between Albinaria and Drosophila. An evolutionary analysis, based on COII sequences, suggests a possible nonmonophyletic origin of molluskan mtDNA. This is supported also by the absence of the ATPase8 gene in the mtDNA of Mytilus and nematodes, while this gene is present in the mtDNA of Albinaria and Cepaea nemoralis and in all other known coelomate metazoan mtDNAs.     

Complete DNA sequence of the mitochondrial genome ofCepaea nemoralis (Gastropoda: Pulmonata)

The nucleotide sequence of a mitochondrial genome of the pulmonate gastropod molluscCepaea nemoralis has been determined. Contained within the 14,100 basepairs (bp) are the two ribosomal RNA genes and 13 protein coding genes typical of metazoan mitochondrial genomes. TheCepaea mtDNA does contain a gene for ATPase subunit 8, like the clausiliid pulmonate,Albinaria, and the chiton,Katharina, but unlike the bivalve mollusc,Mytilus. The mitochondrial genetic code ofCepaea is proposed to be the same as that ofMytilus, Katharina, andDrosophila. Only 14 putative tRNA genes are presented, although there is sufficient unassigned sequence to encode the remainder of the expected total of 22 tRNA genes. These 14 tRNA genes are a mixture of standard cloverleaf structures and nonstandard structures containing TV replacement loops as seen in nematode and mosquito mitochondrial genomes. If the eight unidentified tRNA genes are indeed present, very little unassigned sequence would remain to serve as a control region. Genes are transcribed from both strands of the molecule. Base composition is the least biased for any reported animal mitochondrial genome and is also very little skewed between strands using measures independent of base composition. TheCepaea mitochondrial gene order is quite unlike that of any other reported metazoan mtDNA, with the exception of the recently reported partial sequences ofAlbinaria. No gene bound-aries are shared among all the reported molluscan taxa, demonstrating a complete lack of conservation of mitochondrial gene order across the phylum Mollusca.     

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region- like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region- like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.      

Complete Mitochondrial DNA Sequences of Six Snakes: Phylogenetic Relationships and Molecular Evolution of Genomic Features

Complete mitochondrial DNA (mtDNA) sequences were determined for representative species from six snake families: the acrochordid little file snake, the bold boa constrictor, the cylindrophiid red pipe snake, the viperid himehabu, the pythonid ball python, and the xenopeltid sunbeam snake. Thirteen protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions were identified in these mtDNAs. Duplication of the control region and translocation of the tRNA^Leu gene were two notable features of the snake mtDNAs. The duplicate control regions had nearly identical nucleotide sequences within species but they were divergent among species, suggesting concerted sequence evolution of the two control regions. In addition, the duplicate control regions appear to have facilitated an interchange of some flanking tRNA genes in the viperid lineage. Phylogenetic analyses were conducted using a large number of sites (9570 sites in total) derived from the complete mtDNA sequences. Our data strongly suggested a new phylogenetic relationship among the major families of snakes: ((((Viperidae, Colubridae), Acrochordidae), (((Pythonidae, Xenopeltidae), Cylindrophiidae), Boidae)), Leptotyphlopidae). This conclusion was distinct from a widely accepted view based on morphological characters in denying the sister-group relationship of boids and pythonids, as well as the basal divergence of nonmacrostomatan cylindrophiids. These results imply the significance to reconstruct the snake phylogeny with ample molecular data, such as those from complete mtDNA sequences.[Reviewing Editor: Dr. Bill Ballard]     

Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics

Mitochondrial DNA sequences are often used to construct molecular phylogenetic trees among closely related animals. In order to examine the usefulness of mtDNA sequences for deep-branch phylogenetics, genes in previously reported mtDNA sequences were analyzed among several animals that diverged 20–600 million years ago. Unambiguous alignment was achieved for stem-forming regions of mitochondrial tRNA genes by virtue of their conservative secondary structures. Sequences derived from stem parts of the mitochondrial tRNA genes appeared to accumulate much variation linearly for a long period of time: nearly 100 Myr for transition differences and more than 350 Myr for transversion differences. This characteristic could be attributed, in part, to the structural variability of mitochondrial tRNAs, which have fewer restrictions on their tertiary structure than do nonmitochondrial tRNAs. The tRNA sequence data served to reconstruct a well-established phylogeny of the animals with 100% bootstrap probabilities by both maximum parsimony and neighbor joining methods. By contrast, mitochondrial protein genes coding for cytochrome b and cytochrome oxidase subunit I did not reconstruct the established phylogeny or did so only weakly, although a variety of fractions of the protein gene sequences were subjected to tree-building. This discouraging phylogenetic performance of mitochondrial protein genes, especially with respect to branches originating over 300 Myr ago, was not simply due to high randomness in the data. It may have been due to the relative susceptibility of the protein genes to natural selection as compared with the stem parts of mitochondrial tRNA genes. On the basis of these results, it is proposed that mitochondrial tRNA genes may be useful in resolving deep branches in animal phylogenies with divergences that occurred some hundreds of Myr ago. For this purpose, we designed a set of primers with which mtDNA fragments encompassing clustered tRNA genes were successfully amplified from various vertebrates by the polymerase chain reaction.Abbreviations AA stem amino acid-acceptor stem - AC stem anticodon stem - COI cytochrome oxidase subunit I - cytb cytochrome b - D stem dihydrouridine stem - MP maximum parsimony - mtDNA mitochondrial DNA - Myr million years - NJ neighbor joining - PCR polymerase chain reaction - Ti transition - T stem tC stem - Tv transversion Correspondence to: Y. Kumazawa     

Complete sequence of the mitochondrial genome of Tetrahymena pyriformis and comparison with Paramecium aurelia mitochondrial DNA

We report the complete nucleotide sequence of the Tetrahymena pyriformis mitochondrial genome and a comparison of its gene content and organization with that of Paramecium aurelia mtDNA. T. pyriformis mtDNA is a linear molecule of 47,172 bp (78.7 % A+T) excluding telomeric sequences (identical tandem repeats of 31 bp at each end of the genome). In addition to genes encoding the previously described bipartite small and large subunit rRNAs, the T. pyriformis mitochondrial genome contains 21 protein-coding genes that are clearly homologous to genes of defined function in other mtDNAs, including one (yejR) that specifies a component of a cytochrome c biogenesis pathway. As well, T. pyriformis mtDNA contains 22 open reading frames of unknown function larger than 60 codons, potentially specifying proteins ranging in size from 74 to 1386 amino acid residues. A total of 13 of these open reading frames ("ciliate-specific") are found in P. aurelia mtDNA, whereas the remaining nine appear to be unique to T. pyriformis; however, of the latter, five are positionally equivalent and of similar size in the two ciliate mitochondrial genomes, suggesting they may also be homologous, even though this is not evident from sequence comparisons. Only eight tRNA genes encoding seven distinct tRNAs are found in T. pyriformis mtDNA, formally confirming a long-standing proposal that most T. pyriformis mitochondrial tRNAs are nucleus-encoded species imported from the cytosol. Atypical features of mitochondrial gene organization and expression in T. pyriformis mtDNA include split and rearranged large subunit rRNA genes, as well as a split nad1 gene (encoding subunit 1 of NADH dehydrogenase of respiratory complex I) whose two segments are located on and transcribed from opposite strands, as is also the case in P. aurelia. Gene content and arrangement are very similar in T. pyriformis and P. aurelia mtDNAs, the two differing by a limited number of duplication, inversion and rearrangement events. Phylogenetic analyses using concatenated sequences of several mtDNA-encoded proteins provide high bootstrap support for the monophyly of alveolates (ciliates, dinoflagellates and apicomplexans) and slime molds.     

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

Nearly complete mitochondrial genome of Polyascus gregaria and the phylogenetic relationships among maxillopodans

We determined for the first time the nearly complete mitochondrial genome sequence of the entozoic Polyascus gregaria, a representative of Rhizocephala, Cirripedia. The nearly complete mitogenome was 15, 465 bp in length, consisting of 11 protein-coding genes, two rRNA genes, 22 tRNA genes and one major incomplete noncoding region. In total there are 73 overlapping nucleotides and 17 spacers between genes. All genes sequenced in P. gregaria mtDNA (including RNAs) were encoded on the same strand of the DNA, and the gene arrangement differed from that of other metazoan animals. The mitochondrial genome rearrangements included translocation of at least 8 genes and even inversion of the coding polarity of at least 2 genes. Comparative analysis of the gene orders with other maxillopodan mtDNAs showed that the unique characteristics of the thoracican cirripeds lineage were not observed in this representative of rhizocephalan. Phylogenetic analyses supported a close affinity of Rhizocephala to Thoracica. By adding the mitochondrial genomes from 4 copepods, the reciprocally monophyletic cirripeds and copepods clustered as sister groups, refusing the close relationship between Cirripedia and Remipedia. However, the monophyly of Maxillopoda was not supported in this study.     

Glass sponges and bilaterian animals share derived mitochondrial genomic features: a common ancestry or parallel evolution?

Glass sponges (Hexactinellida) are a group of deep-water benthicanimals that have a unique syncytial organization and possessa characteristic siliceous skeleton. Although hexactinellidsare traditionally grouped with calcareous and demosponges inthe phylum Porifera, the monophyly of sponges and the phylogeneticposition of the Hexactinellida remain contentious. We determinedand analyzed the nearly complete mitochondrial genome sequencesof the hexactinellid sponges Iphiteon panicea and Sympagellanux. Unexpectedly, our analysis revealed several mitochondrialgenomic features shared between glass sponges and bilateriananimals, including an Arg Ser change in the genetic code, acharacteristic secondary structure of one of the serine tRNAs,highly derived tRNA and rRNA genes, and the presence of a singlelarge noncoding region. At the same time, glass sponge mtDNAcontains atp9, a gene previously found only in the mtDNA ofdemosponges (among animals), and encodes a with an atypical A11–U24 pair that is alsofound in demosponges and placozoans. Most of our sequence-basedphylogenetic analyses place Hexactinellida as the sister groupto the Bilateria; however, these results are suspect given acceleratedrates of mitochondrial sequence evolution in these groups. Thus,it remains an open question whether shared mitochondrial genomicfeatures in glass sponges and bilaterian animals reflect theirclose phylogenetic affinity or provide a remarkable exampleof parallel evolution.