Assessing the odd secondary structural properties of nuclear small subunit ribosomal RNA sequences (18S) of the twisted-wing parasites (Insecta: Strepsiptera)

Substantial rearrangements, single nucleotide frameshift deletion and low diversity in mitogenome of Wolbachia-infected strepsipteran endoparasitoid in comparison to its tephritid hosts

Insect mitogenome organisation is highly conserved, yet, some insects, especially with parasitic life cycles, have rearranged mitogenomes. Furthermore, intraspecific mitochondrial diversity can be reduced by fitness-affecting bacterial endosymbionts like Wolbachia due to their maternal coinheritance with mitochondria. We have sequenced mitogenomes of the Wolbachia-infected endoparasitoid Dipterophagus daci (Strepsiptera: Halictophagidae) and four of its 22 known tephritid fruit fly host species using total genomic extracts of parasitised flies collected across > 700 km in Australia. This halictophagid mitogenome revealed extensive rearrangements relative to the four fly mitogenomes which exhibited the ancestral insect mitogenome pattern. Compared to the only four available other strepsipteran mitogenomes, the D. daci mitogenome had additional transpositions of one rRNA and two tRNA genes, and a single nucleotide frameshift deletion in nad5 requiring translational frameshifting or, alternatively, resulting in a large protein truncation. Dipterophagus daci displays an almost completely endoparasitic life cycle when compared to Strepsiptera that have maintained the ancestral state of free-living adults. Our results support the hypothesis that the transition to extreme endoparasitism evolved together with increased levels of mitogenome changes. Furthermore, intraspecific mitogenome diversity was substantially smaller in D. daci than the parasitised flies suggesting Wolbachia reduced mitochondrial diversity because of a role in D. daci fitness.

The SSU rRNA secondary structures of the Plagiorchiida species (Digenea), its applications in systematics and evolutionary inferences

The small subunit ribosomal RNA (SSU rRNA) is widely used phylogenetic marker in broad groups of organisms and its secondary structure increasingly attracts the attention of researchers as supplementary tool in sequence alignment and advanced phylogenetic studies. Its comparative analysis provides a great contribution to evolutionary biology, allowing find out how the SSU rRNA secondary structure originated, developed and evolved. Herein, we provide the first data on the putative SSU rRNA secondary structures of the Plagiorchiida species. The structures were found to be quite conserved across broad range of species studied, well compatible with those of others eukaryotic SSU rRNA and possessed some peculiarities: cross-shaped structure of the ES6b, additional shortened ES6c2 helix, and elongated ES6a helix and h39 + ES9 region. The secondary structures of variable regions ES3 and ES7 appeared to be tissue-specific while ES6 and ES9 were specific at a family level allowing considering them as promising markers for digenean systematics. Their uniqueness more depends on the length than on the nucleotide diversity of primary sequences which evolutionary rates well differ. The findings have important implications for understanding rRNA evolution, developing molecular taxonomy and systematics of Plagiorchiida as well as for constructing new anthelmintic drugs.

The role of floral oils in the evolution of apid bees (Hymenoptera: Apidae)

Most bees collect pollen and nectar for their larvae, while some also collect other resources. We investigated the evolution of floral oil-collecting behaviour in the Apidae and the evolutionary effects of floral oils on host brood cells for cuckoo bees. Focusing on apid bee phylogeny, we reconstructed the evolution of floral oil collection by females, use of floral oils in cell construction and the inclusion of oils in provisioning immatures. Ancestral character reconstruction demonstrates that floral oil-collecting behaviour arose four times independently. We also found that in cuckoo bees, parasitization of oil-collecting bees arose three times (including one secondary return) in Apidae. Except for Ctenoplectrina, oil cuckoo bees are all closely related to each other, forming an independent clade within the Nomadinae. Analysis of evolutionary transition rates indicates that there is a greater tendency for switching from an oil-collecting host to a non-oil-collecting host than the reverse. In apid bees, evolutionary transition rates are higher for switching to cuckoo behaviour from an ancestral lineage in which females collect floral oils than from other pollen-collecting lineages. We conclude that adaptation to oil collection is advantageous for pollen-collecting bees, and that the origin of oil cuckoo bees from non-oil cuckoo bees is constrained.

Strepsiptera, phylogenomics and the long branch attraction problem

