Anatomical features of Leiopelma embryos and larvae: implications for anuran evolution

A controversial issue in anuran systematics is the relationship of Leiopelma to other anurans because recent phylogenetic constructions imply different relationships among the basal frogs. Of particular evolutionary interest is whether early development of Leiopelma resembles an ancestral salamander-like larva, an anuran tadpole, or neither. In the 1950s, Neville G. Stephenson hypothesized that direct development is the primary mode of development in amphibians, based on the fact that Leiopelma spp. lack a free-living (=feeding) larval stage. Although this hypothesis has not been generally accepted, it has not been formally refuted. We review Stephenson's work on Leiopelma and examine the anatomy of embryos/"larvae" of the four extant Leiopelma species for evidence of vestigial larval features that might refute the "direct-developing ancestor" hypothesis. We describe internal oral features in early developmental stages of Leiopelma and compare Leiopelma with a closely related basal anuran, Ascaphus, to assess whether their early developmental stages share any derived features. In Leiopelma hochstetteri, embryos/larvae have open gill slits and some faint rugosities around one gill slit that may be vestiges of gill rakers or filters. They also have more intestinal loops, indicative of an elongated alimentary tract, at earlier rather than late embryonic/larval stages. Collectively, these features support the view that the ancestor of Leiopelma had a free-swimming, free-feeding, aquatic larva. The palatoquadrate of Leiopelma archeyi reorients approximately 40 degrees from a more horizontal to a more vertical position through embryonic/"larval" development. This amount of cranial remodeling is intermediate between that seen in salamanders (17-27 degrees) and that reported for Ascaphus (64 degrees ) and other basal frogs (71-78 degrees) at metamorphosis. We found no internal oral features that Leiopelma shares specifically with Ascaphus. However, Leiopelma embryos have a ventrally positioned mouth and a downturned rostrum, characteristic of Ascaphus and other stream-adapted tadpoles.     

Complete nucleotide sequence of the mitochondrial genome of a Malagasy poison frog Mantella madagascariensis: evolutionary implications on mitochondrial genomes of higher anuran groups

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of a Malagasy poison frog, Mantella madagascariensis (family Mantellidae), and partial sequences of two Mantella (M. baroni and M. bernhardi) and two additional mantellid species (Boophis madagascariensis and Mantidactylus cf. ulcerosus). The M. madagascariensis genome was shown to be the largest (23kbp) of all vertebrate mtDNAs investigated so far. Furthermore, the following unique features were revealed: (1) the positions of some genes and gene regions were rearranged compared to mitochondrial genomes typical for vertebrates and other anuran groups, (2) two distinct genes and a pseudogene corresponding to transfer RNA gene for methionine (tRNA-Met) were encoded, and (3) two control regions with very high sequence homology were present. These features were shared by the two other Mantella species but not the other mantellid species, indicating dynamic genome reorganization in a common ancestor linage before divergence of the Mantella genus. The reorganization pathway could be explained by a model of gene duplication and deletion. Duplication and deletion events also seem to have been responsible for concerted sequence evolution of the control regions in Mantella mt genomes. It is also suggested that the pseudo tRNA-Met gene sustained for a long time in Mantella mt genomes possibly functions as a punctuation marker for NADH dehydrogenase subunit (ND) 2 mRNA processing. Phylogenetic analyses employing a large sequence data set of mt genes supported the monophyly of Mantellidae and Rhacophoridae and other recent phylogenetic views for ranoid frogs. The resultant phylogenetic relationship also suggested parallel occurrence of two tRNA-Met genes, duplicated control regions, and ND5 gene translocation in independent ranoid lineages.     

Complete nucleotide sequence and gene arrangement of the mitochondrial genome of the crab-eating frog Fejervarya cancrivora and evolutionary implications

