On seven undescribed leaf insect species revealed within the recent "Tree of Leaves" (Phasmatodea, Phylliidae)

With the recent advance in molecular phylogenetics focused on the leaf insects (Phasmatodea, Phylliidae), gaps in knowledge are beginning to be filled. Yet, shortcomings are also being highlighted, for instance, the unveiling of numerous undescribed phylliid species. Here, some of these taxa are described, including Phylliumiyadaonsp. nov. from Mindoro Island, Philippines; Phylliumsamarensesp. nov. from Samar Island, Philippines; Phylliumortizisp. nov. from Mindanao Island, Philippines; Pulchriphylliumheraclessp. nov. from Vietnam; Pulchriphylliumdelisleisp. nov. from South Kalimantan, Indonesia; and Pulchriphylliumbhaskaraisp. nov. from Java, Indonesia. Several additional specimens of these species together with a seventh species described herein, Pulchriphylliumanangusp. nov. from southwestern India, were incorporated into a newly constructed phylogenetic tree. Additionally, two taxa that were originally described as species, but in recent decades have been treated as subspecies, are elevated back to species status to reflect their unique morphology and geographic isolation, creating the following new combinations: Pulchriphylliumscythe (Gray, 1843) stat. rev., comb. nov. from Bangladesh and northeastern India, and Pulchriphylliumcrurifolium (Audinet-Serville, 1838) stat. rev., comb. nov. from the Seychelles islands. Lectotype specimens are also designated for Pulchriphylliumscythe (Gray, 1843) stat. rev., comb. nov. and Pulchriphylliumcrurifolium (Audinet-Serville, 1838) stat. rev., comb. nov. from original type material.

Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, non-model insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression.

A comprehensive analysis of the phylogenetic signal in ramp sequences in 211 vertebrates

Ramp sequences increase translational speed and accuracy when rare, slowly-translated codons are found at the beginnings of genes. Here, the results of the first analysis of ramp sequences in a phylogenetic construct are presented. Ramp sequences were compared from 247 vertebrates (114 Mammalian and 133 non-mammalian), where the presence and absence of ramp sequences was analyzed as a binary character in a parsimony and maximum likelihood framework. Additionally, ramp sequences were mapped to the Open Tree of Life synthetic tree to determine the number of parallelisms and reversals that occurred, and those results were compared to random permutations. Parsimony and maximum likelihood analyses of the presence and absence of ramp sequences recovered phylogenies that are highly congruent with established phylogenies. Additionally, 81% of vertebrate mammalian ramps and 81.2% of other vertebrate ramps had less parallelisms and reversals than the mean from 1000 randomly permuted trees. A chi-square analysis of completely orthologous ramp sequences resulted in a p-value < 0.001 as compared to random chance. Ramp sequences recover comparable phylogenies as other phylogenomic methods. Although not all ramp sequences appear to have a phylogenetic signal, more ramp sequences track speciation than expected by random chance. Therefore, ramp sequences may be used in conjunction with other phylogenomic approaches if many orthologs are taken into account. However, phylogenomic methods utilizing few orthologs should be cautious in incorporating ramp sequences because individual ramp sequences may provide conflicting signals.

Codon Pairs are Phylogenetically Conserved: A comprehensive analysis of codon pairing conservation across the Tree of Life

Identical codon pairing and co-tRNA codon pairing increase translational efficiency within genes when two codons that encode the same amino acid are translated by the same tRNA before it diffuses from the ribosome. We examine the phylogenetic signal in both identical and co-tRNA codon pairing across 23 428 species using alignment-free and parsimony methods. We determined that conserved codon pairing typically has a smaller window size than the length of a ribosome, and codon pairing tracks phylogenies across various taxonomic groups. We report a comprehensive analysis of codon pairing, including the extent to which each codon pairs. Our parsimony method generally recovers phylogenies that are more congruent with the established phylogenies than our alignment-free method. However, four of the ten taxonomic groups did not have sufficient orthologous codon pairings and were therefore analyzed using only the alignment-free methods. Since the recovered phylogenies using only codon pairing largely match phylogenies from the Open Tree of Life and the NCBI taxonomy, and are comparable to trees recovered by other algorithms, we propose that codon pairing biases are phylogenetically conserved and should be considered in conjunction with other phylogenomic techniques.

