Comparative analysis of mt LSU rRNA secondary structures of Odonates: structural variability and phylogenetic signal

Characterization of the mitochondrial genome of Arge bella Wei & Du sp. nov. (Hymenoptera: Argidae)

We describe Arge bella Wei & Du sp. nov., a large and beautiful species of Argidae from south China, and report its mitochondrial genome based on high-throughput sequencing data. We present the gene order, nucleotide composition of protein-coding genes (PCGs), and the secondary structures of RNA genes. The nearly complete mitochondrial genome of A. bella has a length of 15,576 bp and a typical set of 37 genes (22 tRNAs, 13 PCGs, and 2 rRNAs). Three tRNAs are rearranged in the A. bella mitochondrial genome as compared to the ancestral type in insects: trnM and trnQ are shuffled, while trnW is translocated from the trnW-trnC-trnY cluster to a location downstream of trnI. All PCGs are initiated by ATN codons, and terminated with TAA, TA or T as stop codons. All tRNAs have a typical cloverleaf secondary structure, except for trnS1. H821 of rrnS and H976 of rrnL are redundant. A phylogenetic analysis based on mitochondrial genome sequences of A. bella, 21 other symphytan species, two apocritan representatives, and four outgroup taxa supports the placement of Argidae as sister to the Pergidae within the symphytan superfamily Tenthredinoidea.

The first mitogenomes of the superfamily Pamphilioidea (Hymenoptera: Symphyta): Mitogenome architecture and phylogenetic inference

The Pamphilioidea represents a small superfamily of the phytophagous suborder Symphyta (Hymenoptera). Here, nearly complete mitochondrial genomes (mitogenomes) of three pamphilioid species: Chinolyda flagellicornis (Pamphiliidae), Megalodontes spiraeae and M. cephalotes (Megalodontesidae) were newly sequenced using next generation sequencing and comparatively analysed with the previously reported symphytan mitogenomes. A positive AT skew (0.013) and a negative GC skew (-0.194) were found in pamphilioid mitogenome, and a deviation from strand asymmetry was also observed in the PCGs encoded on both strands. Several gene rearrangement events were observed in four tRNA gene clusters (WCY, IQM, ARNS1EF and TP clusters), which have not been reported from symphytan mitogenomes to date. As the most parsimonious explanation, compared with the inferred insect ancestral mitogenome architecture, the occurrence of gene rearrangements in pamphilioid mitogenomes requires totally five evolutionary steps, including four transpositions and one inversion. The predicted secondary structures of tRNAs, rrnS and rrnL genes are mostly consistent with reported hymenopteran species. Phylogenetic analyses recovered the monophyly of superfamily Pamphilioidea and indicated the relationship Tenthredinoidea + (Pamphilioidea + (Cephoidea + (Orussoidea + Apocrita))) with strong nodal supports.

Insight into higher-level phylogeny of Neuropterida: Evidence from secondary structures of mitochondrial rRNA genes and mitogenomic data

It is well known that the rRNA structure information is important to assist phylogenetic analysis through identifying homologous positions to improve alignment accuracy. In addition, the secondary structure of some conserved motifs is highly stable among distantly related taxa, which can provide potentially informative characters for estimating phylogeny. In this paper, we applied the high-throughput pooled sequencing approach to the determination of neuropteran mitogenomes. Four complete mitogenome sequences were obtained: Micromus angulatus (Hemerobiidae), Chrysoperla nipponensis (Chrysopidae), Rapisma sp. (Ithonidae), and Thaumatosmylus sp. (Osmylidae). This allowed us to sample more complete mitochondrial RNA gene sequences. Secondary structure diagrams for the complete mitochondrial small and large ribosomal subunit RNA genes of eleven neuropterid species were predicted. Comparative analysis of the secondary structures indicated a closer relationship of Megaloptera and Neuroptera. This result was congruent with the resulting phylogeny inferred from sequence alignments of all 37 mitochondrial genes, namely the hypothesis of (Raphidioptera + (Megaloptera + Neuroptera)).

Inferring phylogenetic trees from the knowledge of rare evolutionary events

