Characterizing regions of ambiguous alignment caused by the expansion and contraction of hairpin-stem loops in ribosomal RNA molecules

A New Parasitic Archamoeba Causing Systemic Granulomatous Disease in Goldfish Extends the Diversity of Pathogenic Endolimax spp

Endolimax is a genus of intestinal amoebae which stands among the least known human protists. Previous studies on amoebic systemic granulomatosis of a marine fish (Solea senegalensis) resulted in the unexpected characterization of a new organism which was related to Endolimax and named E. piscium. The existence of multiple reports of systemic granulomatosis caused presumptively by unidentified amoebae in goldfish lead us to investigate the organism involved in goldfish disease. Analysed goldfish presented small whitish nodules in the kidney, which correspond to chronic granulomatous inflammatory reactions with a ring-layer of amoebae in the periphery. Amoebae were amitochondriate and were located in a parasitophorous vacuole within macrophages, as previous studies on this condition in goldfish and other freshwater fish pointed out. SSU rDNA characterization confirmed a new Endolimax lineage which appears closely related to E. piscium, but the molecular evidence, distinct pathological features and lack of ecological overlapping between the hosts support their assignment to a new species, E. carassius. The results support the existence of a considerable unexplored diversity of Endolimax spp. among fish, and their proper characterization can contribute to an understanding of Archamoebae evolution and pathogenic potential.

Molecular phylogenetic assessment of the tribal classification of Lamiinae (Coleoptera: Cerambycidae)

Lamiinae is the most diverse subfamily of longhorned beetles, with about 20,000 described species classified into 80 tribes. Most of the tribes of Lamiinae were proposed during the 19th century and the suprageneric classification of the subfamily has never been assessed under phylogenetic criteria. In this study, we present the first tribal-level phylogeny of Lamiinae, inferred from 130 terminals (representing 46 tribes, prioritizing generic type species of the tribes) and fragments of two mitochondrial and three nuclear markers (cox1, rrnL, Wg, CPS and LSU; 5,024 aligned positions in total). Analyses were performed under Maximum Likelihood and Bayesian methods based on two datasets: a dataset including all taxa available for the study, and a reduced dataset with 111 terminals where taxa only contributing with mitochondrial markers were excluded from the matrix. The monophyly of Lamiinae was corroborated in three of the four analyses and 11 of the 35 tribes with more than one species represented in the analyses were consistently recovered as monophyletic. However, 15 tribes were not retrieved as monophyletic, requiring a revision of their boundaries: Acanthocinini, Acanthoderini, Agapanthiini, Apomecynini, Desmiphorini, Dorcaschematini, Enicodini, Hemilophini, Monochamini, Onciderini, Parmenini, Phytoeciini, Pogonocherini, Pteropliini and Saperdini. Based on these results, when strong support values for paraphyly were recovered, we argue a number of tribe synonymies, including Moneilemini as synonym of Acanthocinini; Onocephalini of Onciderini; Dorcadionini, Gnomini, Monochamini and Rhodopinini of Lamiini; and Obereini and Phytoeciini of Saperdini. Other taxonomic changes proposed in this study based on the criterion of monophyly and supported by morphological characters include the transfer of Tricondyloides and Stenellipsis to Enicodini, and of Dylobolus stat. rest., which is removed as subgenus of Mecas and restituted as genus, to Hemilophini. Furthermore, our analyses suggest that Ostedes and Neohoplonotus should be removed from Acanthocinini and Parmenini, respectively, and Colobotheini should be redefined to encompass several genera currently placed in Acanthocinini.

Sequence complementarity at the ribosomal Peptidyl Transferase Centre implies self-replicating origin

A feasible scenario for the emergence of life requires the spontaneous materialization and sustainability of a proto-ribosome that could have catalysed the formation of the first peptides. Models of proto-ribosomes were derived from the ribosomal Peptidyl Transferase Centre (PTC) region, but the poor prebiotic copying abilities give rise to the question of their mode of replication. Here, complementarity is demonstrated in bacterial ribosomes, between nucleotides that constitute the two halves of the PTC cavity. The complementarity corroborates the dimeric nature of the proto-ribosome and is likely to underlie the symmetry of the PTC region. Furthermore, it indicates a simple and efficient replication mode; the strand of each monomer could have acted as a template for the synthesis of its counterpart, forming a self-replicating ribozyme.

Size-variable zone in V3 region of 16S rRNA

