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 Sequences of Heptathela hangzhouensis and Ornithoctonus huwena Reveal Unique Gene Arrangements and Atypical tRNAs

We have sequenced the complete mitochondrial genomes of the spiders Heptathela hangzhouensis and Ornithoctonus huwena. Both genomes encode 13 protein-coding genes, 22 tRNA genes, and 2 ribosomal RNA genes. H. hangzhouensis, a species of the suborder Mesothelae and a representative of the most basal clade of Araneae, possesses a gene order identical to that of Limulus polyphemus of Xiphosura. On the other hand, O. huwena, a representative of suborder Opisthothelae, infraorder Mygalomorphae, was found to have seven tRNA genes positioned differently from those of Limulus. The rrnL–trnL1–nad1 arrangement shared by the araneomorph families Salticidae, Nesticidae, and Linyphiidae and the mygalomorph family Theraphosidae is a putative synapomorphy joining the mygalomorph with the araneomorph. Between the two species examined, base compositions also differ significantly. The lengths of most protein-coding genes in H. hangzhouensis and O. huwena mtDNA are either identical to or slightly shorter than their Limulus counterparts. Usage of initiation and termination codons in these protein-coding genes seems to follow patterns conserved among most arthropod and some other metazoan mitochondrial genomes. The sequences of the 3 ends of rrnS and rrnL in the two species are similar to those reported for Limulus, and the entire genes are shortened by about 100–250 nucleotides with respect to Limulus. The lengths of most tRNA genes from the two species are distinctly shorter than those of Limulus and the sequences reveal unusual inferred tRNA secondary structures. Our finding provides new molecular evidence supporting that the suborder Mesothelae is basal to opisthothelids.Reviewing Editor Dr. Rafael Zardoya     

The Cephalopod Loligo bleekeri Mitochondrial Genome: Multiplied Noncoding Regions and Transposition of tRNA Genes

We previously reported the sequence of a 9260-bp fragment of mitochondrial (mt) DNA of the cephalopod Loligo bleekeri [J. Sasuga et al. (1999) J. Mol. Evol. 48:692–702]. To clarify further the characteristics of Loligo mtDNA, we have sequenced an 8148-bp fragment to reveal the complete mt genome sequence. Loligo mtDNA is 17,211 bp long and possesses a standard set of metazoan mt genes. Its gene arrangement is not identical to any other metazoan mt gene arrangement reported so far. Three of the 19 noncoding regions longer than 10 bp are 515, 507, and 509 bp long, and their sequences are nearly identical, suggesting that multiplication of these noncoding regions occurred in an ancestral Loligo mt genome. Comparison of the gene arrangements of Loligo, Katharina tunicata, and Littorina saxatilis mt genomes revealed that 17 tRNA genes of the Loligo mt genome are adjacent to noncoding regions. A majority (15 tRNA genes) of their counterparts is found in two tRNA gene clusters of the Katharina mt genome. Therefore, the Loligo mt genome (17 tRNA genes) may have spread over the genome, and this may have been coupled with the multiplication of the noncoding regions. Maximum likelihood analysis of mt protein genes supports the clade Mollusca + Annelida + Brachiopoda but fails to infer the relationships among Katharina, Loligo, and three gastropod species. Received: 9 May 2001 / Accepted: 3 October 2001     

Genome/transcriptome analysis of the chigger mite Leptotrombidium pallidum,a major vector for scrub typhus,with a special focus on genes more abundantly expressed in larval stage

Leptotrombidium pallidum is the major vector mite for Orientia tsutsugamushi, the causative agent of scrub typhus, in Asian countries, including Korea. Despite its medical importance, L. pallidum has little genetic information available to date. To analyze the L. pallidum genome, we extracted genomic DNA (gDNA) from a single female of a 7-generation inbred L. pallidum colony and amplified the gDNA by whole genome amplification (WGA). The resulting amplified gDNA was used to construct paired-end and mate-pair libraries that were sequenced using Illumina platforms (HiSeq2000 and MiSeq). More than 45 Gb of sequence reads from both paired-end and mate-pair libraries of the WGA gDNA were trimmed and then de novo assembled using CLC Assembly Cell v.4.0 for contig assembly and SSPACE for scaffolding. The assembly generated approximately 6,545 scaffolds with an N50 value of 92,945 and total size of ~ 193 Mb. For gene predictions, the PASA and GeneWise models were used, and ab initio gene predictions were performed independently, resulting in the prediction of 15,842 genes. RNA-Seq expression profiles revealed constitutive expression of 11,572 unique protein-coding genes in larva, 12,364 in protonymph, 12,872 in male adult, and 12,617 in female adult stages. Of the 15,842 predicted genes, 10,885 were commonly expressed through all L. pallidum stages. Genes selectively over-transcribed in the larval stage, which is when host parasitization and disease transmission occur, were further annotated, and their putative roles were discussed.     