Rare events have played an increasing role in molecular phylogenetics as potentially homoplasy-poor characters. In this contribution we analyze the phylogenetic information content from a combinatorial point of view by considering the binary relation on the set of taxa defined by the existence of a single event separating two taxa. We show that the graph-representation of this relation must be a tree. Moreover, we characterize completely the relationship between the tree of such relations and the underlying phylogenetic tree. With directed operations such as tandem-duplication-random-loss events in mind we demonstrate how non-symmetric information constrains the position of the root in the partially reconstructed phylogeny.

Next-Generation Sequencing of Two Mitochondrial Genomes from Family Pompilidae (Hymenoptera: Vespoidea) Reveal Novel Patterns of Gene Arrangement

Animal mitochondrial genomes have provided large and diverse datasets for evolutionary studies. Here, the first two representative mitochondrial genomes from the family Pompilidae (Hymenoptera: Vespoidea) were determined using next-generation sequencing. The sequenced region of these two mitochondrial genomes from the species Auplopus sp. and Agenioideus sp. was 16,746 bp long with an A + T content of 83.12% and 16,596 bp long with an A + T content of 78.64%, respectively. In both species, all of the 37 typical mitochondrial genes were determined. The secondary structure of tRNA genes and rRNA genes were predicted and compared with those of other insects. Atypical trnS1 using abnormal anticodons TCT and lacking D-stem pairings was identified. There were 49 helices belonging to six domains in rrnL and 30 helices belonging to three domains in rrns present. Compared with the ancestral organization, four and two tRNA genes were rearranged in mitochondrial genomes of Auplopus and Agenioideus, respectively. In both species, trnM was shuffled upstream of the trnI-trnQ-trnM cluster, and trnA was translocated from the cluster trnA-trnR-trnN-trnS1-trnE-trnF to the region between nad1 and trnL1, which is novel to the Vespoidea. In Auplopus, the tRNA cluster trnW-trnC-trnY was shuffled to trnW-trnY-trnC. Phylogenetic analysis within Vespoidea revealed that Pompilidae and Mutillidae formed a sister lineage, and then sistered Formicidae. The genomes presented in this study have enriched the knowledge base of molecular markers, which is valuable in respect to studies about the gene rearrangement mechanism, genomic evolutionary processes and phylogeny of Hymenoptera.

Sequence variability in the mitochondrial 12S rRNA and tRNA Val genes of Ixodes scapularis (Acari: Ixodidae) individuals shown previously to be genetically invariant

The DNA sequences of the mitochondrial (mt) 12S rRNA and tRNA(Val) genes were characterized for 82 blacklegged ticks (Ixodes scapularis) that were genetically identical for Domains IV and V of the mt 16S rRNA gene. Thirty-one haplotypes, differed in sequence by 1-9 bp, were detected among the 82 ticks. Most nucleotide alterations in DNA sequence did not affect the stability of the secondary structures of the RNAs. The magnitude of the DNA sequence variation in the mt 12S rRNA and tRNA(Val) genes among blacklegged ticks suggests that this region of the mitochondrial genome has potential as a genetic marker for examining the population genetics and phylogeography of I. scapularis.

Two nearly complete mitogenomes of wheat stem borers, Cephus pygmeus (L.) and Cephus sareptanus Dovnar-Zapolskij (Hymenoptera: Cephidae): an unusual elongation of rrnS gene

Two nearly complete mitochondrial genomes (mitogenomes) of wheat stem borers, Cephus pygmeus and Cephus sareptanus (Hymenoptera: Cephidae), were sequenced, characterised and compared with the previously known mitogenome of Cephus cinctus. The gene orders are mostly conserved, except for translocation of trnM and swapped position of trnI and trnQ. An A+T bias was found, but a deviation from strand asymmetry was also detected on the J strand. All protein coding genes (PCGs) are initiated by ATN codons, except for nad1, nad2 and atp8, and all are terminated with TAA, TA- or T- as a stop codon. The predicted secondary structures of rrnS and rrnL genes are mostly consistent with reported hymenopteran species. However, an unusual elongation in rrnS, not know elsewhere in the order, was discovered in Cephus species. Three autonomous sequences detected in domains I and II are mainly responsible for the length expansions.