The size distribution of complete 16S-rRNA sequences from the SILVA-database and nucleotide shifts that might interfere with the secondary structure of the molecules were evaluated. Overall, 513,309 sequences recorded in SILVA were used to estimate the size of hypervariable regions of the gene. Redundant sequences were treated as a single sequence to achieve a better representation of the molecular diversity. Nucleotides found in each position in 95% of the sequences were considered the consensus sequences for different size-groups (consensus95). The sizes of different regions ranged from 96.7 to 283.1 nucleotides and had similar distribution patterns, except for the V3 region, which exhibited a bimodal distribution composed of 2 main peaks of 161 and 186 nt. The alignment of Consensuses95 of fractions 161 and 186 showed a high degree of similarity and conservation, except for the central positions (gap zone), where the sequence was highly variable and several deletions were observed. Structurally, the gap zone forms the central part of helix 17 (H17), and its extension was directly reflected in the size of this helix. H17 is part of a multihelix conjunction known as the 5-way junction (5 WJ), which is indispensable for 30 S ribosome assembly. However, because a drastic variation in the sequence size of V3 region occurs at a central position in loop H17 without affecting the base of the loop, it has no apparent effect on 5 WJ. Finally, considering that these differences were detected in non-redundant sequences, it can be concluded that this is not an uncommon or isolated event and that the V3 region is possibly more likely to mutate than are other regions.

Origin and higher-level diversification of acariform mites - evidence from nuclear ribosomal genes, extensive taxon sampling, and secondary structure alignment

BACKGROUND

Acariformes is the most species-rich and morphologically diverse radiation of chelicerate arthropods, known from the oldest terrestrial ecosystems. It is also a key lineage in understanding the evolution of this group, with the most vexing question whether mites, or Acari (Parasitiformes and Acariformes) is monophyletic. Previous molecular studies recovered Acari either as monophyletic or non-monophyletic, albeit with a limited taxon sampling. Similarly, relationships between basal acariform groups (include little-known, deep-soil 'endeostigmatan' mites) and major lineages of Acariformes (Sarcoptiformes, Prostigmata) are virtually unknown. We infer phylogeny of chelicerate arthropods, using a large and representative dataset, comprising all main in- and outgroups (228 taxa). Basal diversity of Acariformes is particularly well sampled. With this dataset, we conduct a series of phylogenetically explicit tests of chelicerate and acariform relationships and present a phylogenetic framework for internal relationships of acariform mites.

RESULTS

Our molecular data strongly support a diphyletic Acari, with Acariformes as the sister group to Solifugae (PP =1.0; BP = 100), the so called Poecilophysidea. Among Acariformes, some representatives of the basal group Endeostigmata (mainly deep-soil mites) were recovered as sister-groups to the remaining Acariformes (i. e., Trombidiformes + and most of Sarcoptiformes). Desmonomatan oribatid mites (soil and litter mites) were recovered as the monophyletic sister group of Astigmata (e. g., stored product mites, house dust mites, mange mites, feather and fur mites). Trombidiformes (Sphaerolichida + Prostigmata) is strongly supported (PP =1.0; BP = 98-100). Labidostommatina was inferred as the basal lineage of Prostigmata. Eleutherengona (e. g., spider mites) and Parasitengona (e. g., chiggers, fresh water mites) were recovered as monophyletic. By contrast, Eupodina (e. g., snout mites and relatives) was not. Marine mites (Halacaridae) were traditionally regarded as the sister-group to Bdelloidea (Eupodina), but our analyses show their close relationships to Parasitengona.

CONCLUSIONS

Non-trivial relationships recovered by our analyses with high support (i.e., basal arrangement of endeostigmatid lineages, the position of marine mites, polyphyly of Eupodina) had been proposed by previous underappreciated morphological studies. Thus, we update currently the accepted taxonomic classification to reflect these results: the superfamily Halacaroidea Murray, 1877 is moved from the infraorder Eupodina Krantz, 1978 to Anystina van der Hammen, 1972; and the subfamily Erythracarinae Oudemans, 1936 (formerly in Anystidae Oudemans, 1902) is elevated to family rank, Erythracaridae stat. ressur., leaving Anystidae only with the nominal subfamily. Our study also shows that a clade comprising early derivative Endeostigmata (Alycidae, Nanorchestidae, Nematalycidae, and maybe Alicorhagiidae) should be treated as a taxon with the same rank as Sarcoptiformes and Trombidiformes, and the scope of the superfamily Bdelloidea should be changed. Before turning those findings into nomenclatural changes, however, we consider that our study calls for (i) finding shared apomorphies of the early derivative Endeostigmata clade and the clade including the remaining Acariformes; (ii) a well-supported hypothesis for Alicorhagiidae placement; (iii) sampling the families Proterorhagiidae, Proteonematalycidae and Grandjeanicidae not yet included in molecular analyses; (iv) undertake a denser sampling of clades traditionally placed in Eupodina, Anystina (Trombidiformes) and Palaeosomata (Sarcoptiformes), since consensus networks and Internode certainty (IC) and IC All (ICA) indices indicate high levels of conflict in these tree regions. Our study shows that regions of ambiguous alignment may provide useful phylogenetic signal when secondary structure information is used to guide the alignment procedure and provides an R implementation to the Bayesian Relative Rates test.

