Evolution of the tRNA gene family in mitochondrial genomes of five Meretrix clams (Bivalvia,Veneridae)

In contrast to the extreme conservation of nuclear-encoded tRNAs, organization of the mitochondrial (mt) tRNA gene family in invertebrates is highly dynamic and rapidly evolving. While gene duplication and loss, gene isomerism, recruitment, and rearrangements have occurred sporadically in several invertebrate lineages, little is known regarding the pattern of their evolution. Comparisons of invertebrate mt genomes at a generic level can be extremely helpful in investigating evolutionary patterns of variation, as intermediate stages of the process may be identified. Variation of mitochondrial tRNA organization among Meretrix clams provides good materials to investigate mt tRNA evolution. We characterized the complete mt genome of the lyrate Asiatic hard clam Meretrix lyrata, re-annotated tRNAs of four previously sequenced Meretrix clams, and undertook an intensive comparison of tRNA gene families in these clams. Our results 1) provide evidence that the commonly observed duplication of trnM may have occurred independently in different bivalve lineages and, based on the higher degree of trnM gene similarity, may have occurred more recently than expected; 2) suggest that “horizontal” evolution may have played an important role in tRNA gene family evolution based on frequent gene duplications and gene recruitment events; and 3) reveal the first case of isoacceptor “vertical” tRNA gene recruitment (VTGR) and present the first clear evidence that VTGR allows rapid evolution of tRNAs. We identify the trnS^− UCR gene in Meretrix clams, previously considered missing in this lineage, and speculate that trnS^− UCR lacking the D-arm in both M. lyrata and Meretrix lamarckii may represent the ancestral status. Phylogenetic analysis based on 13 concatenate protein-coding genes provided opportunities to detect rapidly evolved tRNA genes via VTGR and gene isomerism processes. This study suggests that evolution of the mt tRNA gene family in bivalves is more complex than previously thought and that comparison of several congeneric species is a useful strategy in investigating evolutionary patterns and dynamics of mt tRNA genes.     

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 complete mitochondrial genome of Arctic Calanus hyperboreus (Copepoda,Calanoida) reveals characteristic patterns in calanoid mitochondrial genome

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

Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum

Bunostomum trigonocephalum and Bunostomum phlebotomum are blood-feeding hookworms of sheep and cattle, causing considerable economic losses to the live stock industries. Studying genetic variability within and among hookworm populations is critical to addressing epidemiological and ecological questions. Mitochondrial (mt) DNA is known to provide useful markers for investigations of population genetics of hookworms, but mt genome sequence data are scant. In the present study, the complete mitochondrial DNA (mtDNA) sequences of the sheep and goat hookworm B. trigonocephalum were determined for the first time, and the mt genome of B. phlebotomum from yak in China was also sequenced for comparative analyses of their gene contents and genome organizations. The lengths of mt DNA sequences of B. trigonocephalum sheep isolate, B.trigonocephalum goat isolate and B. phlebotomum China yak isolate were 13,764 bp, 13,771 bp and 13,803 bp in size, respectively. The identity of the mt genomes was 99.7% between B. trigonocephalum sheep isolate and B. trigonocephalum goat isolate. The identity of B. phlebotomum China yak isolate mt genomes was 85.3% with B. trigonocephalum sheep isolate, and 85.2% with B. trigonocephalum goat isolate. All the mt genes of the two hookworms were transcribed in the same direction and gene arrangements were consistent with those of the GA3 type, including 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes, but lacking ATP synthetase subunit 8 gene. The mt genomes of B. trigonocephalum and B. phlebotomum were similar to prefer bases A and T, the contents of A + T are 76.5% (sheep isolate), 76.4% (goat isolate) and 76.9% (China yak isolate), respectively. Phylogenetic relationships reconstructed using concatenated amino acid sequences of 12 protein-coding genes with three methods (maximum likelihood, Bayesian inference and neighbor joining) revealed that the B. trigonocephalum and B. phlebotomum represent distinct but closely-related species. These data provide novel and useful genetic markers for studying the systematics, and population genetics of the two ruminant hookworms.     