Codon use and aversion is largely phylogenetically conserved across the tree of life

Using parsimony, we analyzed codon usages across 12,337 species and 25,727 orthologous genes to rank specific genes and codons according to their phylogenetic signal. We examined each codon within each ortholog to determine the codon usage for each species. In total, 890,814 codons were parsimony informative. Next, we compared species that used a codon with species that did not use the codon. We assessed each codon's congruence with species relationships provided in the Open Tree of Life (OTL) and determined the statistical probability of observing these results by random chance. We determined that 25,771 codons had no parallelisms or reversals when mapped to the OTL. Codon usages from orthologous genes spanning many species were 1109× more likely to be congruent with species relationships in the OTL than would be expected by random chance. Using the OTL as a reference, we show that codon usage is phylogenetically conserved within orthologous genes in archaea, bacteria, plants, mammals, and other vertebrates. We also show how to use our provided framework to test different tree hypotheses by confirming the placement of turtles as sister taxa to archosaurs.

CAM: an alignment-free method to recover phylogenies using codon aversion motifs

BACKGROUND

Common phylogenomic approaches for recovering phylogenies are often time-consuming and require annotations for orthologous gene relationships that are not always available. In contrast, alignment-free phylogenomic approaches typically use structure and oligomer frequencies to calculate pairwise distances between species. We have developed an approach to quickly calculate distances between species based on codon aversion.

METHODS

Utilizing a novel alignment-free character state, we present CAM, an alignment-free approach to recover phylogenies by comparing differences in codon aversion motifs (i.e., the set of unused codons within each gene) across all genes within a species. Synonymous codon usage is non-random and differs between organisms, between genes, and even within a single gene, and many genes do not use all possible codons. We report a comprehensive analysis of codon aversion within 229,742,339 genes from 23,428 species across all kingdoms of life, and we provide an alignment-free framework for its use in a phylogenetic construct. For each species, we first construct a set of codon aversion motifs spanning all genes within that species. We define the pairwise distance between two species, A and B, as one minus the number of shared codon aversion motifs divided by the total codon aversion motifs of the species, A or B, containing the fewest motifs. This approach allows us to calculate pairwise distances even when substantial differences in the number of genes or a high rate of divergence between species exists. Finally, we use neighbor-joining to recover phylogenies.

RESULTS

Using the Open Tree of Life and NCBI Taxonomy Database as expected phylogenies, our approach compares well, recovering phylogenies that largely match expected trees and are comparable to trees recovered using maximum likelihood and other alignment-free approaches. Our technique is much faster than maximum likelihood and similar in accuracy to other alignment-free approaches. Therefore, we propose that codon aversion be considered a phylogenetically conserved character that may be used in future phylogenomic studies.

AVAILABILITY

CAM, documentation, and test files are freely available on GitHub at https://github.com/ridgelab/cam.

Missing something? Codon aversion as a new character system in phylogenetics

Although many studies have documented codon usage bias in different species, the importance of codon usage in a phylogenetic framework remains largely unknown. We demonstrate that a phylogenetic signal is present in the codon usage and non-usage biases of 17 717 orthologues evaluated across 72 tetrapod species using a simple parsimony analysis of a binary matrix of codon characters. Phylogenies estimated using stop codons were more congruent with previous hypotheses than phylogenies based on any other single codon or a combination of codons. Although each codon is present in every species, specific genes have different codon preferences and may or may not use every possible codon. This observation allowed us to map the pattern of codon usage and non-usage across the topology. These results suggest that codon usage is phylogenetically conserved across shallow and deep levels within tetrapods.