Phylogeny and niche conservatism in North and Central American triatomine bugs (Hemiptera: Reduviidae: Triatominae), vectors of Chagas' disease

The niche conservatism hypothesis states that related species diverge in niche characteristics at lower rates than expected, given their lineage divergence. Here we analyze whether niche conservatism is a common pattern among vector species (Hemiptera: Reduviidae: Triatominae) of Trypanosoma cruzi that inhabit North and Central America, a highly heterogeneous landmass in terms of environmental gradients. Mitochondrial and nuclear loci were used in a multi-locus phylogenetic framework to reconstruct phylogenetic relationships among species and estimate time of divergence of selected clades to draw biogeographic inferences. Then, we estimated similarity between the ecological niche of sister species and tested the niche conservatism hypothesis using our best estimate of phylogeny. Triatoma is not monophyletic. A primary clade with all North and Central American (NCA) triatomine species from the genera Triatoma, Dipetalogaster, and Panstrongylus, was consistently recovered. Nearctic species within the NCA clade (T. p. protracta, T. r. rubida) diverged during the Pliocene, whereas the Neotropical species (T. phyllosoma, T. longipennis, T. dimidiata complex) are estimated to have diverged more recently, during the Pleistocene. The hypothesis of niche conservatism could not be rejected for any of six sister species pairs. Niche similarity between sister species best fits a retention model. While this framework is used here to infer niche evolution, it has a direct impact on spatial vector dynamics driven by human population movements, expansion of transportation networks and climate change scenarios.

Evolution of the ribosome at atomic resolution

The origins and evolution of the ribosome, 3-4 billion years ago, remain imprinted in the biochemistry of extant life and in the structure of the ribosome. Processes of ribosomal RNA (rRNA) expansion can be "observed" by comparing 3D rRNA structures of bacteria (small), yeast (medium), and metazoans (large). rRNA size correlates well with species complexity. Differences in ribosomes across species reveal that rRNA expansion segments have been added to rRNAs without perturbing the preexisting core. Here we show that rRNA growth occurs by a limited number of processes that include inserting a branch helix onto a preexisting trunk helix and elongation of a helix. rRNA expansions can leave distinctive atomic resolution fingerprints, which we call "insertion fingerprints." Observation of insertion fingerprints in the ribosomal common core allows identification of probable ancestral expansion segments. Conceptually reversing these expansions allows extrapolation backward in time to generate models of primordial ribosomes. The approach presented here provides insight to the structure of pre-last universal common ancestor rRNAs and the subsequent expansions that shaped the peptidyl transferase center and the conserved core. We infer distinct phases of ribosomal evolution through which ribosomal particles evolve, acquiring coding and translocation, and extending and elaborating the exit tunnel.

Equivalent indels--ambiguous functional classes and redundancy in databases

There is considerable interest in studying sequenced variations. However, while the positions of substitutions are uniquely identifiable by sequence alignment, the location of insertions and deletions still poses problems. Each insertion and deletion causes a change of sequence. Yet, due to low complexity or repetitive sequence structures, the same indel can sometimes be annotated in different ways. Two indels which differ in allele sequence and position can be one and the same, i.e. the alternative sequence of the whole chromosome is identical in both cases and, therefore, the two deletions are biologically equivalent. In such a case, it is impossible to identify the exact position of an indel merely based on sequence alignment. Thus, variation entries in a mutation database are not necessarily uniquely defined. We prove the existence of a contiguous region around an indel in which all deletions of the same length are biologically identical. Databases often show only one of several possible locations for a given variation. Furthermore, different data base entries can represent equivalent variation events. We identified 1,045,590 such problematic entries of insertions and deletions out of 5,860,408 indel entries in the current human database of Ensembl. Equivalent indels are found in sequence regions of different functions like exons, introns or 5' and 3' UTRs. One and the same variation can be assigned to several different functional classifications of which only one is correct. We implemented an algorithm that determines for each indel database entry its complete set of equivalent indels which is uniquely characterized by the indel itself and a given interval of the reference sequence.

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.

Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: evidence from case studies in the Metazoa