Molecular cloning and functional analyses of glutathione peroxidase homologous genes from Chlorella sp. NJ-18

The extreme variability of the mitochondrial (mito) genomes of bivalves makes it difficult to understand their evolutionary dynamics, given that species from different families do not share comparable features. We compared the mitogenomes from four Paphia clams (three of them were firstly sequenced) and found that mitogenome reorganization among the four congeneric species is not random but follows phylogenetic trends. Start/stop codon variations are species-correlated rather than gene-correlated, and bear useful phylogenetic information. Unique start/stop codon usage in P. euglypta and A+T content in P. amabilis indicates that these mitogenome-level characters, usually considered to be conservative features in other lineages, may not be phylogenetically evolved, but may have evolved via species-specific mitogenomic maintenance mechanisms. Variable divergence of two trnM genes in different lineages may demonstrate differences in mechanisms by which paralogous trnM genes are maintained. Sequence alignment analysis indicates that the VNTRs in the four mitogenomes have a common origin. The rationale of the subgenus Neotapes Kuroda and Habe, 1971 was supported by evidence from morphological characters, mitogenomic features, as well as phylogenetic analyses using cox1 and rrnS genes. The data suggest that the taxonomic basis of the subgenus should be “smooth surface” but not “undulated lines,” and P. textile should be classified to the Neotapes subgenus.     

Comparison of complete mitochondrial genomes of pine wilt nematode Bursaphelenchus xylophilus and Bursaphelenchus mucronatus (Nematoda: Aphelenchoidea) and development of a molecular tool for species identification

We determined the complete mitochondrial genome sequences for Bursaphelenchus mucronatus, one species of pinewood nematode. The genome is a circular-DNA molecule of 14,583 bp (195 bp smaller than its congener Bursaphelenchus xylophilus) and contains 12 protein-coding genes (lacking atp8), 22 tRNA genes, and 2 rRNA genes encoded in the same direction, consistent with most other nematodes. Based on sequence comparison of mtDNA genomes, we developed a PCR-based molecular assay to differentiate B. xylophilus (highly pathogenic) and B. mucronatus (relatively less virulent) using species-specific primers. The molecular identification system employs multiplex-PCR and is very effective and reliable for discriminating these Bursaphelenchus species, which are economically important, but difficult to distinguish based on morphology. The comparison of the mitochondrial genomes and molecular identification system of the two species of Bursaphelenchus spp. should provide a rich source of genetic information to support the effective control and management (quarantine) of the pine wilt disease caused by pinewood nematodes.     

The complete mitochondrial genome of Pauropus longiramus (Myriapoda: Pauropoda): implications on early diversification of the myriapods revealed from comparative analysis

Myriapods are among the earliest arthropods and may have evolved to become part of the terrestrial biota more than 400 million years ago. A noticeable lack of mitochondrial genome data from Pauropoda hampers phylogenetic and evolutionary studies within the subphylum Myriapoda. We sequenced the first complete mitochondrial genome of a microscopic pauropod, Pauropus longiramus (Arthropoda: Myriapoda), and conducted comprehensive mitogenomic analyses across the Myriapoda. The pauropod mitochondrial genome is a circular molecule of 14,487 bp long and contains the entire set of thirty-seven genes. Frequent intergenic overlaps occurred between adjacent tRNAs, and between tRNA and protein-coding genes. This is the first example of a mitochondrial genome with multiple intergenic overlaps and reveals a strategy for arthropods to effectively compact the mitochondrial genome by overlapping and truncating tRNA genes with neighbor genes, instead of only truncating tRNAs. Phylogenetic analyses based on protein-coding genes provide strong evidence that the sister group of Pauropoda is Symphyla. Additionally, approximately unbiased (AU) tests strongly support the Progoneata and confirm the basal position of Chilopoda in Myriapoda. This study provides an estimation of myriapod origins around 555 Ma (95% CI: 444-704 Ma) and this date is comparable with that of the Cambrian explosion and candidate myriapod-like fossils. A new time-scale suggests that deep radiations during early myriapod diversification occurred at least three times, not once as previously proposed. A Carboniferous origin of pauropods is congruent with the idea that these taxa are derived, rather than basal, progoneatans.     