Sequence and secondary structure of the mitochondrial 16S ribosomal RNA gene of Ixodes scapularis

The complete DNA sequences and secondary structure of the mitochondrial (mt) 16S ribosomal (r) RNA gene were determined for six Ixodes scapularis adults. There were 44 variable nucleotide positions in the 1252 bp sequence alignment. Most (95%) nucleotide alterations did not affect the integrity of the secondary structure of the gene because they either occurred at unpaired positions or represented compensatory changes that maintained the base pairing in helices. A large proportion (75%) of the intraspecific variation in DNA sequence occurred within Domains I, II and VI of the 16S gene. Therefore, several regions within this gene may be highly informative for studies of the population genetics and phylogeography of I. scapularis, a major vector of pathogens of humans and domestic animals in North America.

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

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

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

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

Adipokinetic hormones provide inference for the phylogeny of odonata

Adipokinetic neuropeptides from the corpora cardiaca of 17 species of Odonata encompassing mainly the families Corduliidae and Libellulidae were isolated and structurally elucidated using liquid chromatography coupled with ion trap electrospray ionization mass spectrometry. It became evident that all species of the family Corduliidae studied express the peptide code-named Libau-AKH (pGlu-Val-Asn-Phe-Thr-Pro-Ser-Trp amide), which is also present in all but one libellulid species, Erythemis simplicicollis which expresses Erysi-AKH (pGlu-Leu-Asn-Phe-Thr-Pro-Ser-Trp amide). This divergence from all other Libellulids is due to a nonsynonymous missense single nucleotide polymorphism (SNP) in the nucleotide coding sequence (CDS) of prepro-AKH CDS and supports the polyphyletic nature of Sympetrinae and other subfamilies of libellulids. Despite this exception, these findings then support the hypothesis that Corduliidae and Libellulidae are closely related as stated in most phylogenies. The presence of Anaim-AKH (pGlu-Val-Asn-Phe-Ser-Pro-Ser-Trp amide) in Macromiidae likely distinguishes species in this family from Corduliidae. Current molecular genetic phylogenies and our AKH findings suggest that Syncordulia gracilis, which expresses Anaim-AKH, does not belong in Corduliidae. Evolution of AKHs in anisopteran Odonata are likely due to nucleotide substitution involving nonsynonymous missense SNPs in the CDS of prepro-AKH.

Complete Mitochondrial Genome Sequence of Acrida cinerea (Acrididae: Orthoptera) and Comparative Analysis of Mitochondrial Genomes in Orthoptera

The complete 15,599-bp mitogenome of Acrida cinerea was determined and compared with that of the other 20 orthopterans. It displays characteristic gene content, genome organization, nucleotide composition, and codon usage found in other Caelifera mitogenomes. Comparison of 21 orthopteran sequences revealed that the tRNAs encoded by the H-strand appear more conserved than those by the L-stand. All tRNAs form the typical clover-leaf structure except trnS (agn), and most of the size variation among tRNAs stemmed from the length variation in the arm and loop of TΨC and the loop of DHU. The derived secondary structure models of the rrnS and rrnL from 21 orthoptera species closely resemble those from other insects on CRW except a considerably enlarged loop of helix 1399 of rrnS in Caelifera, which is a potentially autapomorphy of Caelifera. In the A+T-rich region, tandem repeats are not only conserved in the closely related mitogenome but also share some conserved motifs in the same subfamily. A stem-loop structure, 16 bp or longer, is likely to be involved in replication initiation in Caelifera and Grylloidea. A long T-stretch (>17 bp) with conserved stem-loop structure next to rrnS on the H-strand, bounded by a purine at either end, exists in the three species from Tettigoniidae.

The complete mitochondrial genome of Evania appendigaster (Hymenoptera: Evaniidae) has low A+T content and a long intergenic spacer between atp8 and atp6