Insect phylogeny has recently been the focus of renewed interest as advances in sequencing techniques make it possible to rapidly generate large amounts of genomic or transcriptomic data for a species of interest. However, large numbers of markers are not sufficient to guarantee accurate phylogenetic reconstruction, and the choice of the model of sequence evolution as well as adequate taxonomic sampling are as important for phylogenomic studies as they are for single-gene phylogenies. Recently, the sequence of the genome of a strepsipteran has been published and used to place Strepsiptera as sister group to Coleoptera. However, this conclusion relied on a data set that did not include representatives of Neuropterida or of coleopteran lineages formerly proposed to be related to Strepsiptera. Furthermore, it did not use models that are robust against the long branch attraction artifact. Here we have sequenced the transcriptomes of seven key species to complete a data set comprising 36 species to study the higher level phylogeny of insects, with a particular focus on Neuropteroidea (Coleoptera, Strepsiptera, Neuropterida), especially on coleopteran taxa considered as potential close relatives of Strepsiptera. Using models robust against the long branch attraction artifact we find a highly resolved phylogeny that confirms the position of Strepsiptera as a sister group to Coleoptera, rather than as an internal clade of Coleoptera, and sheds new light onto the phylogeny of Neuropteroidea.

A unique box in 28S rRNA is shared by the enigmatic insect order Zoraptera and Dictyoptera

The position of the Zoraptera remains one of the most challenging and uncertain concerns in ordinal-level phylogenies of the insects. Zoraptera have been viewed as having a close relationship with five different groups of Polyneoptera, or as being allied to the Paraneoptera or even Holometabola. Although rDNAs have been widely used in phylogenetic studies of insects, the application of the complete 28S rDNA are still scattered in only a few orders. In this study, a secondary structure model of the complete 28S rRNAs of insects was reconstructed based on all orders of Insecta. It was found that one length-variable region, D3-4, is particularly distinctive. The length and/or sequence of D3-4 is conservative within each order of Polyneoptera, but it can be divided into two types between the different orders of the supercohort, of which the enigmatic order Zoraptera and Dictyoptera share one type, while the remaining orders of Polyneoptera share the other. Additionally, independent evidence from phylogenetic results support the clade (Zoraptera+Dictyoptera) as well. Thus, the similarity of D3-4 between Zoraptera and Dictyoptera can serve as potentially valuable autapomorphy or synapomorphy in phylogeny reconstruction. The clades of (Plecoptera+Dermaptera) and ((Grylloblattodea+Mantophasmatodea)+(Embiodea+Phasmatodea)) were also recovered in the phylogenetic study. In addition, considering the other studies based on rDNAs, this study reached the highest congruence with previous phylogenetic studies of Holometabola based on nuclear protein coding genes or morphology characters. Future comparative studies of secondary structures across deep divergences and additional taxa are likely to reveal conserved patterns, structures and motifs that can provide support for major phylogenetic lineages.

Mermithid nematodes found in adult Anopheles from southeastern Senegal

Background

Over two dozen mermithid nematodes have been described parasitizing mosquitoes worldwide, however, only two species were found in Africa. Mermithid nematodes kill their mosquito host upon emergence, which suggests that they could be developed as biological control agents of mosquitoes. Both Romanomermis culicivorax and Romanomermis iyengari have been reared for mass release to control numerous Anopheles species vector populations, and in one instance this may have led to reduced malaria prevalence in a human population.

Methods

Anopheles mosquitoes were collected during a malaria study in southeastern Senegal. Two different adult blood fed mosquitoes had a single mermithid nematode emerge from their anus while they were being held post-capture. Primers from the 18 S rDNA were developed to sequence nematode DNA and screen mosquitoes for mermithid DNA. 18 S rDNA from the Senegalese mermithid and other mermithid entries in GenBank were used to create a Maximum Parsimony tree of the Mermithidae family.

Results

The mermithid was present in 1.8% (10/551) of the sampled adult Anopheles species in our study area. The mermithid was found in An. gambiae s.s., An. funestus, and An. rufipes from the villages of Ndebou, Boundoucondi, and Damboucoye. Maximum parsimony analysis confirmed that the nematode parasites found in Anopheles were indeed mermithid parasites, and of the mermithid sequences available in GenBank, they are most closely related to Strelkovimermis spiculatus.

Conclusions

To our knowledge, this is the first report of mermithids from adult Anopheles mosquitoes in Senegal. The mermithid appears to infect Anopheles mosquitoes that develop in diverse larval habitats. Although maximum parsimony analysis determined the mermithid was closely related to Strelkovimermis spiculatus, several characteristics of the mermithid were more similar to the Empidomermis genus. Future mermithid isolations will hopefully allow: formal taxonomic identification, laboratory colonization, determination of life history traits and species specificity, and characterize its usefulness as a biological control agent.

A molecular phylogeny of the Chalcidoidea (Hymenoptera)