Gene recruitment--a common mechanism in the evolution of transfer RNA gene families

The evolution of alloacceptor transfer RNAs (tRNAs) has been traditionally thought to occur vertically and reflect the evolution of the genetic code. Yet there have been several indications that a tRNA gene could evolve horizontally, from a copy of an alloacceptor tRNA gene in the same genome. Earlier, we provided the first unambiguous evidence for the occurrence of such "tRNA gene recruitment" in nature--in the mitochondrial (mt) genome of the demosponge Axinella corrugata. Yet the extent and the pattern of this process in the evolution of tRNA gene families remained unclear. Here we analyzed tRNA genes from 21 mt genomes of demosponges as well as nuclear genomes of rhesus macaque, chimpanzee and human. We found four new cases of alloacceptor tRNA gene recruitment in mt genomes and eleven cases in the nuclear genomes. In most of these cases we observed a single nucleotide substitution at the middle position of the anticodon, which resulted in the change of not only the tRNA's amino-acid identity but also the class of the amino-acyl tRNA synthetases (aaRSs) involved in amino-acylation. We hypothesize that the switch to a different class of aaRSs may have prevented the conflict between anticodon and amino-acid identities of recruited tRNAs. Overall our results suggest that gene recruitment is a common phenomenon in tRNA multigene family evolution and should be taken into consideration when tRNA evolutionary history is reconstructed.     

Systematic investigation of Amphioxus (Branchiostoma floridae) microRNAs

Ascaris lumbricoides and Ascaris suum are parasitic nematodes living in the small intestine of humans and pigs, and can cause the disease ascariasis. For long, there has been controversy as to whether the two ascaridoid taxa represent the same species due to their significant resemblances in morphology. However, the complete mitochondrial (mt) genome data have been lacking for A. lumbricoides in spite of human and animal health significance and socio-economic impact globally of these parasites. In the present study, we sequenced the complete mt genomes of A. lumbricoides and A. suum (China isolate), which was 14,303 bp and 14,311 bp in size, respectively. The identity of the mt genomes was 98.1% between A. lumbricoides and A. suum (China isolate), and 98.5% between A. suum (China isolate) and A. suum (USA isolate). Both genomes are circular, and consist of 36 genes, including 12 genes for proteins, 2 genes for rRNA and 22 genes for tRNA, which are consistent with that of all other species of ascaridoid studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T (71.7% for A. lumbricoides and 71.8% for A. suum). The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. Phylogenetic analyses of A. lumbricoides and A. suum using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (Bayesian analysis, maximum likelihood and maximum parsimony) all clustered in a clade with high statistical support, indicating that A. lumbricoides and A. suum was very closely related. These mt genome data and the results provide some additional genetic evidence that A. lumbricoides and A. suum may represent the same species. The mt genome data presented in this study are also useful novel markers for studying the molecular epidemiology and population genetics of Ascaris.     

Chimeric mitochondrial minichromosomes of the human body louse, Pediculus humanus: evidence for homologous and non-homologous recombination