Background

Failure to account for covariation patterns in helical regions of ribosomal RNA (rRNA) genes has the potential to misdirect the estimation of the phylogenetic signal of the data. Furthermore, the extremes of length variation among taxa, combined with regional substitution rate variation can mislead the alignment of rRNA sequences and thus distort subsequent tree reconstructions. However, recent developments in phylogenetic methodology now allow a comprehensive integration of secondary structures in alignment and tree reconstruction analyses based on rRNA sequences, which has been shown to correct some of these problems. Here, we explore the potentials of RNA substitution models and the interactions of specific model setups with the inherent pattern of covariation in rRNA stems and substitution rate variation among loop regions.

Results

We found an explicit impact of RNA substitution models on tree reconstruction analyses. The application of specific RNA models in tree reconstructions is hampered by interaction between the appropriate modelling of covarying sites in stem regions, and excessive homoplasy in some loop regions. RNA models often failed to recover reasonable trees when single-stranded regions are excessively homoplastic, because these regions contribute a greater proportion of the data when covarying sites are essentially downweighted. In this context, the RNA6A model outperformed all other models, including the more parametrized RNA7 and RNA16 models.

Conclusions

Our results depict a trade-off between increased accuracy in estimation of interdependencies in helical regions with the risk of magnifying positions lacking phylogenetic signal. We can therefore conclude that caution is warranted when applying rRNA covariation models, and suggest that loop regions be independently screened for phylogenetic signal, and eliminated when they are indistinguishable from random noise. In addition to covariation and homoplasy, other factors, like non-stationarity of substitution rates and base compositional heterogeneity, can disrupt the signal of ribosomal RNA data. All these factors dictate sophisticated estimation of evolutionary pattern in rRNA data, just as other molecular data require similarly complicated (but different) corrections.

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 position of the acariform mites: sensitivity to homology assessment under total evidence

BACKGROUND

Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei.

RESULTS

The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infra-ordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes.

CONCLUSIONS

Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

A multi-locus phylogeny of Nectogalini shrews and influences of the paleoclimate on speciation and evolution

Nectogaline shrews are a major component of the small mammalian fauna of Europe and Asia, and are notable for their diverse ecology, including utilization of aquatic habitats. So far, molecular phylogenetic analyses including nectogaline species have been unable to infer a well-resolved, well-supported phylogeny, thus limiting the power of comparative evolutionary and ecological analyses of the group. Here, we employ Bayesian phylogenetic analyses of eight mitochondrial and three nuclear genes to infer the phylogenetic relationships of nectogaline shrews. We subsequently use this phylogeny to assess the genetic diversity within the genus Episoriculus, and determine whether adaptation to aquatic habitats evolved independently multiple times. Moreover, we both analyze the fossil record and employ Bayesian relaxed clock divergence dating analyses of DNA to assess the impact of historical global climate change on the biogeography of Nectogalini. We infer strong support for the polyphyly of the genus Episoriculus. We also find strong evidence that the ability to heavily utilize aquatic habitats evolved independently in both Neomys and Chimarrogale+Nectogale lineages. Our Bayesian molecular divergence analysis suggests that the early history of Nectogalini is characterized by a rapid radiation at the Miocene/Pliocene boundary, thus potentially explaining the lack of resolution at the base of the tree. Finally, we find evidence that nectogalines once inhabited northern latitudes, but the global cooling and desiccating events at the Miocene/Pliocene and Pliocene/Pleistocene boundaries and Pleistocene glaciation resulted in the migration of most Nectogalini lineages to their present day southern distribution.

Phylogeny of gammaproteobacteria

