The complete mitochondrial genome of the nudibranch Roboastra europaea (Mollusca: Gastropoda) supports the monophyly of opisthobranchs

The complete nucleotide sequence (14,472 bp) of the mitochondrial genome of the nudibranch Roboastra europaea (Gastropoda: Opisthobranchia) was determined. This highly compact mitochondrial genome is nearly identical in gene organization to that found in opisthobranchs and pulmonates (Euthyneura) but not to that in prosobranchs (a paraphyletic group including the most basal lineages of gastropods). The newly determined mitochondrial genome differs only in the relative position of the trnC gene when compared with the mitochondrial genome of Pupa strigosa, the only opisthobranch mitochondrial genome sequenced so far. Pupa and Roboastra represent the most basal and derived lineages of opisthobranchs, respectively, and their mitochondrial genomes are more similar in sequence when compared with those of pulmonates. All phylogenetic analyses (maximum parsimony, minimum evolution, maximum likelihood, and Bayesian) based on the deduced amino acid sequences of all mitochondrial protein-coding genes supported the monophyly of opisthobranchs. These results are in agreement with the classical view that recognizes Opisthobranchia as a natural group and contradict recent phylogenetic studies of the group based on shorter sequence data sets. The monophyly of opisthobranchs was further confirmed when a fragment of 2,500 nucleotides including the mitochondrial cox1, rrnL, nad6, and nad5 genes was analyzed in several species representing five different orders of opisthobranchs with all common methods of phylogenetic inference. Within opisthobranchs, the polyphyly of cephalaspideans and the monophyly of nudibranchs were recovered. The evolution of mitochondrial tRNA rearrangements was analyzed using the cox1+rrnL+nad6+nad5 gene phylogeny. The relative position of the trnP gene between the trnA and nad6 genes was found to be a synapomorphy of opisthobranchs that supports their monophyly.     

The Australian fresh water isopod (Phreatoicidea: Isopoda) allows insights into the early mitogenomic evolution of isopods

The complete mitochondrial (mt) genome sequence of the Australian fresh water isopod Eophreatoicus sp.-14 has been determined. The new species is a member of the taxon Phreatoicidea, a clade of particular interest, as it is often regarded as the sister group to all other Isopoda. Although the overall genome organization of Eophreatoicus sp.-14 conforms to the typical state of Metazoa—it is a circular ring of DNA hosting the usual 37 genes and one major non-coding region—it bears a number of derived characters that fall within the scope of “genome morphology”. Earlier studies have indicated that the isopod mitochondrial gene order is not as conserved as that of other crustaceans. Indeed, the mt genome of Eophreatoicus sp.-14 shows an inversion of seven genes (including cox1), which is as far as we know unique. Even more interesting is the derived arrangement of nad1, trnL(CUN), rrnS, control region, cob, trnT, nad5 and trnF that is shared by nearly all available isopod mt genomes. A striking feature is the close proximity of the rearranged genes to the mt control region. Inferable gene translocation events are, however, more suitable to trace the evolution of mt genomes. Genes like nad1/trnL(CUN) and nad5/trnF, which retained their adjacent position after being rearranged, were most likely translocated together. A very good example for the need to understand the mechanisms of translocations is the remolding of trnL(UUR) to trnL(CUN). Both tRNA genes are adjacent and have a high sequence similarity, probably the result of a gene duplication and subsequent anticodon mutation. Modified secondary structures were found in three tRNAs of Eophreatoicus sp.-14, which are all characterized by the loss of the DHU-arm. This is common to crustaceans for tRNA Serine(AGY), while the arm-loss in tRNA Cysteine within Malacostraca is only shared by other isopods. Modification of the third tRNA, Isoleucine, is not known from any other related species. Nucleotide frequencies of genes have been found to be indirectly correlated to the orientation of the mitochondrial replication process. In Eophreatoicus sp.-14 and in other Isopoda the associated nucleotide bias is inversed to the state of other Malacostraca. This is a strong indication for an inversion of the control region that most likely evolved in the isopod ancestor.     

The complete mitochondrial genome of the relict frog Leiopelma archeyi: insights into the root of the frog Tree of Life