The complete nucleotide sequence of the mitochondrial (mt) genome of the crab-eating frog, Fejervarya cancrivora Gravenhorst (Amphibia: Anura: Ranidae), was determined. The mt genome is 17,843 bp long and contains 13 protein-coding (ATP6, ATP8, COI-III, ND1-6 and 4L, and Cyt b) and two ribosomal RNA (12S and 16SrRNA) genes. Although metazoan mt genomes typically encode 22 transfer RNA genes (tRNAs), the F. cancrivora mtDNA contains 23 tRNAs due to the presence of an extra copy of tRNA(Met). A major noncoding region and a prominent intergenic spacer corresponding to the control region and light-strand replication origin were also found. To confirm the phylogenetic position of F. cancrivora, we compared the gene arrangement with that of other anurans and performed phylogenetic analyses based on mt genomic data. The genome organization of F. cancrivora mtDNA differs from that of typical vertebrates and neobatrachian frogs but is identical with that of F. limnocharis, suggesting that the unique gene arrangement occurred in the common ancestor of the genus. Phylogenetic analyses supported the monophyly of the Fejervarya species used here as well as the dicroglossini clade. Although the family Ranidae as previously recognized (= Ranidae, Discoglossidae, and some other natatanuran families; sensu Frost et al., 2006) is shown as a clade in the maximum parsimony analysis, the maximum likelihood and the Bayesian analyses suggest the paraphyly of the Ranidae with respect to the families, Mantellidae and Rhacophoridae. Three-tandem duplications of gene regions followed by subsequent deletions of supernumerary genes were proposed to explain the evolution of the extra tRNA(Met) and translocation of ND5 from the original neobatrachian gene order.     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.     

Highly conserved gene arrangement of the mitochondrial genomes of 23 Plasmodium species

Mitochondrial (mt) genomes from diverse phylogenetic groups vary considerably in size, structure and organization. The genus Plasmodium, the causative agent of malaria, has the smallest mt genome in the form of a tandemly repeated, linear element of 6 kb. The Plasmodium mt genome encodes only three protein genes (cox1, cox3 and cob) and large- and small-subunit ribosomal RNA (rRNA) genes, which are highly fragmented with 19 identified rRNA pieces. The complete mt genome sequences of 21 Plasmodium species have been published but a thorough investigation of the arrangement of rRNA gene fragments has been undertaken for only Plasmodium falciparum, the human malaria parasite. In this study, we determined the arrangement of mt rRNA gene fragments in 23 Plasmodium species, including two newly determined mt genome sequences from P. gallinaceum and P. vinckei vinckei, as well as Leucocytozoon caulleryi, an outgroup of Plasmodium. Comparative analysis reveals complete conservation of the arrangement of rRNA gene fragments in the mt genomes of all the 23 Plasmodium species and L. caulleryi. Surveys for a new rRNA gene fragment using hidden Markov models enriched with recent mt genome sequences led us to suggest the mtR-26 sequence as a novel candidate LSU rRNA fragment in the mt genomes of the 24 species. Additionally, we found 22-25 bp-inverted repeat sequences, which may be involved in the generation of lineage-specific mt genome arrangements after divergence from a common ancestor of the genera Eimeria and Plasmodium/Leucocytozoon.     

节肢动物线粒体基因组与系统发生重建

对mt基因组的比较研究是探讨节肢动物系统发生的有效手段之一。基因的排列和DNA序列可以为重建节肢动物的系统发生提供有用的信息。目前,已测定mt基因组全序列的节肢动物已增加到44种。归纳、总结了节肢动物mt基因组的基本特征、基因顺序、基因重排的发生和机制等。简要评述基于mt基因组的节肢动物系统发生研究。     

The Complete Mitochondrial Genome of Galba pervia (Gastropoda: Mollusca), an Intermediate Host Snail of Fasciola spp