Mitochondrial gene sequences alone or combined with ITS region sequences provide firm molecular criteria for the classification of Lecanicillium species

The recent revision of Verticillium sect. Prostrata led to the introduction of the genus Lecanicillium, which comprises the majority of the entomopathogenic strains. Sixty-five strains previously classified as Verticillium lecanii or Verticillium sp. from different geographical regions and hosts were examined and their phylogenetic relationships were determined using sequences from three mitochondrial (mt) genes [the small rRNA subunit (rns), the NADH dehydrogenase subunits 1 (nad1) and 3 (nad3)] and the ITS region. In general, single gene phylogenetic trees differentiated and placed the strains examined in well-supported (by BS analysis) groups of L. lecanii, L. longisporum, L. muscarium, and L. nodulosum, although in some cases a few uncertainties still remained. nad1 was the most informative single gene in phylogenetic analyses and was also found to contain group I introns with putative open reading frames (ORFs) encoding for GIY–YIG endonucleases. The combined use of mt gene sequences resolved taxonomic uncertainties arisen from ITS analysis and, alone or in combination with ITS sequences, helped in placing uncharacterised Verticillium lecanii and Verticillium sp. firmly into Lecanicillium species. Combined gene data from all the mt genes and all the mt genes and the ITS region together, were very similar. Furthermore, a relaxed correlation with host specificity—at least for Homoptera—was indicated for the rns and the combined mt gene sequences. Thus, the usefulness of mt gene sequences as a convenient molecular tool in phylogenetic studies of entomopathogenic fungi was demonstrated.     

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.     

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

Background

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

Results

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

Conclusions

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Electronic supplementary material

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

Cloning and characterization of the Chlamydomonas moewusii mitochondrial genome.

Summary We report that the mitochondrial genome of Chlamydomonas moewusii has a 22 kb circular map and thus contrasts with the mitochondrial genome of Chlamydomonas reinhardtii, which is linear and about 6 kb shorter. Overlapping restriction fragments spanning over 90% of the C. moewusii mitochondrial DNA (mtDNA) were identified in a clone bank constructed using a Sau3AI partial digest of a C. moewusii DNA fraction enriched for mtDNA by preparative CsCI density gradient centrifugation. Overlapping Sau3AI clones were identified by a chromosome walk initiated with a clone of C. moewusii mtDNA. The mtDNA map was completed by Southern blot analysis of the C. moewusii mtDNA fraction using isolated mtDNA clones. Regions that hybridized to C. reinhardtii or wheat mitochondrial gene probes for subunit I of cytochrome oxidase (cox1), apocytochrome b (cob), three subunits of NADH dehydrogenase (nadl, nad2 and nad5) and the small and the large ribosomal RNAs (rrnS and rrnL, respectively) were localized on the C. moewusii mtDNA map by Southern blot analysis. The results show that the order of genes in the mitochondrial genome of C. moewusii is completely rearranged relative to that of C. reinhardtii.     

Comparative Analyses of Coding and Noncoding DNA Regions Indicate that Acropora (Anthozoa: Scleractina) Possesses a Similar Evolutionary Tempo of Nuclear vs. Mitochondrial Genomes as in Plants

Evidence suggests that the mitochondrial (mt)DNA of anthozoans is evolving at a slower tempo than their nuclear DNA; however, parallel surveys of nuclear and mitochondrial variations and calibrated rates of both synonymous and nonsynonymous substitutions across taxa are needed in order to support this scenario. We examined species of the scleractinian coral genus Acropora, including previously unstudied species, for molecular variations in protein-coding genes and noncoding regions of both nuclear and mt genomes. DNA sequences of a calmodulin (CaM)-encoding gene region containing three exons, two introns and a 411-bp mt intergenic spacer (IGS) spanning the cytochrome b (cytb) and NADH 2 genes, were obtained from 49 Acropora species. The molecular evolutionary rates of coding and noncoding regions in nuclear and mt genomes were compared in conjunction with published data, including mt cytochrome b, the control region, and nuclear Pax-C introns. Direct sequencing of the mtIGS revealed an average interspecific variation comparable to that seen in published data for mt cytb. The average interspecific variation of the nuclear genome was two to five times greater than that of the mt genome. Based on the calibration of the closure of Panama Isthmus (3.0 mya) and closure of the Tethy Seaway (12 mya), synonymous substitution rates ranged from 0.367% to 1.467% Ma⁻¹ for nuclear CaM, which is about 4.8 times faster than those of mt cytb (0.076–0.303% Ma⁻¹). This is similar to the findings in plant genomes that the nuclear genome is evolving at least five times faster than those of mitochondrial counterparts. I-Ping Chen and Chung-Yu Tang, co-first author (equal contribution)     

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

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