Determining the root of the anuran Tree of Life is still a contentious and open question in frog systematics. Two genera with disjunct distributions have been traditionally considered the most basal among extant frogs: Leiopelma, which is endemic to New Zealand, and Ascaphus, which lives in North America. However, their specific phylogenetic position is rather elusive because each genus shows many autapomorphies, and together they retain many symplesiomorphic characters. Therefore, several alternative hypotheses have been proposed regarding the relative phylogenetic position of both Leiopelma and Ascaphus. In order to distinguish among these competing phylogenetic hypotheses, we sequenced the complete mitochondrial (mt) genome of Leiopelma archeyi and used it along with previously reported frog mt genomes (including that of Ascaphus truei) to infer a robust phylogeny of major anuran lineages. The reconstructed maximum likelihood and Bayesian inference phylogenies recovered identical topology, which supports the sister group relationship of Ascaphus and Leiopelma, and the placement of this clade at the base of the anuran tree. Interestingly, the mt genome of L. archeyi displays a novel gene arrangement in frog mt genomes affecting the relative position of cytochrome b, trnT, NADH dehydrogenase subunit 6, trnE, and trnP genes. The tandem duplication-random loss model of gene order change explains the origin of this novel frog mt genome arrangement, which is convergent with others reported in some fishes and salamanders. These results, together with comparative data for other available vertebrate mt genomes, provide evidence that the 5' end of the control region is a hot spot for gene order rearrangement.     

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.     

The first complete mitogenome of the South China deep‐sea giant isopod Bathynomus sp. (Crustacea: Isopoda: Cirolanidae) allows insights into the early mitogenomic evolution of isopods

In this study, the complete mitochondrial (mt) genome sequence of the South China deep‐sea giant isopod Bathynomus sp. was determined, and this study is the first to explore in detail the mt genome of a deep‐sea member of the order Isopoda. This species belongs to the genus Bathynomus, the members of which are saprophagous residents of the deep‐sea benthic environment; based on their large size, Bathynomus is included in the “supergiant group” of isopods. The mt genome of Bathynomus sp. is 14,965 bp in length and consists of 13 protein‐coding genes, two ribosomal RNA genes, only 18 transfer RNA genes, and a noncoding control region 362 bp in length, which is the smallest control region discovered in Isopoda to date. Although the overall genome organization is typical for metazoans, the mt genome of Bathynomus sp. shows a number of derived characters, such as an inversion of 10 genes when compared to the pancrustacean ground pattern. Rearrangements in some genes (e.g., cob, trnT, nad5, and trnF) are shared by nearly all isopod mt genomes analyzed thus far, and when compared to the putative isopod ground pattern, five rearrangements were found in Bathynomus sp. Two tRNAs exhibit modified secondary structures: The TΨC arm is absent from trnQ, and trnC lacks the DHU. Within the class Malacostraca, trnC arm loss is only found in other isopods. Phylogenetic analysis revealed that Bathynomus sp. (Cymothoida) and Sphaeroma serratum (Sphaeromatidea) form a single clade, although it is unclear whether Cymothoida is monophyletic or paraphyletic. Moreover, the evolutionary rate of Bathynomus sp. (dN/dS [nonsynonymous mutational rate/synonymous mutational rate] = 0.0705) is the slowest measured to date among Cymothoida, which may be associated with its relatively constant deep‐sea environment. Overall, our results may provide useful information for understanding the evolution of deep‐sea Isopoda species.     

Mitogenomic phylogeny of Trochoidea (Gastropoda: Vetigastropoda): New insights from increased complete genomes