New record of a phoretic flea associated with earwigs (Dermaptera, Arixeniidae) and a redescription of the bat flea Lagaropsylla signata (Siphonaptera, Ischnopsyllidae)

Lagaropsyllasignata (Wahlgren, 1903), previously known only from the Island of Java, Indonesia is redescribed and reported for the first time in Deer Cave, Gunung Mulu National Park, Sarawak, Malaysia (west coast of Borneo). Many were found clinging to the earwig Arixeniaesau Jordan, 1909. A similar account of a phoretic flea (Lagaropsyllaturba Smit, 1958) on the same species of cave-dwelling earwig has been reported in peninsular Malaysia in a well-documented association with the hairless naked bulldog bat, Cheiromelestorquatus Horsfield, 1824. The association of Lagaropsyllasignata with Arixeniaesau is parallel to the evolution and co-existence with bats in Deer Cave just as in the case of Lagaropsyllaturba, Arixeniaesau, and Cheiromelestorquatus. The evidence suggests that Lagaropsyllaturba and Lagaropsyllasignata are obligate phoretic parasites whose survival depends on Arixeniaesau to access a bat host. Arixeniaesau is reported for the first time in Deer Cave and the occurrence of Lagaropsyllasignata on the island of Borneo represented a new record, previously being found only on the island of Java. Images of Lagaropsyllasignata attached to Arixeniaesau are provided. Xeniariajacobsoni (Burr, 1912), often associated with Arixeniaesau in other geographical areas, was not present in the material examined from Deer Cave. The natural history of the earwig genera Arixenia Jordan, 1909 and Xeniaria Maa, 1974 are discussed and summarized relative to their associations with phoretic fleas and their bat hosts.

Reassessing the phylogenetic position of the epizoic earwigs (Insecta: Dermaptera)

Dermaptera is a relatively small order of free-living insects that typically feed on detritus and other plant material. However, two earwig lineages - Arixeniidae and Hemimeridae - are epizoic on Cheiromeles bats and Beamys and Cricetomys rats respectively. Both of these epizoic families are comprised of viviparous species. The monophyly of these epizoic lineages and their placement within dermapteran phylogeny has remained unclear. A phylogenetic analyses was performed on a diverse sample of 47 earwig taxa for five loci (18S rDNA, 28S rDNA, COI, Histone 3, and Tubulin Alpha I). Our results support two independent origins of the epizoic lifestyle within Dermaptera, with Hemimeridae and Arixeniidae each derived from a different lineage of Spongiphoridae. Our analyses places Marava, a genus of spongiphorids that includes free-living but viviparous earwigs, as sister group to Arixeniidae, suggesting that viviparity evolved prior to the shift to the epizoic lifestyle. Additionally, our results support the monophyly of Forficulidae and Chelisochidae and the paraphyly of Labiduridae, Pygidicranidae, Spongiphoridae, and Anisolabididae.

On the probability of dinosaur fleas

Recently, a set of publications described flea fossils from Jurassic and Early Cretaceous geological strata in northeastern China, which were suggested to have parasitized feathered dinosaurs, pterosaurs, and early birds or mammals. In support of these fossils being fleas, a recent publication in BMC Evolutionary Biology described the extended abdomen of a female fossil specimen as due to blood feeding.We here comment on these findings, and conclude that the current interpretation of the evolutionary trajectory and ecology of these putative dinosaur fleas is based on appeal to probability, rather than evidence. Hence, their taxonomic positioning as fleas, or stem fleas, as well as their ecological classification as ectoparasites and blood feeders is not supported by currently available data.

Reinventing the leaf: multiple origins of leaf-like wings in katydids (Orthoptera : Tettigoniidae)

