Molecular coevolution among cryptically simple expansion segments of eukaryotic 26S/28S rRNAs

The set of "expansion segments" of any eukaryotic 26S/28S ribosomal RNA (rRNA) gene is responsible for the bulk of the difference in length between the prokaryotic 23S rRNA gene and the eukaryotic 26S/28S rRNA gene. The expansion segments are also responsible for interspecific fluctuations in length during eukaryotic evolution. They show a consistent bias in base composition in any species; for example, they are AT rich in Drosophila melanogaster and GC rich in vertebrate species. Dot-matrix comparisons of sets of expansion segments reveal high similarities between members of a set within any 28S rRNA gene of a species, in contrast to the little or spurious similarity that exists between sets of expansion segments from distantly related species. Similarities among members of a set of expansion segments within any 28S rRNA gene cannot be accounted for by their base-compositional bias alone. In contrast, no significant similarity exists within a set of "core" segments (regions between expansion segments) of any 28S rRNA gene, although core segments are conserved between species. The set of expansion segments of a 26S/28S gene is coevolving as a unit in each species, at the same time as the family of 28S rRNA genes, as a whole, is undergoing continual homogenization, making all sets of expansion segments from all ribosomal DNA (rDNA) arrays in a species similar in sequence. Analysis of DNA simplicity of 26S/28S rRNA genes shows a direct correlation between significantly high relative simplicity factors (RSFs) and sequence similarity among a set of expansion segments. A similar correlation exists between RSF values, overall rDNA lengths, and the lengths of individual expansion segments. Such correlations suggest that most length fluctuations reflect the gain and loss of simple sequence motifs by slippage-like mechanisms. We discuss the molecular coevolution of expansion segments, which takes place against a background of slippage-like and unequal crossing-over mechanisms of turnover that are responsible for the accumulation of interspecific differences in rDNA sequences.      

The sequence of 28S ribosomal RNA varies within and between human cell lines.

The primary structure of 28S ribosomal RNA constitutes a conserved core which is similar among most 23S-like rRNAs and expansion segments which occur at specific positions in the sequence. The expansion segments account for most of the size difference between prokaryotic (archaeal and eubacterial) and eukaryotic rRNAs and they exhibit a sequence variation which is unique among rRNAs. We have investigated the sequence variation of one of the expansion segments, V8, by sequencing a total of 111 V8 segments from 9 different human cell lines and tissues and have found 35 different variants. The variation occur mainly at two 'hot spots' which are separated by 170 nucleotides in the primary sequence but are neighbours in the secondary structure. The sequence of V8 segments varies both within and between human cell lines and tissues. The implications for the evolution of the eukaryotic 28S rRNA are discussed together with possible functions of the expansion segments. We also present a secondary structure model for the V8 segment based on comparative sequence analysis and chemical and enzymatic foot printing.     

Structural alignment of 18S and 28S rDNA sequences provides insights into phylogeny of Phytophaga (Coleoptera: Curculionoidea and Chrysomeloidea)

We performed a comparative study of partial rDNA sequences from a variety of Coleoptera taxa to construct an annotated alignment based on secondary structure information, which in turn, provides improved rRNA structure models useful for phylogenetic reconstruction. Subsequent phylogenetic analysis was performed to test monophyly and interfamilial relationships of the megadiverse plant feeding beetle group known as ‘Phytophaga’ (Curculionoidea and Chrysomeloidea), as well as to discover their closest relatives among the Cucujiformia. Parsimony and Bayesian analyses were performed based on the structural alignment of segments of 18S rRNA (variable regions V4‐V5, V7‐V9) and 28S rRNA (expansion segment D2). A total of 104 terminal taxa of Coleoptera were included: 96 species of Cucujiformia beetles, representing the families and most ‘subfamilies’ of weevils and chrysomeloids (Phytophaga), as well as several families of Cleroidea, Tenebrionoidea and Cucujoidea, and eight outgroups from three other polyphagan series: Scarabaeiformia, Elateriformia and Bostrichiformia. The results from the different methods of analysis agree — recovering the monophyly of the ‘Phytophaga’, including Curculionoidea and Chrysomeloidea as sister groups. The curculionoid and chrysomeloid phylogeny recovered from the aligned 18S and 28S rDNA segments, which is independent of morphological data, is in agreement with recent hypotheses or concepts based on morphological evidence, particularly with respect to familial relationships. Our results provide clues about the evolutionary origin of the phytophagan beetles within the megaclade Cucujiformia, suggesting that the sister group of ‘Curculionoidea + Chrysomeloidea’ is a clade of the ‘Cucujoidea’, represented in this study by species in Boganiidae, Erotylidae, Nitidulidae, Cucujidae and Silvanidae. The Coccinellidae and Endomychidae are not grouped with the latter, and the remaining terminal taxa are nested in Tenebrionoidea and Cleroidea. We propose that the combination of structurally aligned ribosomal RNA gene regions 18S (V4‐V5, V7‐V9) and 28S (D2) are useful in testing monophyly and resolving relationships among beetle superfamilies and families.     

