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.     

Detection of two different mitochondrial genomes in a gall wasp species,Andricus mairei (Hymenoptera: Cynipoidea: Cynipidae)

Mitochondrial genomes (mitogenomes) have been used widely in comparative and evolutionary genomics, molecular evolution, phylogenetics, and population genetics, but very limited information is available for the family Cynipidae. In this report, we described the mitogenome of Andricus mairei. The mitogenome of A. mairei was 16,514 bp in length and contained a typical set of 37 genes. Two control regions (CRs) were detected, one being a partial reverse repeat of the other. In a comparison with the putative ancestral mitogenome, gene rearrangements were found in transfer RNA (tRNA) genes, protein-coding genes and ribosomal RNA (rRNA) genes. Consistent with other Cynipidae species, the gene rearrangement of A. mairei had four obvious characteristics: trnE and trnF had inverted and swapped positions; rrnL and rrnS genes had moved into the cob–nad1 junction; a novel tRNA gene cluster trnL1–trnI–trnL2–trnW–trnM–trnQ had been formed between nad1 and nad2; and trnV had inverted and moved to the nad2–cox1 gene junction. Furthermore, A. mairei had two types of mitochondrial circular DNA molecules. Type II differed from type I in an inverted rearrangement of a large fragment of 3349 nucleotides, including two CRs and two rRNA genes.     

The complete mitochondrial genome of Neobenedenia melleni (Platyhelminthes: Monogenea): mitochondrial gene content,arrangement and composition compared with two Benedenia species

The complete mitochondrial (mt) genome sequences of Neobenedenia melleni were determined and compared with those of Benedenia seriolae and B. hoshinai. This circular genome comprises 13,270 bp and includes all 36 typical mt genes found in flatworms. Total AT content of N. melleni is 75.9 %. ATG is the most common start codon, while nad4L is initiated by GTG. All protein-coding genes are predicted to terminate with TAG and TAA. N. melleni has the trnR with a TCG anticodon, which is the same to B. seriolae but different from B. hoshinai (ACG). The mt gene arrangement of N. melleni is similar to that of B. seriolae and B. hoshinai with the exception of three translocations (trnF, trnT and trnG). The overlapped region between nad4L and nad4 was found in the N. melleni mt genome, which was also reported for the published Gyrodactylus species, but it was not found in those of B. seriolae and B. hoshinai, which are non-coding regions instead. The present study provides useful molecular characters for species or strain identification and systematic studies of this parasite.     

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.     

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.     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

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

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

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

The first complete mitochondrial genome of a parasitic isopod supports Epicaridea Latreille, 1825 as a suborder and reveals the less conservative genome of isopods

The complete mitochondrial genome sequence of the holoparasitic isopod Gyge ovalis (Shiino, 1939) has been determined. The mitogenome is 14,268 bp in length and contains 34 genes: 13 protein-coding genes, two ribosomal RNA, 19 tRNA and a control region. Three tRNA genes (trnE, trnI and trnS1) are missing. Most of the tRNA genes show secondary structures which derive from the usual cloverleaf pattern except for trnC which is characterised by the loss of the DHU-arm. Compared to the isopod ground pattern and Eurydice pulchra Leach, 1815 (suborder Cymothoida Wägele, 1989), the genome of G. ovalis shows few differences, with changes only around the control region. However, the genome of G. ovalis is very different from that of non-cymothoidan isopods and reveals that the gene order evolution in isopods is less conservative compared to other crustaceans. Phylogenic trees were constructed using maxiumum likelihood and Bayesian inference analyses based on 13 protein-coding genes. The results do not support the placement of G. ovalis with E. pulchra and Bathynomus sp. in the same suborder; rather, G. ovalis appears to have a closer relationship to Ligia oceanica (Linnaeus, 1767), but this result suggests a need for more data and further analysis. Nevertheless, these results cast doubt that Epicaridea Latreille, 1825 can be placed as an infraorder within the suborder Cymothoida, and Epicaridea appears to also deserve subordinal rank. Further development of robust phylogenetic relationships across Isopoda Latreille, 1817 will require more genetic data from a greater diversity of taxa belonging to all isopod suborders.     

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

Complete mitochondrial genome of Bactrocera ritsemai (Insecta: Tephritidae) and phylogenetic relationship with its congeners and related tephritid taxa