Complete mitochondrial (mt) genomes and the gene rearrangements are increasingly used as molecular markers for investigating phylogenetic relationships. Contributing to the complete mt genomes of Gastropoda, especially Pulmonata, we determined the mt genome of the freshwater snail Galba pervia, which is an important intermediate host for Fasciola spp. in China. The complete mt genome of G. pervia is 13,768 bp in length. Its genome is circular, and consists of 37 genes, including 13 genes for proteins, 2 genes for rRNA, 22 genes for tRNA. The mt gene order of G. pervia showed novel arrangement (tRNA-His, tRNA-Gly and tRNA-Tyr change positions and directions) when compared with mt genomes of Pulmonata species sequenced to date, indicating divergence among different species within the Pulmonata. A total of 3655 amino acids were deduced to encode 13 protein genes. The most frequently used amino acid is Leu (15.05%), followed by Phe (11.24%), Ser (10.76%) and IIe (8.346%). Phylogenetic analyses using the concatenated amino acid sequences of the 13 protein-coding genes, with three different computational algorithms (maximum parsimony, maximum likelihood and Bayesian analysis), all revealed that the families Lymnaeidae and Planorbidae are closely related two snail families, consistent with previous classifications based on morphological and molecular studies. The complete mt genome sequence of G. pervia showed a novel gene arrangement and it represents the first sequenced high quality mt genome of the family Lymnaeidae. These novel mtDNA data provide additional genetic markers for studying the epidemiology, population genetics and phylogeographics of freshwater snails, as well as for understanding interplay between the intermediate snail hosts and the intra-mollusca stages of Fasciola spp..     

The mitochondrial genomes of Amphiascoides atopus and Schizopera knabeni (Harpacticoida: Miraciidae) reveal similarities between the copepod orders Harpacticoida and Poecilostomatoida

Members of subclass Copepoda are abundant, diverse, and—as a result of their variety of ecological roles in marine and freshwater environments—important, but their phylogenetic interrelationships are unclear. Recent studies of arthropods have used gene arrangements in the mitochondrial (mt) genome to infer phylogenies, but for copepods, only seven complete mt genomes have been published. These data revealed several within-order and few among-order similarities. To increase the data available for comparisons, we sequenced the complete mt genome (13,831 base pairs) of Amphiascoides atopus and 10,649 base pairs of the mt genome of Schizopera knabeni (both in the family Miraciidae of the order Harpacticoida). Comparison of our data to those for Tigriopus japonicus (family Harpacticidae, order Harpacticoida) revealed similarities in gene arrangement among these three species that were consistent with those found within and among families of other copepod orders. Comparison of the mt genomes of our species with those known from other copepod orders revealed the arrangement of mt genes of our Harpacticoida species to be more similar to that of Sinergasilus polycolpus (order Poecilostomatoida) than to that of T. japonicus. The similarities between S. polycolpus and our species are the first to be noted across the boundaries of copepod orders and support the possibility that mt-gene arrangement might be used to infer copepod phylogenies. We also found that our two species had extremely truncated transfer RNAs and that gene overlaps occurred much more frequently than has been reported for other copepod mt genomes.     

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.     

Complete nucleotide sequence of the mitochondrial genome of Schlegel's tree frog Rhacophorus schlegelii (family Rhacophoridae): duplicated control regions and gene rearrangements

The complete nucleotide sequence (21,359 bp) of the mitochondrial DNA of the rhacophorid frog Rhacophorus schlegelii was determined. The gene content, nucleotide composition, and codon usage of this genome corresponded to those typical of vertebrates. However, the Rh. schlegelii genome was unusually large due to the inclusion of two control regions and the accumulation of lengthy repetitive sequences in these regions. The two control regions had 97% sequence similarity over 1,510 bp, suggesting the occurrence of concerted sequence evolution. Comparison of the gene organizations among anuran species revealed that the mitochondrial gene arrangement of Rh. schlegelii diverged from that of typical vertebrates but was similar to that of Buergeria buergeri. The positions of the tRNA-Leu(CUN) and tRNA-Thr genes were exchanged between Rh. schlegelii and B. buergeri. Based on parsimonious consideration and the basal phylogenetic position of B. buergeri, these genes seemed to have been rearranged in an ancestral lineage leading to Rh. schlegelii.     

A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.     