Cotranscription and processing of 23S, 4.5S and 5S rRNA in chloroplasts from Zea mays

The termini of rRNA processing intermediates and of mature rRNA species encoded by the 3' terminal region of 23S rDNA, by 4.5S rDNA, by the 5' terminal region of 5S rDNA and by the 23S/4.5S/5S intergenic regions from Zea mays chloroplast DNA were determined by using total RNA isolated from maize chloroplasts and 32P-labelled rDNA restriction fragments of these regions for nuclease S1 and primer extension mapping. Several processing sites detectable by both 3' and 5' terminally labelled probes could be identified and correlated to the secondary structure for the 23S/4.5S intergenic region. The complete 4.5S/5S intergenic region can be reverse transcribed and a common processing site for maturation of 4.5S and 5S rRNA close to the 3' end of 4.5S rRNA was detected. It is therefore concluded that 23S, 4.5S and 5S rRNA are cotranscribed.     

Probing the structure of mouse Ehrlich ascites cell 5.8S, 18S and 28S ribosomal RNA in situ.

The secondary structure of mouse Ehrlich ascites 18S, 5.8S and 28S ribosomal RNA in situ was investigated by chemical modification using dimethyl sulphate and 1-cyclohexyl-3-(morpholinoethyl) carbodiimide metho-p-toluene sulphonate. These reagents specifically modify unpaired bases in the RNA. The reactive bases were localized by primer extension followed by gel electrophoresis. The three rRNA species were equally accessible for modification i.e. approximately 10% of the nucleotides were reactive. The experimental data support the theoretical secondary structure models proposed for 18S and 5.8/28S rRNA as almost all modified bases were located in putative single-strand regions of the rRNAs or in helical regions that could be expected to undergo dynamic breathing. However, deviations from the suggested models were found in both 18S and 28S rRNA. In 18S rRNA some putative helices in the 5'-domain were extensively modified by the single-strand specific reagents as was one of the suggested helices in domain III of 28S rRNA. Of the four eukaryote specific expansion segments present in mouse Ehrlich ascites cell 28S rRNA, segments I and III were only partly available for modification while segments II and IV showed average to high modification.     

Molecular characterization of gap region in 28S rRNA molecules in brine shrimp <Emphasis Type="Italic">Artemia parthenogenetica</Emphasis> and planarian <Emphasis Type="Italic">Dugesia japonica</Emphasis>

In most insects and some other protostomes, a small stretch of nucleotides can be removed from mature 28S rRNA molecules, which could create two 28S rRNA subunits (28Sα and 28Sβ). Thus, during electrophoresis, the rRNA profiles of these organisms may differ significantly from the standard benchmark since the two subunits co-migrate with the 18S rRNA. To understand the structure and mechanism of the atypical 28S rRNA molecule, partial fragments of 28Sα and 28Sβ in brine shrimp Artemia parthenogenetica and planarian Dugesia japonica were cloned using a modified technology based on terminal transferase. Alignment with the corresponding sequences of 28S rDNAs indicates that there are 41 nucleotides in A. parthenogenetica and 42 nucleotides in D. japonica absent from the mature rRNAs. The AU content of the gap sequences of D. japonica and A. parthenogenetica is high. Both the gaps may form stem-loop structure. In D. japonica a UAAU cleavage signal is identified in the loop, but it is absent in A. parthenogenetica. Thus, it is proposed that the gap processing of 28S rRNA was a late enzyme-dependent cleavage event in the rRNA maturational process based on the AU rich gap sequence and the formation of the stem-loop structure to expose the processing segment, while the deletion of the gap region would not affect the structure and function of the 28S rRNA molecule.     

Primary and secondary structure of rat 28 S ribosomal RNA.

The primary structure of rat (Rattus norvegicus) 28 S rRNA is determined inferred from the sequence of cloned rDNA fragments. The rat 28 S rRNA contains 4802 nucleotides and has an estimated relative molecular mass (Mr, Na-salt) of 1.66 X 10(6). Several regions of high sequence homology with S. cerevisiae 25 S rRNA are present. These regions can be folded in characteristic base-paired structures homologous to those proposed for Saccharomyces and E. coli. The excess of about 1400 nucleotides in the rat 28 S rRNA (as compared to Saccharomyces 25 S rRNA) is accounted for mainly by the presence of eight distinct G+C-rich segments of different length inserted within the regions of high sequence homology. The G+C content of the four insertions, containing more than 200 nucleotides, is in the range of 78 to 85 percent. All G+C-rich segments appear to form strongly base-paired structures. The two largest G+C-rich segments (about 760 and 560 nucleotides, respectively) are located near the 5'-end and in the middle of the 28 S rRNA molecule. These two segments can be folded into long base-paired structures, corresponding to the ones observed previously by electron microscopy of partly denatured 28 S rRNA molecules.     