Insects have developed incredible means to avoid detection by predators. At least five insect orders have species that resemble leaves. Katydids (Orthoptera : Tettigoniidae) are the most diverse and wide-ranging of the leaf-like insects. At least 14 of the 20 extant katydid subfamilies contain species with leaf-like wings. Although it is undisputed that many katydids resemble leaves, methods for delineating the leaf-like from non-leaf-like forms have varied by author and in many cases are not explicitly stated. We provide a simple ratio method that can be used to differentiate the leaf-like and non-leaf-like forms. Geometric morphometrics were used to validate the ratio method. Leaf-like wings have been independently derived in at least 15 katydid lineages. Relative rates of speciation were found to be greater in the non-leaf-like forms, suggesting that leaf-like wings within Tettigoniidae are not a driver of diversification. Likewise, throughout Tettigoniidae, selection seems to be favouring the transition away from leaf-like wings. However, within the large Phaneropterinae subclade, relative speciation and transition rates between the leaf-like and non-leaf-like forms do not differ significantly.

Phylogenetic analyses of termite post-embryonic sequences illuminate caste and developmental pathway evolution

Termites are highly eusocial insects with a caste polyphenism (i.e., discontinuous morphological differences between castes) and elaborated behaviors. While the developmental pathways leading to caste occurrence are well-known in many species, the evolutionary origin of these pathways is still obscure. Recent molecular phylogenetic studies suggest multiple independent origins of sterile castes in termites, reviving a 30 years old debate. We demonstrate here that diploid sterile castes ("true" workers) evolved several times independently in this group and that this caste was lost at least once in a lineage with developmentally more flexible workers called pseudergates or "false" workers. We also infer that flexibility in post-embryonic development was acquired multiple times independently during termite evolution. We suggest that focusing on detailed developmental pathways in phylogenetic analyses is essential for elucidating the origin of caste polyphenism in termites.

Rampant nuclear insertion of mtDNA across diverse lineages within Orthoptera (Insecta)

Nuclear mitochondrial pseudogenes (numts) are non-functional fragments of mtDNA inserted into the nuclear genome. Numts are prevalent across eukaryotes and a positive correlation is known to exist between the number of numts and the genome size. Most numt surveys have relied on model organisms with fully sequenced nuclear genomes, but such analyses have limited utilities for making a generalization about the patterns of numt accumulation for any given clade. Among insects, the order Orthoptera is known to have the largest nuclear genome and it is also reported to include several species with a large number of numts. In this study, we use Orthoptera as a case study to document the diversity and abundance of numts by generating numts of three mitochondrial loci across 28 orthopteran families, representing the phylogenetic diversity of the order. We discover that numts are rampant in all lineages, but there is no discernable and consistent pattern of numt accumulation among different lineages. Likewise, we do not find any evidence that a certain mitochondrial gene is more prone to nuclear insertion than others. We also find that numt insertion must have occurred continuously and frequently throughout the diversification of Orthoptera. Although most numts are the result of recent nuclear insertion, we find evidence of very ancient numt insertion shared by highly divergent families dating back to the Jurassic period. Finally, we discuss several factors contributing to the extreme prevalence of numts in Orthoptera and highlight the importance of exploring the utility of numts in evolutionary studies.

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.

Mitochondrial genome deletions and minicircles are common in lice (Insecta: Phthiraptera)

Background

The gene composition, gene order and structure of the mitochondrial genome are remarkably stable across bilaterian animals. Lice (Insecta: Phthiraptera) are a major exception to this genomic stability in that the canonical single chromosome with 37 genes found in almost all other bilaterians has been lost in multiple lineages in favour of multiple, minicircular chromosomes with less than 37 genes on each chromosome.

Results

