Complete mitochondrial genomes of two cockroaches, Blattella germanica and Periplaneta americana, and the phylogenetic position of termites

Primitive new termites (Blattodea, Termitoidae) in Cretaceous amber from Myanmar

Mastotermitidae, the first-diverging extant family of termites, has only one relic extant species; however, this family had greater richness during the Mesozoic and Cenozoic eras. Fossil termites from the Cretaceous provide information on the early evolution of termites and the transition between extinct families. Herein, two new Mastotermitidae species found in upper Cretaceous (Cenomanian) Kachin amber are reported. One is a female imago described as Angustitermesreflexusgen. et sp. nov. and assigned to the subfamily Mastotermitinae. The other is Mastotermesreticulatussp. nov., which is described from an isolated forewing. With the comparison especially of the antenna and venation, these new mastotermitids further increase our knowledge of the diversity and morphology of Mastotermitidae during the Mesozoic.

Complete mitochondrial genome of Episymploce splendens (Blattodea: Ectobiidae): A large intergenic spacer and lacking of two tRNA genes

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

Genome size evaluations in cockroaches: new entries

In this paper, we report genome size (GS) values for nine cockroaches (order Blattodea, families Blattidae, Blaberidae and Ectobiidae, ex Blattelidae), three of which are original additions to the ten already present in the GS database: the death’s head roach (Blaberus craniifer), the Surinam cockroach (Pycnoscelus surinamensis) and the Madeira cockroach (Leucophaea maderae). Regarding the American cockroach (Periplaneta americana), the GS database contains two contrasting values (2.72 vs 3.41 pg); likely, the 2.72 pg value is the correct one as it is strikingly similar to our sperm DNA content evaluation (2.80 ± 0.11 pg). Also, we suggest halving the published GS of the Argentine cockroach Blaptica dubia and the spotted cockroach (the gray cockroach) Nauphoeta cinerea discussing i) the occurrence of a correlation between increasing 2n chromosome number and GS within the order Blattodea; and ii) the possible occurrence of a polyploidization phenomenon doubling a basic GS of 0.58 pg of some termite families (superfamily Blattoidea, epifamily Termitoidae).

Increasing 28 mitogenomes of Ephemeroptera, Odonata and Plecoptera support the Chiastomyaria hypothesis with three different outgroup combinations

BACKGROUND

The phylogenetic relationships of Odonata (dragonflies and damselflies) and Ephemeroptera (mayflies) remain unresolved. Different researchers have supported one of three hypotheses (Palaeoptera, Chiastomyaria or Metapterygota) based on data from different morphological characters and molecular markers, sometimes even re-assessing the same transcriptomes or mitochondrial genomes. The appropriate choice of outgroups and more taxon sampling is thought to eliminate artificial phylogenetic relationships and obtain an accurate phylogeny. Hence, in the current study, we sequenced 28 mt genomes from Ephemeroptera, Odonata and Plecoptera to further investigate phylogenetic relationships, the probability of each of the three hypotheses, and to examine mt gene arrangements in these species. We selected three different combinations of outgroups to analyze how outgroup choice affected the phylogenetic relationships of Odonata and Ephemeroptera.

METHODS

Mitochondrial genomes from 28 species of mayflies, dragonflies, damselflies and stoneflies were sequenced. We used Bayesian inference (BI) and Maximum likelihood (ML) analyses for each dataset to reconstruct an accurate phylogeny of these winged insect orders. The effect of outgroup choice was assessed by separate analyses using three outgroups combinations: (a) four bristletails and three silverfish as outgroups, (b) five bristletails and three silverfish as outgroups, or (c) five diplurans as outgroups.

RESULTS