Evolution of the secondary structures and compensatory mutations of the ribosomal RNAs of Drosophila melanogaster

This paper examines the effects of DNA sequence evolution on RNA secondary structures and compensatory mutations. Models of the secondary structures of Drosophila melanogaster 18S ribosomal RNA (rRNA) and of the complex between 2S, 5.8S, and 28S rRNAs have been drawn on the basis of comparative and energetic criteria. The overall AU richness of the D. melanogaster rRNAs allows the resolution of some ambiguities in the structures of both large rRNAs. Comparison of the sequence of expansion segment V2 in D. melanogaster 18S rRNA with the same region in three other Drosophila species and the tsetse fly (Glossina morsitans morsitans) allows us to distinguish between two models for the secondary structure of this region. The secondary structures of the expansion segments of D. melanogaster 28S rRNA conform to a general pattern for all eukaryotes, despite having highly divergent sequences between D. melanogaster and vertebrates. The 70 novel compensatory mutations identified in the 28S rRNA show a strong (70%) bias toward A-U base pairs, suggesting that a process of biased mutation and/or biased fixation of A and T point mutations or AT-rich slippage-generated motifs has occurred during the evolution of D. melanogaster rDNA. This process has not occurred throughout the D. melanogaster genome. The processes by which compensatory pairs of mutations are generated and spread are discussed, and a model is suggested by which a second mutation is more likely to occur in a unit with a first mutation as such a unit begins to spread through the family and concomitantly through the population. Alternatively, mechanisms of proofreading in stem-loop structures at the DNA level, or between RNA and DNA, might be involved. The apparent tolerance of noncompensatory mutations in some stems which are otherwise strongly supported by comparative criteria within D. melanogaster 28S rRNA must be borne in mind when compensatory mutations are used as a criterion in secondary-structure modeling. Noncompensatory mutation may extend to the production of unstable structures where a stem is stabilized by RNA- protein or additional RNA-RNA interactions in the mature ribosome. Of motifs suggested to be involved in rRNA processing, one (CGAAAG) is strongly overrepresented in the 28S rRNA sequence. The data are discussed both in the context of the forces involved with the evolution of multigene families and in the context of molecular coevolution in the rDNA family in particular.      

A 12 S precursor to 5.8 S rRNA associated with rat liver nucleolar 28 S rRNA

The pre-rRNA and rRNA components of rat and mouse liver nucleolar RNA were analysed. It was shown that upon denaturation, part of the 32 S pre-rRNA is converted into 28 S rRNA and 12 S RNA. The 12 S RNA from mouse (Mr, 0.36 X 10(6)) is larger than the one from rat (Mr, 0.32 X 10(6). The 12 S RNA chain is intact and resists denaturation treatment. The non-covalent binding of this RNA with nucleolar 28 S rRNA is stronger than that of 5.8 S rRNA with 28 S rRNA. Hybridization with a rat internal-transcribed spacer rDNA fragment identifies 12 S RNA as corresponding to the 5'-end non-conserved segment of 32 S pre-rRNA, including 5.8 S rRNA. The significance of the formation of a 12 S precursor to 5.8 S rRNA in the biogenesis of ribosomes in mammalian cells is discussed.     

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 10 kb Drosophila virilis 28S rDNA intervening sequence is flanked by a direct repeat of 14 base pairs of coding sequence

Most repeat units of rDNA in Drosophila virilis are interrupted in the 28S rRNA coding region by an intervening sequence about 10 kb in length; uninterrupted repeats have a length of about 11 kb. We have sequenced the coding/intervening sequence junctions and flanking regions in two independent clones of interrupted rDNA, and the corresponding 28S rRNA coding region in a clone of uninterrupted rDNA. The intervening sequence is terminated at both ends by a direct repeat of a fourteen nucleotide sequence that is present once in the corresponding region of an intact gene. This is a phenomenon associated with transposable elements in other eukaryotes and in prokaryotes, and the Drosophila rDNA intervening sequence is discussed in this context. We have compared more than 200 nucleotides of the D. virilis 28S rRNA gene with sequences of homologous regions of rDNA in Tetrahymena pigmentosa (Wild and Sommer, 1980) and Xenopus laevis (Gourse and Gerbi, 1980): There is 93% sequence homology among the diverse species, so that the rDNA region in question (about two-thirds of the way into the 28S rRNA coding sequence) has been very highly conserved in eukaryote evolution. The intervening sequence in T. pigmentosa is at a site 79 nucleotides upstream from the insertion site of the Drosophila intervening sequence.     

Sequence requirements for maturation of the 5' terminus of human 18 S rRNA in vitro.

Creation of the mature 5' terminus of human 18 S rRNA in vitro occurs via a two-step processing reaction. In the first step, an endonucleolytic activity found in HeLa cell nucleolar extract cleaves an rRNA precursor spanning the external transcribed spacer-18 S boundary at a position 3 bases upstream from the mature 18 S terminus leaving 2',3'-cyclic phosphate, 5' hydroxyl termini. In the second step, a nucleolytic activity(s) found in HeLa cell cytoplasmic extract removes the 3 extra bases and creates the authentic 5'-phosphorylated terminus of 18 S rRNA. Here we have examined the sequence requirements for the trimming reaction. The trimming activity(s), in addition to requiring a 5' hydroxyl terminus, prefers the naturally occurring adenosine as the 5'-terminal base. By a combination of deletion, site-directed mutagenesis, and chemical modification interference approaches we have also identified a region of 18 S rRNA spanning bases +6 to +25 (with respect to the mature 5' end) which comprises a critical recognition sequence for the trimming activity(s).     