Minicircular mt genomes are found in six of the ten louse species examined to date and three types of minicircles were identified: heteroplasmic minicircles which coexist with full sized mt genomes (type 1); multigene chromosomes with short, simple control regions, we infer that the genome consists of several such chromosomes (type 2); and multiple, single to three gene chromosomes with large, complex control regions (type 3). Mapping minicircle types onto a phylogenetic tree of lice fails to show a pattern of their occurrence consistent with an evolutionary series of minicircle types. Analysis of the nuclear-encoded, mitochondrially-targetted genes inferred from the body louse, Pediculus, suggests that the loss of mitochondrial single-stranded binding protein (mtSSB) may be responsible for the presence of minicircles in at least species with the most derived type 3 minicircles (Pediculus, Damalinia).

Conclusions

Minicircular mt genomes are common in lice and appear to have arisen multiple times within the group. Life history adaptive explanations which attribute minicircular mt genomes in lice to the adoption of blood-feeding in the Anoplura are not supported by this expanded data set as minicircles are found in multiple non-blood feeding louse groups but are not found in the blood-feeding genus Heterodoxus. In contrast, a mechanist explanation based on the loss of mtSSB suggests that minicircles may be selectively favoured due to the incapacity of the mt replisome to synthesize long replicative products without mtSSB and thus the loss of this gene lead to the formation of minicircles in lice.

Intron retention in the Drosophila melanogaster Rieske Iron Sulphur Protein gene generated a new protein

Genomes can encode a variety of proteins with unrelated architectures and activities. It is known that protein-coding genes of de novo origin have significantly contributed to this diversity. However, the molecular mechanisms and evolutionary processes behind these originations are still poorly understood. Here we show that the last 102 codons of a novel gene, Noble, assembled directly from non-coding DNA following an intronic deletion that induced alternative intron retention at the Drosophila melanogaster Rieske Iron Sulphur Protein (RFeSP) locus. A systematic analysis of the evolutionary processes behind the origin of Noble showed that its emergence was strongly biased by natural selection on and around the RFeSP locus. Noble mRNA is shown to encode a bona fide protein that lacks an iron sulphur domain and localizes to mitochondria. Together, these results demonstrate the generation of a novel protein at a naturally selected site.

Identification of a developmentally and hormonally regulated Delta-Class glutathione S-transferase in rice moth Corcyra cephalonica

Glutathione S-transferases (GSTs) are a large family of multifunctional enzymes, known for their role in cellular detoxification. Here we report a cytosolic GST with optimal activity at alkaline pH (8.3) from the visceral fat body of late-last instar (LLI) larvae of a lepidopteran insect rice moth Corcyra cephalonica. All previously known GSTs are active between pH 6.0 to 6.5. Purification and characterization revealed the Corcyra cephalonica GST (CcGST) as a 23-kDa protein. HPLC and 2D analysis showed a single isoform of the protein in the LLI visceral fat body. Degenerate primer based method identified a 701-nucleotide cDNA and the longest open reading frame contained 216 amino acids. Multiple sequence and structural alignment showed close similarity with delta-class GSTs. CcGST is present mainly in the fat body with highest activity at the late-last instar larval stage. Juvenile hormone (JH) negatively inhibits the CcGST activity both ex vivo and in vivo. We speculate that high expression and activity of CcGST in the fat body of the late-last instar larvae, when endogenous JH titer is low may have role in the insect post-embryonic development unrelated to their previously known function.

Pupal deposition and ecology of bat flies (Diptera: Streblidae): Trichobius sp. (caecus group) in a Mexican cave habitat

We studied the deposition of pupae of the winged bat fly Trichobius sp. (caecus group; Diptera), an ectoparasite of Natalus stramineus (Chiroptera, Natalidae), in a natural cave in Tamaulipas, Mexico. For the first time, we show a strong spatial segregation of populations of a streblid bat fly at different stages of development. Using molecular techniques we were able to match developmental stages to adults. Only 5 pupae were present in the main bat roosts. The overwhelming majority occurred exclusively in the bat flyway passages at a considerable distance from roosting bats. Pupal density corresponded positively with the average flight height of bats in the cave passage. Taken together, observations suggest that these ectoparasites must actively seek out their hosts by moving onto passing or roosting bats. The scarceness of pupae in the main roost may be dictated by environmental constraints for their development. The estimated population of viable pupae far exceeds the population of imagoes on the bats, and predation on adults by spiders is common.

