Nucleotide sequence and presumed secondary structure of the 28S rRNA of pea aphid: implication for diversification of insect rRNA

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

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

Informative characteristics of 12 divergent domains in complete large subunit rDNA sequences from the harmful dinoflagellate genus, Alexandrium (Dinophyceae)

The genus Alexandrium includes organisms of interest, both for the study of dinoflagellate evolution and for their impacts as toxic algae affecting human health and fisheries. Only partial large subunit (LSU) rDNA sequences of Alexandrium and other dinoflagellates are available, although they contain much genetic information. Here, we report complete LSU rDNA sequences from 11 strains of Alexandrium, including Alexandrium affine, Alexandrium catenella, Alexandrium fundyense, Alexandrium minutum, and Alexandrium tamarense, and discuss their segmented domains and structure. Putative LSU rRNA coding regions were recorded to be around 3,400 bp long. Their GC content (about 43.7%) is among the lowest when compared with other organisms. Furthermore, no AT-rich regions were found in Alexandrium LSU rDNA, although a low GC content was recorded within the LSU rDNA. No intron-like sequences were found. The secondary structure of the LSU rDNA and parsimony analyses showed that most variation in LSU rDNA is found in the divergent (D) domains with the D2 region being the most informative. This high D domain variability can allow members of the diverse Alexandrium genus to be categorized at the species level. In addition, phylogenetic analysis of the alveolate group using the complete LSU sequences strongly supported previous findings that the dinoflagellates and apicomplexans form a clade.

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.

Heterologous rRNA gene expression: internal fragmentation of Sciara coprophila 28S rRNA within microinjected Xenopus laevis oocytes

Species-specific pre-rRNA processing variations may result in fragmented 18S, 5.8S and 28S rRNAs. Some insect 5.8S and 28S rRNAs are further cleaved, creating within a 'hidden break' or 'gap'. We investigated the specificity of the processing mechanism by microinjecting Sciara coprophila (fungus fly) rDNA into Xenopus laevis oocytes to examine insect rRNA maturation within a cell where endogenous rRNAs are not cleaved at homologous sites. Results confirm insect rDNA transcription and pre-28S rRNA fragmentation, demonstrating that fly-specific processing machinery is not required. Instead, oocytes may provide required accessory factors, suggesting that the insect gap processing mechanism is served by an evolutionarily conserved apparatus. Alternatively, these results may suggest that processing in some lineages is an autocatalytic property of the rRNA.

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.

X-linked heterochromatin distribution in the holocentric chromosomes of the green apple aphid Aphis pomi

Chromatin organization in the holocentric chromosomes of the green apple aphid Aphis pomi has been investigated at a cytological level after C-banding, NOR, Giemsa, fluorochrome staining and fluorescent in situ hybridization (FISH). C-banding technique showed that heterochromatic bands are exclusively located on X chromosomes. This data represents a peculiar feature that clearly contradicts the equilocal distribution of heterochromatin typical of monocentric chromosomes. Moreover, silver staining and FISH carried out with a 28S rDNA probe localized rDNA genes on one telomere of each X chromosome; CMA3 staining reveals that these silver positive telomeres are the only GC-rich regions among A. pomi heterochromatin, whereas all other C-positive bands are DAPI positive thus containing AT-rich DNA.

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

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

The structure of a single unit of ribosomal RNA gene (rDNA) including intergenic subrepeats in the Australian bulldog ant Myrmecia croslandi (Hymenoptera: Formicidae)

A complete single unit of a ribosomal RNA gene (rDNA) of M. croslandi was sequenced. The ends of the 18S, 5.8S and 28S rRNA genes were determined by using the sequences of D. melanogaster rDNAs as references. Each of the tandemly repeated rDNA units consists of coding and non-coding regions whose arrangement is the same as that of D. melanogaster rDNA. The intergenic spacer (IGS) contains, as in other species, a region with subrepeats, of which the sequences are different from those previously reported in other insect species. The length of IGSs was estimated to be 7-12 kb by genomic Southern hybridization, showing that an rDNA repeating unit of M. croslandi is 14-19 kb-long. The sequences of the coding regions are highly conserved, whereas IGS and ITS (internal transcribed spacer) sequences are not. We obtained clones with insertions of various sizes of R2 elements, the target sequence of which was found in the 28S rRNA coding region. A short segment in the IGS that follows the 3' end of the 28S rRNA gene was predicted to form a secondary structure with long stems.

The role of rDNA genes in X chromosome association in the aphid Acyrthosiphon pisum

Silver staining of mitotic metaphases of the aphid A. pisum reveals the presence of argentophilic bridges connecting the two X chromosomes. The presence of nucleolar material connecting sex chromosomes seems to be quite a common phenomenon in organisms belonging to very different phyla, and suggests a role of nucleolar proteins in chromosome association and disjunction. In somatic cells of A. pisum, bridges connecting X chromosomes are detectable not only after silver staining but also after CMA3 staining. This finding suggests that GC rich DNA is involved in this type of association. Molecular analysis of rDNA intergenic spacers shows several 247 bp repeats containing short sequences having a high level of homology with the chi sequence of Escherichia coli and with the consensus core region of human hypervariable minisatellites. Moreover, each 247 bp repeat presents a perfect copy of a promoter sequence for polymerase I. These aphid repeats show structural homologies with a 240 bp repeat, which is considered to be responsible for sex chromosome pairing in Drosophila, not only in view of their common presence within rDNA spacers but also for their length and structure. The presence of chi sequences in the IGS of A. pisum, by promoting unequal crossing-over between rDNA genes, could thus give rise to the nucleolar organizing region (NOR) heteromorphism described in different aphid species. Although X pairing at NORs is fundamental in aphid male determination, the presence of heteromorphism of rDNA genes does not inhibit male determination in the A. pisum clone utilized for our experiments.

Molecular phylogenetics at the Felsenstein zone: approaching the Strepsiptera problem using 5.8S and 28S rDNA sequences

Recent efforts to reconstruct the phylogenetic position of the insect order Strepsiptera have elicited a major controversy in molecular phylogenetics. We sequenced the 5.8S rDNA and major parts of the 28S rDNA 5' region of the strepsipteran species Stylops melittae. Their evolutionary dynamics were analyzed together with previously published insect rDNA sequences to identify tree estimation bias risks and to explore additional sources of phylogenetic information. Several major secondary structure changes were found as being autapomorphic for the Diptera, the Strepsiptera, or the Archaeognatha. Besides elevated substitution rates a significant AT bias was present in dipteran and strepsipteran 28S rDNA which, however, was restricted to stem sites in the Diptera while also affecting single-stranded sites in the Strepsiptera. When dipteran taxa were excluded from tree estimation all methods consistently supported the placement of Strepsiptera to within the Holometabola. When dipteran taxa were included maximum likelihood continued to favor a sister-group relationship of Strepsiptera with Mecoptera while remaining methods strongly supported a sister-group relationship with Diptera. Parametric bootstrap analysis revealed maximum likelihood as a consistent estimator if rate heterogeneity across sites was taken into account. Though the position of Strepsiptera within Holometabola remains elusive, we conclude that the evolution of dipteran and strepsipteran rDNA involved similar yet independent changes of substitution parameters.