Bactrocera ritsemai is a dacine fruit fly found in Indonesia. We report here the complete mitogenome of this fruit fly from Lombok, Indonesia determined by Illumina MiSeq sequencing and its phylogenetic relationship with its congeners and related tephritid taxa. The whole mitogenome of B. ritsemai had a total length of 15,927 bp, comprising 37 genes – 13 protein-coding genes (PCGs), 2 ribosomal ribonucleic acid (rRNA) and 22 transfer ribonucleic acid (tRNA) genes – and a control region (D-loop). Of the PCGs, 6 (atp6, cob, cox2, cox3, nad4, nad4l) had ATG start codon, 4 (nad2, nad3, nad5, nad6) had ATT, and one each had ATA (nad1), GTG (atp8) and TCG (cox1). Seven PCGs (atp6, atp8, cox2, cox3, nad2, nad4l, nad6) had TAA stop codon, 3 (cob, nad3, nad4) had TAG, and 3 had incomplete stop codon (cox1 – TA; nad1, nad5 – T). The TΨC-loop of tRNA was absent in trnF while trnS1 lacked the DHU-loop. Phylogenetic analysis based on 15 mt-genes (13 PCGs + 2 rRNA genes) indicated B. ritsemai forming a sister group with B. umbrosa and the subgenus Bactrocera was monophyletic. The genera Bactrocera and Zeugodacus were monophyletic while the subfamilies Dacinae and Tephritinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

Involvement of two differenturf-s related mitochondrial sequences in the molecular evolution of the CMS-specificS-Pcf locus in petunia

In petunia, a mitochondrial (mt) locus,S-Pcf, has been found to be strongly associated with cytoplasmic male sterility (CMS). TheS-Pcf locus consists of three open reading frames (ORF) that are co-transcribed. The first ORF,Pcf, contains parts of theatp9 andcoxII genes and an unidentified reading frame,urf-s. The second and third ORFs contain NADH dehydrogenase subunit 3 (nad3) and ribosomal protein S12 (rps12) sequences, respectively. Thenad3 andrps12 sequences included in theS-Pcf locus are identical to the corresponding sequences on the mt genome of fertile petunia. In both CMS and fertile petunia, only a single copy ofnad3 andrps12 has been detected on the physical map of the main mt genome. The origin of theurf-s sequence and the molecular events leading to the formation of the chimericS-Pcf locus are not known. This paper presents evidence indicating that two different mt sequences, related tourf-s and found in fertile petunia lines (orf-h and Rf-1), might have been involved in the molecular evolution of theS-Pcf locus. Southern analysis of mtDNA derived from both fertile and sterile petunia plants suggests that one of theseurf-s related sequences (showing 100% homology tourf-s and termedorf-h) is located on a sublimon. An additional, low-homologyurf-s related sequence (Rf-1) is shown to be located on the main mt genome 5′ to thenad3 gene. It is, thus, suggested that the sequence of events leading to the generation of theS-Pcf locus might have involved introduction of theorf-h sequence, via homologous recombination, into the main mt genome 5′ tonad3 at the region where the Rf-1 sequence is located.     

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: An ultra-large metazoan mitochondrial genome

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8–10 length variants. The motifs of the HTR arrays are about 353–362 bp and the number of repeats ranges from 1 to 11.     

The complete nucleotide sequence of the mitochondrial genome of sugar beet (Beta vulgaris L.) reveals a novel gene for tRNA(Cys)(GCA)

We determined the complete nucleotide sequence of the mitochondrial genome of an angiosperm, sugar beet (Beta vulgaris cv TK81-O). The 368 799 bp genome contains 29 protein, five rRNA and 25 tRNA genes, most of which are also shared by the mitochondrial genome of Arabidopsis thaliana, the only other completely sequenced angiosperm mitochondrial genome. However, four genes identified here (namely rps13, trnF-GAA, ccb577 and trnC2-GCA) are missing in Arabidopsis mitochondria. In addition, four genes found in Arabidopsis (ccb228, rpl2, rpl16 and trnY2-GUA) are entirely absent in sugar beet or present only in severely truncated form. Introns, duplicated sequences, additional reading frames and inserted foreign sequences (chloroplast, nuclear and plasmid DNA sequences) contribute significantly to the overall size of the sugar beet mitochondrial genome. Nevertheless, 55.6% of the genome has no obvious features of information. We identified a novel tRNA^Cys gene (trnC2-GCA) which shows no sequence homology with any tRNA^Cys genes reported so far in higher plants. Intriguingly, this tRNA gene is actually transcribed into a mature tRNA, whereas the native tRNA^Cys gene (trnC1-GCA) is most likely a pseudogene.     

The mitochondrial genome of <Emphasis Type="Italic">Euphausia superba</Emphasis> (Prydz Bay) (Crustacea: Malacostraca: Euphausiacea) reveals a novel gene arrangement and potential molecular markers

Euphausiid krill are dominant organisms in the zooplankton population and play a central role in marine ecosystems. In this paper, we described the gene organization, gene rearrangement and codon usage in the mitochondrial genome of Euphausia superba Dana 1852 (sampling from Prydz Bay, PB). The mitochondrial genome of E. superba is more than 15,498 bp in length (partial non-coding region was not determined). Translocation of four tRNAs (trnL ₁ , trnL ₂ , trnW and trnI) and duplication of one tRNA (trnN) were founded in the mitochondrial genome of E. superba when comparing its genome with the pancrustacean ground pattern. To investigate the phylogenetic relationship within Malacostraca, phylogenetic trees based on currently available malacostracan mitochondrial genomes were built with the maximum likelihood and the Bayesian models. All analyses based on nucleotide and amino acid data strongly support the monophyly of Stomatopoda, Penaeidae, Caridea, and Brachyura, which is consistent with previous research. However, the taxonomic position of Euphausiacea within Malacostraca is unstable. From comparing the mitochondrial genome between E. superba (PB) and E. superba (sampling from Weddell Sea, WS), we found that nad2 gene contains maximal variation with 61 segregating sites, following by nad5 gene which has 12 segregating sites. Thus, nad2 and nad5 genes may be used as potential molecular markers to study the inherit diversity among different E. superba groups, which would be helpful to the exploitation and management of E. superba resources.     