The apocritan Hymenoptera show extraordinary features in mitochondrial genomes, but no complete sequence has been reported for the basal lineage, Evanioidea. Here, we sequenced the complete mitochondrial genome of Evania appendigaster. This genome is 17,817 bp long; with low A+T content, 77.8%, compared with other hymenopteran species. Four tRNA genes were rearranged, among which remote inversion is the dominant gene rearrangement event. Gene shuffling is caused by tandem duplication-random loss while remote inversion is best explained by recombination. The start codon of nad1 was found as TTG, which might be common across Hymenoptera. trnS2 and trnK use abnormal anticodons TCT and TTT, respectively, and the D-stem pairings in trnS2 are absent. The secondary structure of two rRNA genes are predicted and compared with those in other insects. Five long intergenic spacers were present, including a long intergenic spacer between atp8 and atp6, where these two genes overlap in the previously reported animal genomes. A conserved motif was found between trnS1 and nad1, which is proposed to be associated with mtTERM. The A+T-rich region is 2,325 bp long, among the longest in insects, and contains a tandem repeat region.

The complete mitochondrial genome of Diadegma semiclausum (Hymenoptera: Ichneumonidae) indicates extensive independent evolutionary events

Few complete mitochondrial genomes representing limited families in the order Hymenoptera have been sequenced. Here, we sequenced the complete mitochondrial genome of Diadegma semiclausum (Hymenoptera: Ichneumonidae). This genome is 18 728 bp long, the second largest hexapod mitochondrial genome sequenced in its entirety and that with the highest A+T content at 87.4%. Four tRNAs are rearranged compared with the ancestral arrangement. Gene rearrangement mechanisms are different among all three rearranged regions. Six tRNAs have a large variable loop, which is not found in other metazoan mitochondrial genomes. trnS(AGY) uses the abnormal anticodon TCT but trnK uses the normal CTT. The A+T-rich region is very long (2161 bp). An extremely A+T-rich (99.1%) 1515 bp tandem repeat region with three types of repeat elements is located between cox1 and cox2, and the most likely ancestral element originated from the 3′ end of cox1. Independent tandem duplications followed by mutation–insertion–deletion is the best model to explain the formation of this region. These results indicate that independent evolutionary events occurred extensively, such as gene rearrangement events, gene rearrangement mechanisms, derivation of tRNA variable loops, and tandem repeat region evolutionary processes, all of which likely contribute to the diversified features of hymenopteran mitochondrial genomes.

Phylogeny of the higher Libelluloidea (Anisoptera: Odonata): An exploration of the most speciose superfamily of dragonflies