Among these sequenced mitogenomes we found the gene arrangement IMQM in Heptageniidae (Ephemeroptera), and an inverted and translocated tRNA-Ile between the 12S RNA gene and the control region in Ephemerellidae (Ephemeroptera). The IMQM gene arrangement in Heptageniidae (Ephemeroptera) can be explained via the tandem-duplication and random loss model, and the transposition and inversion of tRNA-Ile genes in Ephemerellidae can be explained through the recombination and tandem duplication-random loss (TDRL) model. Our phylogenetic analysis strongly supported the Chiastomyaria hypothesis in three different outgroup combinations in BI analyses. The results also show that suitable outgroups are very important to determining phylogenetic relationships in the rapid evolution of insects especially among Ephemeroptera and Odonata. The mt genome is a suitable marker to investigate the phylogeny of inter-order and inter-family relationships of insects but outgroup choice is very important for deriving these relationships among winged insects. Hence, we must carefully choose the correct outgroup in order to discuss the relationships of Ephemeroptera and Odonata.

Phylosymbiosis across Deeply Diverging Lineages of Omnivorous Cockroaches (Order Blattodea)

The gut microbiome is shaped by both host diet and host phylogeny. However, separating the relative influence of these two factors over long periods of evolutionary time is often difficult. We conducted a 16S rRNA gene amplicon-based survey of the gut microbiome from 237 individuals and 19 species of omnivorous cockroaches from the order Blattodea. The order Blattodea represents an ancient lineage of insects that emerged over 300 million years ago, have a diverse gut microbiota, and have a typically gregarious lifestyle. All cockroaches shared a broadly similar gut microbiota, with 66 microbial families present across all species and 13 present in every individual examined. Although our network analysis of the cockroach gut microbiome showed a large amount of connectivity, we demonstrated that gut microbiota cluster strongly by host species. We conducted follow-up tests to determine if cockroaches exhibit phylosymbiosis, or the tendency of host-associated microbial communities to parallel the phylogeny of related host species. Across the full data set, gut microbial community similarity was not found to correlate with host phylogenetic distance. However, a weak but significant phylosymbiotic signature was observed using the matching cluster metric, which allows for localized changes within a phylogenetic tree that are more likely to occur over long evolutionary distances. This finding suggests that host phylogeny plays a large role in structuring the cockroach gut microbiome over shorter evolutionary distances and a weak but significant role in shaping the gut microbiome over extended periods of evolutionary time.IMPORTANCE The gut microbiome plays a key role in host health. Therefore, it is important to understand the evolution of the gut microbiota and how it impacts, and is impacted by, host evolution. In this study, we explore the relationship between host phylogeny and gut microbiome composition in omnivorous, gregarious cockroaches within the Blattodea order, an ancient lineage that spans 300 million years of evolutionary divergence. We demonstrate a strong relationship between host species identity and gut microbiome composition and found a weaker but significant role for host phylogeny in determining microbiome similarity over extended periods of evolutionary time. This study advances our understanding of the role of host phylogeny in shaping the gut microbiome over different evolutionary distances.

Genetic Diversity of Blattella germanica Isolates from Central China based on Mitochondrial Genes

Background:

Blattella germanica is a widespread urban invader insect that can spread numerous types of human pathogens, including bacteria, fungi, and protozoa. Despite the medical significance of B. germanica, the genetic diversity of this species has not been investigated across its wide geographical distribution in China.

Objective:

In this study, the genetic variation of B. germanica was evaluated in central China.

Methods:

Fragments of the mitochondrial cytochrome c oxidase subunit I (COI) gene and the 16S rRNA gene were amplified in 36 B. germanica isolates from 7 regions. The sequence data for COI and 16S rRNA genes were analyzed using bioinformatics methods.

Results:

In total, 13 haplotypes were found among the concatenated sequences. Each sampled population, and the total population, had high haplotype diversity (Hd) that was accompanied by low nucleotide diversity (Pi). Molecular genetic variation analysis indicated that 84.33% of the genetic variation derived from intra-region sequences. Phylogenetic analysis indicated that the B. germanica isolates from central China should be classified as a single population. Demographic analysis rejected the hypothesis of sudden population expansion of the B. germanica population.

Conclusion:

The 36 isolates of B. germanica sampled in this study had high genetic variation and belonged to the same species. They should be classified as a single population. The mismatch distribution analysis and BSP analysis did not support a demographic population expansion of the B. germanica population, which provided useful knowledge for monitoring changes in parasite populations for future control strategies.

Sex loss in insects: causes of asexuality and consequences for genomes