Increased mitochondrial (mt) genomes can provide more sets of genome‐level characteristics for resolving deeper phylogeny. Limited information with respect to the Trochoidea mitochondrial genome organization is available; besides, monophyly and internal relationships of the superfamily still remain a matter of discussion. To resolve the monophyly and internal phylogenetic controversies of Trochoidea and expand our understanding for mt genomic characteristic evolution among Trochoidea, the phylogenetic trees were reconstructed using 13 newly sequenced complete mt genomes and 35 genomes from GenBank, and both the maximum likelihood and Bayesian inference analyses were highly supported. Vetigastropoda phylogenetic analyses recovered the monophyly of Trochoidea. Trochoidea phylogenetic analyses and genetic distances supported the non‐monophyly of Tegulidae and Tegula, indicating that the taxonomic status of several genera (Rochia, Tectus and Cittarium) should be revised and Tegula, Omphalius and Chlorostoma should be placed as a same genus. The close affinity between Tectus virgatus and Rochia was also revealed. Three‐nucleotide insertion in nad1, nine‐nucleotide insertion and six‐nucleotide deletion in nad5 are detected in Tegulidae, Tectus and Rochia, respectively. Gene orders within Trochoidea are stable, with gene rearrangements exclusive to tRNA genes observed. Homoplasious convergences because of trnT rearrangement display translocation in Turbinidae and reversion in Trochidae and Calliostomatida. For trnE and trnG, we identify 11 arrangement types, suggesting that the gene rearrangement history needs to be further evaluated. Our study emphasizes the importance of mt genomes in resolving phylogenetic relationships within Trochoidea. In addition, the mt genomic characters would contribute new insights into the classification of Trochoidea.     

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.     

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% A+T). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.     

The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms

Complete mitochondrial genome sequences for the schistosomes Schistosoma haematobium and Schistosoma. spindale have been characterized. S. haematobium is the causative agent of urinary schistosomiasis in humans and S. spindale uses ruminants as its definitive host; both are transmitted by freshwater snail intermediate hosts. Results confirm a major gene order rearrangement among schistosomes in all traditional Schistosoma species groups other than Schistosoma japonicum; i.e., species groups S. mansoni, S. haematobium, and S. indicum. These data lend support to the 'out of Asia' (East and Southeast Asia) hypothesis for Schistosoma. The gene order change involves translocation of atp6-nad2-trnA and a rearrangement of nad3-nad1 relative to other parasitic flatworm mt genomes so far sequenced. Gene order and tRNA secondary structure changes (loss and acquisition of the DHU and/or TPsiC arms of trnC, trnF, and trnR) between mitochondrial genomes of these and other (digenean and cestode) flatworms were inferred by character mapping onto a phylogeny estimated from nuclear small subunit rRNA gene sequences of these same species, in order to find additional rare genomic changes suitable as synapomorphies. Denser and wider taxon sampling of mt genomes across the Platyhelminthes will validate these putative characters.     

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 genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

Background

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

Principal Findings

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

Conclusions/Significance

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.     

土蜂科线粒体基因组序列测定和分析

【目的】本研究旨在通过针尾部(Aculeata)昆虫线粒体基因组系统发育分析认知土蜂科(Scoliidae)的单系性及系统发育位置。【方法】利用Illumina Hiseq2500二代测序技术测序土蜂科3属5种的线粒体基因组,并进行注释和分析;基于针尾部昆虫36个线粒体基因组13个蛋白质编码基因(protein-coding genes, PCGs)和2个rRNA基因序列采用最大似然法(maximum likelihood, ML)和贝叶斯法(Bayesian inference, BI)法构建系统发育树。【结果】新测序的土蜂科5个线粒体基因组为五带波壁土蜂Colpa quinquecincta线粒体基因组(GenBank登录号：OM103696),齿石波壁土蜂Colpa tartara线粒体基因组(GenBank登录号：OM103697),厚大长腹土蜂Megacampsomeris grossa线粒体基因组(GenBank登录号：OM103796),台湾大长腹土蜂Megacampsomeris formosensis线粒体基因组(GenBank登录号：OM142776)和斯式土蜂Sc...     

Molecular variability in the Celleporella hyalina (Bryozoa; Cheilostomata) species complex: evidence for cryptic speciation from complete mitochondrial genomes