Chalcidoidea (Hymenoptera) are extremely diverse with more than 23,000 species described and over 500,000 species estimated to exist. This is the first comprehensive phylogenetic analysis of the superfamily based on a molecular analysis of 18S and 28S ribosomal gene regions for 19 families, 72 subfamilies, 343 genera and 649 species. The 56 outgroups are comprised of Ceraphronoidea and most proctotrupomorph families, including Mymarommatidae. Data alignment and the impact of ambiguous regions are explored using a secondary structure analysis and automated (MAFFT) alignments of the core and pairing regions and regions of ambiguous alignment. Both likelihood and parsimony approaches are used to analyze the data. Overall there is no impact of alignment method, and few but substantial differences between likelihood and parsimony approaches. Monophyly of Chalcidoidea and a sister group relationship between Mymaridae and the remaining Chalcidoidea is strongly supported in all analyses. Either Mymarommatoidea or Diaprioidea are the sister group of Chalcidoidea depending on the analysis. Likelihood analyses place Rotoitidae as the sister group of the remaining Chalcidoidea after Mymaridae, whereas parsimony nests them within Chalcidoidea. Some traditional family groups are supported as monophyletic (Agaonidae, Eucharitidae, Encyrtidae, Eulophidae, Leucospidae, Mymaridae, Ormyridae, Signiphoridae, Tanaostigmatidae and Trichogrammatidae). Several other families are paraphyletic (Perilampidae) or polyphyletic (Aphelinidae, Chalcididae, Eupelmidae, Eurytomidae, Pteromalidae, Tetracampidae and Torymidae). Evolutionary scenarios discussed for Chalcidoidea include the evolution of phytophagy, egg parasitism, sternorrhynchan parasitism, hypermetamorphic development and heteronomy.

The first molecular phylogeny of Strepsiptera (Insecta) reveals an early burst of molecular evolution correlated with the transition to endoparasitism

A comprehensive model of evolution requires an understanding of the relationship between selection at the molecular and phenotypic level. We investigate this in Strepsiptera, an order of endoparasitic insects whose evolutionary biology is poorly studied. We present the first molecular phylogeny of Strepsiptera, and use this as a framework to investigate the association between parasitism and molecular evolution. We find evidence of a significant burst in the rate of molecular evolution in the early history of Strepsiptera. The evolution of morphological traits linked to parasitism is significantly correlated with the pattern in molecular rate. The correlated burst in genotypic-phenotypic evolution precedes the main phase of strepsipteran diversification, which is characterised by the return to a low and even molecular rate, and a period of relative morphological stability. These findings suggest that the transition to endoparasitism led to relaxation of selective constraint in the strepsipteran genome. Our results indicate that a parasitic lifestyle can affect the rate of molecular evolution, although other causal life-history traits correlated with parasitism may also play an important role.

Large-scale analysis of phylogenetic search behavior

Phylogenetic analysis is used in all branches of biology with applications ranging from studies on the origin of human populations to investigations of the transmission patterns of HIV. Most phylogenetic analyses rely on effective heuristics for obtaining accurate trees. However, relatively little work has been done to analyze quantitatively the behavior of phylogenetic heuristics in tree space. A better understanding of local search behavior can facilitate the design of better heuristics, which ultimately lead to more accurate depictions of the true evolutionary relationships. In this paper, we present new and novel insights into local search behavior for maximum parsimony on three biological datasets consisting of 44, 60, and 174 taxa. By analyzing all trees from search, we find that, as the search algorithm climbs the hill to local optima, the trees in the neighborhood surrounding the current solution improve as well. Furthermore, the search is quite robust to a small number of randomly selected neighbors. Thus, our work shows how to gain insights into the behavior of local search algorithm by exploring a large diverse collection of trees.

Origin of clothing lice indicates early clothing use by anatomically modern humans in Africa

Clothing use is an important modern behavior that contributed to the successful expansion of humans into higher latitudes and cold climates. Previous research suggests that clothing use originated anywhere between 40,000 and 3 Ma, though there is little direct archaeological, fossil, or genetic evidence to support more specific estimates. Since clothing lice evolved from head louse ancestors once humans adopted clothing, dating the emergence of clothing lice may provide more specific estimates of the origin of clothing use. Here, we use a Bayesian coalescent modeling approach to estimate that clothing lice diverged from head louse ancestors at least by 83,000 and possibly as early as 170,000 years ago. Our analysis suggests that the use of clothing likely originated with anatomically modern humans in Africa and reinforces a broad trend of modern human developments in Africa during the Middle to Late Pleistocene.

