The mitochondrial genome of the water spider Argyroneta aquatica (Araneae: Cybaeidae)

We sequenced nearly the entire mitochondrial genome of Argyroneta aquatica, a wholly underwater‐living spider, thereby enhancing the available genomic information for Arachnida. The confirmed sequences contained the complete set of known genes present in other metazoan mitochondrial genomes. However, the mitochondrial gene order of A. aquatica was distinctly different from that of the most distant Chelicerata Limulus polyphemus (Xiphosura), probably because of a series of gene translocations and/or inversions. Comparison of arachnid mitochondrial gene orders for the purpose of phylogenetic inference is only minimally useful, but provides a strong signal in closely related lineages. To test the basal relationships and the evolutionary pattern of tRNA gene rearrangements among Arachnida, phylogenetic analyses using amino acid sequences of the 13 protein‐coding genes were performed. An interesting feature, the five 135‐bp tandem repeats and two 363‐bp tandem repeats, was identified in the putative control region. Although control region tandem repeats have been reported in many other arachnid and metazoan species, this is the first time it has been described in spiders.     

Contrasting intra versus interspecies DNA sequence variation for representatives of the Batrachospermales (Rhodophyta): Insights from a DNA barcoding approach

Sixty‐five accessions of the species‐rich freshwater red algal order Batrachospermales were characterized through DNA sequencing of two regions: the mitochondrial cox1 gene (664 bp), which is proposed as the DNA barcode for red algae, and the UPA (universal plastid amplicon) marker (370 bp), which has been recently identified as a universally amplifying region of the plastid genome. upgma phenograms of both markers were consistent in their species‐level relationships, although levels of sequence divergence were very different. Intraspecific variation of morphologically identified accessions for the cox1 gene ranged from 0 to 67 bp (divergences were highest for the two taxa with the greatest number of accessions; Batrachospermum helminthosum and Batrachospermum macrosporum); while in contrast, the more conserved universal plastid amplicon exhibited much lower intraspecific variation (generally 0–3 bp). Comparisons to previously published mitochondrial cox2–3 spacer sequences for B. helminthosum indicated that the cox1 gene and cox2–3 spacer were characterized by similar levels of sequence divergence, and phylogeographic patterns based on these two markers were consistent. The two taxa represented by the largest numbers of specimens (B. helminthosum and B. macrosporum) have cox1 intraspecific divergence values that are substantially higher than previously reported, but no morphological differences can be discerned at this time among the intraspecific groups revealed in the analyses. DNA barcode data, which are based on a short fragment of an organellar genome, need to be interpreted in conjunction with other taxonomic characters, and additional batrachospermalean taxa need to be analyzed in detail to be able to draw generalities regarding intraspecific variation in this order. Nevertheless, these analyses reveal a number of batrachospermalean taxa worthy of more detailed DNA barcode study, and it is predicted that such research will have a substantial effect on the taxonomy of species within the Batrachospermales in the future.     

Mitochondrial Genome of Prasinophyte Alga Pyramimonas parkeae

Prasinophytes are a paraphyletic assemblage of nine heterogeneous lineages in the Chlorophyta clade of Archaeplastida. Until now, seven complete mitochondrial genomes have been sequenced from four prasinophyte lineages. Here, we report the mitochondrial genome of Pyramimonas parkeae, the first representative of the prasinophyte clade I. The circular‐mapping molecule is 43,294 bp long, AT rich (68.8%), very compact and it comprises two 6,671 bp long inverted repeat regions. The gene content is slightly smaller than the gene‐richest prasinophyte mitochondrial genomes. The single identified intron is located in the cytochrome c oxidase subunit 1 gene (cox1). Interestingly, two exons of cox1 are encoded on the same strand of DNA in the reverse order and the mature mRNA is formed by trans‐splicing. The phylogenetic analysis using the data set of 6,037 positions assembled from 34 mtDNA‐encoded proteins of 48 green algae and plants is not in compliance with the branching order of prasinophyte clades revealed on the basis of 18S rRNA genes and cpDNA‐encoded proteins. However, the phylogenetic analyses based on all three genomic elements support the sister position of prasinophyte clades Pyramimonadales and Mamiellales.     

Phylogeny and species delineation in the marine diatom Pseudo‐nitzschia (Bacillariophyta) using cox1, LSU,and ITS2 rRNA genes: A perspective in character evolution