The First Mitochondrial Genomes of Antlion (Neuroptera: Myrmeleontidae) and Split-footed Lacewing (Neuroptera: Nymphidae), with Phylogenetic Implications of Myrmeleontiformia

In the holometabolous insect order Neuroptera (lacewings), the cosmopolitan Myrmeleontidae (antlions) are the most species-rich family, while the closely related Nymphidae (split-footed lacewings) are a small endemic family from the Australian-Malesian region. Both families belong to the suborder Myrmeleontiformia, within which controversial hypotheses on the interfamilial phylogenetic relationships exist. Herein, we describe the complete mitochondrial (mt) genomes of an antlion (Myrmeleon immanis Walker, 1853) and a split-footed lacewing (Nymphes myrmeleonoides Leach, 1814), representing the first mt genomes for both families. These mt genomes are relatively small (respectively composed of 15,799 and 15,713 bp) compared to other lacewing mt genomes, and comprise 37 genes (13 protein coding genes, 22 tRNA genes and two rRNA genes). The arrangement of these two mt genomes is the same as in most derived Neuroptera mt genomes previously sequenced, specifically with a translocation of trnC. The start codons of all PCGs are started by ATN, with an exception of cox1, which is ACG in the M. immanis mt genome and TCG in N. myrmeleonoides. All tRNA genes have a typical clover-leaf structure of mitochondrial tRNA, with the exception of trnS1(AGN). The secondary structures of rrnL and rrnS are similar with those proposed insects and the domain I contains nine helices rather than eight helices, which is common within Neuroptera. A phylogenetic analysis based on the mt genomic data for all Neuropterida sequenced thus far, supports the monophyly of Myrmeleontiformia and the sister relationship between Ascalaphidae and Myrmeleontidae.     

The rearranged mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae), a parasitoid wasp of mantis

The complete sequence of the mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae) is described. The mitogenome was 15,845 bp in size, and contained typical sets of mitochondrial genes. The base composition of the Podagrion sp. mitogenome was also biased toward A + T bases (81.8%). The mitochondrial genome of Podagrion sp. has a weak AT skew (0.07) and a strong GC skew (?0.26). Podagrion sp. exhibits a novel rearrangement compared with the ancestral order, including six protein-coding genes (nad3, cox3, atp6, atp8, cox2 and cox1), which have inverted to the minor strand from the major strand. The A + T-rich region of Podagrion sp., which is located between trnN and trnI, have five tandem repeats. The apomorphic rearrangements, including the conserved block “cox3-atp6-atp8-cox2-cox1-nad5-nad4-nad4l-nad6-cob” and the special locations of trnV and trnA, were mapped onto the phylogeny of Proctotrupomorpha.     

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

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

Complete mitochondrial genome of Dacus conopsoides (Insecta: Tephritidae) with tRNA gene duplication and molecular phylogeny of Dacini tribe

To date there is only a single report on the complete mitochondrial genome of the Dacus fruit flies. We report here the whole mitogenome of Dacus conopsoides with first report of tRNA gene duplication in tephritid fruit flies determined using next-generation sequencing and discuss the molecular phylogeny of Dacini tribe. It had a total length of 15,852 bp, comprising 13 protein coding genes, 2 rRNA genes, 23 tRNA genes, and a non-coding region (A + T-rich control region). The 65-bp trnF gene was duplicated, and the 68-bp trnE gene was partially duplicated resulting in a 31-bp pseudogene. The cloverleaf structure for trnN, trnH, and trnF lacked the TΨC-loop, while trnS lacked the D-stem. The start codons for the protein coding genes included 6 ATG, 3 ATC, 2 ATA, and 1 each of ATT and TCG. Seven PCGs had TAA stop codon, two had TAG and four had incomplete T stop codon. Molecular phylogeny based on 15 mt-genes (13 PCGs +2 rRNA genes) and 30 taxa of Tephritidae indicated D. conopsoides forming a monophyletic sister group with D. longicornis supported by high bootstrap value. The lineage containing also the monophyletic genus Zeugodacus. The Dacini and Ceratitidini tribes of the subfamily Dacinae were monophyletic but the subfamilies Dacinae and Trypetinae were paraphyletic. A broader taxa sampling of the Tephritidae is needed to better elucidate the phylogenetics and systematics of the tribes and subfamilies of tephritid fruit flies.     

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