Phylogeny and evolution of male genitalia within the praying mantis genus Tenodera (Mantodea:Mantidae)

The mantis genus Tenodera is composed of several species distributed across Africa, Asia and Australasia, along with recent human introductions to North America. Species of the genus are morphologically similar and utilise equivalent habitats across their distribution. Relationships among these species and the morphological characters used to diagnose them have never been formally tested, leaving authors to disagree as to the species composition of Tenodera. With DNA sequence data from five molecular loci and morphological characters from male genitalia, we reconstructed the phylogeny of Tenodera using multiple optimality criteria. All included species were found to be monophyletic in analyses of the combined data. Tenodera sinensis and T. bokiana were both supported as distinct species recovered in separate clades, resolving confusion as to their placement and classification. Our analysis identified a previously undescribed species of Tenodera collected in India, recovered as sister to T. aridifolia and T. sinensis, and exhibiting distinct male genital morphology. In light of the phylogeny, we characterise for the first time, and investigate the evolution of, the male genitalia, which allowed us to discover several transitions in structural forms. We also consider the connection of these transitions to sexual cannibalism and how this behaviour may have led to rapid evolution of the male genitalia.

A comparative analysis of mitochondrial genomes in Coleoptera (Arthropoda: Insecta) and genome descriptions of six new beetles

Coleoptera is the most diverse group of insects with over 360,000 described species divided into four suborders: Adephaga, Archostemata, Myxophaga, and Polyphaga. In this study, we present six new complete mitochondrial genome (mtgenome) descriptions, including a representative of each suborder, and analyze the evolution of mtgenomes from a comparative framework using all available coleopteran mtgenomes. We propose a modification of atypical cox1 start codons based on sequence alignment to better reflect the conservation observed across species as well as findings of TTG start codons in other genes. We also analyze tRNA-Ser(AGN) anticodons, usually GCU in arthropods, and report a conserved UCU anticodon as a possible synapomorphy across Polyphaga. We further analyze the secondary structure of tRNA-Ser(AGN) and present a consensus structure and an updated covariance model that allows tRNAscan-SE (via the COVE software package) to locate and fold these atypical tRNAs with much greater consistency. We also report secondary structure predictions for both rRNA genes based on conserved stems. All six species of beetle have the same gene order as the ancestral insect. We report noncoding DNA regions, including a small gap region of about 20 bp between tRNA-Ser(UCN) and nad1 that is present in all six genomes, and present results of a base composition analysis.

Mitochondrial genome organization and phylogeny of two vespid wasps

We sequenced the entire mitochondrial genome of Abispa ephippium (Hymenoptera: Vespoidea: Vespidae: Eumeninae) and most of the mitochondrial genome of Polistes humilis synoecus (Hymenoptera: Vespoidea: Vespidae: Polistinae). The arrangement of genes differed between the two genomes and also differed slightly from that inferred to be ancestral for the Hymenoptera. The genome organization for both vespids is different from that of all other mitochondrial genomes previously reported. A number of tRNA gene rearrangements were identified that represent potential synapomorphies for a subset of the Vespidae. Analysis of all available hymenopteran mitochondrial genome sequences recovered an uncontroversial phylogeny, one consistent with analyses of other types of data.

Head morphology of Caurinus (Boreidae, Mecoptera) and its phylogenetic implications

