Comparative analysis of secondary structure of insect mitochondrial small subunit ribosomal RNA using maximum weighted matching

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

Comparative and phylogenomic analyses of mitochondrial genomes in Coccinellidae (Coleoptera: Coccinelloidea)

The Coccinellidae are one of the most familiar beetle families, the ladybirds. Despite the great ecological and economic significance, the phylogenetic relationships of Coccinellidae remain poorly understood. One of the reasons is that the sequenced mitogenomes available for this family are very limited. We sequenced complete or nearly complete mitogenomes from seven species of the tribe Coccinellini with next-generation sequencing. All species have the same gene content and gene order as the putatively ancestral insect mitogenome. A large intergenic spacer region (> 890 bp) was found located between trnI and trnQ. The potential for using secondary structures of the large and small ribosomal subunits for phylogenetic reconstruction was predicted. The phylogenetic relationships were explored through comparative analyses across more than 30 coccinellid species. We performed phylogenetic analyses with both concatenation methods (Maximum Likelihood and Bayesian Inference) and multispecies coalescent method (ASTRAL). Phylogenetic results strongly supported the monophyly of Coccinellidae. Within Coccinellidae, the Epilachnini and the Coccinellini including Halyziini were monophyletic, while the Scymnini and Coccidulini were non-monophyletic.

Characteristic and variability of five complete aphid mitochondrial genomes: Aphis fabae mordvilkoi, Aphis craccivora, Myzus persicae, Therioaphis tenera and Appendiseta robiniae (Hemiptera; Sternorrhyncha; Aphididae)

The complete mitochondrial genomes of aphids Aphis fabae mordvilkoi, A. craccivora, Myzus persicae from Aphidinae as well as Therioaphis tenera and Appendiseta robiniae from Calaphidinae were sequenced and compared with the genomes of other aphid species. A. fabae mordvilkoi, Th. tenera and A. robiniae mitogenomes were sequenced and analyzed for the first time. The annotation of A. craccivora and M. persicae were corrected compared to what was previously published. According to our data there is no translocation of tRNA-Tyr gene in A. craccivora mitogenome and this aphid species has an ancestral type of mitochondrial gene order. A + T content in all 5 mitogenomes was higher than 80%. A + T content in the Th. tenera CR was 59.5% which is untypically low. CRs of all 5 studied mitogenomes had 2 conserved motifs at their ends and extended G + C rich region. A. craccivora, M. persicae and Th. tenera had large tandem repeats inside the CRs. Detailed molecular analysis of all 5 aphid mitochondrial genomes showed the importance of a deep understanding of the molecular organization of all the functional regions of the mitochondrial DNA, which helps to avoid mistakes during genome annotation.

Comparative Mitogenomic Analysis of Species Representing Six Subfamilies in the Family Tenebrionidae

To better understand the architecture and evolution of the mitochondrial genome (mitogenome), mitogenomes of ten specimens representing six subfamilies in Tenebrionidae were selected, and comparative analysis of these mitogenomes was carried out in this study. Ten mitogenomes in this family share a similar gene composition, gene order, nucleotide composition, and codon usage. In addition, our results show that nucleotide bias was strongly influenced by the preference of codon usage for A/T rich codons which significantly correlated with the G + C content of protein coding genes (PCGs). Evolutionary rate analyses reveal that all PCGs have been subjected to a purifying selection, whereas 13 PCGs displayed different evolution rates, among which ATPase subunit 8 (ATP8) showed the highest evolutionary rate. We inferred the secondary structure for all RNA genes of Tenebrio molitor (Te2) and used this as the basis for comparison with the same genes from other Tenebrionidae mitogenomes. Some conserved helices (stems) and loops of RNA structures were found in different domains of ribosomal RNAs (rRNAs) and the cloverleaf structure of transfer RNAs (tRNAs). With regard to the AT-rich region, we analyzed tandem repeat sequences located in this region and identified some essential elements including T stretches, the consensus motif at the flanking regions of T stretch, and the secondary structure formed by the motif at the 3′ end of T stretch in major strand, which are highly conserved in these species. Furthermore, phylogenetic analyses using mitogenomic data strongly support the relationships among six subfamilies: ((Tenebrionidae incertae sedis + (Diaperinae + Tenebrioninae)) + (Pimeliinae + Lagriinae)), which is consistent with phylogenetic results based on morphological traits.