Analyses of the mitochondrial cox1, the nuclear‐encoded large subunit (LSU), and the internal transcribed spacer 2 (ITS2) RNA coding region of Pseudo‐nitzschia revealed that the P. pseudodelicatissima complex can be phylogenetically grouped into three distinct clades (Groups I–III), while the P. delicatissima complex forms another distinct clade (Group IV) in both the LSU and ITS2 phylogenetic trees. It was elucidated that comprehensive taxon sampling (sampling of sequences), selection of appropriate target genes and outgroup, and alignment strategies influenced the phylogenetic accuracy. Based on the genetic divergence, ITS2 resulted in the most resolved trees, followed by cox1 and LSU. The morphological characters available for Pseudo‐nitzschia, although limited in number, were overall in agreement with the phylogenies when mapped onto the ITS2 tree. Information on the presence/absence of a central nodule, number of rows of poroids in each stria, and of sectors dividing the poroids mapped onto the ITS2 tree revealed the evolution of the recently diverged species. The morphologically based species complexes showed evolutionary relevance in agreement with molecular phylogeny inferred from ITS2 sequence–structure data. The data set of the hypervariable region of ITS2 improved the phylogenetic inference compared to the cox1 and LSU data sets. The taxonomic status of P. cuspidata and P. pseudodelicatissima requires further elucidation.     

Phylogeny of Rhigonematomorpha based on the complete mitochondrial genome of Rhigonema thysanophora (Nematoda: Chromadorea)

We determined the complete mitochondrial genome sequence of Rhigonema thysanophora, the first representative of Rhigonematomorpha, and used this sequence along with 57 other nematode species for phylogenetic analyses. The R. thysanophora mtDNA is 15 015 bp and identical to all other chromadorean nematode mtDNAs published to date in that it contains 36 genes (lacking atp8) encoded in the same direction. Phylogenetic analyses of nucleotide and amino acid sequence data for the 12 protein‐coding genes recovered Rhigonematomorpha as the sister group to the heterakoid species, Ascaridia columbae (Ascaridomorpha). The organization of R. thysanophora mtDNA resembles the most common pattern for the Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha clade in gene order, but with some substantial gene rearrangements. This similarity in gene order is in agreement with the sequence‐based analyses that indicate a close relationship between Rhigonematomorpha and Rhabditomorpha+Ascaridomorpha+Diplogasteromorpha. These results are consistent with certain analyses of nuclear SSU rDNA for R. thysanophora and some earlier classification systems that asserted phylogenetic affinity between Rhigonematomorpha and Ascaridomorpha, but inconsistent with morphology‐based phylogenetic hypotheses that suggested a close (taxonomic) relationship between rhigonematomorphs and oxyuridomorphs (pinworms). These observations must be tempered by noting that few rhigonematomorph species have been sequenced and included in phylogenetic analyses, and preliminary studies based on SSU rDNA suggest the group is not monophyletic. Additional mitochondrial genome sequences of rhigonematids are needed to characterize their phylogenetic relationships within Chromadorea, and to increase understanding of mitochondrial genome evolution.     

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

Complete mitochondrial genome of the black‐tailed hornet,Vespa ducalis (Hymenoptera: Vespidae): Genomic comparisons in Vespoidea

We sequenced the complete mitochondrial genome (mitogenome) of the black‐tailed hornet, Vespa ducalis (Hymenoptera: Vespidae). The genome was 15,779‐bp long and contained typical sets of genes [13 protein‐coding genes (PCGs), 22 tRNAs, and 2 rRNAs]. The V. ducalis A + T‐rich region was 166‐bp long and was the shortest of all sequenced Vespoidea genomes, including Vespa. The genome was highly biased toward A/T nucleotides—80.1 % in the whole genome, 77.8 % in PCGs, 83.4–85.6 % in RNAs, and 92.8 % in the A + T‐rich region. These values are well within the typical range for genes and regions of Vespoidea mitogenomes. Start and stop codons in several Vespa species—including V. ducalis—were diversified, despite these species belonging to the same genus. In comparison with the ancestral mitogenomes, Vespa mitogenomes—including that of V. ducalis—showed substantial gene rearrangement; however, we detected no gene rearrangement among Vespa species. We conducted phylogenetic reconstruction based on concatenated sequences of 13 PCGs and two rRNAs (12,755 bp ) in available species of Vespoidea—21 species in six subfamilies in two families (Vespidae and Formicidae). The Bayesian inference and maximum likelihood (ML) methods revealed that each family formed strong monophyletic groups [Bayesian posterior probability (BPP) = 1; ML, 100 %]. Moreover, V. ducalis and V. mandarinia formed a strong sister group (BPP = 1; ML, 94 %).     