Phylogenetic relationships among insect orders based on three nuclear protein-coding gene sequences

Many attempts to resolve the phylogenetic relationships of higher groups of insects have been made based on both morphological and molecular evidence; nonetheless, most of the interordinal relationships of insects remain unclear or are controversial. As a new approach, in this study we sequenced three nuclear genes encoding the catalytic subunit of DNA polymerase delta and the two largest subunits of RNA polymerase II from all insect orders. The predicted amino acid sequences (In total, approx. 3500 amino acid sites) of these proteins were subjected to phylogenetic analyses based on the maximum likelihood and Bayesian analysis methods with various models. The resulting trees strongly support the monophyly of Palaeoptera, Neoptera, Polyneoptera, and Holometabola, while within Polyneoptera, the groupings of Isoptera/"Blattaria"/Mantodea (Superorder Dictyoptera), Dictyoptera/Zoraptera, Dermaptera/Plecoptera, Mantophasmatodea/Grylloblattodea, and Embioptera/Phasmatodea are supported. Although Paraneoptera is not supported as a monophyletic group, the grouping of Phthiraptera/Psocoptera is robustly supported. The interordinal relationships within Holometabola are well resolved and strongly supported that the order Hymenoptera is the sister lineage to all other holometabolous insects. The other orders of Holometabola are separated into two large groups, and the interordinal relationships of each group are (((Siphonaptera, Mecoptera), Diptera), (Trichoptera, Lepidoptera)) and ((Coleoptera, Strepsiptera), (Neuroptera, Raphidioptera, Megaloptera)). The sister relationship between Strepsiptera and Diptera are significantly rejected by all the statistical tests (AU, KH and wSH), while the affinity between Hymenoptera and Mecopterida are significantly rejected only by AU and KH tests. Our results show that the use of amino acid sequences of these three nuclear genes is an effective approach for resolving the relationships of higher groups of insects.

Phylogeny of the millipede order Spirobolida (Arthropoda: Diplopoda: Helminthomorpha)

This study examines relationships within the millipede order Spirobolida using an exemplar approach, sampling within families to maximize geographical and morphological diversity; due to lack of available material, Allopocockiidae and Hoffmanobolidae were not included in analyses. The focus of this study was to test monophyly of the order, the suborders, and the families of Spirobolida and to propose interfamilial relationships using morphological and molecular data in a total-evidence approach. Both maximum-parsimony analyses and Bayesian inference were employed to analyse two datasets consisting of combined morphological and molecular data, one aligned using progressive alignment methods and the second aligned by secondary structure models. Rhinocricidae was recovered sister to all remaining spirobolidan millipedes and is elevated to suborder status as suborder Rhinocricidea. Trigoniulidea was recovered as monophyletic as was Spirobolidea excluding Rhinocricidae; Spirobolidea is redefined to reflect this change. All previously recognized families were recovered, with the exception of Spirobolidae; in all instances, this family was paraphyletic or part of a polytomy that lacked sufficient resolution to assess its monophyly. The results reaffirm much of the existing taxonomic foundation within Spirobolida. This study provides the first phylogenetic test of higher-level relationships within Spirobolida and will serve as a foundation for future work in this group at finer levels. © The Willi Hennig Society 2010.

Positional and neighboring base pair effects on the thermodynamic stability of RNA single mismatches

Many naturally occurring RNA structures contain single mismatches, many of which occur near the ends of helices. However, previous thermodynamic studies have focused their efforts on thermodynamically characterizing centrally placed single mismatches. Additionally, algorithms currently used to predict secondary structure from sequence are based on two assumptions for predicting the stability of RNA duplexes containing this motif. It has been assumed that the thermodynamic contribution of small RNA motifs is independent of both its position in the duplex and the identity of the non-nearest neighbors. Thermodynamically characterizing single mismatches three nucleotides from both the 3' and 5' ends (i.e., off-center) of an RNA duplex and comparing these results to those of the same single mismatch-nearest neighbor combination centrally located have allowed for the investigation of these effects. The thermodynamic contributions of 13 single mismatch-nearest neighbor combinations are reported, but only nine combinations are studied at all three duplex positions and are used to determine trends and patterns. In general, the 5'- and 3'-shifted single mismatches are relatively similar, on average, and more favorable in free energy than centrally placed single mismatches. However, close examination and comparison shows there are several associated idiosyncrasies with these identified general trends. These peculiarities may be due, in part, to the identities of the single mismatch, the nearest neighbors, and the non-nearest neighbors, along with the effects of the single mismatch position in the duplex. The prediction algorithm recently proposed by Davis and Znosko [Davis, A. R., and Znosko, B. M. (2008) Biochemistry 47, 10178-10187] is used to predict the thermodynamic parameters of single mismatch contribution, and those values are compared to the measured values presented here. This comparison suggests the proposed model is a good approximation but could be improved by the addition of parameters that account for positional and/or non-nearest neighbor effects. However, more data are required to improve our understanding of these effects and to accurately account for them.