Boasting a staggering diversity of reproductive strategies, insects provide attractive models for the comparative study of the causes and consequences of transitions to asexuality. We provide an overview of some contemporary studies of reproductive systems in insects and compile an initial database of asexual insect genome resources. Insect systems have already yielded some important insights into various mechanisms by which sex is lost, including genetic, endosymbiont-mediated, and hybridization. Studies of mutation and substitution after loss of sex provide the strongest empirical support for hypothesized effects of asexuality, whereas there is mixed evidence for ecological hypotheses such as increased parasite load and altered niche breadth in asexuals. Most hypotheses have been explored in a select few taxa (e.g. stick insects, aphids), such that much of the great taxonomic breadth of insects remain understudied. Given the variation in the proximate causes of asexuality in insects, we argue for expanding the taxonomic breadth of study systems. Despite some challenges for investigating sex in insects, the increasing cost-effectiveness of genomic sequencing makes data generation for closely-related asexual and sexual lineages increasingly feasible.

The complete mitochondrial genome of red-clawed crab Chiromantes haematochir (Sesarmidae: Grapsidae)

The complete mitogenome of Chiromantes haematocheir is 15 899 bp in length contains the typical set of 37 genes, including 13 protein-coding genes (PCGs) (ATP6, ATP8, COI-III, ND1-6, ND4L, and Cyt b), two rRNAs (12S rRNA and 16S rRNA), 22 tRNAs, and a putative CR (D-loop). All PCGs start with an ATG codon except ND1 start with ATT and ND3 start with ATA. TAA is the most frequent stop codon, although COI end with TA-, and COII stop with the single nucleotide T-. The new mtDNA sequence contains 12S rRNA and 16S rRNA of rRNAs, separated with tRNA^val. All tRNAs possess the typical clover leaf secondary structure except for tRNA^Ser(AGN) and tRNA^Leu(CUN), which lacks a dihydroxyuridine (DHU) arm. The CR is 826 bp in length, located between 12S rRNA and tRNA^Glx. The phylogenetic trees support Chiromantes haematochir has close relative with Sesarmops sinensis and Sesarmane glectum.

The Complete Mitochondrial Genome of Coptotermes 'suzhouensis' (syn. Coptotermes formosanus) (Isoptera: Rhinotermitidae) and Molecular Phylogeny Analysis

Coptotermes suzhouensis (Isoptera: Rhinotermitidae) is a significant subterranean termite pest of wooden structures and is widely distributed in southeastern China. The complete mitochondrial DNA sequence of C. suzhouensis was analyzed in this study. The mitogenome was a circular molecule of 15,764 bp in length, which contained 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A+T-rich region with a gene arrangement typical of Isoptera mitogenomes. All PCGs were initiated by ATN codons and terminated by complete termination codons (TAA), except COX2, ND5, and Cytb, which ended with an incomplete termination codon T. All tRNAs displayed a typical clover-leaf structure, except for tRNASer(AGN), which did not contain the stem-loop structure in the DHU arm. The A+T content (69.23%) of the A+T-rich region (949 bp) was higher than that of the entire mitogenome (65.60%), and two different sets of repeat units (A+B) were distributed in this region. Comparison of complete mitogenome sequences with those of Coptotermes formosanus indicated that the two taxa have very high genetic similarity. Forty-one representative termite species were used to construct phylogenetic trees by maximum likelihood, maximum parsimony, and Bayesian inference methods. The phylogenetic analyses also strongly supported (BPP, MLBP, and MPBP = 100%) that all C. suzhouensis and C. formosanus samples gathered into one clade with genetic distances between 0.000 and 0.002. This study provides molecular evidence for a more robust phylogenetic position of C. suzhouensis and inferrs that C. suzhouensis was the synonymy of C. formosanus.

The Complete Mitogenome of the Wood-Feeding Cockroach Cryptocercus meridianus (Blattodea: Cryptocercidae) and Its Phylogenetic Relationship among Cockroach Families