Complete mitochodrial genome of the crab spider Ebrechtella tricuspidata (Araneae: Thomisidae): A novel tRNA rearrangement and phylogenetic implications for Araneae

Mitochondrial genomes are widely used for phylogenetic and phylogeographic analyses among arthropods, but there is a lack of sufficient mitochondrial genome sequence data for spiders. Herein, we sequenced and characterized the complete mitochondrial genome of a crab spider Ebrechtella tricuspidata (Araneae: Thomisidae). The circular mitochondrial genome is 14,352 bp long, including a standard set of 37 genes and an A + T-rich region. Nucleotide composition is highly biased toward A + T nucleotides (77.3%). A novel gene order rearrangement is detected by a tRNA (trnL1) translocation. Tandem repeats are not identified in the A + T-rich region. Most of the tRNAs are greatly reduced in size and cannot be folded into typical cloverleaf-shaped secondary structures. The phylogenetic analysis confirms that the mitochondrial genome sequences are useful in resolving higher-level relationship of Araneae. Overall, our data present in this study will elevate our knowledge on the architecture and evolution of spider mitochondrial genome.     

The complete chloroplast genome of the spring ephemeral plant Alyssum desertorum and its implications for the phylogenetic position of the tribe Alysseae within the Brassicaceae

Alyssum desertorum (Alysseae, Brassicaceae) is an annual spring ephemeral plant whose life cycle is only 2–3 months. It typically has high photosynthetic capacity and a high growth rate. However, little was known about the chloroplast (cp) genome structure of this species. Furthermore, the phylogenetic position of the tribe Alysseae relative to other tribes in the Brassicaceae has not been established and there appear to be inconsistences between different DNA markers. This study is the first report on a cp genome of the genus Alyssum and discusses the phylogenetic relationships of the tribe Alysseae relative to other tribes in the family. The complete cp genome of A. desertorum was 151 677 bp in size and is thus the smallest cp genome of Brassicaceae sequenced to date. The genome includes a large single‐copy region of 81 551 bp, a small single‐copy region of 17 804 bp, and two inverted repeats of 26 161 bp each. The genome contains 132 genes, including 86 protein‐coding genes (PCGs), 38 tRNA genes and 8 rRNA genes. A total of 16 genes contained introns, including 10 PCGs and 6 tRNA genes; the ycf3 and clpP genes contained two introns, and the remaining genes each contained one. Compared to the cp genomes of 21 other Brassicaceae species, the cp genome of Alyssum desertorum was the smallest, as due to variation in gene content and gene length, such as a lack of the rps16 gene and the deletion of some coding genes. Additionally, deletions of introns and intergenic spacers were observed, but their total length was not significantly shorter than those of other taxa. Phylogenetic analysis at the tribal level based on a cp genome dataset revealed that the tribe Alysseae is an early‐diverging lineage that is sister to other species within subclade B of clade II.     

Presence of nuclear copies of mitochondrial origin (NUMTs) in two related species of stingless bee genus Melipona (Hymenoptera: Meliponini)

Transferred copies of mitochondrial DNA (mtDNA) into the nuclear genome (numts) have been reported in several Hymenoptera species, even at a high density in the honey bee nuclear genome. The accidental amplification of numts in phylogenetic studies focused on mtDNA highlights the importance of a correct determination of numts and their related mtDNA sequences. We report here the presence of numts derived from a mitochondrial rDNA 16S gene in the genome of the stingless bee species Melipona colimana and M. fasciata (tribe Meliponini) from Western Mexico. PCR products were cloned in both species obtaining thirty paralogous numts. Numts were identified by the presence of insertions and deletions and the disruption of the 16S secondary structure. Further phylogenetic analyses including alternative mitochondrial cox1 and nuclear ITS1 genes have revealed the presence of another numt (cox1) in the nuclear genome of these two species, and place both as sister lineages within the subgenus Michmelia. This is one of the first studies reporting the presence of numts in Meliponini species, and supports previous studies suggesting frequent transfer of mtDNA to the nuclear genome in Hymenoptera.     

The complete mitochondrial genome of the Sichuan taimen (Hucho bleekeri): repetitive sequences in the control region and phylogenetic implications for Salmonidae