The phylogeny of the large bacterial class Gammaproteobacteria has been difficult to resolve. Here we apply a telescoping multiprotein approach to the problem for 104 diverse gammaproteobacterial genomes, based on a set of 356 protein families for the whole class and even larger sets for each of four cohesive subregions of the tree. Although the deepest divergences were resistant to full resolution, some surprising patterns were strongly supported. A representative of the Acidithiobacillales routinely appeared among the outgroup members, suggesting that in conflict with rRNA-based phylogenies this order does not belong to Gammaproteobacteria; instead, it (and, independently, "Mariprofundus") diverged after the establishment of the Alphaproteobacteria yet before the betaproteobacteria/gammaproteobacteria split. None of the orders Alteromonadales, Pseudomonadales, or Oceanospirillales were monophyletic; we obtained strong support for clades that contain some but exclude other members of all three orders. Extreme amino acid bias in the highly A+T-rich genome of Candidatus Carsonella prevented its reliable placement within Gammaproteobacteria, and high bias caused artifacts that limited the resolution of the relationships of other insect endosymbionts, which appear to have had multiple origins, although the unbiased genome of the endosymbiont Sodalis acted as an attractor for them. Instability was observed for the root of the Enterobacteriales, with nearly equal subsets of the protein families favoring one or the other of two alternative root positions; the nematode symbiont Photorhabdus was identified as a disruptor whose omission helped stabilize the Enterobacteriales root.

Marvaldi A. Systematic revision, cladistics and biogeography of the genus Neogutierrezia Martínez (Coleoptera:Scarabaeidae) and its phylogenetic placement in Rutelinae based on structural alignment of 28S rDNA sequences

The Argentinean endemic genus Neogutierrezia Martínez, 1953 (Scarabaeidae : Rutelinae) is revised and seven new species are described: N. bicolor Ocampo & Ruiz-Manzanos, sp. nov., N. chelii Ocampo & Ruiz-Manzanos, sp. nov., N. galileoi Ocampo & Ruiz-Manzanos, sp. nov., N. lagosae Ocampo & Ruiz-Manzanos, sp. nov., N. payuniensis Ocampo & Ruiz-Manzanos, sp. nov., N. scutata Ocampo & Ruiz-Manzanos, sp. nov. and N. variabilis Ocampo & Ruiz-Manzanos, sp. nov. Species status is proposed for N. affinis Martínez stat. nov., which was formerly considered as subspecies of N. mirabilis Martínez. The genus Neogutierrezia now includes 10 species distributed in the Monte biogeographic province in Argentina. In order to clarify the systematic placement of the genus Neogutierrezia, a molecular phylogenetic analysis was performed using structurally aligned 28S rDNA sequences (expansion segments D2 and D3) from 23 taxa in Scarabaeoidea, including two representative species of Neogutierrezia. This is the first report of an annotated secondary structure alignment of the D2 and D3 segments of 28S rRNA that spans a wide sample of scarabaeoids, providing a useful homology template for further phylogenetic reconstruction in these and closely related beetles. Results of the molecular parsimony analysis strongly indicate that the genus is closely related to members of the Rutelinae (Scarabaeidae), and thus Neogutierrezia Martínez is transferred from Melolonthinae: Pachydemini to Rutelinae, new placement. A morphological cladistic analysis of the genus was also undertaken, including all the 10 known species in the genus plus two outgroup taxa in Rutelinae, and based on 53 adult characters. The most-parsimonious cladogram provides evidence for the monophyly of the genus, which shows three main clades, distributed in Central Monte and Southern Monte. The adult morphology of the 10 species is described and a key is provided, along with illustrations of the diagnostic characters. The biogeography of species in the genus is discussed.

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.

A Monte Carlo approach successfully identifies randomness in multiple sequence alignments: a more objective means of data exclusion

Random similarity of sequences or sequence sections can impede phylogenetic analyses or the identification of gene homologies. Additionally, randomly similar sequences or ambiguously aligned sequence sections can negatively interfere with the estimation of substitution model parameters. Phylogenomic studies have shown that biases in model estimation and tree reconstructions do not disappear even with large data sets. In fact, these biases can become pronounced with more data. It is therefore important to identify possible random similarity within sequence alignments in advance of model estimation and tree reconstructions. Different approaches have been already suggested to identify and treat problematic alignment sections. We propose an alternative method that can identify random similarity within multiple sequence alignments (MSAs) based on Monte Carlo resampling within a sliding window. The method infers similarity profiles from pairwise sequence comparisons and subsequently calculates a consensus profile. This consensus profile represents a summary of all calculated single similarity profiles. In consequence, consensus profiles identify dominating patterns of nonrandom similarity or randomness within sections of MSAs. We show that the approach clearly identifies randomness in simulated and real data. After the exclusion of putative random sections, node support drastically improves in tree reconstructions of both data. It thus appears to be a powerful tool to identify possible biases of tree reconstructions or gene identification. The method is currently restricted to nucleotide data but will be extended to protein data in the near future.