In this study, the complete mitochondrial genome of Cryptocercus meridianus was sequenced. The circular mitochondrial genome is 15,322 bp in size and contains 13 protein-coding genes, two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA genes, and one D-loop region. We compare the mitogenome of C. meridianus with that of C. relictus and C. kyebangensis. The base composition of the whole genome was 45.20%, 9.74%, 16.06%, and 29.00% for A, G, C, and T, respectively; it shows a high AT content (74.2%), similar to the mitogenomes of C. relictus and C. kyebangensis. The protein-coding genes are initiated with typical mitochondrial start codons except for cox1 with TTG. The gene order of the C. meridianus mitogenome differs from the typical insect pattern for the translocation of tRNA-Ser^AGN, while the mitogenomes of the other two Cryptocercus species, C. relictus and C. kyebangensis, are consistent with the typical insect pattern. There are two very long non-coding intergenic regions lying on both sides of the rearranged gene tRNA-Ser^AGN. The phylogenetic relationships were constructed based on the nucleotide sequence of 13 protein-coding genes and two ribosomal RNA genes. The mitogenome of C. meridianus is the first representative of the order Blattodea that demonstrates rearrangement, and it will contribute to the further study of the phylogeny and evolution of the genus Cryptocercus and related taxa.

The comparative mitogenomics and phylogenetics of the two grouse-grasshoppers (Insecta, Orthoptera, Tetrigoidea)

OBJECTIVE

This study aimed to reveal the mitochondrial genomes (mtgenomes) of Tetrix japonica and Alulatettix yunnanensis, and the phylogenetics of Orthoptera species.

METHODS

The mtgenomes of A. yunnanensis and T. japonica were firstly sequenced and assembled through partial sequences amplification, and then the genome organization and gene arrangement were analyzed. Based on nucleotide/amino acid sequences of 13 protein-coding genes and whole mtgenomes, phylogenetic trees were established on 37 Orthoptera species and 5 outgroups, respectively.

RESULTS

Except for a regulation region (A+T rich region), a total of 37 genes were found in mtgenomes of T. japonica and A. yunnanensis, including 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes, which exhibited similar characters with other Orthoptera species. Phylogenetic tree based on 13 concatenated protein-coding nucleotide sequences were considered to be more suitable for phylogenetic reconstruction of Orthoptera species than amino acid sequences and mtgenomes. The phylogenetic relationships of Caelifera species were Acridoidea and Pamphagoidea > Pyrgomorphoidea > Pneumoroidea > Eumastacoidea > Tetrigoidea > Tridactyloidea. Besides, a sister-group relationship between Tettigonioidea and Rhaphidophoroidea was revealed in Ensifera.

CONCLUSION

Concatenated protein-coding nucleotide sequences of 13 genes were suitable for reconstruction of phylogenetic relationship in orthopteroid species. Tridactyloidea was a sister group of Tetrigoidea in Caelifera, and Rhaphidophoroidea was a sister group of Tettigonioidea in Ensifera.

Complete mitochondrial genomes of two blattid cockroaches, Periplaneta australasiae and Neostylopyga rhombifolia, and phylogenetic relationships within the Blattaria

Complete mitochondrial genomes (mitogenomes) of two cockroach species, Periplaneta australasiae and Neostylopyga rhombifolia, 15,605 bp and 15,711 bp in length, respectively, were determined. As reported for other cockroach mitogenomes, the two mitogenomes possessed typical ancestral insect mitogenome gene composition and arrangement. Only several small intergenic spacers were found: one, which was common in all sequenced cockroach mitogenomes except for the genus Cryptocercus, was between tRNA-Ser (UCN) and ND1 and contained a 7bp highly conserved motif (WACTTAA). Three different types of short tandem repeats in the N. rhombifolia control region (CR) were observed. The homologous alignments of these tandem repeats with other six cockroach mitogenome CRs revealed a low similarity. Three conserved sequence blocks (CSB) were detected in both cockroach mitochondrial CRs. CSB1 was specific for blattinine mitogenomes and was highly conserved with 95% similarity, speculating that this block was a possible molecular synapomorphy for this subfamily. CSB3 located nearby downstream of CSB1 and has more variations within blattinine mitogenomes compared with CSB1. The CSB3 was capable of forming stable stem-loop structure with a small T-stretch in the loop portion. We assessed the influence of four datasets and two inference methods on topology within Orthopteroidea. All genes excluding the third codon positions of PCGs could generate more stable topology, and higher posterior probabilities than bootstrap values were presented at some branch nodes. The phylogenetic analysis with different datasets and analytical methods supported the monophyly of Dictyoptera and supported strongly the proposal that Isoptera should be classified as a family (Termitidae) of the Blattaria. Specifically, Shelfordella lateralis was inserted in the clade Periplaneta. Considering the K2P genetic distance, morphological characters, and the phylogenetic trees, we suggested that S. lateralis should be placed in the genus Periplaneta.