The bryozoan Celleporella has been shown to be composed of multiple, often cryptic, lineages. We sequenced two complete mitochondrial (mt) genomes of the Celleporella hyalina species complex from Wales, UK and Norway (i) to determine genetic divergence at the complete mt genome level, and (ii) to design new molecular markers for examining the interrelationships amongst the major lineages. In addressing (i), we estimated genetic divergence at three levels: (a) nucleotide diversity (π), (b) genome size, and (c) gene order. Genes nad4L, nad6, and atp8 showed the highest levels of divergence, and rrnL, rrnS, and cox1 showed the lowest levels. Inter-genome nucleotide divergence of protein-coding and ribosomal RNA genes, measured as π, was 0.21. The two genomes differed substantially in size, with the Norwegian genome being 2,573 base pairs (bp) longer than the Welsh genome, 17,265 and 14,692 bp, respectively. This difference in size is attributable to long non-coding regions present in the Norwegian genome. Both genomes exhibit similar gene orders, except for the translocation of one transfer RNA (trnA). Considering the high nucleotide diversity, genome size difference and change in gene order, these mt genomes are considered sufficiently divergent to have originated from two distinct species. In addressing (ii) we designed PCR primers that flank the most conserved regions of the genome: 1,300 bp of cox1 and a contiguous 2,000 bp fragment of rrnL + rrnS. The primers have yielded products for tissue from Wales, Norway, New Zealand, Alaska and Chile and should provide useful tools in establishing species- and population-level diversity within the Celleporella complex.     

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

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

New features of Asian Crassostrea oyster mitochondrial genomes: A novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.     

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: a case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another 11 species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.     

Mitochondrial genome reorganization characterizes various lineages of mesostigmatid mites (Acari: Parasitiformes)

Mesostigmata is an extremely diverse group of mites with more than 11,000 described species in 109 families. The complete mitochondrial (mt) genomes of five species of mesostigmatid mites from three families (Varroidae, Ologamasidae, Phytoseiidae) have been reported previously; all of them are rearranged or highly rearranged in gene order. However, it is unclear when mt genome reorganization occurred and how common it is in mesostigmatid mites. We sequenced the mt genomes of ten species of mesostigmatid mites from five more families (Blattisociidae, Diplogyniidae, Laelapidae, Macrochelidae, Parasitidae). We found that species in the families Diplogyniidae and Parasitidae have retained the ancestral mt genome organization of arthropods, which is in stark contrast to the highly rearranged mt genomes in the Phytoseiidae species. As in the Varroidae and Ologamasidae species, the mt genomes of the Blattisociidae, Macrochelidae and Laelapidae species are also rearranged but are less rearranged than in the Phytoseiidae species. Each of the six mesostigmatid families that have rearranged mt genomes is characterized by unique gene order not seen in other mesostigmatid families. Furthermore, the mt genome organization also differs among three genera of the Phytoseiidae, between two genera of the Laelapidae, and among three Macrocheles species of the Macrochelidae. Our results indicate that: (a) the most recent common ancestor of mesostigmatid mites likely retained the ancestral mt genome organization of arthropods; and (b) mt genome organization characterizes various lineages of mesostigmatid mites and provides a valuable source of information for understanding their phylogeny and evolution.     

Complete mtDNA sequence of the North American freshwater mussel,Lampsilis ornata (Unionidae): an examination of the evolution and phylogenetic utility of mitochondrial genome organization in Bivalvia (Mollusca)

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. Two factors inhibiting our understanding the evolution of gene rearrangement in bivalves are inadequate taxonomic sampling and failure to examine gene order in a phylogenetic framework. Here, we report the first complete nucleotide sequence (16,060 bp) of the mitochondrial (mt) genome of a North American freshwater bivalve, Lampsilis ornata (Mollusca: Paleoheterodonta: Unionidae). Gene order and mt genome content is examined in a comparative phylogenetic framework for Lampsilis and five other bivalves, representing five families. Mitochondrial genome content is shown to vary by gene duplication and loss among taxa and between male and female mitotypes within a species. Although mt gene arrangement is highly variable among bivalves, when optimized on an independently derived phylogenetic hypothesis, it allows for the reconstruction of ancestral gene order states and indicates the potential phylogenetic utility of the data. However, the interpretation of reconstructed ancestral gene order states must take in to account both the accuracy of the phylogenetic estimation and the probability of character state change across the topology, such as the presence/absence of atp8 in bivalve lineages. We discuss what role, if any, doubly uniparental inheritance (DUI) and recombination between sexual mitotypes may play in influencing gene rearrangement of the mt genome in some bivalve lineages.