Ribosomal protein genes of holometabolan insects reject the Halteria, instead revealing a close affinity of Strepsiptera with Coleoptera

The phylogenetic relationships among holometabolan insect orders remain poorly known, despite a wealth of previous studies. In particular, past attempts to clarify the sister-group of the enigmatic order Strepsiptera with rRNA genes have led to intense debate about long-branch attraction (the 'Strepsiptera problem'), without resolving the taxonomic question at hand. Here, we appealed to alternative nuclear sequences of 27 ribosomal proteins (RPs) to generate a data matrix of 10,731 nucleotides for 22 holometabolan taxa, including two strepsipteran species. Phylogenetic relationships among holometabolan insects were analyzed under several nucleotide-coding schemes to explore differences in signal and systematic biases. Saturation and compositional bias particularly affected third positions, which greatly differed in AT content (18-72%). Such confounding factors were best reduced by R-Y coding and removal of third codon positions, resulting in more strongly supported topologies, whereas amino acid coding gave poor resolution. The placement of Strepsiptera with Coleoptera (the Coleopterida) was recovered under most coding schemes and analytical methods, if often with modest support and ambiguity. In contrast, an alternative sister-group with Diptera (the Halteria) was only found in one analysis using parsimony, and weakly supported. The topologies here generally support a Coleoptera+Strepsiptera as sister-group to Mecopterida (Siphonaptera+Mecoptera+Diptera+Lepidoptera+Trichoptera), while Hymenoptera were always recovered as sister-group to the remaining Holometabola.

The mitochondrial genome of the 'twisted-wing parasite' Mengenilla australiensis (Insecta, Strepsiptera): a comparative study

BACKGROUND

Strepsiptera are an unusual group of sexually dimorphic, entomophagous parasitoids whose evolutionary origins remain elusive. The lineage leading to Mengenilla australiensis (Family Mengenillidae) is the sister group to all remaining extant strepsipterans. It is unique in that members of this family have retained a less derived condition, where females are free-living from pupation onwards, and are structurally much less simplified. We sequenced almost the entire mitochondrial genome of M. australiensis as an important comparative data point to the already available genome of its distant relative Xenos vesparum (Family Xenidae). This study represents the first in-depth comparative mitochondrial genomic analysis of Strepsiptera.

RESULTS

The partial genome of M. australiensis is presented as a 13421 bp fragment, across which all 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes and 18 transfer RNA (tRNA) sequences are identified. Two tRNA translocations disrupt an otherwise ancestral insect mitochondrial genome order. A+T content is measured at 84.3%, C-content is also very skewed. Compared with M. australiensis, codon bias in X. vesparum is more balanced. Interestingly, the size of the protein coding genome is truncated in both strepsipterans, especially in X. vesparum which, uniquely, has 4.3% fewer amino acids than the average holometabolan complement. A revised assessment of mitochondrial rRNA secondary structure based on comparative structural considerations is presented for M. australiensis and X. vesparum.

CONCLUSIONS

The mitochondrial genome of X. vesparum has undergone a series of alterations which are probably related to an extremely derived lifestyle. Although M. australiensis shares some of these attributes; it has retained greater signal from the hypothetical most recent common ancestor (MRCA) of Strepsiptera, inviting the possibility that a shift in the mitochondrial selective environment might be related to the specialization accompanying the evolution of a small, morphologically simplified completely host-dependent lifestyle. These results provide useful insights into the nature of the evolutionary transitions that accompanied the emergence of Strepsiptera, but we emphasize the need for adequate sampling across the order in future investigations concerning the extraordinary developmental and evolutionary origins of this group.

Single-copy nuclear genes resolve the phylogeny of the holometabolous insects

BACKGROUND

Evolutionary relationships among the 11 extant orders of insects that undergo complete metamorphosis, called Holometabola, remain either unresolved or contentious, but are extremely important as a context for accurate comparative biology of insect model organisms. The most phylogenetically enigmatic holometabolan insects are Strepsiptera or twisted wing parasites, whose evolutionary relationship to any other insect order is unconfirmed. They have been controversially proposed as the closest relatives of the flies, based on rDNA, and a possible homeotic transformation in the common ancestor of both groups that would make the reduced forewings of Strepsiptera homologous to the reduced hindwings of Diptera. Here we present evidence from nucleotide sequences of six single-copy nuclear protein coding genes used to reconstruct phylogenetic relationships and estimate evolutionary divergence times for all holometabolan orders.