Comparative mitogenomic analysis of the superfamily Pentatomoidea (Insecta: Hemiptera: Heteroptera) and phylogenetic implications

BACKGROUND

Insect mitochondrial genomes (mitogenomes) are the most extensively used genetic marker for evolutionary and population genetics studies of insects. The Pentatomoidea superfamily is economically important and the largest superfamily within Pentatomomorpha with over 7,000 species. To better understand the diversity and evolution of pentatomoid species, we sequenced and annotated the mitogenomes of Eurydema gebleri and Rubiconia intermedia, and present the first comparative analysis of the 11 pentatomoid mitogenomes that have been sequenced to date.

RESULTS

We obtained the complete mitogenome of Eurydema gebleri (16,005 bp) and a nearly complete mitogenome of Rubiconia intermedia (14,967 bp). Our results show that gene content, gene arrangement, base composition, codon usage, and mitochondrial transcription termination factor sequences are highly conserved in pentatomoid species, especially for species in the same family. Evolutionary rate analyses of protein-coding genes reveal that the highest and lowest rates are found in atp8 and cox1 and distinctive evolutionary patterns are significantly correlated with the G + C content of genes. We inferred the secondary structures for two rRNA genes for eleven pentatomoid species, and identify some conserved motifs of RNA structures in Pentatomidea. All tRNA genes in pentatomoid mitogenomes have a canonical cloverleaf secondary structure, except for two tRNAs (trnS1 and trnV) which appear to lack the dihydrouridine arm. Regions that are A + T-rich have several distinct characteristics (e.g. size variation and abundant tandem repeats), and have potential as species or population level molecular markers. Phylogenetic analyses based on mitogenomic data strongly support the monophyly of Pentatomoidea, and the estimated phylogenetic relationships are: (Urostylididae + (Plataspidae + (Pentatomidae + (Cydnidae + (Dinidoridae + Tessaratomidae))))).

CONCLUSIONS

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the superfamily Pentatomoidea. Mitogenomes can be effectively used to resolve phylogenetic relationships of pentatomomorphan insects at various taxonomic levels. Sequencing more mitogenomes at various taxonomic levels, particularly from closely related species, will improve the annotation accuracy of mitochondrial genes, as well as greatly enhance our understanding of mitogenomic evolution and phylogenetic relationships in pentatomoids.

The Complete Mitochondrial Genome of Corizus tetraspilus (Hemiptera: Rhopalidae) and Phylogenetic Analysis of Pentatomomorpha

Insect mitochondrial genome (mitogenome) are the most extensively used genetic information for molecular evolution, phylogenetics and population genetics. Pentatomomorpha (>14,000 species) is the second largest infraorder of Heteroptera and of great economic importance. To better understand the diversity and phylogeny within Pentatomomorpha, we sequenced and annotated the complete mitogenome of Corizus tetraspilus (Hemiptera: Rhopalidae), an important pest of alfalfa in China. We analyzed the main features of the C. tetraspilus mitogenome, and provided a comparative analysis with four other Coreoidea species. Our results reveal that gene content, gene arrangement, nucleotide composition, codon usage, rRNA structures and sequences of mitochondrial transcription termination factor are conserved in Coreoidea. Comparative analysis shows that different protein-coding genes have been subject to different evolutionary rates correlated with the G+C content. All the transfer RNA genes found in Coreoidea have the typical clover leaf secondary structure, except for trnS1 (AGN) which lacks the dihydrouridine (DHU) arm and possesses a unusual anticodon stem (9 bp vs. the normal 5 bp). The control regions (CRs) among Coreoidea are highly variable in size, of which the CR of C. tetraspilus is the smallest (440 bp), making the C. tetraspilus mitogenome the smallest (14,989 bp) within all completely sequenced Coreoidea mitogenomes. No conserved motifs are found in the CRs of Coreoidea. In addition, the A+T content (60.68%) of the CR of C. tetraspilus is much lower than that of the entire mitogenome (74.88%), and is lowest among Coreoidea. Phylogenetic analyses based on mitogenomic data support the monophyly of each superfamily within Pentatomomorpha, and recognize a phylogenetic relationship of (Aradoidea + (Pentatomoidea + (Lygaeoidea + (Pyrrhocoroidea + Coreoidea)))).