Gastropod evolutionary rates and phylogenetic relationships assessed using partial 28S rDNA and histone H3 sequences

Sequence data for two segments of 28S and Histone H3 from 36 gastropod taxa, a chiton, two bivalves and Nautilus are used to test recently published morphology‐based phylogenetic hypotheses of gastropod relationships. Statistical results suggest that the accuracy of the available hypotheses could be improved. The data support the monophyly of the Patellogastropoda (true limpets), Euthyneura and the ‘higher’ vetigastropods and the polyphyly of the ‘Cocculiniformia’. The division of the gastropods into two major clades (Eogastropoda and Orthogastropoda) as has been proposed on morphological grounds is not supported, and neither the Caenogastropoda nor Heterobranchia is well supported. Within the Euthyneura, opisthobranchs are paraphyletic with respect to the pulmonates. The hot vent taxon, Depressigyra, groups with the lower vetigastropod Pleurotomaria in some analyses. Much of the variability in the 28S rDNA segments lies in discrete areas of the sequence. Forone of the segments, corresponding to positions 691–942 of the mosquito Aedes albopictus 28S sequence, the variable regions represent known expansion regions (D4 and D5). For the other segment, corresponding to positions 2259–2538 of the A. albopictus sequence, the variable area, which is found in the patellogastropods, vetigastropods and Nautilus, represents an unreported expansion region. The data show marked variability in the rate of evolution in both segments of the 28S rDNA, whether or not the expansion regions are included. The variability is largely clade specific. Rates are high in the patellogastropods, vetigastropods, the lower heterobranch ‘Heterostropha’ (Cornirostra and Philippea), Depressigyra and the deep sea cocculinid limpet Coccopigya and substantially lower in other taxa. Rate variation in the histone H3 data is less extreme. The correlation between evolutionary rates in the two 28S rDNA segments is very high, andis also significant for the the pairing of each of the 28S rDNA segments with H3. The rate variability may be due to differential selection but no causative factor has been identified. The histone H3 data have high codon usage bias. For all amino acids encoded by multiple codons, at least some triplets occur at a frequency of less than a quarter of their expected usage. For all three‐, four‐and sixfold degenerate amino acids, the most abundant triplet occurs at least twice as frequently as expected. Despite the usage bias, there is a large amount of apparent homoplasy in synonymous alternatives at both the first and third codon positions.     

28 S ribosomal RNA in vertebrates. Locations of large-scale features revealed by electron microscopy in relation to other features of the sequences.

The 28 S rRNA from several vertebrate species, when examined by electron microscopy, is seen to contain regions of extensive secondary structure, as first reported for HeLa-cell 28 S rRNA by Wellauer & Dawid [(1973) Proc. Natl. Acad. Sci. U.S.A. 70, 2827-2831]. Here we correlate the locations of these regions, determined from the electron-microscopic data, with the primary structure of 28 S rRNA from human, mouse and Xenopus laevis determined by sequence analysis of rDNA. The secondary-structure features observed by electron microscopy correspond closely to phylogenetically variable G + C-rich regions that largely comprise the eukaryotic expansion segments in these three species. In most if not all cases the features can be identified with long G + C-rich helices deduced from sequence data. Correlations are given between the locations of the secondary-structure features and several 'landmark' restriction sites in 28 S rDNA. By correlating the locations of the rRNA methyl groups reported elsewhere [Maden (1988) J. Mol. Biol. 201, 289-314] with the present findings it is concluded that the rRNA secondary-structure features revealed by electron microscopy are largely or wholly unmethylated.