Evolution of the parasitic wasp subfamily Rogadinae (Braconidae): phylogeny and evolution of lepidopteran host ranges and mummy characteristics

BACKGROUND

The braconid subfamily Rogadinae is a large, cosmopolitan group of endoparasitoid wasps characterised by 'mummifying' their lepidopteran host larvae, from which the adult subsequently emerges. Rogadines attack a variety of both macro- and microlepidopteran taxa, although the speciose genus Aleiodes almost exclusively attacks macrolepidopterans. Here, we investigate the phylogenetic history of the Rogadinae, revise their higher-level classification and assess the evolution of their host ranges and mummy types. We also assess the divergence times within the subfamily and discuss the reasons for the extraordinary evolutionary diversification of Aleiodes.

RESULTS

Our Bayesian analyses weakly support the monophyly of the subfamily. A clade comprising all Aleiodes species and some other taxa is not nested within the tribe Rogadini as previously supposed, but instead is recovered as sister to the Yeliconini, with the remaining Rogadini genera being recovered as sister to the Stiropiini. The Rogadinae is estimated to have originated during the mid to late Eocene, 36.1-51.62 MYA. Molecular dating gives a more recent origin for the Aleiodes clade (17.98-41.76 MYA) compared to the origins proposed for two of its principal lepidopteran host groups (Noctuidae: 60.7-113.4 MYA; Geometridae 48-62 MYA). The Bayesian ancestral reconstruction of the emergence habits from the mummified hosts weakly recovered an anterior emergence as the ancestral condition for the subfamily. Producing a hard mummy has evolved at various times independently, though most of the species with this biology belong to the Aleiodes clade.

CONCLUSION

Based on our results, we erect the tribe Aleiodini nov. to include Aleiodes and Heterogamus stat. rev. Cordylorhogas, Pholichora and Hemigyroneuron are synonymised with Aleiodes. The molecular dating of clades and the ancestral reconstruction of host ranges support the hypothesis that radiation within Aleiodes s. s. was due to host recruitment leading to host range expansion followed by speciation, and not to parasitoid-host coevolution. Within the Rogadinae, variation in the site of emergence from the mummified host probably evolved as a consequence of the mummy's site and mode of formation, and the extent of mummy tanning/hardness to the degree of protection needed in relation to the cost of providing it.

Comparative study of the evolution of nuclear ribosomal spacers incorporating secondary structure analyzes within Dodonaeoideae, Hippocastanoideae and Xanthoceroideae (Sapindaceae)

Ribosomal DNA internal transcribed spacers (ITS) and partial external transcribed spacers (ETSf) are popularly used to infer evolutionary hypotheses. However, there is generally little consideration given to the secondary structures of these small RNA molecules and their potential effects on sequence alignment and phylogenetic analyzes. Intergeneric relationships amongst three of the four major lineages in the Sapindaceae, the Dodonaeoideae, Hippcastanoideae and Xanthoceroideae were assessed by firstly, generating secondary structure predictions for ITS and partial ETSf sequences, and then these predictions were used to assist alignment of the sequences. Secondly, the alignment was analyzed using RNA specific models of sequence evolution that account for the variation in nucleotide evolution in the independent loops and covariating stems regions of the ribosomal spacers. These models and phylogeny drawn from these analyzes were compared with that from analyzes using 'traditional' 4-state models and previous plastid analyzes. These analyzes identified that paired-site models developed to deal specifically with stem structures in RNA encoding sequences more appropriately account for the evolutionary history of the sequences than traditional 4-state substitution models.

Comparative analysis of sequences and secondary structures of the rRNA internal transcribed spacer 2 (ITS2) in pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae): potential use of slippage-derived sequences in molecular systematics

A comparative analysis of ITS2 sequences and secondary structures in 89 species of pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae) was performed. The ITS2 folding pattern was highly conserved and comparable with the general model proposed for eukaryotes. Simple sequence repeats (SSRs) were responsible for most of the observed nucleotide variability (approximately 1-3%) and length variation (359-459bp). When plotted on secondary structures, SSRs mapped in expansion segments positioned at the apices of three ITS2 helices ('A', 'B' and 'D1') and appeared to have evolved under mechanisms of compensatory slippage. Homologies among SSRs nucleotides could not be unambiguously assigned, and thus were not useful to resolve phylogeny. However, slippage-derived motifs provided some preliminary genetic support for newly proposed taxonomic arrangements of several genera and subgenera of Meligethinae, corroborating existing morphological and ecological datasets.

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.