The complete mitochondrial DNA genome of the Sichuan taimen (Hucho bleekeri) was determined by the long and accurate polymerase chain reaction (LA-PCR) and primer walking sequence method. The entire mitochondrial genome of this species is 16,997 bp in length, making it the longest among the completely sequenced Salmonidae mitochondrial genomes. It consists of two ribosomal RNA (rRNA) genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes, and one control region (CR). The gene arrangement, nucleotide composition, and codon usage pattern of the mitochondrial genome are similar to those of other teleosts. A T-type mononucleotide microsatellite and an 82 bp tandem repeat were identified in the control region, which were almost identical among the three H. bleekeri individuals examined. Both phylogenetic analyses based on 12 concatenated protein-coding genes of the heavy strand and on just the control region show that H. bleekeri is a basal species in Salmoninae. In addition, Salmo, Salvelinus and Oncorhynchus all represent monophyletic groups, respectively. All freshwater species occupied basal phylogenetic positions, and also possessed various tandem repeats in their mitochondrial control regions. These results support established phylogenetic relationships among genera in Salmonidae based on morphological and molecular analyses, and are consistent with the hypothesis that Salmonidae evolved from freshwater species.     

Complete mitochondrial genome of the darkling beetle Gonocephalum outreyi (Coleoptera: Tenebrionidae) with phylogenetic implications

The complete mitochondrial genome (mitogenome) of Gonocephalum outreyi was determined by using next-generation sequencing approach. The full length of this mitogenome is 15,836?bp, which consists of 37 typical metazoan mitochondrial genes with an identical genome organization to ancestral insects. The majority of the protein-coding genes begin with the codon ATN, except for cox1 and cox2 with AAT and AAA, respectively. To elucidate the phylogenetic position of G. outreyi, we used various sequence coding schemes for protein-coding genes and the combined nucleotide sequences of all mitochondrial genes for tree building under the Bayesian and Maximum Likelihood inferences. The phylogenetic results consistently supported G. outreyi as a member of the family Tenebrionidae. The monophyly of both Tenebrionoidea and Tenebrionidae were strongly supported. The Scraptiidae and Melandryidae were recovered to be non-monophyletic in regards to the Osphya. Within Tenebrionidae, the subfamilies Diaperinae and Tenebrioninae were found to be non-monophyletic.     

Complete mitochondrial genome of Tubulipora flabellaris (Bryozoa: Stenolaemata): the first representative from the class Stenolaemata with unique gene order

Mitochondrial genomes play a significant role in the reconstruction of phylogenetic relationships within metazoans. There are still many controversies concerning the phylogenetic position of the phylum Bryozoa. In this research, we have finished the complete mitochondrial genome of one bryozoan (Tubulipora flabellaris), which is the first representative from the class Stenolaemata. The complete mitochondrial genome of T. flabellaris is 13,763 bp in length and contains 36 genes, which lacks the atp8 gene in contrast to the typical metazoan mitochondrial genomes. Gene arrangement comparisons indicate that the mitochondrial genome of T. flabellaris has unique gene order when compared with other metazoans. The four known bryozoans complete mitochondrial genomes also have very different gene arrangements, indicates that bryozoan mitochondrial genomes have experienced drastic rearrangements. To investigate the phylogenetic relationship of Bryozoa, phylogenetic analyses based on amino acid sequences of 11 protein coding genes (excluding atp6 and atp8) from 26 metazoan complete mitochondrial genomes were made utilizing Maximum Likelihood (ML) and Bayesian methods, respectively. The results indicate the monopoly of Lophotrochozoa and a close relationship between Chaetognatha and Bryozoa. However, more evidences are needed to clarify the relationship between two groups. Lophophorate appeared to be polyphyletic according to our analyses. Meanwhile, neither analysis supports close relationship between Branchiopod and Phoronida. Four bryozoans form a clade and the relationship among them is T. flabellaris + (F. hispida + (B. neritina + W. subtorquata)), which is in coincidence with traditional classification system.     

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.     

Analysis of the mitochondrial genome of Schistocerca gregaria gregaria (Orthoptera: Acrididae)

In the present study, we present the full sequence of the mitochondrial genome of the African desert locust Schistocerca gregaria gregaria. The size of 15625 bp reported matches very well with mitochondrial genomes of other Orthopteriodea. The mitochondrial genome comprises 13 protein‐coding genes, two ribosomal RNAs and 22 t‐RNAs with two t‐RNA (trnD and trnK) rearrangements that are typical for the taxon Caelifera. We compared the sequence with 12 mitochondrial genes of Schistocerca gregaria flaviventris and Schistocerca americana and used some of these data to construct phylogenetic trees, which confirm the close relationship between the two subspecies S. g. flaviventris and S. g. gregaria. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 296–305.     

Phylogeny of Veneridae (Bivalvia) based on mitochondrial genomes