Mitochondrial Genome of <Emphasis Type="Italic">Ciona savignyi</Emphasis> (Urochordata,Ascidiacea, Enterogona): Comparison of Gene Arrangement and tRNA Genes with <Emphasis Type="Italic">Halocynthia roretzi</Emphasis> Mitochondrial Genome

The complete nucleotide sequence of the urochordate Ciona savignyi (Ascidiacea, Enterogona) mitochondrial (mt) genome (14,737 bp) was determined. The Ciona mt genome does not encode a gene for ATP synthetase subunit 8 but encodes an additional tRNA^Gly gene (anticodon UCU), as is the case in another urochordate, Halocynthia roretzi (Ascidiacea, Pleurogona), mt genome. In addition, the Ciona mt genome encodes two tRNA^Met genes; anticodon CAT and anticodon TAT. The tRNA^Cys gene is thought to lack base pairs at the D-stem. Thus, the Ciona mt genome encodes 12 protein, 2 rRNA, and 24 tRNA genes. The gene arrangement of the Ciona mt genome differs greatly from those of any other metazoan mt genomes reported to date. Only three gene boundaries are shared between the Halocynthia and the Ciona mt genomes. Molecular phylogenetic analyses based on amino acid sequences of mt protein genes failed to demonstrate the monophyly of the chordates.     

The mitochondrial genome of the brachiopod Laqueus rubellus

The complete nucleotide sequence of the 14,017-bp mitochondrial (mt) genome of the articulate brachiopod Laqueus rubellus is presented. Being one of the smallest of known mt genomes, it has an extremely compact gene organization. While the same 13 polypeptides, two rRNAs, and 22 tRNAs are encoded as in most other animal mtDNAs, lengthy noncoding regions are absent, with the longest apparent intergenic sequence being 54 bp in length. Gene-end sequence overlaps are prevalent, and several stop codons are abbreviated. The genes are generally shorter, and three of the protein-coding genes are the shortest among known homologues. All of the tRNA genes indicate size reduction in either or both of the putative TPsiC and DHU arms compared with standard tRNAs. Possession of a TV (TPsiC arm-variable loop) replacement loop is inferred for tRNA(R) and tRNA(L-tag). The DHU arm appears to be unpaired not only in tRNA(S-tct) and tRNA(S-tga), but also in tRNA(C), tRNA(I), and tRNA(T), a novel condition. All the genes are encoded in the same DNA strand, which has a base composition rich in thymine and guanine. The genome has an overall gene arrangement drastically different from that of any other organisms so far reported, but contains several short segments, composed of 2-3 genes, which are found in other mt genomes. Combined cooccurrence of such gene assortments indicates that the Laqueus mt genome is similar to the annelid Lumbricus, the mollusc Katharina, and the octocoral Sarcophyton mt genomes, each with statistical significance. Widely accepted schemes of metazoan phylogeny suggest that the similarity with the octocoral could have arisen through a process of convergent evolution, while it appears likely that the similarities with the annelid and the mollusc reflect phylogenetic relationships.     

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.     

The complete mitochondrial genome of Arctic Calanus hyperboreus (Copepoda,Calanoida) reveals characteristic patterns in calanoid mitochondrial genome

Copepoda is the most diverse and abundant group of crustaceans, but its phylogenetic relationships are ambiguous. Mitochondrial (mt) genomes are useful for studying evolutionary history, but only six complete Copepoda mt genomes have been made available and these have extremely rearranged genome structures. This study determined the mt genome of Calanus hyperboreus, making it the first reported Arctic copepod mt genome and the first complete mt genome of a calanoid copepod. The mt genome of C. hyperboreus is 17,910 bp in length and it contains the entire set of 37 mt genes, including 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. It has a very unusual gene structure, including the longest control region reported for a crustacean, a large tRNA gene cluster, and reversed GC skews in 11 out of 13 protein-coding genes (84.6%). Despite the unusual features, comparing this genome to published copepod genomes revealed retained pan-crustacean features, as well as a conserved calanoid-specific pattern. Our data provide a foundation for exploring the calanoid pattern and the mechanisms of mt gene rearrangement in the evolutionary history of the copepod mt genome.     