后生动物线粒体基因组:起源、大小和基因排列进化

由于受到强烈的进化约束,后生动物线粒体基因组在大小和基因含量上一直保持稳定,相比之下核基因组则发生了巨大的改变。后生动物线粒体基因组结构的可塑性在一定程度上归功于可能由tRNA基因介导的基因重排事件,虽然亲缘关系密切的物种间也可能出现基因重排,但同门内的线粒体基因组仍趋向于具有类似的结构特征。我们对后生动物线粒体基因组的起源、大小和基因排列进化方面的特点进行了介绍。     

Mitogenomic analysis of Montipora cactus and Anacropora matthai (cnidaria; scleractinia; acroporidae) indicates an unequal rate of mitochondrial evolution among Acroporidae corals

The complete nucleotide sequence of the mitochondrial (mt) genome was determined for specimens of the coral species Montipora cactus (Bernard 1897) and Anacropora matthai (Pillai 1973), representing two morphologically distinct genera of the family Acroporidae. These sequences were compared with the published mt genome sequence for the confamilial species, Acropora tenuis (Dana 1846). The size of the mt genome was 17,887 bp and 17,888 bp for M. cactus and A. matthai. Gene content and organization was found to be very similar among the three Acroporidae mt genomes with a group I intron occurring in the NADH dehyrogenase 5 (nad5) gene. The intergenic regions were also similar in length among the three corals. The control region located between the small ribosomal RNA (ms) and the cytochrome oxidase 3 (cox3) gene was significantly smaller in M. cactus and A. matthai (both 627 bp) than in A. tenuis (1086 bp). Only one set of repeated sequences was identified at the 3′-end of the control regions in M. cactus and A. matthai. A lack of the abundant repetitive elements which have been reported for A. tenuis, accounts for the relatively short control regions in M. cactus and A. matthai. Pairwise distances and relative rate analyses of 13 protein coding genes, the group I intron and the largest intergenic region, igr3, revealed significant differences in the rate of molecular evolution of the mt genome among the three species, with an extremely slow rate being seen between Montipora and Anacropora. It is concluded that rapid mt genome evolution is taking place in genus Acropora relative to the confamilial genera Montipora and Anacropora although all are within the relatively slow range thought to be typical of Anthozoa.     

Gene expression suggests conserved mechanisms patterning the heads of insects and myriapods

Segmentation, i.e. the subdivision of the body into serially homologous units, is one of the hallmarks of the arthropods. Arthropod segmentation is best understood in the fly Drosophila melanogaster. But different from the situation in most arthropods in this species all segments are formed from the early blastoderm (so called long-germ developmental mode). In most other arthropods only the anterior segments are formed in a similar way (so called short-germ developmental mode). Posterior segments are added one at a time or in pairs of two from a posterior segment addition zone. The segmentation mechanisms are not universally conserved among arthropods and only little is known about the genetic patterning of the anterior segments. Here we present the expression patterns of the insect head patterning gene orthologs hunchback (hb), orthodenticle (otd), buttonhead-like (btdl), collier (col), cap-n-collar (cnc) and crocodile (croc), and the trunk gap gene Krüppel (Kr) in the myriapod Glomeris marginata. Conserved expression of these genes in insects and a myriapod suggests that the anterior segmentation system may be conserved in at least these two classes of arthropods. This finding implies that the anterior patterning mechanism already existed in the last common ancestor of insects and myriapods.     

Phylogenomic Analysis of the PEBP Gene Family in Cereals

The TFL1 and FT genes, which are key genes in the control of flowering time in Arabidopsis thaliana, belong to a small multigene family characterized by a specific phosphatidylethanolamine-binding protein domain, termed the PEBP gene family. Several PEBP genes are found in dicots and monocots, and act on the control of flowering time. We investigated the evolution of the PEBP gene family in cereals. First, taking advantage of the complete rice genome sequence and EST databases, we found 19 PEBP genes in this species, 6 of which were not previously described. Ten genes correspond to five pairs of paralogs mapped on known duplicated regions of the rice genome. Phylogenetic analysis of Arabidopsis and rice genes indicates that the PEBP gene family consists of three main homology classes (the so-called TFL1-LIKE, MFT-LIKE, and FT-LIKE subfamilies), in which gene duplication and/or loss occurred independently in Arabidopsis and rice. Second, phylogenetic analyses of genomic and EST sequences from five cereal species indicate that the three subfamilies of PEBP genes have been conserved in cereals. The tree structure suggests that the ancestral grass genome had at least two MFT-like genes, two TFL1-like genes, and eight FT-like genes. A phylogenomic approach leads to some hypotheses about conservation of gene function within the subfamilies. [Reviewing Editor: Dr. Yves Van de Peer]     

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.     

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.     

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.     

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.