R, Quicke DLJ. Molecular phylogeny and historical biogeography of the cosmopolitan parasitic wasp subfamily Doryctinae (Hymenoptera:Braconidae)

The phylogenetic relationships among representatives of 64 genera of the cosmopolitan parasitic wasps of the subfamily Doryctinae were investigated based on nuclear 28S ribosomal (r) DNA (~650 bp of the D2–3 region) and cytochrome c oxidase I (COI) mitochondrial (mt) DNA (603 bp) sequence data. The molecular dating of selected clades and the biogeography of the subfamily were also inferred. The partitioned Bayesian analyses did not recover a monophyletic Doryctinae, though the relationships involved were only weakly supported. Strong evidence was found for rejecting the monophylies of both Doryctes Haliday, 1836 and Spathius Nees, 1818. Our results also support the recognition of the Rhaconotini as a valid tribe. A dispersal–vicariance analysis showed a strong geographical signal for the taxa included, with molecular dating estimates for the origin of Doryctinae and its subsequent radiation both occurring during the late Paleocene–early Eocene. The divergence time estimates suggest that diversification in the subfamily could have in part occurred as a result of continental break-up events that took place in the southern hemisphere, though more recent dispersal events account for the current distribution of several widespread taxa.

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.

Nematode small subunit phylogeny correlates with alignment parameters

The number of nuclear small subunit (SSU) ribosomal RNA (rRNA) sequences for Nematoda has increased dramatically in recent years, and although their use in constructing phylogenies has also increased, relatively little attention has been given to their alignment. Here we examined the sensitivity of the nematode SSU data set to different alignment parameters and to the removal of alignment ambiguous regions. Ten alignments were created with CLUSTAL W using different sets of alignment parameters (10 full alignments), and each alignment was examined by eye and alignment ambiguous regions were removed (creating 10 reduced alignments). These alignment ambiguous regions were analyzed as a third type of data set, culled alignments. Maximum parsimony, neighbor-joining, and parsimony bootstrap analyses were performed. The resulting phylogenies were compared to each other by the symmetric difference distance tree comparison metric (SymD). The correlation of the phylogenies with the alignment parameters was tested by comparing matrices from SymD with corresponding matrices of Manhattan distances representing the alignment parameters. Differences among individual parsimony trees from the full alignments were frequently correlated with the differences among alignment parameters (580/1000 tests), as were trees from the culled alignments (403/1000 tests). Differences among individual parsimony trees from the reduced alignments were less frequently correlated with the differences among alignment parameters (230/1000 tests). Differences among majority-rule consensus trees (50%) from the parsimony analysis of the full alignments were significantly correlated with the differences among alignment parameters, whereas consensus trees from the reduced and culled analyses were not correlated with the alignment parameters. These patterns of correlation confirm that choice of alignment parameters has the potential to bias the resultant phylogenies for the nematode SSU data set, and suggest that the removal of alignment ambiguous regions reduces this effect. Finally, we discuss the implications of conservative phylogenetic hypotheses for Nematoda produced by exploring alignment space and removing alignment ambiguous regions for SSU rDNA.

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.

Multigene phylogeny reveals eusociality evolved twice in vespid wasps

Eusocial wasps of the family Vespidae are thought to have derived their social behavior from a common ancestor that had a rudimentary caste-containing social system. In support of this behavioral scenario, the leading phylogenetic hypothesis of Vespidae places the eusocial wasps (subfamilies Stenogastrinae, Polistinae, and Vespinae) as a derived monophyletic clade, thus implying a single origin of eusocial behavior. This perspective has shaped the investigation and interpretation of vespid social evolution for more than two decades. Here we report a phylogeny of Vespidae based on data from four nuclear gene fragments (18S and 28S ribosomal DNA, abdominal-A and RNA polymerase II) and representatives from all six extant subfamilies. In contrast to the current phylogenetic perspective, our results indicate two independent origins of vespid eusociality, once in the clade Polistinae+Vespinae and once in the Stenogastrinae. The stenogastrines appear as an early diverging clade distantly related to the vespines and polistines and thus evolved their distinctive form of social behavior from a different ancestor than that of Polistinae+Vespinae. These results support earlier views based on life history and behavior and have important implications for interpreting transitional stages in vespid social evolution.