RESULTS

Our results strongly support Hymenoptera as the earliest branching holometabolan lineage, the monophyly of the extant orders, including the fleas, and traditionally recognized groupings of Neuropteroidea and Mecopterida. Most significantly, we find strong support for a close relationship between Coleoptera (beetles) and Strepsiptera, a previously proposed, but analytically controversial relationship. Exploratory analyses reveal that this relationship cannot be explained by long-branch attraction or other systematic biases. Bayesian divergence times analysis, with reference to specific fossil constraints, places the origin of Holometabola in the Carboniferous (355 Ma), a date significantly older than previous paleontological and morphological phylogenetic reconstructions. The origin and diversification of most extant insect orders began in the Triassic, but flourished in the Jurassic, with multiple adaptive radiations producing the astounding diversity of insect species for which these groups are so well known.

CONCLUSION

These findings provide the most complete evolutionary framework for future comparative studies on holometabolous model organisms and contribute strong evidence for the resolution of the 'Strepsiptera problem', a long-standing and hotly debated issue in insect phylogenetics.

Using tree diversity to compare phylogenetic heuristics

BACKGROUND

Evolutionary trees are family trees that represent the relationships between a group of organisms. Phylogenetic heuristics are used to search stochastically for the best-scoring trees in tree space. Given that better tree scores are believed to be better approximations of the true phylogeny, traditional evaluation techniques have used tree scores to determine the heuristics that find the best scores in the fastest time. We develop new techniques to evaluate phylogenetic heuristics based on both tree scores and topologies to compare Pauprat and Rec-I-DCM3, two popular Maximum Parsimony search algorithms.

RESULTS

Our results show that although Pauprat and Rec-I-DCM3 find the trees with the same best scores, topologically these trees are quite different. Furthermore, the Rec-I-DCM3 trees cluster distinctly from the Pauprat trees. In addition to our heatmap visualizations of using parsimony scores and the Robinson-Foulds distance to compare best-scoring trees found by the two heuristics, we also develop entropy-based methods to show the diversity of the trees found. Overall, Pauprat identifies more diverse trees than Rec-I-DCM3.

CONCLUSION

Overall, our work shows that there is value to comparing heuristics beyond the parsimony scores that they find. Pauprat is a slower heuristic than Rec-I-DCM3. However, our work shows that there is tremendous value in using Pauprat to reconstruct trees-especially since it finds identical scoring but topologically distinct trees. Hence, instead of discounting Pauprat, effort should go in improving its implementation. Ultimately, improved performance measures lead to better phylogenetic heuristics and will result in better approximations of the true evolutionary history of the organisms of interest.

Host-parasitoid associations in Strepsiptera

Strepsiptera are obligate endoparasitoids that exhibit extreme sexual dimorphism and parasitize seven orders and 33 families of Insecta. The adult males and the first instar larvae in the Mengenillidia and Stylopidia are free-living, whereas the adult females in Mengenillidia are free-living but in the suborder Stylopidia they remain endoparasitic in the host. Parasitism occurs at the host larval/nymphal stage and continues in a mobile host until that host's adult stage. The life of the host is lengthened to allow the male strepsipteran to complete maturation and the viviparous female to release the first instar larvae when the next generation of the host's larvae/nymphs has been produced. The ability of strepsipterans to parasitize a wide range of hosts, in spite of being endoparasitoids, is perhaps due to their unique immune avoidance system. Aspects of virulence, heterotrophic heteronomy in the family Myrmecolacidae, cryptic species, genomics, immune response, and behavior of stylopized hosts are discussed in this chapter.

Phylogenetic comparison of local length plasticity of the small subunit of nuclear rDNAs among all Hexapoda orders and the impact of hyper-length-variation on alignment

The SSU nrDNA (18S), is one of the most frequently sequenced molecular markers in phylogenetic studies. However, the length-hyper-variation at multiple positions of this gene can affect the accuracy of alignment greatly and this length variation makes alignment across arthropod orders a serious problem. The analyses of Hexapoda phylogeny is such a case. A more clear recognition of the distribution of the length-variable-regions is needed. In this study, the secondary structure of some length-variable-regions in the SSU nrRNA of Arthropoda was adjusted by the principle of co-variation. It is found that the extent of plasticity of some length-variable-region can extraordinarily be higher than 600 bases in hexapods. And the numbers of hyper length-variable-regions are largest in Strepsiptera and Sternorrhyncha (Hemiptera). Our study shows that some length-variable-regions can serve as synapomorphies for some groups. The phylogenetic comparison also suggested that the expansion of a lateral bulge could be the origin of a helix.

Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura)