Genetic code evolution reveals the neutral emergence of mutational robustness, and information as an evolutionary constraint

The standard genetic code (SGC) is central to molecular biology and its origin and evolution is a fundamental problem in evolutionary biology, the elucidation of which promises to reveal much about the origins of life. In addition, we propose that study of its origin can also reveal some fundamental and generalizable insights into mechanisms of molecular evolution, utilizing concepts from complexity theory. The first is that beneficial traits may arise by non-adaptive processes, via a process of "neutral emergence". The structure of the SGC is optimized for the property of error minimization, which reduces the deleterious impact of point mutations. Via simulation, it can be shown that genetic codes with error minimization superior to the SGC can emerge in a neutral fashion simply by a process of genetic code expansion via tRNA and aminoacyl-tRNA synthetase duplication, whereby similar amino acids are added to codons related to that of the parent amino acid. This process of neutral emergence has implications beyond that of the genetic code, as it suggests that not all beneficial traits have arisen by the direct action of natural selection; we term these "pseudaptations", and discuss a range of potential examples. Secondly, consideration of genetic code deviations (codon reassignments) reveals that these are mostly associated with a reduction in proteome size. This code malleability implies the existence of a proteomic constraint on the genetic code, proportional to the size of the proteome (P), and that its reduction in size leads to an "unfreezing" of the codon - amino acid mapping that defines the genetic code, consistent with Crick's Frozen Accident theory. The concept of a proteomic constraint may be extended to propose a general informational constraint on genetic fidelity, which may be used to explain variously, differences in mutation rates in genomes with differing proteome sizes, differences in DNA repair capacity and genome GC content between organisms, a selective pressure in the evolution of sexual reproduction, and differences in translational fidelity. Lastly, the utility of the concept of an informational constraint to other diverse fields of research is explored.

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

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

The complete mitochondrial genome of the stalk-eyed bug Chauliops fallax Scott, and the monophyly of Malcidae (Hemiptera: Heteroptera)

Chauliops fallax Scott, 1874 (Hemiptera: Heteroptera: Malcidae: Chauliopinae) is one of the most destructive insect pests of soybean and rice fields in Asia. Here we sequenced the complete mitochondrial genome of this pest. This genome is 15,739 bp long, with an A+T content of 73.7%, containing 37 typical animal mitochondrial genes and a control region. All genes were arranged in the same order as most of other Heteroptera. A remarkable strand bias was found for all nine protein coding genes (PCGs) encoded by the majority strand were positive AT-skew and negative GC-skew, whereas the reverse were found in the remaining four PCGs encoded by the minority strand and two rRNA genes. The models of secondary structures for the two rRNA genes of sequenced true bugs and Lygaeoidea were predicted. 16S rRNA consisted of six domains (domain III is absent as in other known arthropod mitochondrial genomes) and 45 helices, while three domains and 27 helices for 12S rRNA. The control region consists of five subregions: a microsatellite-like region, a tandem repeats region and other three motifs. The unusual intergenic spacer between tRNA-H and ND4 only found in the species of Lygaeoidea, not in other heteropteran species, may be the synapomorphy of this superfamily. Phylogenetic analyses were carried out based on all the 13 PCGs showed that Chauliopinae was the sister group of Malcinae and the monophyly of Lygaeoidea.

The complete mitochondrial genome of the damsel bug Alloeorhynchus bakeri (Hemiptera: Nabidae)

The complete sequence of the mitochondrial DNA (mtDNA) of the damsel bug, Alloeorhynchus bakeri, has been completed and annotated in this study. It represents the first sequenced mitochondrial genome of heteropteran family Nabidae. The circular genome is 15, 851 bp in length with an A+T content of 73.5%, contains the typical 37 genes that are arranged in the same order as that of the putative ancestor of hexapods. Nucleotide composition and codon usage are similar to other known heteropteran mitochondrial genomes. All protein-coding genes (PCGs) use standard initiation codons (methionine and isoleucine), except COI, which started with TTG. Canonical TAA and TAG termination codons are found in eight protein-coding genes, the remaining five (COI, COII, COIII, ND5, ND1) have incomplete termination codons (T or TA). PCGs of two strands present opposite CG skew which is also reflected by the nucleotide composition and codon usage. All tRNAs have the typical clover-leaf structure, except the dihydrouridine (DHU) arm of tRNA(Ser (AGN))which forms a simple loop as known in many other metazoa. Secondary structure models of the ribosomal RNA genes of A. bakeri are presented, similar to those proposed for other insect orders. There are six domains and 45 helices and three domains and 27 helices in the secondary structures of rrnL and rrnS, respectively. The major non-coding region (also called control region) between the small ribosomal subunit and the tRNA(Ile )gene includes two special regions. The first region includes four 133 bp tandem repeat units plus a partial copy of the repeat (28 bp of the beginning), and the second region at the end of control region contains 4 potential stem-loop structures. Finally, PCGs sequences were used to perform a phylogenetic study. Both maximum likelihood and Bayesian inference analyses highly support Nabidae as the sister group to Anthocoridae and Miridae.

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

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