Veneridae is one of the most diverse families of bivalve molluscs. However, their phylogenetic relationships among subfamilies have been debated for years. To explore phylogenetic relationships of Veneridae, we sequenced 13 complete mitochondrial genome sequences from eight subfamilies and compared with available complete mitochondrial genome of other Veneridae taxa (18 previously reported sequences). Phylogenetic analyses using probabilistic methods recovered two highly supported clades. In addition, the protein‐coding gene order revealed a highly conserved pattern among the same subclade lineages. According to our molecular analyses, Tapetinae should be recognized as a valid subfamily, but the genera formed para‐polyphyletic clades. Chioninae was recovered not monophyletic that differs from a previously molecular phylogeny. Furthermore, the reconstructed chronogram calibrated with fossils recovered the Veneridae have originated during the early Permian (about 290 million years ago). Noticeably, programmed frameshift was found in the nad4 gene of Leukoma jedoensis, Anomalodiscus squamosus and Antigona lamellaris and cob gene of L. jedoensis. This is the first time that the presence of the programmed frameshift has been found in the protein‐coding genes of Heterodonta species. Our results improved the phylogenetic resolution within Veneridae, and a more taxonomic sampling analysis of the subfamily Chioninae is supposed to construct.     

A Comparative Characterization of the Mitochondrial Genomes of Paramoeba aparasomata and Neoparamoeba pemaquidensis (Amoebozoa,Paramoebidae)

Marine amebae of the genus Paramoeba (Amoebozoa, Dactylopodida) normally contain a eukaryotic endosymbiont known as Perkinsela‐like organism (PLO). This is one of the characters to distinguish the genera Neoparamoeba and Paramoeba from other Dactylopodida. It is known that the PLO may be lost, but PLO‐free strains of paramoebians were never available for molecular studies. Recently, we have described the first species of the genus Paramoeba which has no parasome—Paramoeba aparasomata. In this study, we present a mitochondrial genome of this species, compare it with that of Neoparamoeba pemaquidensis, and analyze the evolutionary dynamics of gene sequences and gene order rearrangements between these species. The mitochondrial genome of P. aparasomata is 46,254 bp long and contains a set of 31 protein‐coding genes, 19 tRNAs, two rRNA genes, and 7 open reading frames. Our results suggest that these two mitochondrial genomes within the genus Paramoeba have rather similar organization and gene order, base composition, codon usage, the composition and structure of noncoding, and overlapping regions.     

Extensive structural variations between mitochondrial genomes of CMS and normal peppers (Capsicum annuum L.) revealed by complete nucleotide sequencing

Background

Cytoplasmic male sterility (CMS) is an inability to produce functional pollen that is caused by mutation of the mitochondrial genome. Comparative analyses of mitochondrial genomes of lines with and without CMS in several species have revealed structural differences between genomes, including extensive rearrangements caused by recombination. However, the mitochondrial genome structure and the DNA rearrangements that may be related to CMS have not been characterized in Capsicum spp.

Results

We obtained the complete mitochondrial genome sequences of the pepper CMS line FS4401 (507,452 bp) and the fertile line Jeju (511,530 bp). Comparative analysis between mitochondrial genomes of peppers and tobacco that are included in Solanaceae revealed extensive DNA rearrangements and poor conservation in non-coding DNA. In comparison between pepper lines, FS4401 and Jeju mitochondrial DNAs contained the same complement of protein coding genes except for one additional copy of an atp6 gene (ψatp6-2) in FS4401. In terms of genome structure, we found eighteen syntenic blocks in the two mitochondrial genomes, which have been rearranged in each genome. By contrast, sequences between syntenic blocks, which were specific to each line, accounted for 30,380 and 17,847 bp in FS4401 and Jeju, respectively. The previously-reported CMS candidate genes, orf507 and ψatp6-2, were located on the edges of the largest sequence segments that were specific to FS4401. In this region, large number of small sequence segments which were absent or found on different locations in Jeju mitochondrial genome were combined together. The incorporation of repeats and overlapping of connected sequence segments by a few nucleotides implied that extensive rearrangements by homologous recombination might be involved in evolution of this region. Further analysis using mtDNA pairs from other plant species revealed common features of DNA regions around CMS-associated genes.

Conclusions

Although large portion of sequence context was shared by mitochondrial genomes of CMS and male-fertile pepper lines, extensive genome rearrangements were detected. CMS candidate genes located on the edges of highly-rearranged CMS-specific DNA regions and near to repeat sequences. These characteristics were detected among CMS-associated genes in other species, implying a common mechanism might be involved in the evolution of CMS-associated genes.

Electronic supplementary material

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

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.