Cospeciation between hosts and parasites offers a unique opportunity to use information from parasites to infer events in host evolutionary history. Although lice (Insecta: Phthiraptera) are known to cospeciate with their hosts and have frequently served as important markers to infer host evolutionary history, most molecular studies are based on only one or two markers. Resulting phylogenies may, therefore, represent gene histories (rather than species histories), and analyses of multiple molecular markers are needed to increase confidence in the results of phylogenetic analyses. Herein, we phylogenetically examine nine molecular markers in primate sucking lice (Phthiraptera: Anoplura) and we use these markers to estimate divergence times among louse lineages. Individual and combined analyses of these nine markers are, for the most part, congruent, supporting relationships hypothesized in previous studies. Only one marker, the nuclear protein-coding gene Histone 3, has a significantly different tree topology compared to the other markers. The disparate evolutionary history of this marker, however, has no significant effect on topology or nodal support in the combined phylogenetic analyses. Therefore, phylogenetic results from the combined data set likely represent a solid hypothesis of species relationships. Additionally, we find that simultaneous use of multiple markers and calibration points provides the most reliable estimates of louse divergence times, in agreement with previous studies estimating divergences among species. Estimates of phylogenies and divergence times also allow us to verify the results of [Reed, D.L., Light, J.E., Allen, J.M., Kirchman, J.J., 2007. Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biol. 5, 7.]; there was probable contact between gorilla and archaic hominids roughly 3 Ma resulting in a host switch of Pthirus lice from gorillas to archaic hominids. Thus, these results provide further evidence that data from cospeciating organisms can yield important information about the evolutionary history of their hosts.

What's in a name: the taxonomic status of human head and body lice

Human head lice (Anoplura: Pediculidae: Pediculus) are pandemic, parasitizing countless school children worldwide due to the evolution of insecticide resistance, and human body (clothing) lice are responsible for the deaths of millions as a result of vectoring several deadly bacterial pathogens. Despite the obvious impact these lice have had on their human hosts, it is unclear whether head and body lice represent two morphological forms of a single species or two distinct species. To assess the taxonomic status of head and body lice, we provide a synthesis of publicly available molecular data in GenBank, and we compare phylogenetic and population genetic methods using the most diverse geographic and molecular sampling presently available. Our analyses find reticulated networks, gene flow, and a lack of reciprocal monophyly, all of which indicate that head and body lice do not represent genetically distinct evolutionary units. Based on these findings, as well as inconsistencies of morphological, behavioral, and ecological variability between head and body lice, we contend that no known species concept would recognize these louse morphotypes as separate species. We recommend recognizing head and body lice as morphotypes of a single species, Pediculus humanus, until compelling new data and analyses (preferably analyses of fast evolving nuclear markers in a coalescent framework) indicate otherwise.

Large expansion segments in 18S rDNA support a new sponge clade (Class Demospongiae, Order Haplosclerida)

Newly emerging molecular phylogenetic hypotheses involving the sponge Order Haplosclerida (Class Demospongiae) are far removed from traditional views on their classification using morphology. In the new grouping of marine haplosclerid taxa by molecular data all members of one highly supported clade were found to have three large indels in the 18S rRNA gene. These indels were not found in this gene in other marine haplosclerids or in any other demosponges analysed. These indels were found in the variable V4 and V7 region of the gene, had high GC contents and formed stable double stranded helices in the 18S rRNA secondary structure. These indels are very important synapomorphies, provide high support for an alternative taxonomic scheme and could help resolve the phylogeny of this order in conjunction with other phylogenetically informative characters.

Molecular evolution of rDNA in early diverging Metazoa: first comparative analysis and phylogenetic application of complete SSU rRNA secondary structures in Porifera

Background

The cytoplasmic ribosomal small subunit (SSU, 18S) ribosomal RNA (rRNA) is the most frequently-used gene for molecular phylogenetic studies. However, information regarding its secondary structure is neglected in most phylogenetic analyses. Incorporation of this information is essential in order to apply specific rRNA evolutionary models to overcome the problem of co-evolution of paired sites, which violates the basic assumption of the independent evolution of sites made by most phylogenetic methods. Information about secondary structure also supports the process of aligning rRNA sequences across taxa. Both aspects have been shown to increase the accuracy of phylogenetic reconstructions within various taxa.