Mitochondrial rRNA secondary structures and genome arrangements distinguish chelicerates: comparisons with a harvestman (Arachnida: Opiliones: Phalangium opilio)

Arachnids are a highly diverse group of arthropods, and many of the mitochondrial genomes that have been sequenced from arachnids possess unusual features in their inferred gene structures and genome organization. The first complete sequence of a mitochondrial genome from the arachnid order Opiliones (harvestmen) is presented here. Secondary structures of the two mitochondrial ribosomal subunits of Phalangium opilio are inferred and compared to mitochondrial rRNA structures of a hexapod and a chelicerate. The large subunit rRNA of P. opilio is found to have more helices conserved than in other arthropods, while the small subunit rRNA shows a complexity similar to that of other arthropods. These comparisons suggest that a reduction in rRNA complexity occurred in Pancrustacea after the divergence of Pancrustacea and Chelicerata from a common ancestor. The gene arrangement of the mitochondrial genome of P. opilio is compared with the gene order of taxa from all seven other orders of arachnids for which representative mitochondrial genomes have been sequenced. Taxa from five of these seven orders possess gene arrangements identical to that of Limulus polyphemus, and P. opilio is found to have a similar arrangement. However, in P. opilio, some genes near the putative control region are rearranged, with the suite of genes encoding tRNA(Gln), the control region, and tRNA(Ile) located downstream of the two ribosomal RNA genes, and upstream of where they are typically located in chelicerates. The genome encodes only 21 of the typical 22 mitochondrial tRNA genes and lacks the gene for tRNA(Leu(CUN)). The protein-coding genes in the mitochondrial genome of P. opilio show a significantly decreased use of codons recognized by tRNA(Leu(CUN)), likely due to selection to utilize the more specific tRNA(Leu(UUR)) anticodon. The gene arrangement and lack of a tRNA(Leu(CUN)) gene in P. opilio is most parsimoniously explained by the occurrence of at least two translocation events, one of which probably destroyed the function of the tRNA(Leu(CUN)) gene. Phylogenetic relationships among the major orders of arachnids are inferred, using all 13 mt protein-coding genes, and gene rearrangements are mapped onto the phylogeny. The phylogenetic analyses are unable to resolve the placement of P. opilio but are generally consistent with an early divergence of members of the Dromopoda (harvestmen, scorpions, and solifuges) from the Micruran arachnids (spiders, whip spiders, vinegaroons, ricinuleids, and mites). However, unlike some morphologically based phylogenetic analyses, the existence of a clade of Dromopoda is not supported. While data on genome arrangement and gene loss do not provide further information to help resolve relationships among the arachnid orders, they distinguish some groups of arachnids, distinguish chelicerates from other arthropods, and further clarify the ancestral gene order of this diverse group of arthropods.

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

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

Defining the limits of taxonomic conservatism in host-plant use for phytophagous insects: molecular systematics and evolution of host-plant associations in the seed-beetle genus Bruchus Linnaeus (Coleoptera: Chrysomelidae: Bruchinae)

