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 Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

Complete mitochondrial genome sequence of the Arctic gammarid,Onisimus nanseni (Crustacea; Amphipoda): Novel gene structures and unusual control region features

To analyze the mitogenome of the amphipod Onisimus nanseni, we amplified the complete mitogenome of O. nanseni using long-PCR and genome walking techniques. The mitogenome of O. nanseni is circular and contains all the typical mt genes (2 rRNAs, 22 tRNAs, and 13 protein-coding genes). It has two peculiar non-coding regions of 148 bp and 194 bp. The latter can be involved in replication and termination processes. The total length of the pooled protein-coding, rRNA, and tRNA genes is shorter than those of other crustaceans. In addition, the intergenic spacers of the O. nanseni mitogenome are considerably shorter in length than those of other crustaceans. Fourteen adjacent genes overlap, resulting in a compact mitogenomic structure. In the O. nanseni mitogenome, the AT composition is elevated, particularly in the control regions (78.9% AT), as has been demonstrated for two other amphipods. The tRNA order is highly rearranged compared to other arthropod mitogenomes, but the order of protein-coding genes and rRNAs is largely conserved. The gene cluster between the CO1 and CO3 genes is completely conserved among all amphipods compared. This provides insights into the evolution and gene structures of crustacean mitochondrial genomes, particularly in amphipods.     

The Complete Maternally and Paternally Inherited Mitochondrial Genomes of the Endangered Freshwater Mussel Solenaia carinatus (Bivalvia: Unionidae) and Implications for Unionidae Taxonomy

Doubly uniparental inheritance (DUI) is an exception to the typical maternal inheritance of mitochondrial (mt) DNA in Metazoa, and found only in some bivalves. In species with DUI, there are two highly divergent gender-associated mt genomes: maternal (F) and paternal (M), which transmit independently and show different tissue localization. Solenaia carinatus is an endangered freshwater mussel species exclusive to Poyang Lake basin, China. Anthropogenic events in the watershed greatly threaten the survival of this species. Nevertheless, the taxonomy of S. carinatus based on shell morphology is confusing, and the subfamilial placement of the genus Solenaia remains unclear. In order to clarify the taxonomic status and discuss the phylogenetic implications of family Unionidae, the entire F and M mt genomes of S. carinatus were sequenced and compared with the mt genomes of diverse freshwater mussel species. The complete F and M mt genomes of S. carinatus are 16716 bp and 17102 bp in size, respectively. The F and M mt genomes of S. carinatus diverge by about 40% in nucleotide sequence and 48% in amino acid sequence. Compared to F counterparts, the M genome shows a more compact structure. Different gene arrangements are found in these two gender-associated mt genomes. Among these, the F genome cox2-rrnS gene order is considered to be a genome-level synapomorphy for female lineage of the subfamily Gonideinae. From maternal and paternal mtDNA perspectives, the phylogenetic analyses of Unionoida indicate that S. carinatus belongs to Gonideinae. The F and M clades in freshwater mussels are reciprocal monophyly. The phylogenetic trees advocate the classification of sampled Unionidae species into four subfamilies: Gonideinae, Ambleminae, Anodontinae, and Unioninae, which is supported by the morphological characteristics of glochidia.     

Comparative Mitogenomics of the Assassin Bug Genus Peirates (Hemiptera: Reduviidae: Peiratinae) Reveal Conserved Mitochondrial Genome Organization of P. atromaculatus,P. fulvescens and P. turpis

In this study, we sequenced four new mitochondrial genomes and presented comparative mitogenomic analyses of five species in the genus Peirates (Hemiptera: Reduviidae). Mitochondrial genomes of these five assassin bugs had a typical set of 37 genes and retained the ancestral gene arrangement of insects. The A+T content, AT- and GC-skews were similar to the common base composition biases of insect mtDNA. Genomic size ranges from 15,702 bp to 16,314 bp and most of the size variation was due to length and copy number of the repeat unit in the putative control region. All of the control region sequences included large tandem repeats present in two or more copies. Our result revealed similarity in mitochondrial genomes of P. atromaculatus, P. fulvescens and P. turpis, as well as the highly conserved genomic-level characteristics of these three species, e.g., the same start and stop codons of protein-coding genes, conserved secondary structure of tRNAs, identical location and length of non-coding and overlapping regions, and conservation of structural elements and tandem repeat unit in control region. Phylogenetic analyses also supported a close relationship between P. atromaculatus, P. fulvescens and P. turpis, which might be recently diverged species. The present study indicates that mitochondrial genome has important implications on phylogenetics, population genetics and speciation in the genus Peirates.     

Decoding and analysis of organelle genomes of Indian tea (Camellia assamica) for phylogenetic confirmation