Although libelluloid dragonflies are diverse, numerous, and commonly observed and studied, their phylogenetic history is uncertain. Over 150 years of taxonomic study of Libelluloidea Rambur, 1842, beginning with Hagen (1840), [Rambur, M.P., 1842. Neuropteres. Histoire naturelle des Insectes, Paris, pp. 534; Hagen, H., 1840. Synonymia Libellularum Europaearum. Dissertation inaugularis quam consensu et auctoritate gratiosi medicorum ordinis in academia albertina ad summos in medicina et chirurgia honores.] and Selys (1850), [de Selys Longchamps, E., 1850. Revue des Odonates ou Libellules d'Europe [avec la collaboration de H.A. Hagen]. Muquardt, Bruxelles; Leipzig, 1-408.], has failed to produce a consensus about family and subfamily relationships. The present study provides a well-substantiated phylogeny of the Libelluloidea generated from gene fragments of two independent genes, the 16S and 28S ribosomal RNA (rRNA), and using models that take into account non-independence of correlated rRNA sites. Ninety-three ingroup taxa and six outgroup taxa were amplified for the 28S fragment; 78 ingroup taxa and five outgroup taxa were amplified for the 16S fragment. Bayesian, likelihood and parsimony analyses of the combined data produce well-resolved phylogenetic hypotheses and several previously suggested monophyletic groups were supported by each analysis. Macromiinae, Corduliidae s. s., and Libellulidae are each monophyletic. The corduliid (s.l.) subfamilies Synthemistinae, Gomphomacromiinae, and Idionychinae form a monophyletic group, separate from the Corduliinae. Libellulidae comprises three previously accepted subfamilies (Urothemistinae, a very restricted Tetrathemistinae, and a modified Libellulinae) and five additional consistently recovered groups. None of the other previously proposed subfamilies are supported. Bayesian analyses run with an additional 71 sequences obtained from GenBank did not alter our conclusions. The evolution of adult and larval morphological characters is discussed here to suggest areas for future focus. This study shows the inherent problems in using poorly defined and sometimes inaccurately scored characters, basing groups on symplesiomorphies, and failure to recognize the widespread effects of character correlation and convergence, especially in aspects of wing venation.

Phylogeny and Bayesian divergence time estimations of small-headed flies (Diptera: Acroceridae) using multiple molecular markers

The first formal analysis of phylogenetic relationships among small-headed flies (Acroceridae) is presented based on DNA sequence data from two ribosomal (16S and 28S) and two protein-encoding genes: carbomoylphosphate synthase (CPS) domain of CAD (i.e., rudimentary locus) and cytochrome oxidase I (COI). DNA sequences from 40 species in 22 genera of Acroceridae (representing all three subfamilies) were compared with outgroup exemplars from Nemestrinidae, Stratiomyidae, Tabanidae, and Xylophagidae. Parsimony and Bayesian simultaneous analyses of the full data set recover a well-resolved and strongly supported hypothesis of phylogenetic relationships for major lineages within the family. Molecular evidence supports the monophyly of traditionally recognised subfamilies Philopotinae and Panopinae, but Acrocerinae are polyphyletic. Panopinae, sometimes considered "primitive" based on morphology and host-use, are always placed in a more derived position in the current study. Furthermore, these data support emerging morphological evidence that the type genus Acrocera Meigen, and its sister genus Sphaerops, are atypical acrocerids, comprising a sister lineage to all other Acroceridae. Based on the phylogeny generated in the simultaneous analysis, historical divergence times were estimated using Bayesian methodology constrained with fossil data. These estimates indicate Acroceridae likely evolved during the late Triassic but did not diversify greatly until the Cretaceous.

Bayesian analysis of new and old malaria parasite DNA sequence data demonstrates the need for more phylogenetic signal to clarify the descent of Plasmodium falciparum

Molecular systematic studies published during the last 15 years to clarify the phylogenetic relationships among the malaria parasites have led to two major hypotheses on the descent of Plasmodium falciparum: One supports an avian origin as a result of a relatively recent host switch, and another one favours the evolutionary development of P. falciparum together with its human host from primate ancestors. In this paper, we present phylogenetic analyses of three different Plasmodium genes, the nuclear 18 small sub-unit (SSU) ribosomal ribonucleic acid (rRNA), the mitochondrial cytochrome b (cyt b) and the plastid caseinolytic protease C (ClpC) gene, using numerous haemosporidian parasite DNA sequences obtained from the GenBank as well as several new sequences for major malaria parasites including the avian one Plasmodium cathemerium, which has never been considered in molecular phylogenetic analyses before. Most modern and sophisticated DNA substitution models based on Bayesian inference analysis were applied to estimate the cyt b and ClpC phylogenetic trees, whereas the 18 SSU rRNA gene was examined with regards to its secondary structure using PHASE software. Our results indicate that the data presently available are generally neither sufficient in number nor in information to solve the problem of the phylogenetic origin of P. falciparum.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

The mitochondrial genome of the entomophagous endoparasite Xenos vesparum (Insecta: Strepsiptera)

In this study, the nearly complete sequence (14,519 bp) of the mitochondrial DNA (mtDNA) of the entomophagous endoparasite Xenos vesparum (Insecta: Strepsiptera) is described. All protein coding genes (PCGs) are in the arrangement known to be ancestral for insects, but three tRNA genes (trnA, trnS(gcu), and trnL(uag)) have transposed to derived positions and there are three tandem copies of trnH, each of which is potentially functional. All of these rearrangements except for that of trnL(uag) is within the short span between nad3 and nad4 and there are numerous blocks of unassignable sequence in this region, perhaps as remnants of larger scale predisposing rearrangements. X. vesparum mtDNA nucleotide composition is strongly biased toward A and T, as is typical for insect mtDNAs. There is also a significant strand skew in the distribution of these nucleotides, with the J-strand being richer in A than T and in C than G, and the N-strand showing an opposite skew for complementary pairs of nucleotides. The hypothetical secondary structure of the LSU rRNA has also been reconstructed, obtaining a structural model similar to that of other insects.

The mitochondrial genomes of Campodea fragilis and Campodea lubbocki (Hexapoda: Diplura): High genetic divergence in a morphologically uniform taxon

Complete mitochondrial genome sequences are presented from two dipluran hexapods (i.e., a group of "primarily wingless insects") of the genus Campodea and compared to those of other arthropods. Their gene order is the same as in most other hexapods and crustaceans. Structural changes have occurred in tRNA-C, tRNA-R, tRNA-S1 and tRNA-S2 as well as in both ribosomal RNAs. These mtDNAs have striking biases in nucleotide and amino acid composition. Although the two Campodea species are morphologically highly similar, their genetic divergence is larger than expected, suggesting a long evolutionary history, perhaps under stable ecological conditions.

Higher phylogeny of zygaenid moths (Insecta: Lepidoptera) inferred from nuclear and mitochondrial sequence data and the evolution of larval cuticular cavities for chemical defence

Zygaenid moths are capable of releasing hydrogen cyanide in their defense by enzymatic break-down of cyanoglucosides, but only larvae of chalcosiine and zygaenine moths store cyanogenic compounds in cuticular cavities and thus are able to discharge defense droplets, which effectively deter potential predators. A previously proposed phylogeny of Zygaenidae hypothesized a sister group relationship of chalcosiine and zygaenine moths because of their similar larval defense system. Not all chalcosiine taxa possess cuticular cavities, however, and a comparable defense mechanism has been reported in larvae of the zygaenoid family Heterogynidae. Considering sequence data of seven molecular loci, the present study estimates the posterior probability of phylogenetic hypotheses explaining the occurrence of larval cuticular cavities. The molecular data confirm the previous exclusion of Himantopteridae from Zygaenidae and suggest their close affinity to Somabrachyidae. The sequence data also corroborate the recently proposed exclusion of the Phaudinae from the Zygaenidae, because this subfamily is recovered in a reasonably well supported species cluster consisting of members of the families Lacturidae, Limacodidae, Himantopteridae, and Somabrachyidae. We consequently agree to raise Phaudinae to family rank. Within Zygaenidae, the subfamilies Callizygaeninae, Chalcosiinae, and Procridinae most likely constitute a monophyletic group, which is sister to the Zygaeninae. Our results imply that cuticular cavities were probably present in the larvae of the most recent common ancestor of Zygaenidae. Heterogynidae cannot be confirmed as sister taxon to this family, but appear at the very first split of the Zygaenoidea, although with poor support. The specific pattern of taxa in the molecular phylogeny showing larval cuticular cavities opens the possibility that these structures could have been already present in the most recent common ancestor of the Zygaenoidea.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

A Hexapod nuclear SSU rRNA secondary-structure model and catalog of taxon-specific structural variation

RNA molecules and in particular the nuclear SSU RNA play an important role in molecular systematics. With the advent of increasingly parameterized substitution models in systematic research, the incorporation of secondary-structure information became a realistic option compensating interdependence of character variation. As a prerequisite, consensus structures of eukaryotic SSU RNA molecules have become available through extensive comparative analyses and crystallographic studies. Despite extensive research in hexapod phylogenetics, consensus SSU RNA secondary structures focusing on hexapods have not yet been explored. In this study, we compiled a representative hexapod SSU data set of 261 sequences and inferred a specific consensus SSU secondary-structure model. Our search for conserved structural motives relied on a combined approach of thermodynamic and covariation analyses. The hexapod consensus-structure model deviates from the canonical eukaryotic model in a number of helices. Additionally, in several helices the hexapod sequences did not support a single consensus structure. We provide consensus structures of these sections of single less-inclusive taxa, thus facilitating the adaptation of the consensus hexapod model to less-inclusive phylogenetic questions. The secondary-structure catalog will foster the application of RNA structure models in phylogenetic analyses using the SSU rRNA molecule, and it will improve the realism of substitution models and the reliability of reconstructions based on rRNA sequences.

Predicted Secondary Structure for 28S and 18S rRNA from Ichneumonoidea (Insecta: Hymenoptera: Apocrita): Impact on Sequence Alignment and Phylogeny Estimation

We utilize the secondary structural properties of the 28S rRNA D2-D10 expansion segments to hypothesize a multiple sequence alignment for major lineages of the hymenopteran superfamily Ichneumonoidea (Braconidae, Ichneumonidae). The alignment consists of 290 sequences (originally analyzed in Belshaw and Quicke, Syst Biol 51:450-477, 2002) and provides the first global alignment template for this diverse group of insects. Predicted structures for these expansion segments as well as for over half of the 18S rRNA are given, with highly variable regions characterized and isolated within conserved structures. We demonstrate several pitfalls of optimization alignment and illustrate how these are potentially addressed with structure-based alignments. Our global alignment is presented online at (http://hymenoptera.tamu.edu/rna) with summary statistics, such as basepair frequency tables, along with novel tools for parsing structure-based alignments into input files for most commonly used phylogenetic software. These resources will be valuable for hymenopteran systematists, as well as researchers utilizing rRNA sequences for phylogeny estimation in any taxon. We explore the phylogenetic utility of our structure-based alignment by examining a subset of the data under a variety of optimality criteria using results from Belshaw and Quicke (2002) as a benchmark.

A Secondary Structural Model of the 28S rRNA Expansion Segments D2 and D3 for Chalcidoid Wasps (Hymenoptera: Chalcidoidea)

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.

The adipokinetic hormones of Odonata: a phylogenetic approach

Adipokinetic neuropeptides from the corpora cardiaca of the major families of all three suborders of the Odonata were identified by one or more of the following methods: (1) Isolation of the peptides from a methanolic extract of the corpora cardiaca by liquid chromatography, peak monitoring by fluorescence of the Trp residue and comparison of the retention time with those of known synthetic peptides of Odonata. (2) Hyperlipaemic bioassays of the HPLC-generated fractions either in Locusta migratoria or, in a few cases, in Anax imperator or Orthetrum julia. (3) Sequencing of the isolated, bioactive HPLAC fraction by Edman degradation. (4) Mass spectrometric measurement of the isolated, bioactive fraction. Sequence assignment revealed that the investigated Odonata species always contain only one adipokinetic peptide. This is always an octapeptide. The suborder Zygoptera contains the peptide code-named Psein-AKH, the Anisozygoptera and the families Aeshnidae, Cordulegastridae and Macromiidae of the Anisoptera contain Anaim-AKH, whereas Gomphidae, Corduliidae (with the exception of Syncordulia gracilis) and Libellulidae contain Libau-AKH; one species of Libellulidae has Erysi-AKH, a very conservative modification of Libau-AKH (one point mutation). When these structural data are interpreted in conjunction with existing phylogenies of Odonata, they support the following: (1) Zygoptera are monophyletic and not paraphyletic. (2) Anisozygoptera and Anisoptera are sister groups and contain the ancestral Anaim-AKH which is independently and convergently mutated to Libau-AKH in Gomphidae and Libellulidae. (3) The Corduliidae are of special interest. Only Corduliidae sensu stricto appear to contain Libau-AKH, other species placed into this family by most authorities contain the ancestral Anaim-AKH. Possibly, assignments of AKHs can untangle the paraphyly of this family.

A secondary structural model of the 28S rRNA expansion segments D2 and D3 from rootworms and related leaf beetles (Coleoptera: Chrysomelidae; Galerucinae)

We analysed the secondary structure of two expansion segments (D2, D3) of the 28S rRNA gene from 229 leaf beetles (Coleoptera: Chrysomelidae), the majority of which are in the subfamily Galerucinae. The sequences were compared in a multiple sequence alignment, with secondary structure inferred primarily from the compensatory base changes in the conserved helices of the rRNA molecules. This comparative approach yielded thirty helices comprised of base pairs with positional covariation. Based on these leaf beetle sequences, we report an annotated secondary structural model for the D2 and D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related beetle taxa. This predicted structure, consisting of seven major compound helices, is mostly consistent with previously proposed models for the D2 and D3 expansion segments in insects. Despite a lack of conservation in the primary structure of these regions of insect 28S rRNA, the evolution of the secondary structure of these seven major motifs may be informative above the nucleotide level for higher-order phylogeny reconstruction of major insect lineages.