In this study, we have investigated the limits of taxonomic conservatism in host-plant use in the seed-beetle genus Bruchus. To reconstruct the insect phylogeny, parsimony and multiple partitioned Bayesian inference analyses were conducted on a combined data set of four genes. Permutation tests and both global and local maximum-likelihood optimizations of host preferences at distinct taxonomic levels revealed that host-fidelity is still discernible beyond the host-plant tribe level, suggesting the existence of more important than previously thought evolutionary constraints, which are further discussed in details. Our tree topologies are also mostly consistent with extant taxonomic groups. Through the analysis of this empirical data set we also provide meaningful insights on two methodological issues. First, Bayesian inference analyses suggest that partitioning by using codon positions greatly increase the accuracy of phylogenetical reconstructions. Regarding reconstruction of ancestral character states through maximum likelihood, the present study also highlights the usefulness of local optimizations. The issue of over-parameterization is also addressed, as the optimizations with the most parameter-rich models have returned the most counterintuitive results.

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.

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

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

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

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

The mitochondrial genome of the Japanese freshwater crab, Geothelphusa dehaani (Crustacea: Brachyura): Evidence for its evolution via gene duplication

The translocation of tRNA genes has occurred frequently among multiple independent arthropod lineages. However, the mechanisms of such gene rearrangement are still unclear. In this study we sequenced the nearly complete mitochondrial genome of the Japanese freshwater crab, Geothelphusa dehaani (Decapoda, Brachyura, Potamidae). Its size is relatively large (app. 24 kb) for higher animal mtDNA and is due to the presence of many intergenic noncoding sequences (IGNs), one of which is a large repeat region (>7.0 kb). Excluding the repeat region and the putative control region (514 bp), the total size of IGNs (1953 bp) is still the largest reported in arthropods. In comparison to the proposed ancestral arrangement for the insect-crustacean clade, G. dehaani displays changes in the positional order for three tRNA genes (trnQ, trnL(uur) and trnH). At the putative former and current locations for trnQ and trnH, there are IGNs with high similarities to corresponding sequences in the ancestral arrangement. This suggests the occurrence of partial duplication, followed by a loss of function for one copy of these two duplicated genes in G. dehaani. Furthermore, three trnLs (one trnL(uur) and two trnL(cun)s) are found in G. dehaani; their sequences are highly similar to each other (70.8-82.5%) and similar to the trnL(cun) of other arthropods. This suggests that duplication of trnL(cun) occurred at least twice, and the current trnL(uur) of G. dehaani arose by a subsequent anticodon alternation. These observations provide direct evidence for the "duplication-random loss" model for mtDNA gene rearrangement. Comparisons with other available Brachyuran data(http://www.ncbi.nlm.nih.gov/genomes/ORGANELLES/6657.html) suggest that two of the three tRNA translocations are unique in G. dehaani. For G. dehaani, where duplications have occurred but deletion of extra sequences is incomplete, the IGNs found at the ancestral or trans-locations are likely to be footprints of relatively recent mitochondrial genome evolution.

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.

Secondary structure, high variability and conserved motifs for domain III of 12S rRNA in the Arthropleona (Hexapoda; Collembola)

The robustness of different hypothetical structural models for ribosomal RNAs can be established by comparing sequences from widely different taxonomic groups. Secondary structure models of the two mitochondrial ribosomal subunits have been proposed for a large number of pterygote insects, but little information is available for the more ancient entognath hexapods. We have investigated the most frequently sequenced of the four domains of the mitochondrial small subunit rRNA in twenty-two collembolan species in the taxon Arthropleona. We provide secondary structure models of this third domain for representative species, and a consensus based on all species studied. This consensus is in partial agreement with previous models. High levels of variation, in terms of length and of primary sequence, are present in the peripheral portions of the domain. Some of the structural elements differ more among genera of springtails than among insect orders. Conversely, a limited number of structural differences compared with other taxa have been detected in the core regions.

Widespread polyphyly among Alopiinae snail genera: when phylogeny mirrors biogeography more closely than morphology

Consider a group of species that is evenly divided by an easily identifiable complex morphological character. Most biologists would assume that this character should provide better phylogenetic information than, say, the spatial distribution of these species over a fairly continuous 500-km radius area. Paradoxically, this is not the case among terrestrial snail genera in the clausiliid subfamily Alopiinae. Phylogenetic analysis using the nuclear markers ITS1/ITS2 and mitochondrial markers COI/12S reveals widespread homoplasy in the clausilial apparatus (a complex aperture-closing mechanism), and concomitant extensive polyphyly among Carinigera, Isabellaria, and Sericata. In contrast, phylogenetic relationships as revealed by molecular data are closely congruent with biogeography at a relatively small scale. A combination of extremely low vagility and extremely high morphological convergence has conspired to produce this unexpected result. Implications as to the function of the clausilial apparatus are discussed.

ILM: a web server for predicting RNA secondary structures with pseudoknots

The ILM web server provides a web interface to two algorithms, iterated loop matching and maximum weighted matching, for efficiently predicting RNA secondary structures with pseudoknots. The algorithms can utilize either thermodynamic or comparative information or both, and thus can work on both aligned and individual sequences. Predicted secondary structures are presented in several formats compatible with a variety of existing visualization tools. The service can be accessed at http://cic.cs.wustl.edu/RNA/.

Secondary structure and molecular evolution of the mitochondrial small subunit ribosomal RNA in Agaricales (Euagarics clade, Homobasidiomycota)

The complete sequences and secondary structures of the mitochondrial small subunit (SSU) ribosomal RNAs of both mostly cultivated mushrooms Agaricus bisporus (1930 nt) and Lentinula edodes (2164 nt) were achieved. These secondary structures and that of Schizophyllum commune (1872 nt) were compared to that previously established for Agrocybe aegerita. The four structures are near the model established for Archae, Bacteria, plastids, and mitochondria; particularly the helices 23 and 37, described as specific to bacteria, are present. Within the four Agaricales (Homobasidiomycota), the SSU-rRNA "core" is conserved in size (966 to 1009 nt) with the exception of an unusual extension of 40 nt in the H17 helix of S. commune. The four core sequences possess 76% of conserved positions and a cluster of C in their 3' end, which could constitute a signal involved in the RNA maturation process. Among the nine putative variable domains, three (V3, V5, V7) do not show significant length variations and possess similar percentages of conserved positions (69%) than the core. The other six variable domains show important length variations, due to independent large size inserted/deleted sequences, and higher rates of nucleotide substitutions than the core (only 31% of conserved positions between the four species). Interestingly, the inserted/deleted sequences are located in few preferential sites (hot spots for insertion/deletion) where they seem to arise or disappear haphazardly during evolution. These sites are located on the surface of the tertiary structure of the 30S ribosomal subunit, at the beginning of hairpin loops; the insertions lead to a lengthening of existing hairpins or to branching loops bearing up to five additional helices.

Linking coevolutionary history to ecological process: doves and lice

Many host-specific parasites are restricted to a limited range of host species by ecological barriers that impede dispersal and successful establishment. In some cases, microevolutionary differentiation is apparent on top of host specificity, as evidenced by significant parasite population genetic structure among host populations. Ecological barriers responsible for specificity and genetic structure can, in principle, reinforce macroevolutionary processes that generate congruent host-parasite phylogenies. However, few studies have explored both the micro- and macroevolutionary ramifications of close association in a single host-parasite system. Here we compare the macroevolutionary histories of two genera of feather lice (Phthiraptera: Ischnocera) that both parasitize New World pigeons and doves (Aves: Columbiformes). Earlier work has shown that dove body lice (genus Physconelloides) are more host specific and have greater population genetic structure than dove wing lice (Columbicola). We reconstructed phylogenies for representatives of the two genera of lice and their hosts, using nuclear and mitochondrial DNA sequences. The phylogenies were well resolved and generally well supported. We compared the phylogenies of body lice and wing lice to the host phylogeny using reconciliation analyses. We found that dove body lice show strong evidence of cospeciation whereas dove wing lice do not. Although the ecology of body and wing lice is very similar, differences in their dispersal ability may underlie these joint differences in host specificity, population genetic structure, and coevolutionary history.

Louse (Insecta: Phthiraptera) mitochondrial 12S rRNA secondary structure is highly variable

Lice are ectoparasitic insects hosted by birds and mammals. Mitochondrial 12S rRNA sequences obtained from lice show considerable length variation and are very difficult to align. We show that the louse 12S rRNA domain III secondary structure displays considerable variation compared to other insects, in both the shape and number of stems and loops. Phylogenetic trees constructed from tree edit distances between louse 12S rRNA structures do not closely resemble trees constructed from sequence data, suggesting that at least some of this structural variation has arisen independently in different louse lineages. Taken together with previous work on mitochondrial gene order and elevated rates of substitution in louse mitochondrial sequences, the structural variation in louse 12S rRNA confirms the highly distinctive nature of molecular evolution in these insects.