Species Delimitation and Phylogenetic Relationships in Ectobiid Cockroaches (Dictyoptera, Blattodea) from China

We collected Ectobiidae cockroach specimens from 44 locations in the south of the Yangtze valley. We obtained 297 COI sequences specimens and carried out phylogenetic and divergence dating analyses, as well as species delimitation analysis using a General Mixed Yule Coalescent (GMYC) framework. The intraspecific and interspecific sequence divergence in Ectobiidae cockroaches ranged from 0.0 to 7.0% and 4.6 to 30.8%, respectively. GMYC analysis resulted in 53 (confidence interval: 37-65) entities (likelihood ratio = 103.63) including 14 downloaded species. The COI GMYC groups partly corresponded to the ectobiid species and 52 ectobiid species were delimited successfully based on the combination of GMYC result with morphological information. We used the molecular data and 6 cockroach fossil calibrations to obtain a preliminary estimate of the timescale of ectobiid evolution. The major subfamilies in the group were found to have diverged between ~125-110 Ma, and morphospecies pairs were found to have diverged ~10 or more Ma.

De Novo Transcriptome Analysis and Detection of Antimicrobial Peptides of the American Cockroach Periplaneta americana (Linnaeus)

Cockroaches are surrogate hosts for microbes that cause many human diseases. In spite of their generally destructive nature, cockroaches have recently been found to harbor potentially beneficial and medically useful substances such as drugs and allergens. However, genomic information for the American cockroach (Periplaneta americana) is currently unavailable; therefore, transcriptome and gene expression profiling is needed as an important resource to better understand the fundamental biological mechanisms of this species, which would be particularly useful for the selection of novel antimicrobial peptides. Thus, we performed de novo transcriptome analysis of P. americana that were or were not immunized with Escherichia coli. Using an Illumina HiSeq sequencer, we generated a total of 9.5 Gb of sequences, which were assembled into 85,984 contigs and functionally annotated using Basic Local Alignment Search Tool (BLAST), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) database terms. Finally, using an in silico antimicrobial peptide prediction method, 86 antimicrobial peptide candidates were predicted from the transcriptome, and 21 of these peptides were experimentally validated for their antimicrobial activity against yeast and gram positive and -negative bacteria by a radial diffusion assay. Notably, 11 peptides showed strong antimicrobial activities against these organisms and displayed little or no cytotoxic effects in the hemolysis and cell viability assay. This work provides prerequisite baseline data for the identification and development of novel antimicrobial peptides, which is expected to provide a better understanding of the phenomenon of innate immunity in similar species.

Complete mitochondrial genome of Polyphaga plancyi (Blattaria: Polyphagidae)

The complete mitochondrial genome of Polyphaga plancyi is a 15,547 bp circular molecule, which contains 37 typical mitochondrial genes (13 protein-coding genes, 2 rRNAs and 22 tRNAs) and a 854 bp D-loop. Its gene arrangement pattern is identical with typical other cockroaches. All protein-coding genes start with an ATG codon except COI, ND3, ND5, ND4, ND4L ND6 and ND1. TAA is the most frequent stop codon, and TAG, GCA, TA- and T- are also occurred very common. The mtDNA sequence contains 12S rRNA and 16S rRNA of rRNA. Except for tRNA^Ser(AGY) and tRNA^Leu(CUN) without the dihydrouridine (DHU) arm, all tRNAs could be folded into canonical cloverleaf secondary structures. The phylogenetic trees from the ML and BI analyses based on the complete mtDNA of nine cockroaches’ species, share similar topologies and high node support values. Polyphaga plancy has close relative with Eupolyphaga sinensis.