External and internal head structures of Caurinus dectes were examined and described in detail. The features are compared to conditions found in other groups of Antliophora. Caurinus is obviously crucial for the reconstruction of the mecopteran and antliophoran groundplan. It displays a remarkable series of plesiomorphic character states such as a complete clypeolabral suture, the presence of M. hypopharyngomandibularis (M. 13) and M. frontohypopharyngalis (M. 41), a subdivided clypeus, a short head without rostrum, a dorsal tentorial arm attached to the head capsule, the absence of a cranial dilator of the antenna, and large mandibles with a well developed apical tooth, two distinct subapical teeth, and a basal molar part. The first three plesiomorphic features render potential autapomorphies of Mecoptera in the traditional sense invalid. Autapomorphies of Caurinus are the distinctly flattened labrum, the absence of the labroepipharyngeal muscle, the very large size of M. 13, the strongly enlarged penultimate palpomeres, the partition of M. 41, the very strongly developed precerebral sucking chamber, strongly curved optic lobes, the presence of a large protocerebral extension in the genal region and deep posterior excavations of the protocerebrum. The maxillolabial plate, the absence of cardines as separate structures, the reduction of ocelli, and the origin of maxillary palp muscles on a median ridge or area of the maxillolabial plate are likely autapomorphies of Boreidae. Another potential autapomorphy of the family is the presence of longitudinal furrows on the mandibles. However, they are absent in Boreus. The thick strongly sclerotised, median ridge of the maxillolabial plate, the missing retractibility of the prementum, the absence of extrinsic labial muscles, and the presence of a median ridge on the prepharyngeal roof suggest a clade Boreus+Hesperoboreus. The origin of extrinsic maxillary muscles from the clypeus has probably evolved independently in Boreus and Hesperoboreus, and in Panorpa, respectively. The absence of M. craniolacinialis and the presence of a row of several subapical mandibular teeth are autapomorphies of Boreus. The presence of a specific intrinsic muscle of the salivary duct and a membranous galea enclosing the labrum and mandibular base are derived features shared by Boreidae and Pistillifera (galea absent in Nannochorista, Siphonaptera and Diptera). The loss of M. frontolabralis (M. 8) is a potential apomorphy of Mecoptera incl. Siphonaptera. A sister group relationship between Boreidae and Siphonaptera is not supported by characters of the adult head. Head structures of Siphonaptera are extremely modified in correlation with ectoparasitic habits.

The phylogeny of termites (Dictyoptera: Isoptera) based on mitochondrial and nuclear markers: Implications for the evolution of the worker and pseudergate castes, and foraging behaviors

A phylogenetic hypothesis of termite relationships was inferred from DNA sequence data. Seven gene fragments (12S rDNA, 16S rDNA, 18S rDNA, 28S rDNA, cytochrome oxidase I, cytochrome oxidase II and cytochrome b) were sequenced for 40 termite exemplars, representing all termite families and 14 outgroups. Termites were found to be monophyletic with Mastotermes darwiniensis (Mastotermitidae) as sister group to the remainder of the termites. In this remainder, the family Kalotermitidae was sister group to other families. The families Kalotermitidae, Hodotermitidae and Termitidae were retrieved as monophyletic whereas the Termopsidae and Rhinotermitidae appeared paraphyletic. All of these results were very stable and supported with high bootstrap and Bremer values. The evolution of worker caste and foraging behavior were discussed according to the phylogenetic hypothesis. Our analyses suggested that both true workers and pseudergates ("false workers") were the result of at least two different origins. Our data support a traditional hypothesis of foraging behavior, in which the evolutionary transition from a one-piece type to a separate life type occurred through an intermediate behavioral form.

The complete mitochondrial genome of the tobacco hornworm, Manduca sexta, (Insecta: Lepidoptera: Sphingidae), and an examination of mitochondrial gene variability within butterflies and moths