Experimental infection of self-cured Leiopelma archeyi with the amphibian chytrid Batrachochytrium dendrobatidis

The susceptibility of Archey's frog Leiopelma archeyi to Batrachochytrium dendrobatidis (Bd) is unknown, although one large population is thought to have declined sharply due to chytridiomycosis. As primary infection experiments were not permitted in this endangered New Zealand species, 6 wild-caught L. archeyi that naturally cleared infections with Bd while in captivity were exposed again to Bd to assess their immunity. These frogs were from an infected population at Whareorino, which has no known declines. All 6 L. archeyi became reinfected at low intensities, but rapidly self cured, most by 2 wk. Six Litoria ewingii were used as positive controls and developed heavier infections and clinical signs by 3 wk, demonstrating that the zoospore inoculum was virulent. Six negative controls of each species remained uninfected and healthy. Our results show that L. archeyi that have self cured have resistance to chytridiomycosis when exposed. The pattern is consistent with innate or acquired immunity to Bd, and immunological studies are needed to confirm this.     

The Mitochondrial Genome of a Deep-Sea Bamboo Coral (Cnidaria,Anthozoa, Octocorallia,Isididae): Genome Structure and Putative Origins of Replication Are Not Conserved Among Octocorals

Octocoral mitochondrial (mt) DNA is subject to an exceptionally low rate of substitution, and it has been suggested that mt genome content and structure are conserved across the subclass, an observation that has been supported for most octocorallian families by phylogenetic analyses using PCR products spanning gene boundaries. However, failure to recover amplification products spanning the nad4L-msh1 gene junction in species from the family Isididae (bamboo corals) prompted us to sequence the complete mt genome of a deep-sea bamboo coral (undescribed species). Compared to the "typical" octocoral mt genome, which has 12 genes transcribed on one strand and 5 genes on the opposite (cox2, atp8, atp6, cox3, trnM), in the bamboo coral genome a contiguous string of 5 genes (msh1, rnl, nad2, nad5, nad4) has undergone an inversion, likely in a single event. Analyses of strand-specific compositional asymmetry suggest that (i) the light-strand origin of replication was also inverted and is adjacent to nad4, and (ii) the orientation of the heavy-strand origin of replication (OriH) has reversed relative to that of previously known octocoral mt genomes. Comparative analyses suggest that intramitochondrial recombination and errors in replication at OriH may be responsible for changes in gene order in octocorals and hexacorals, respectively. Using primers flanking the regions at either end of the inverted set of five genes, we examined closely related taxa and determined that the novel gene order is restricted to the deep-sea subfamily Keratoisidinae; however, we found no evidence for strand-specific mutational biases that may influence phylogenetic analyses that include this subfamily of bamboo corals.     

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.     

啮总目昆虫的线粒体基因组多样性及系统发育研究进展

啮总目包括啮虫目（皮虱和书虱）和虱目（羽虱和吸虱）,是农业和医学等领域具有重要经济意义和研究价值的类群,目前已鉴定和描述的物种超过10 000个。啮总目昆虫线粒体基因组的变异性在昆虫各类群中最为剧烈,这些变异包括基因组的结构、基因排序、基因含量和链上分布等诸多方面。本文全面分析和总结了啮总目昆虫裂化线粒体基因组的进化属性,并结合两侧对称动物线粒体基因组的裂化特征重构了线粒体基因组环裂化的过程。引入“线粒体基因组核型”的概念来描述动物线粒体基因组丰富的变异程度。动物线粒体的染色体有减小的趋势,而线粒体基因组的裂化正是体现这种趋势的一种重要策略。同时,总结和探讨了目前具有争议的啮总目主要类群间的系统发育关系。本综述为啮总目昆虫线粒体基因组学、啮总目系统发生关系以及两侧对称动物线粒体基因组进化模式的研究提供一个新的视角。