Stochastic effects in adaptive reconstruction of body damage: implied the creativity of natural selection

After an injury occurs, mechanical/biochemical loads on muscles influence the composition and structure of recovering muscles; this effect likely occurs in other tissues, cells and biological molecules as well owing to the similarity, interassociation and interaction among biochemical reactions and molecules. The 'damage and reconstruction' model provides an explanation for how an ideal cytoarchitecture is created by reducing components not suitable for bearing loads; in this model, adaptive changes are induced by promoting the stochasticity of biochemical reactions. Biochemical and mechanical loads can direct the stochasticity of biochemical reactions, which can in turn induce cellular changes. Thus, mechanical and biochemical loads, under natural selection pressure, modify the direction of cell- and tissue-level changes and guide the formation of new structures and traits, thereby influencing microevolution. In summary, the 'damage and reconstruction' model accounts for the role of natural selection in the formation of new organisms, helps explain punctuated equilibrium, and illustrates how macroevolution arises from microevolution.

The evolutionary history of termites as inferred from 66 mitochondrial genomes

Termites have colonized many habitats and are among the most abundant animals in tropical ecosystems, which they modify considerably through their actions. The timing of their rise in abundance and of the dispersal events that gave rise to modern termite lineages is not well understood. To shed light on termite origins and diversification, we sequenced the mitochondrial genome of 48 termite species and combined them with 18 previously sequenced termite mitochondrial genomes for phylogenetic and molecular clock analyses using multiple fossil calibrations. The 66 genomes represent most major clades of termites. Unlike previous phylogenetic studies based on fewer molecular data, our phylogenetic tree is fully resolved for the lower termites. The phylogenetic positions of Macrotermitinae and Apicotermitinae are also resolved as the basal groups in the higher termites, but in the crown termitid groups, including Termitinae + Syntermitinae + Nasutitermitinae + Cubitermitinae, the position of some nodes remains uncertain. Our molecular clock tree indicates that the lineages leading to termites and Cryptocercus roaches diverged 170 Ma (153-196 Ma 95% confidence interval [CI]), that modern Termitidae arose 54 Ma (46-66 Ma 95% CI), and that the crown termitid group arose 40 Ma (35-49 Ma 95% CI). This indicates that the distribution of basal termite clades was influenced by the final stages of the breakup of Pangaea. Our inference of ancestral geographic ranges shows that the Termitidae, which includes more than 75% of extant termite species, most likely originated in Africa or Asia, and acquired their pantropical distribution after a series of dispersal and subsequent diversification events.

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.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

Terrestrial isopods (Crustacea: Isopoda: Oniscidea) in termite nests (Blattodea: Termitidae) in a cocoa plantation in Brazil

Various animals live in termite nests as secondary occupants. Among them are terrestrial isopods. We found five species of terrestrial isopods (Oniscidea) in 17 of 34 termite nests of Nasutitermes in a cocoa plantation near Ilhéus, Bahia, Brazil (in 2 aging and 15 abandoned nests). The 14 active, 1 aging and 2 abandoned nests bore no isopods. The oniscideans were Atlantoscia rubromarginata, an undescribed species of Atlantoscia (Philosciidae), Pudeoniscus birabeni, P. obscurus (Pudeoniscidae) and Neotroponiscus carolii (Bathytropidae). Of the 146 individuals collected, 122 were A. rubromarginata. Due to the absence of isopods in active nests, these terrestrial isopods are considered termitariophilous.

Mitochondrial genome of the Emberiza rustica (Emberizidae: Emberiza)

Emberiza rustica, least concern species (IUCN), is a passerine bird in the bunting family with wide geographical range. The complete mitochondrial genome of E. rustica (16,798 bp in length) had been analyzed for building the database. Similar to the typical mtDNA of vertebrates, it contained 37 genes (13 protein-coding genes, 2 rRNA genes and 22 tRNA genes) and a non-coding region (D-loop). All the protein-coding genes in E. rustica were distributed on the H-strand, except for the ND6 subunit gene and 10 tRNA genes which were encoded on the L-strand.