The mitochondrial (mt) genome of the human body louse, Pediculus humanus, consists of 18 minichromosomes. Each minichromosome is 3 to 4 kb long and has 1 to 3 genes. There is unequivocal evidence for recombination between different mt minichromosomes in P. humanus. It is not known, however, how these minichromosomes recombine. Here, we report the discovery of eight chimeric mt minichromosomes in P. humanus. We classify these chimeric mt minichromosomes into two groups: Group I and Group II. Group I chimeric minichromosomes contain parts of two different protein-coding genes that are from different minichromosomes. The two parts of protein-coding genes in each Group I chimeric minichromosome are joined at a microhomologous nucleotide sequence; microhomologous nucleotide sequences are hallmarks of non-homologous recombination. Group II chimeric minichromosomes contain all of the genes and the non-coding regions of two different minichromosomes. The conserved sequence blocks in the non-coding regions of Group II chimeric minichromosomes resemble the "recombination repeats" in the non-coding regions of the mt genomes of higher plants. These repeats are essential to homologous recombination in higher plants. Our analyses of the nucleotide sequences of chimeric mt minichromosomes indicate both homologous and non-homologous recombination between minichromosomes in the mitochondria of the human body louse.     

The complete mitochondrial genome of the Chinese hook snout carp Opsariichthys bidens (Actinopterygii: Cypriniformes) and an alternative pattern of mitogenomic evolution in vertebrate

The complete mitochondrial genome sequence of the Chinese hook snout carp, Opsariichthys bidens, was newly determined using the long and accurate polymerase chain reaction method. The 16,611-nucleotide mitogenome contains 13 protein-coding genes, two rRNA genes (12S, 16S), 22 tRNA genes, and a noncoding control region. We use these data and homologous sequence data from multiple other ostariophysan fishes in a phylogenetic evaluation to test hypothesis pertaining to codon usage pattern of O. bidens mitochondrial protein genes as well as to re-examine the ostariophysan phylogeny. The mitochondrial genome of O. bidens reveals an alternative pattern of vertebrate mitochondrial evolution. For the mitochondrial protein genes of O. bidens, the most frequently used codon generally ends with either A or C, with C preferred over A for most fourfold degenerate codon families; the relative synonymous codon usage of G-ending codons is greatly elevated in all categories. The codon usage pattern of O. bidens mitochondrial protein genes is remarkably different from the general pattern found previously in the relatively closely related zebrafish and most other vertebrate mitochondria. Nucleotide bias at third codon positions is the main cause of codon bias in the mitochondrial protein genes of O. bidens, as it is biased particularly in favor of C over A. Bayesian analysis of 12 concatenated mitochondrial protein sequences for O. bidens and 46 other teleostean taxa supports the monophyly of Cypriniformes and Otophysi and results in a robust estimate of the otophysan phylogeny.     

Beyond barcoding: A mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.     

The complete mitochondrial genome of Taeniogonalos taihorina (Bischoff) (Hymenoptera: Trigonalyidae) reveals a novel gene rearrangement pattern in the Hymenoptera