The NCBI database has >15 chloroplast (cp) genome sequences available for different Camellia species but none for C. assamica. There is no report of any mitochondrial (mt) genome in the Camellia genus or Theaceae family. With the strong believes that these organelle genomes can play a great tool for taxonomic and phylogenetic analysis, we successfully assembled and analyzed cp and mt genome of C. assamica. We assembled the complete mt genome of C. assamica in a single circular contig of 707,441 bp length comprising of a total of 66 annotated genes, including 35 protein-coding genes, 29 tRNAs and two rRNAs. The first ever cp genome of C. assamica resulted in a circular contig of 157,353 bp length with a typical quadripartite structure. Phylogenetic analysis based on these organelle genomes showed that C. assamica was closely related to C. sinensis and C. leptophylla. It also supports Caryophyllales as Superasterids.     

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.     

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

Background

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

Principal Findings

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

Conclusions/Significance

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

Fragmented mitochondrial genomes are present in both major clades of the blood-sucking lice (suborder Anoplura): evidence from two Hoplopleura rodent lice (family Hoplopleuridae)

Background

The suborder Anoplura contains 540 species of blood-sucking lice that parasitize over 840 species of eutherian mammals. Fragmented mitochondrial (mt) genomes have been found in the lice of humans, pigs, horses and rats from four families: Pediculidae, Pthiridae, Haematopinidae and Polyplacidae. These lice, eight species in total, are from the same major clade of the Anoplura. The mt genomes of these lice consist of 9–20 minichromosomes; each minichromosome is 1.5–4 kb in size and has 1–8 genes. To understand mt genome fragmentation in the other major clade of the Anoplura, we sequenced the mt genomes of two species of rodent lice in the genus Hoplopleura (family Hoplopleuridae).

Results

We identified 28 mt genes on 10 minichromosomes in the mouse louse, Ho. akanezumi; each minichromosome is 1.7–2.7 kb long and has 1–6 genes. We identified 34 mt genes on 11 minichromosomes in the rat louse, Ho. kitti; each minichromosome is 1.8–2.8 kb long and has 1–5 genes. Ho. akanezumi also has a chimeric minichromosome with parts of two rRNA genes and a full-length tRNA gene for tyrosine. These two rodent lice share the same pattern for the distribution of all of the protein-coding and rRNA genes but differ in tRNA gene content and gene arrangement in four minichromosomes. Like the four genera of blood-sucking lice that have been investigated in previous studies, the Hoplopleura species have four minichromosomes that are only found in this genus.

Conclusions

Our results indicate that fragmented mt genomes were present in the most recent common ancestor of the two major clades of the blood-sucking lice, which lived ~75 million years ago. Intra-genus variation in the pattern of mt genome fragmentation is common in the blood-sucking lice (suborder Anoplura) and genus-specific minichromosomes are potential synapomorphies. Future studies should expand into more species, genera and families of blood-sucking lice to explore further the phylogenetic utility of the novel features associated with fragmented mt genomes.

Electronic supplementary material

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

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.     

Sexual Size Dimorphism in the Color Pattern Elements of Two Mimetic Heliconius Butterflies

Sexual dichromatism and sexual dimorphism of body size are reasonably well studied in butterflies. Sexual size dimorphism of color pattern elements, however, is much less explored. The object of this study is Heliconius, a genus of butterflies well known for the coevolution between mate color preferences and mimicry. Given the sexual role of wing coloration, we investigated the existence of sexual size dimorphism in the wing color elements of a mimetic pair—Heliconius erato phyllis Fabricius and Heliconius besckei Ménétriés—and analyzed the allometric patterns of these traits. Correlation between size of elements in the dorsal and ventral wing surfaces were also estimated. In both species, three out of four elements were larger in males, but the non-dimorphic element was not the same. With regard to the allometric patterns, our most important finding was that smaller males of one species have proportionally larger yellow bars. This is the first study specifically concerning quantitative sexual dimorphism in the coloration of this well-known genus of butterflies and it opens new prospects to investigate sex-related natural selection and sexual selection of color pattern elements.     

The First Mitochondrial Genome for the Fishfly Subfamily Chauliodinae and Implications for the Higher Phylogeny of Megaloptera