Here, we explore SSU rRNA secondary structures from the three extant classes of Phylum Porifera (Grant, 1836), a pivotal, but largely unresolved taxon of early branching Metazoa. This is the first phylogenetic study of poriferan SSU rRNA data to date that includes detailed comparative secondary structure information for all three sponge classes.

Results

We found base compositional and structural differences in SSU rRNA among Demospongiae, Hexactinellida (glass sponges) and Calcarea (calcareous sponges). We showed that analyses of primary rRNA sequences, including secondary structure-specific evolutionary models, in combination with reconstruction of the evolution of unusual structural features, reveal a substantial amount of additional information. Of special note was the finding that the gene tree topologies of marine haplosclerid demosponges, which are inconsistent with the current morphology-based classification, are supported by our reconstructed evolution of secondary structure features. Therefore, these features can provide alternative support for sequence-based topologies and give insights into the evolution of the molecule itself. To encourage and facilitate the application of rRNA models in phylogenetics of early metazoans, we present 52 SSU rRNA secondary structures over the taxonomic range of Porifera in a database, along with some basic tools for relevant format-conversion.

Conclusion

We demonstrated that sophisticated secondary structure analyses can increase the potential phylogenetic information of already available rDNA sequences currently accessible in databases and conclude that the importance of SSU rRNA secondary structure information for phylogenetic reconstruction is still generally underestimated, at least among certain early branching metazoans.

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.

The secondary structure of the unusually long 18S ribosomal RNA of the myxozoan Sphaerospora truttae and structural evolutionary trends in the Myxozoa

The nearly complete 18S rRNA sequence of the myxozoan parasite Sphaerospora truttae shows an extraordinary length (2,552bp) in comparison with other myxozoans and with metazoans in general (average 1,800-1,900bp). The sequence shows nucleotide insertions in most variable regions of the 18S rRNA (V2, V4, V5 and V7), with especially large expansion segments in V4 and V7. In the myxozoans, nucleotide insertions and specific secondary structures in these regions of the gene were found to be strongly related to large scale phylogenetic clustering and thus with the invertebrate host type. Whereas expansion segments were generally found to be absent in the malacasporeans and the clade of primary marine myxozoan species, they occur in all taxa of the clade containing freshwater species, where they showed a consistent secondary structure throughout. The longest expansion segments occur in S. truttae, Sphaerospora elegans and Leptotheca ranae, which represent a clade that has emerged after the malacosporeans and before the radiation of all other myxozoan genera. These three species demonstrate structural links to the malacosporeans as well as other unique features. A smaller number of nucleotide insertions in different subhelices and specific secondary structures appear to have evolved independently in two marine genera, i.e. Ceratomyxa and Parvicapsula. The secondary structural elements of V4 and V7 of the myxozoan 18S rRNAs were found to be highly informative and revealed evolutionary trends of various regions of the gene hitherto unknown, since previous analyses have been based on primary sequence data excluding these regions. Furthermore, the unique features of the V4 region in S. truttae allowed for the design of a highly specific PCR assay for this species.

Molecular phylogeny and biogeography of an ancient Holarctic lineage of mygalomorph spiders (Araneae: Antrodiaetidae: Antrodiaetus)

The mygalomorph spider genera Antrodiaetus and Atypoides (Antrodiaetidae) belong to an ancient lineage that has persisted since at least the Cretaceous. These spiders display a classic disjunct Holarctic distribution with species in the eastern Palaearctic plus the western and eastern Nearctic. Prior phylogenetic analyses of this group have been proposed on the basis of morphology, but lack strong support and independent corroboration. Here we present the first phylogenetic analysis of species-level relationships based on molecular data obtained from the mitochondrial (cytochrome c oxidase subunit I) and nuclear (18S and 28S rRNA) genomes. Analyses corroborate earlier findings that Atypoides forms a paraphyletic grade with respect to Antrodiaetus, and consequently, that genus is formally synonymized under Antrodiaetus. In addition, our results support the relatively early divergence of Antrodiaetus roretzi. Antrodiaetus pacificus is "paraphyletic" with respect to the A. lincolnianus group and is likely an assemblage of numerous species. The final topology based on a combined molecular dataset, in conjunction with two different molecular dating techniques (penalized likelihood plus a Bayesian approach) and ancestral distribution reconstructions, was used to infer the historical biogeography of these spiders. Trans-Beringian and trans-Atlantic routes appear to account for the present-day distribution of Antrodiaetus in Japan and North America. Future studies on Antrodiaetus phylogeny will be used to address questions regarding morphological stasis and the evolution of quantitative morphological characters.

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.