Phylogeny of parasitiform mites (Acari) based on rRNA

Acari (mites and ticks) form one the most diverse lineages of arthropods, but basal relationships in the group are still poorly understood. The current study addresses this issue for one of its two main lineages, the order Parasitiformes. Relationships are examined at the subordinal and infraordinal level using complete 18S and partial 28S nuclear rRNA sequence data. Most currently recognized lineages are recovered with good support, suggesting that nuclear rRNA, and specifically 18S rRNA, is very well suited for analyzing relationships at this level in this lineage. These results were found despite quite variable rates of sequence evolution, with rates "ratcheting up" from relatively low in most non-mite arachnid lineages, to intermediate in Pseudoscorpiones, the mite order Acariformes, and the parasitiform suborders Opilioacarida, Holothyrida, and Ixodida, to high in the parasitiform suborder Mesostigmata. The most species rich mesostigmatid infraorder, Dermanyssina, shows huge distances to the outgroups, but remarkably low within-group divergence in nuclear rRNA. This suggests the possibility of a relatively recent origin of this lineage.

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.

Systematics of the cyclostome subfamilies of braconid parasitic wasps (Hymenoptera: Ichneumonoidea): A simultaneous molecular and morphological Bayesian approach

Phylogenetic relationships among 95 genera collectively representing 17 of the 18 currently recognized cyclostome braconid wasp subfamilies were investigated based on DNA sequence fragments of the mitochondrial COI and the nuclear 28S rDNA genes, in addition to morphological data. The treatment of sequence length variation of the 28S partition was explored by either excluding ambiguously aligned regions and indel information (28SN) or recoding them (28SA) using the 'fragment-level' alignment method with a modified coding approach. Bayesian MCMC analyses were performed for the separate and combined data sets and a maximum parsimony analysis was also carried out for the simultaneous molecular and morphological data sets. There was a significant incongruence between the two genes and between 28S and morphology, but not for morphology and COI. Different analyses with the 28SA data matrix resulted in topologies that were generally similar to the ones from the 28SN matrix; however, the former topologies recovered a higher number of significantly supported clades and had a higher mean Bayesian posterior probability, thus supporting the inclusion of information from ambiguously aligned regions and indel events in phylogenetic analyses where possible. Based on the significantly supported clades obtained from the simultaneous molecular and morphological analyses, we propose that a total of 17 subfamilies should be recognized within the cyclostome group. The subfamilial placements of several problematic cyclostome genera were also established.

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

We report the entire sequence (2864 nts) and secondary structure of the nuclear small subunit ribosomal RNA (SSU rRNA) gene (18S) from the twisted-wing parasite Caenocholax fenyesi texensis Kathirithamby & Johnston (Strepsiptera: Myrmecolacidae). The majority of the base pairings in this structural model map on to the SSU rRNA secondary and tertiary helices that were previously predicted with comparative analysis. These regions of the core rRNA were unambiguously aligned across all Arthropoda. In contrast, many of the variable regions, as previously characterized in other insect taxa, had very large insertions in C. f. texensis. The helical base pairs in these regions were predicted with a comparative analysis of a multiple sequence alignment (that contains C. f. texensis and 174 published arthropod 18S rRNA sequences, including eleven strepsipterans) and thermodynamic-based algorithms. Analysis of our structural alignment revealed four unusual insertions in the core rRNA structure that are unique to animal 18S rRNA and in general agreement with previously proposed insertion sites for strepsipterans. One curious result is the presence of a large insertion within a hairpin loop of a highly conserved pseudoknot helix in variable region 4. Despite the extraordinary variability in sequence length and composition, this insertion contains the conserved sequences 5'-AUUGGCUUAAA-3' and 5'-GAC-3' that immediately flank a putative helix at the 5'- and 3'-ends, respectively. The longer sequence has the potential to form a nine base pair helix with a sequence in the variable region 2, consistent with a recent study proposing this tertiary interaction. Our analysis of a larger set of arthropod 18S rRNA sequences has revealed possible errors in some of the previously published strepsipteran 18S rRNA sequences. Thus we find no support for the previously recovered heterogeneity in the 18S molecules of strepsipterans. Our findings lend insight to the evolution of RNA structure and function and the impact large insertions pose on genome size. We also provide a novel alignment template that will improve the phylogenetic placement of the Strepsiptera among other insect taxa.

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.