Megaloptera are a basal holometabolous insect order with larvae exclusively predacious and aquatic. The evolutionary history of Megaloptera attracts great interest because of its antiquity and important systematic status in Holometabola. However, due to the difficulties identifying morphological apomorphies for the group, controversial hypotheses on the monophyly and higher phylogeny of Megaloptera have been proposed. Herein, we describe the complete mitochondrial (mt) genome of a fishfly species, Neochauliodes punctatolosus Liu & Yang, 2006, representing the first mt genome of the subfamily Chauliodinae. A phylogenomic analysis was carried out based on the mt genomic sequences of 13 mt protein-coding genes (PCGs) and two rRNA genes of nine Neuropterida species, comprising all three orders of Neuropterida and all families and subfamilies of Megaloptera. Both maximum likelihood and Bayesian inference analyses highly support the monophyly of Megaloptera, which was recovered as the sister of Neuroptera. Within Megaloptera, the sister relationship between Corydalinae and Chauliodinae was corroborated. The divergence time estimation suggests that stem lineage of Neuropterida and Coleoptera separated in the Early Permian. The interordinal divergence within Neuropterida might have occurred in the Late Permian.     

Defense mechanisms against herbivory in Picea: sequence evolution and expression regulation of gene family members in the phenylpropanoid pathway

Background

The insect order Neuroptera encompasses more than 5,700 described species. To date, only three neuropteran mitochondrial genomes have been fully and one partly sequenced. Current knowledge on neuropteran mitochondrial genomes is limited, and new data are strongly required. In the present work, the mitochondrial genome of the ascalaphid owlfly Libelloides macaronius is described and compared with the known neuropterid mitochondrial genomes: Megaloptera, Neuroptera and Raphidioptera. These analyses are further extended to other endopterygotan orders.

Results

The mitochondrial genome of L. macaronius is a circular molecule 15,890 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. The gene order of this newly sequenced genome is unique among Neuroptera and differs from the ancestral type of insects in the translocation of trnC. The L. macaronius genome shows the lowest A+T content (74.50%) among known neuropterid genomes. Protein-coding genes possess the typical mitochondrial start codons, except for cox1, which has an unusual ACG. Comparisons among endopterygotan mitochondrial genomes showed that A+T content and AT/GC-skews exhibit a broad range of variation among 84 analyzed taxa. Comparative analyses showed that neuropterid mitochondrial protein-coding genes experienced complex evolutionary histories, involving features ranging from codon usage to rate of substitution, that make them potential markers for population genetics/phylogenetics studies at different taxonomic ranks. The 22 tRNAs show variable substitution patterns in Neuropterida, with higher sequence conservation in genes located on the α strand. Inferred secondary structures for neuropterid rrnS and rrnL genes largely agree with those known for other insects. For the first time, a model is provided for domain I of an insect rrnL. The control region in Neuropterida, as in other insects, is fast-evolving genomic region, characterized by AT-rich motifs.

Conclusions

The new genome shares many features with known neuropteran genomes but differs in its low A+T content. Comparative analysis of neuropterid mitochondrial genes showed that they experienced distinct evolutionary patterns. Both tRNA families and ribosomal RNAs show composite substitution pathways. The neuropterid mitochondrial genome is characterized by a complex evolutionary history.     

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.     

Mitochondrial genomic comparisons of the subterranean termites from the Genus Reticulitermes (Insecta: Isoptera: Rhinotermitidae).

Termites of the genus Reticulitermes are some of the most significant pests of structural timber and tree farming in the northern hemisphere, causing losses in the billions of dollars annually because of direct damage and termite control costs. This group has been frequently targeted for population genetic, phylogenetic, and species limit studies, most of which use mitochondrial (mt) genes; however, only a small fraction of the genome has been sequenced. The entire mt genome was sequenced for the eastern North American members of Reticulitermes: R. flavipes, R. santonensis, R. virginicus, and R. hageni. The mt genome has the same gene content and organization as that found in most insect species; however, the nucleotide composition and skew are highly biased (AT% low, strong A- and C-skew). Both the protein-coding and transfer RNA genes show high absolute levels of nucleotide substitution, suggesting that the high rates of mutation within Reticulitermes inferred from analyses of single mt genes are a general characteristic of the entire mt genome. The AT-rich or control region has a remarkable structure not previously observed in insect mt genomes. The majority of the control region is made up of 2 sets of repeat units, typically with 2 full and 1 partial copies of both the A (or small; 186 bp) and B (or large; 552 bp) repeats. The partial repeat units overlap by 36 bp. The size, location, and degree of overlap for the partial repeat units correspond to highly conserved stem/loop structures within the repeat units, suggesting that these structures are involved in the replication-mediated processes that govern repeat-unit evolution within mt genomes. Finally, molecular variation within the mt gene regions was compared with previous regions used in molecular diagnostics or phylogenetics of Reticulitermes. High numbers of single nucleotide polymorphisms were found in each of the mt genes, and some of the highest variability was found in gene regions that have not previously been investigated in this group. The whole mt genome sequence can thus be used to predict useful regions for future investigation.     

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.     

The 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     

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.