The entire mitochondrial genome of the tobacco hornworm, Manduca sexta (Lepidoptera: Spinghidae) was sequenced -- a circular molecular 15516 bp in size. The arrangement of the protein coding genes (PCGs) was the same as that found in the ancestral insect, however Manduca possessed the derived tRNA arrangement of CR-M-I-Q which has been found in all Lepidoptera sequenced to date. Additionally, Manduca, like all lepidopteran mt genomes, has numerous large intergenic spacer regions and microsatellite-like repeat regions. Nucleotide composition is highly A+T biased, and the lepidopterans have the second most biased nucleotide composition of the insect orders after Hymenoptera. Secondary structural features of the PCGs identified in other Lepidoptera were present but highly modified by the presence of microsatellite-like repeat regions which may significantly alter their function in the post-transcriptional modification of pre-mRNAs. Secondary structure models of the ribosomal RNA genes of Manduca are presented and are similar to those proposed for other insect orders. Conserved regions were identified within non-translated spacer regions which correspond to sites for the origin and termination of replication and transcription. Comparisons of gene variability across the order suggest that the mitochondrial genes most frequently used in phylogenetic analysis of the Lepidoptera, cox1 and cox2, are amongst the least variable genes in the genome and phylogenetic resolution could be improved by using alternative, higher variability genes such as nad2, nad3, nad4 and nad5.

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.

The complete mitochondrial genome sequence of the Mormon cricket (Anabrus simplex: Tettigoniidae: Orthoptera) and an analysis of control region variability

The Anabrus simplex is a swarming plague orthopteran found in western North America. The genome is 15 766 bp in length and genome organization follows the ancestral insect gene arrangement. atp6 lacked any readily identifiable stop codon. Examination of mRNA secondary structure for this gene suggested a stem/loop-mediated mRNA post-transcriptional processing to liberate a mature atp6 mRNA with a complete stop codon produced by polyadenylation. Comparison of similar protein with protein gene boundaries in other insect species reveal a general mechanism for mRNA excision and provide further supporting evidence for post-transcriptional mRNA processing in mitochondrial genomes. The A + T-rich region, or control region, was sequenced for 55 A. simplex individuals from 12 different populations. Variance studies between these individuals show that the A + T-rich region contains significant phylogenetic signal to be used in population studies.

Phylogeny and biogeography of ice crawlers (Insecta: Grylloblattodea) based on six molecular loci: designating conservation status for Grylloblattodea species

Ice crawlers (Insecta: Grylloblattodea) are rarely encountered insects that consist of five genera representing 26 species from North America and Asia. Asian grylloblattids are the most diverse, but North American ice crawlers (genus Grylloblatta) are known for their adaptation to cold conditions. Phylogenetic relationships among grylloblattid species and genera are not known. Six genes were sampled in 35 individuals for 18S rRNA, 28S rRNA, histone 3, 12S rRNA, 16S rRNA, and cytochrome oxidase II from 21 populations of Grylloblatta, three populations from Japan (genus Galloisiana), and three populations from Russia (genus Grylloblattina). Phylogenetic analysis of these data with two mantophasmid outgroups in POY supported monophyletic genera, with Grylloblatta as sister to Grylloblattina. Grylloblatta was shown to contain two major lineages: a clade in Northern California and Oregon and a clade in Washington and Oregon. One new species and six candidate species are proposed. IUCN Red List Conservation Criteria were implemented to designate conservation status for each lineage.

Identification and phylogenetic analysis of Arsenophonus- and Photorhabdus-type bacteria from adult Hippoboscidae and Streblidae (Hippoboscoidea)

This is the first report of Arsenophonus- and Photorhabdus-type bacteria from Streblidae (bat flies) and Hippoboscidae (louse flies, keds). Strains were detected by means of polymerase chain reaction of 16S rDNA, and phylogenetic analysis determined the relationship of the obtained sequences to previously reported sequences in GenBank. Phylogenetic analysis by means of maximum parsimony revealed that all isolated Arsenophonus spp. 16S rDNA sequences formed a monophyletic sub-clade within other insect Arsenophonus spp., while the Photorhabdus spp. sequences are part of a monophyletic clade including Photorhabdus spp., Xenorhabdus spp. and Proteus spp.