The family Trigonalyidae is considered to be one of the most basal lineages in the suborder Apocrita of Hymenoptera. Here, we determine the first complete mitochondrial genome of the Trigonalyidae, from the species Taeniogonalos taihorina (Bischoff, 1914). This mitochondrial genome is 15,927 bp long, with a high A + T-content of 84.60%. It contains all of the 37 typical animal mitochondrial genes and an A + T-rich region. The orders and directions of all genes are different from those of previously reported hymenopteran mitochondrial genomes. Eight tRNA genes, three protein-coding genes and the A + T-rich region were rearranged, with the dominant gene rearrangement events being translocation and local inversion. The arrangements of three tRNA clusters, trnY–trnM–trnI–trnQ, trnW–trnL2–trnC, and trnH–trnA–trnR–trnN–trnS–trnE–trnF, and the position of the cox1 gene, are novel to the Hymenoptera, even the insects. Six long intergenic spacers are present in the genome. The secondary structures of the RNA genes are normal, except for trnS2, in which the D-stem pairing is absent.     

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A + T content of the complete mitogenome is 80.83%. The AT skew value ([A − T] / [A + T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of “ATAGA” + 19-bp poly T. It also contains a microsatellite-like (TA)₂₆, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.     

Complete nucleotide sequence and gene organization of the mitochondrial genome of Paa spinosa (Anura: Ranoidae)

The mt genome of Paa spinosa (Anura: Ranoidae) is a circular molecule of 18,012 bp in length, containing 38 genes (including an extra copy of tRNA-Met gene). This mt genome is characterized by three distinctive features: a cluster of rearranged tRNA genes (LTPF tRNA gene cluster), a tandem duplication of tRNA-Met gene (Met1 and Met2), and distinct repeat regions at both 5′ and 3′-sides in the control region. Comparing the locations and the sequences of all tRNA-Met genes among Ranoidae, and constructing NJ tree of the nucleotide of those tRNA-Met genes, we suggested a tandem duplication of tRNA-Met gene can be regarded as a synapomorphy of Dicroglossinae. To further investigate the phylogenetic relationships of anurans, phylogenetic analyses (BI, ML and MP) based on the nucleotide dataset and the corresponding amino acid dataset of 11 protein-coding genes (except ND5 and ATP8) arrived at the similar topology.     

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.     

The complete mitochondrial genome of the sycamore lace bug Corythucha ciliata (Hemiptera: Tingidae)

The complete mitochondrial genome of the sycamore lace bug, Corythucha ciliata, was sequenced in this study. It represents the first sequenced mitogenome of family Tingidae in Heteroptera. The mitogenome of C. ciliata is 15,257 bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a large non-coding region. Gene arrangement, nucleotide content, codon usage, and amino acid composition and asymmetry indicate a high degree of conservation with six other species of Cimicomorpha. The 13 PCGs initiated with ATN as the start codon and terminated with TAA, TA or T as stop codon. The evolutionary rate of each PCG was different, among which ATP8 showed the highest rate while ATP6 indicated the lowest rate. The 22 tRNAs genes apparently fold into a typical cloverleaf structure; however, the anticodon (TTC) of trnSer (AGN) differs from other Heteropteran insects. Secondary structure modeling of rRNA genes revealed similarity to other insects, except for two incomplete helices (H1648 and H2735) in lrRNA. The predicted secondary structure of lrRNA indicates 45 helices in six domains, whereas srRNA has 27 helices in three domains. Three potential stem–loops and two tandem repeats (–TCTAAT–) were identified in the A+T-rich region. Phylogenetic analysis indicated that C. ciliata is a sister group to other Heteroptera species based on analysis of the 13 PCGs.     

The complete nucleotide sequence of the mitochondrial genome of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae)

The complete mitochondrial genome of the oriental fruit fly Bactrocera dorsalis s.s. has been sequenced, and is here described and compared with the homologous sequences of Bactrocera oleae and Ceratitis capitata. The genome is a circular molecule of 15,915 bp, and encodes the set of 37 genes generally found in animal mitochondrial genomes. The structure and organization of the molecule is typical and similar to the two closely related species B. oleae and C. capitata, although it presents an interesting case of putative intra-molecular recombination. The relevance of the growing comparative dataset of tephritid complete mitochondrial genomes is discussed in relation to the possibility to develop robust assays for species discrimination in quarantine and agricultural monitoring practices, as well as basic phylogeography/population genetic studies.     

The phylogeny of Ephemeroptera in Pterygota revealed by the mitochondrial genome of Siphluriscus chinensis (Hexapoda: Insecta)

The mayfly species Siphluriscus chinensis (Siphluriscidae) has valuable structures useful for phylogeny reconstruction, given its putative basal position within the Ephemeroptera. Here its nearly complete mitochondrial genome is sequenced. We built phylogenetic trees through multiple analytical strategies with some other insect mitogenomes. Structurally, the obtained mitochondrial genome of S. chinensis is 16,616 bp in length, ¹ containing 37 genes and an extra trnK-like (trnK2 (AAA)) gene. The 12 PCGs start with typical ATN codons, except the nad1 gene which starts with an unnormalized TTG. Like other known mayfly mitogenomes, the strand bias has negative AT-skew and negative GC-skew. Phylogenetically, our topologies suggest that Odonata is the basally diverged clade in Pterygota; Ephemeroptera is the sister group of the Neoptera; and S. chinensis is indeed the most basal mayfly branch.