Comparisons of large subunit rRNAs reveal some eukaryote-specific elements of secondary structure

The story of rRNA expansion segments: Finding functionality amidst diversity

Expansion segments (ESs) are multinucleotide insertions present across phyla at specific conserved positions in eukaryotic rRNAs. ESs are generally absent in bacterial rRNAs with some exceptions, while the archaeal rRNAs have microexpansions at regions that coincide with those of eukaryotic ESs. Although there is an increasing prominence of ribosomes, especially the ribosomal proteins, in fine-tuning gene expression through translation regulation, the role of rRNA ESs is relatively underexplored. While rRNAs have been established as the major catalytic hub in ribosome function, the presence of ESs widens their scope as a species-specific regulatory hub of protein synthesis. In this comprehensive review, we have elaborately discussed the current understanding of the functional aspects of rRNA ESs of cytoplasmic eukaryotic ribosomes and discuss their past, present, and future. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems Translation > Ribosome Structure/Function Translation > Regulation.

Supersized Ribosomal RNA Expansion Segments in Asgard Archaea

The ribosome’s common core, comprised of ribosomal RNA (rRNA) and universal ribosomal proteins, connects all life back to a common ancestor and serves as a window to relationships among organisms. The rRNA of the common core is similar to rRNA of extant bacteria. In eukaryotes, the rRNA of the common core is decorated by expansion segments (ESs) that vastly increase its size. Supersized ESs have not been observed previously in Archaea, and the origin of eukaryotic ESs remains enigmatic. We discovered that the large ribosomal subunit (LSU) rRNA of two Asgard phyla, Lokiarchaeota and Heimdallarchaeota, considered to be the closest modern archaeal cell lineages to Eukarya, bridge the gap in size between prokaryotic and eukaryotic LSU rRNAs. The elongated LSU rRNAs in Lokiarchaeota and Heimdallarchaeota stem from two supersized ESs, called ES9 and ES39. We applied chemical footprinting experiments to study the structure of Lokiarchaeota ES39. Furthermore, we used covariation and sequence analysis to study the evolution of Asgard ES39s and ES9s. By defining the common eukaryotic ES39 signature fold, we found that Asgard ES39s have more and longer helices than eukaryotic ES39s. Although Asgard ES39s have sequences and structures distinct from eukaryotic ES39s, we found overall conservation of a three-way junction across the Asgard species that matches eukaryotic ES39 topology, a result consistent with the accretion model of ribosomal evolution.

Root of the Tree: The Significance, Evolution, and Origins of the Ribosome

The ribosome is an ancient molecular fossil that provides a telescope to the origins of life. Made from RNA and protein, the ribosome translates mRNA to coded protein in all living systems. Universality, economy, centrality and antiquity are ingrained in translation. The translation machinery dominates the set of genes that are shared as orthologues across the tree of life. The lineage of the translation system defines the universal tree of life. The function of a ribosome is to build ribosomes; to accomplish this task, ribosomes make ribosomal proteins, polymerases, enzymes, and signaling proteins. Every coded protein ever produced by life on Earth has passed through the exit tunnel, which is the birth canal of biology. During the root phase of the tree of life, before the last common ancestor of life (LUCA), exit tunnel evolution is dominant and unremitting. Protein folding coevolved with evolution of the exit tunnel. The ribosome shows that protein folding initiated with intrinsic disorder, supported through a short, primitive exit tunnel. Folding progressed to thermodynamically stable β-structures and then to kinetically trapped α-structures. The latter were enabled by a long, mature exit tunnel that partially offset the general thermodynamic tendency of all polypeptides to form β-sheets. RNA chaperoned the evolution of protein folding from the very beginning. The universal common core of the ribosome, with a mass of nearly 2 million Daltons, was finalized by LUCA. The ribosome entered stasis after LUCA and remained in that state for billions of years. Bacterial ribosomes never left stasis. Archaeal ribosomes have remained near stasis, except for the superphylum Asgard, which has accreted rRNA post LUCA. Eukaryotic ribosomes in some lineages appear to be logarithmically accreting rRNA over the last billion years. Ribosomal expansion in Asgard and Eukarya has been incremental and iterative, without substantial remodeling of pre-existing basal structures. The ribosome preserves information on its history.

Translation: The Universal Structural Core of Life

The Universal Gene Set of Life (UGSL) is common to genomes of all extant organisms. The UGSL is small, consisting of <100 genes, and is dominated by genes encoding the translation system. Here we extend the search for biological universality to three dimensions. We characterize and quantitate the universality of structure of macromolecules that are common to all of life. We determine that around 90% of prokaryotic ribosomal RNA (rRNA) forms a common core, which is the structural and functional foundation of rRNAs of all cytoplasmic ribosomes. We have established a database, which we call the Sparse and Efficient Representation of the Extant Biology (the SEREB database). This database contains complete and cross-validated rRNA sequences of species chosen, as far as possible, to sparsely and efficiently sample all known phyla. Atomic-resolution structures of ribosomes provide data for structural comparison and validation of sequence-based models. We developed a similarity statistic called pairing adjusted sequence entropy, which characterizes paired nucleotides by their adherence to covariation and unpaired nucleotides by conventional conservation of identity. For canonically paired nucleotides the unit of structure is the nucleotide pair. For unpaired nucleotides, the unit of structure is the nucleotide. By quantitatively defining the common core of rRNA, we systematize the conservation and divergence of the translational system across the tree of life, and can begin to understand the unique evolutionary pressures that cause its universality. We explore the relationship between ribosomal size and diversity, geological time, and organismal complexity.

Expression of distinct maternal and somatic 5.8S, 18S, and 28S rRNA types during zebrafish development

There is mounting evidence that the ribosome is not a static translation machinery, but a cell-specific, adaptive system. Ribosomal variations have mostly been studied at the protein level, even though the essential transcriptional functions are primarily performed by rRNAs. At the RNA level, oocyte-specific 5S rRNAs are long known for Xenopus. Recently, we described for zebrafish a similar system in which the sole maternal-type 5S rRNA present in eggs is replaced completely during embryonic development by a somatic-type. Here, we report the discovery of an analogous system for the 45S rDNA elements: 5.8S, 18S, and 28S. The maternal-type 5.8S, 18S, and 28S rRNA sequences differ substantially from those of the somatic-type, plus the maternal-type rRNAs are also replaced by the somatic-type rRNAs during embryogenesis. We discuss the structural and functional implications of the observed sequence differences with respect to the translational functions of the 5.8S, 18S, and 28S rRNA elements. Finally, in silico evidence suggests that expansion segments (ES) in 18S rRNA, previously implicated in ribosome-mRNA interaction, may have a preference for interacting with specific mRNA genes. Taken together, our findings indicate that two distinct types of ribosomes exist in zebrafish during development, each likely conducting the translation machinery in a unique way.

Eukaryote-specific rRNA expansion segments function in ribosome biogenesis

The secondary structure of ribosomal RNA (rRNA) is largely conserved across all kingdoms of life. However, eukaryotes have evolved extra blocks of rRNA sequences, relative to those of prokaryotes, called expansion segments (ES). A thorough characterization of the potential roles of ES remains to be done, possibly because of limitations in the availability of robust systems to study rRNA mutants. We sought to systematically investigate the potential functions, if any, of the ES in 25S rRNA of Saccharomyces cerevisiae by deletion mutagenesis. We deleted 14 of the 16 different eukaryote-specific ES in yeast 25S rRNA individually and assayed their phenotypes. Our results show that all but two of the ES tested are necessary for optimal growth and are required for production of 25S rRNA, suggesting that ES play roles in ribosome biogenesis. Further, we classified expansion segments into groups that participate in early nucleolar, middle, and late nucleoplasmic steps of ribosome biogenesis, by assaying their pre-rRNA processing phenotypes. This study is the first of its kind to systematically identify the functions of eukaryote-specific expansion segments by showing that they play roles in specific steps of ribosome biogenesis. The catalog of phenotypes we identified, combined with previous investigations of the roles ribosomal proteins in large subunit biogenesis, leads us to infer that assembling ribosomes are composed of distinct RNA and protein structural neighborhood clusters that participate in specific steps of ribosome biogenesis.

History of the ribosome and the origin of translation

We present a molecular-level model for the origin and evolution of the translation system, using a 3D comparative method. In this model, the ribosome evolved by accretion, recursively adding expansion segments, iteratively growing, subsuming, and freezing the rRNA. Functions of expansion segments in the ancestral ribosome are assigned by correspondence with their functions in the extant ribosome. The model explains the evolution of the large ribosomal subunit, the small ribosomal subunit, tRNA, and mRNA. Prokaryotic ribosomes evolved in six phases, sequentially acquiring capabilities for RNA folding, catalysis, subunit association, correlated evolution, decoding, energy-driven translocation, and surface proteinization. Two additional phases exclusive to eukaryotes led to tentacle-like rRNA expansions. In this model, ribosomal proteinization was a driving force for the broad adoption of proteins in other biological processes. The exit tunnel was clearly a central theme of all phases of ribosomal evolution and was continuously extended and rigidified. In the primitive noncoding ribosome, proto-mRNA and the small ribosomal subunit acted as cofactors, positioning the activated ends of tRNAs within the peptidyl transferase center. This association linked the evolution of the large and small ribosomal subunits, proto-mRNA, and tRNA.

The ribosome challenge to the RNA world

An RNA World that predated the modern world of polypeptide and polynucleotide is one of the most widely accepted models in origin of life research. In this model, the translation system shepherded the RNA World into the extant biology of DNA, RNA, and protein. Here, we examine the RNA World Hypothesis in the context of increasingly detailed information available about the origins, evolution, functions, and mechanisms of the translation system. We conclude that the translation system presents critical challenges to RNA World Hypotheses. Firstly, a timeline of the RNA World is problematic when the ribosome is incorporated. The mechanism of peptidyl transfer of the ribosome appears distinct from evolved enzymes, signaling origins in a chemical rather than biological milieu. Secondly, we have no evidence that the basic biochemical toolset of life is subject to substantive change by Darwinian evolution, as required for the transition from the RNA world to extant biology. Thirdly, we do not see specific evidence for biological takeover of ribozyme function by protein enzymes. Finally, we can find no basis for preservation of the ribosome as ribozyme or the universality of translation, if it were the case that other information transducing ribozymes, such as ribozyme polymerases, were replaced by protein analogs and erased from the phylogenetic record. We suggest that an updated model of the RNA World should address the current state of knowledge of the translation system.

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.

Evolution of helix formation in the ribosomal Internal Transcribed Spacer 2 (ITS2) and its significance for RNA secondary structures

Helices are the most common elements of RNA secondary structure. Despite intensive investigations of various types of RNAs, the evolutionary history of the formation of new helices (novel helical structures) remains largely elusive. Here, by studying the nuclear ribosomal Internal Transcribed Spacer 2 (ITS2), a fast-evolving part of the eukaryotic nuclear ribosomal operon, we identify two possible types of helix formation: one type is "dichotomous helix formation"--transition from one large helix to two smaller helices by invagination of the apical part of a helix, which significantly changes the shape of the original secondary structure but does not increase its complexity (i.e., the total length of the RNA). An alternative type is "lateral helix formation"--origin of an extra helical region by the extension of a bulge loop or a spacer in a multi-helix loop of the original helix, which does not disrupt the pre-existing structure but increases RNA size. Moreover, we present examples from the RNA sequence literature indicating that both types of helix formation may have implications for RNA evolution beyond ITS2.

The GC content of LSU rRNA evolves across topological and functional regions of the ribosome in all three domains of life

Large-subunit rRNA is the ribozyme that catalyzes protein synthesis by translation, and many of its features vary along a deep-to-superficial gradient. By measuring the G+C proportions in this rRNA in all three domains of life (60 bacteria, 379 eukaryote, and 23 archaean sequences), we tested whether the proportion of GC nucleotides varies along this in-out gradient. The rRNA regions used were several zones identified by Bokov and Steinberg (2009) as being arranged from deep to superficial within the LSU. To the Bokov-Steinberg zones, we added the most superficial zone of all, the divergent domains (expansion segments), which are greatly enlarged in eukaryotes. Regression lines constructed from the hundreds of species of organisms revealed the expected in-out gradient, showing that species with high %GC (or high %AT) in their rRNA distribute more of these abundant nucleotides into the peripheral zones. This could be explained by the evolutionary rates of replacement of all nucleotides (A, C, G, T), because these latter rates are fastest at the periphery and slowest near the conserved core. As an overall explanation, we propose that when extrinsic factors (whole-genome nucleotide composition, or environmental temperature) demand the percentage of GC in the rRNA of a species be high or low, then the deep-lying zones are buffered against GC variation because they are the slowest to evolve. The deep, conserved zones are also the most involved in translation, hinting that stabilizing selection there prevents a high GC variability that would diminish LSU rRNA's core functions. We found only a few domain-specific trends in rRNA-GC distribution, which relate to many Archaea living at high temperatures or to the highly complex genes and adaptations of Eukaryota. Use of rRNA sequences in molecular phylogenetic studies, for reconstructing the relationships of organisms across the tree of life, requires accurate models of how rRNA evolves. The demonstration that GC distributes in regular patterns across rRNA regions can improve these tree-reconstruction models in the future and should yield phylogenies of greater accuracy.

Identification, molecular characterization, and evolution of group I introns at the expansion segment D11 of 28S rDNA in Rhizoctonia species

The nuclear ribosomal DNA of Rhizoctonia species is polymorphic in terms of the nucleotide composition and length. Insertions of 349-410 nucleotides in length with characteristics of group I introns were detected at a single insertion point at the expansion segment D11 of 28S rDNA in 12 out of 64 isolates. Eleven corresponded to Rhizoctonia solani (teleomorph: Thanatephorous) and one (AG-Q) to Rhizoctonia spp. (teleomorph: Ceratobasidium). Sequence data showed that all but AG-Q contained conserved DNA catalytic core regions (P, Q, R, and S) essential for selfsplicing. The predicted secondary structure revealed that base-paired helices corresponded to subgroup IC1. Isolates from same anastomosis group and even subgroups within R. solani were variable with regard to possession of introns. Phylogenetic analyses indicated that introns were vertically transmitted. Unfortunately, sequence data from the conserved region from all 64 isolates were not useful for delimiting species. Analyses with IC1 introns at same insertion point, of both Ascomycota and Basidiomycota indicated the possibility of horizontal transfer at this site. The present study uncovered new questions on evolutionary pattern of change of these introns within Rhizoctonia species.

Molecular phylogeny of pearl oysters and their relatives (Mollusca, Bivalvia, Pterioidea)

BACKGROUND

The superfamily Pterioidea is a morphologically and ecologically diverse lineage of epifaunal marine bivalves distributed throughout the tropical and subtropical continental shelf regions. This group includes commercially important pearl culture species and model organisms used for medical studies of biomineralization. Recent morphological treatment of selected pterioideans and molecular phylogenetic analyses of higher-level relationships in Bivalvia have challenged the traditional view that pterioidean families are monophyletic. This issue is examined here in light of molecular data sets composed of DNA sequences for nuclear and mitochondrial loci, and a published character data set of anatomical and shell morphological characters.

RESULTS

The present study is the first comprehensive species-level analysis of the Pterioidea to produce a well-resolved, robust phylogenetic hypothesis for nearly all extant taxa. The data were analyzed for potential biases due to taxon and character sampling, and idiosyncracies of different molecular evolutionary processes. The congruence and contribution of different partitions were quantified, and the sensitivity of clade stability to alignment parameters was explored.

CONCLUSIONS

Four primary conclusions were reached: (1) the results strongly supported the monophyly of the Pterioidea; (2) none of the previously defined families (except for the monotypic Pulvinitidae) were monophyletic; (3) the arrangement of the genera was novel and unanticipated, however strongly supported and robust to changes in alignment parameters; and (4) optimizing key morphological characters onto topologies derived from the analysis of molecular data revealed many instances of homoplasy and uncovered synapomorphies for major nodes. Additionally, a complete species-level sampling of the genus Pinctada provided further insights into the on-going controversy regarding the taxonomic identity of major pearl culture species.

Comprehensive molecular phylogenetic analysis and evolution of the genus Phyllactinia (Ascomycota: Erysiphales) and its allied genera

Phyllactinia is a unique genus within the Erysiphales (Ascomycota) having a partly endo-parasitic nature of the mycelium within the host plant tissues. We constructed phylogenetic trees for the genus Phyllactinia and its allied genera based on a total of 120 nucleotide sequences of the 28S rDNA and ITS regions to discuss their phylogenetic relationships with special references to host plants, biogeography, evolutionary dating, and taxonomy. The analysis of the Erysiphales confirmed the monophyly of the endo-parasitic genera, i.e. Leveillula, Phyllactinia, and Pleochaeta. Phyllactinia specimens used in this study were divided into six distinctive groups and three subgroups. Interestingly, Leveillula, an obligately endo-parasitic genus of the Erysiphales, grouped together with Phyllactinia, although this was not significantly supported by the Kishino-Hasegawa and Shimodaira-Hasegawa tests. This suggests that the evolution within this group of fungi occurred from partial endo-parasitism to obligate endo-parasitism. The host range of Phyllactinia is mostly confined to woody plants, especially deciduous trees. Betulaceae, Fagaceae, Ulmaceae, Moraceae, and Rosaceae may have close connections to the divergence of the groups and subgroups of Phyllactinia concerned. Most of these plant families are known as major members of the boreotropical flora of the Tertiary, which suggests an early Tertiary origin of this genus. A comparison of the phylogenies of hosts and parasites revealed that host range expansion at higher taxonomic levels (higher than family level) is independent of the phylogeny of plants. Conversely, host range expansions in lower taxonomic levels (infrafamilial or infrageneric) tend to occur within a single family or genus. An estimation of the evolutionary timing using a molecular clock approach suggested that Phyllactinia split from Pleochaeta about 60 M years ago (Ma) in the early Tertiary and divergence of the six major clades of Phyllactinia occurred between 5 and 40 Ma during the Oligocene and Miocene. Divergence within the major clades and within Leveillula occurred maybe from more than 5 Ma onwards during the Pliocene and Quaternary. This is the first comprehensive phylogenetic study of Phyllactinia and other endo-parasitic genera of the Erysiphales.

Tetrahymena ORC contains a ribosomal RNA fragment that participates in rDNA origin recognition

The Tetrahymena thermophila ribosomal DNA (rDNA) replicon contains dispersed cis-acting replication determinants, including reiterated type I elements that associate with sequence-specific, single-stranded binding factors, TIF1 through TIF4. Here, we show that TIF4, previously implicated in cell cycle-controlled DNA replication and rDNA gene amplification, is the T. thermophila origin recognition complex (TtORC). We further demonstrate that TtORC contains an integral RNA subunit that participates in rDNA origin recognition. Remarkably, this RNA, designated 26T, spans the terminal 282 nts of 26S ribosomal RNA. 26T RNA exhibits extensive complementarity to the type I element T-rich strand and binds the rDNA origin in vivo. Mutations that disrupt predicted interactions between 26T RNA and its complementary rDNA target change the in vitro binding specificity of ORC and diminish in vivo rDNA origin utilization. These findings reveal a role for ribosomal RNA in chromosome biology and define a new mechanism for targeting ORC to replication initiation sites.

Identification of conserved secondary structures and expansion segments in enod40 RNAs reveals new enod40 homologues in plants

enod40 is a plant gene that participates in the regulation of symbiotic interaction between leguminous plants and bacteria or fungi. Furthermore, it has been suggested to play a general role in non-symbiotic plant development. Although enod40 seems to have multiple functions, being present in many land plants, the molecular mechanisms of its activity are unclear; they may be determined though, by short peptides and/or RNA structures encoded in the enod40 genes. We utilized conserved RNA structures in enod40 sequences to search nucleotide sequence databases and identified a number of new enod40 homologues in plant species that belong to known, but also, to yet unknown enod40-containing plant families. RNA secondary structure predictions and comparative sequence analysis of enod40 RNAs allowed us to determine the most conserved structural features, present in all known enod40 genes. Remarkably, the topology and evolution of one of the conserved structural domains are similar to those of the expansion segments found in structural RNAs such as rRNAs, RNase P and SRP RNAs. Surprisingly, the enod40 RNA structural elements are much more stronger conserved than the encoded peptides. This finding suggests that some general functions of enod40 gene could be determined by the encoded RNA structure, whereas short peptides may be responsible for more diverse functions found only in certain plant families.

Application of the secondary structure model of rRNA for phylogeny: D2-D3 expansion segments of the LSU gene of plant-parasitic nematodes from the family Hoplolaimidae Filipjev, 1934

Knowledge of rRNA structure is increasingly important to assist phylogenetic analysis through reconstructing optimal alignment, utilizing molecule features as an additional source of data and refining appropriate models of evolution of the molecule. We describe a procedure of optimization for alignment and a new coding method for nucleotide sequence data using secondary structure models of the D2 and D3 expansion fragments of the LSU-rRNA gene reconstructed for fifteen nematode species of the agriculturally important and diverse family Hoplolaimidae, order Tylenchida. Using secondary structure information we converted the original sequence data into twenty-eight symbol codes and submitted the transformed data to maximum parsimony analysis. We also applied the original sequence data set for Bayesian inference. This used the doublet model with sixteen states of nucleotide doublets for the stem region and the standard model of DNA substitution with four nucleotide states for loops and bulges. By this approach, we demonstrate that using structural information for phylogenetic analyses led to trees with lower resolved relationships between clades and likely eliminated some artefactual support for misinterpreted relationships, such as paraphyly of Helicotylenchus or Rotylenchus. This study as well as future phylogenetic analyses is herein supported by the development of an on-line database, NEMrRNA, for rRNA molecules in a structural format for nematodes. We also have developed a new computer program, RNAstat, for calculation of nucleotide statistics designed and proposed for phylogenetic studies.

Analysis of the secondary structure of expansion segment 39 in ribosomes from fungi, plants and mammals

The structure of expansion segment 39, ES39, in eukaryotic 23 S-like ribosomal RNA was analysed using a combination of chemical and enzymic reagents. Ribosomes were isolated from yeast, wheat, mouse, rat and rabbit, five organisms representing three different eukaryotic kingdoms. The isolated ribosomes were treated with structure-sensitive chemical and enzymic reagents and the modification patterns analysed by primer extension and gel electrophoresis on an ABI 377 automated DNA sequencer. The expansion segment was relatively accessible to modification by both enzymic and chemical probes, suggesting that ES39 was exposed on the surface of the ribosomes. The collected modification data were used in secondary structure modelling of the expansion segment. Despite considerable variation in both sequence and length between organisms from different kingdoms, the structure analysis of the expansion segment gave rise to structural fingerprints that allowed identification of homologous structures in ES39 from fungi, plants and mammals. The homologous structures formed an initial helix and an invariant hairpin connected to the initial helix via a long single-stranded loop. The remaining part of the ES39 sequences accounted for most of the length variation seen between the analysed species. This part could form additional, albeit less similar, hairpins. A comparison of ES39 sequences from other fungi, plants and mammals showed that identical structures could be formed in these organisms.

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.

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.

First complete large subunit ribosomal RNA sequence and secondary structure for a parasitic nematode: phylogenetic and diagnostic implications

The complete sequence and secondary structure of the large subunit of nuclear ribosomal RNA(LSUrRNA) were determined for the parasitic nematode Labiostrongylus bipapillosus (order Strongylida). Its LSU rRNA sequence was shorter (by 18 bp) than that of the free-living nematode, Caenorhabditis elegans (order Rhabditida), the only other species within the Nematoda for which a complete LSU rRNA sequence has been determined. Interspecific differences in sequence were greater in the 12 D domains compared with the core segments, with the secondary structure being maintained by partial or complete compensatory base pair changes. The magnitude of interspecific sequence difference in each D domain (except for D6 and D12) was similar, suggesting that several domains contain informative genetic markers for phylogenetic studies of the phylum Nematoda at different taxonomic levels. The LSU rRNA may also provide species-specific markers for the identification of some bursate nematodes of veterinary and medical importance.

A DNA-based method for identification of krill species and its application to analysing the diet of marine vertebrate predators

Accurate identification of species that are consumed by vertebrate predators is necessary for understanding marine food webs. Morphological methods for identifying prey components after consumption often fail to make accurate identifications of invertebrates because prey morphology becomes damaged during capture, ingestion and digestion. Another disadvantage of morphological methods for prey identification is that they often involve sampling procedures that are disruptive for the predator, such as stomach flushing or lethal collection. We have developed a DNA-based method for identifying species of krill (Crustacea: Malacostraca), an enormously abundant group of invertebrates that are directly consumed by many groups of marine vertebrates. The DNA-based approach allows identification of krill species present in samples of vertebrate stomach contents, vomit, and, more importantly, faeces. Utilizing samples of faeces from vertebrate predators minimizes the impact of dietary studies on the subject animals. We demonstrate our method first on samples of Adelie penguin (Pygoscelis adeliae) stomach contents, where DNA-based species identification can be confirmed by prey morphology. We then apply the method to faeces of Adelie penguins and to faeces of the endangered pygmy blue whale (Balaenoptera musculus brevicauda). In each of these cases, krill species consumed by the predators could be identified from their DNA present in faeces or stomach contents.

Ribin, a protein encoded by a message complementary to rRNA, modulates ribosomal transcription and cell proliferation

The control of rRNA transcription, tightly coupled to the cell cycle and growth state of the cell, is a key process for understanding the mechanisms that drive cell proliferation. Here we describe a novel protein, ribin, found in rodents, that binds to the rRNA promoter and stimulates its activity. The protein also interacts with the basal rRNA transcription factor UBF. The open reading frame encoding ribin is 96% complementary to a central region of the large rRNA. This demonstrates that ribosomal DNA-related sequences in higher eukaryotes can be expressed as protein-coding messages. Ribin contains two predicted nuclear localization sequence elements, and green fluorescent protein-ribin fusion proteins localize in the nucleus. Cell lines overexpressing ribin exhibit enhanced rRNA transcription and faster growth. Furthermore, these cells significantly overcome the suppression of rRNA synthesis caused by serum deprivation. On the other hand, the endogenous ribin level correlates positively with the amount of serum in the medium. The data show that ribin is a limiting stimulatory factor for rRNA synthesis in vivo and suggest its involvement in the pathway that adapts ribosomal transcription and cell proliferation to physiological changes.

Evolution of transcript structure and base composition of rDNA expansion segment D3 in ticks

Four to thirty-two copies of the rDNA 28S gene expansion segment D3 and flanking H14 stem were sequenced in six species of ticks (Ixodes: Ixodidae: Acari). Sequence match among species varied from 66% to 97%. Sequence length averaged 130 bases in I. persulcatus across eight Eurasian sites and averaged 186 bases in five other species across 19 Eurasian and North American sites. The difference in length represents one or more deletions totalling about 60 bases that correspond to stems S3 or S4 of the folded transcript. The typical transcript conformation was observed as one possible low energy structure in the five species of longer D3. The structure entails a basal loop with four stem/loop structures, S1-S4 (moving 5' to 3') atop stem H14. A secondary structure lacking S4 but possessing all other putative standard features of the D3 transcript is possible with the shorter I. persulcatus sequences. Interspecific sequence differences occur at higher frequency in loops and bulges vs. complementary pairing regions of stems. Insertion/deletion events (indels) and base substitutions accounted equally for sequence differences. Indels are flanked by similar sequences, suggesting that they occur by slippage during replication. The D3 of Ixodes species is composed of a degenerate set of subrepeats. Thus, unequal exchange among subrepeats may have caused the reduction in length of the I. persulcatus D3. Compensatory base substitution and compensatory insertion/deletion events are indicated by the failure of mutations to affect secondary structure. Transversions accounted for 64% of sequence differences and were biased toward the gain of G and U and the loss of A and C. This bias could re-establish intramolecular base pairing when disrupted by insertions or deletions that shift one side of a stem relative to the other. The distribution of sequence differences, biased substitution, and conservation of transcript conformation in D3 suggest selective constraint.

Proposed secondary structure of eukaryote specific expansion segment 15 in 28S rRNA from mice, rats, and rabbits

The expansion segments in eukaryotic ribosomal RNAs are additional RNA sequences not found in the RNA core common to both prokaryotes and eukaryotes. These regions show large species-dependent variations in sequence and size. This makes it difficult to create secondary structure models for the expansion segments exclusively based on phylogenetic sequence comparison. Here we have used a combination of experimental data and computational methods to generate secondary structure models for expansion segment 15 in 28S rRNA in mice, rats, and rabbits. The experimental data were collected using the structure sensitive reagents DMS, CMCT, kethoxal, micrococcal nuclease, RNase T(1), RNase CL3, RNase V(1), and lead(II) acetate. ES15 was folded with the computer program RNAStructure 3.5 using modification data and phylogenetic similarities between different ES15 sequences. This program uses energy minimization to find the most stable secondary structure of an RNA sequence. The presented secondary structure models include several common structural motifs, but they also have characteristics unique to each organism. Overall, the secondary structure models showed indications of an energetically stable but dynamic structure, easily accessible from the solution by the modification reagents, suggesting that the expansion segment is located on the ribosomal surface.

Evolution of mitochondrial and ribosomal gene sequences in anophelinae (Diptera: Culicidae): implications for phylogeny reconstruction

In this study, two mitochondrial genes, cyt b and ND5, and the D2 expansion segment of the 28S nuclear ribosomal gene were used to reconstruct a phylogeny of the mosquito subfamily Anophelinae. The ingroup consisted of all three genera of Anophelinae and five of six subgenera of Anopheles. Six genera of Culicinae were used as the outgroup. Extreme conservation at the protein level coupled with rapid saturation of synonymous positions probably accounted for the lack of meaningful phylogenetic signal in the cyt b gene. In contrast, abundant variation at all codon positions of the ND5 gene allowed recovery of the basal and most of the recent relationships. Phylogenetic analysis of D2 produced results consistent with those of ND5. Combined analysis indicated well-supported monophyletic Anophelinae (with Chagasia basal), Anopheles + Bironella, and subgeneric clades within the genus Anopheles. Moreover, subgenera Nyssorhynchus and Kerteszia were supported as a monophyletic lineage. The Kishino-Hasegawa test could not reject the monophyly of Anopheles, whereas the recently proposed hypothesis of close affinity of Bironella to the subgenus Anopheles was rejected by the analyses of ND5 and combined data sets. The lack of resolution of Bironella and Anopheles clades, or basal relationships among subgeneric clades within Anopheles, suggests their rapid diversification. Recovery of relationships consistent with morphology and previous molecular studies provides evidence of substantial phylogenetic signal in D2 and ND5 genes at levels of divergence from closely related species to subfamily in mosquitoes.

Evolutionary assembly of the conifer fauna: distinguishing ancient from recent associations in bark beetles

Several shifts from ancestral conifer feeding to angiosperm feeding have been implicated in the unparalleled diversification of beetle species. The single largest angiosperm-feeding beetle clade occurs in the weevils, and comprises the family Curculionidae and relatives. Most authorities confidently place the bark beetles (Scolytidae) within this radiation of angiosperm feeders. However, some clues indicate that the association between conifers and some scolytids, particularly in the tribe Tomicini, is a very ancient one. For instance, several fragments of Gondwanaland (South America, New Caledonia, Australia and New Guinea) harbour endemic Tomicini specialized on members of the formerly widespread and abundant conifer family Araucariaceae. As a first step towards resolving this seeming paradox, we present a phylogenetic analysis of the beetle family Scolytidae with particularly intensive sampling of conifer-feeding Tomicini and allies. We sequenced and analysed elongation factor 1alpha and nuclear rDNAs 18S and 28S for 45 taxa, using members of the weevil family Cossoninae as an out-group. Our results indicate that conifer feeding is the ancestral host association of scolytids, and that the most basal lineages of scolytids feed on Aramucaria. If scolytids are indeed nested within a great angiosperm-feeding clade, as many authorities have held, then a reversion to conifer feeding in ancestral scolytids appears to have occurred in the Mesozoic, when Araucaria still formed a major component of the woody flora.

Caspase-dependent cleavage of nucleic acids

Autoimmune diseases are frequently characterized by the presence of autoantibodies directed against nucleic acid-protein complexes present in the nucleus of the cell. The mechanisms by which these autoantigenic molecules escape immunological tolerance are largely unknown, although a number of recent observations suggest that modified self-proteins generated during apoptosis may play an important role in the development of autoimmunity. It has been hypothesized that the recognition of these modified self-proteins by the immune system may promote autoantibody production. While apoptosis is specifically characterized by posttranslational modification of proteins, recent findings also show that nucleic acids are modified. This review summarizes the specific cleavages of some of these key nucleic acids, i.e. chromosomal DNA, ribosomal RNA and small structural RNAs (U1 snRNA, Y RNA), in apoptotic cells.

Common secondary structures for the second internal transcribed spacer pre-rRNA of two subfamilies of trichostrongylid nematodes

Sequences of the second internal transcribed spacer ribosomal DNA for the parasitic trichostrongylid nematodes Trichostrongylus probolurus, Trichostrongylus rugatus and Camelostrongylus mentulatus were compared with previously published sequences for five other species within the genus Trichostrongylus. The secondary structures of the second internal transcribed spacer pre-rRNA for these nematodes were predicted using an energy minimisation method. The results indicate that a common secondary structure of the second internal transcribed spacer of these nematodes is maintained despite distinct differences in primary sequence between species. Sequence differences among Trichostrongylus species ranged from 1.3 to 7.6%, but each species differed by 22-26% in sequence when compared with C. mentulatus which belongs to a different subfamily.

Microsporidian Encephalitozoon cuniculi, a unicellular eukaryote with an unusual chromosomal dispersion of ribosomal genes and a LSU rRNA reduced to the universal core

Microsporidia are eukaryotic parasites lacking mitochondria, the ribosomes of which present prokaryote-like features. In order to better understand the structural evolution of rRNA molecules in microsporidia, the 5S and rDNA genes were investigated in Encephalitozoon cuniculi . The genes are not in close proximity. Non-tandemly arranged rDNA units are on every one of the 11 chromosomes. Such a dispersion is also shown in two other Encephalitozoon species. Sequencing of the 5S rRNA coding region reveals a 120 nt long RNA which folds according to the eukaryotic consensus structural shape. In contrast, the LSU rRNA molecule is greatly reduced in length (2487 nt). This dramatic shortening is essentially due to truncation of divergent domains, most of them being removed. Most variable stems of the conserved core are also deleted, reducing the LSU rRNA to only those structural features preserved in all living cells. This suggests that the E.cuniculi LSU rRNA performs only the basic mechanisms of translation. LSU rRNA phylogenetic analysis with the BASEML program favours a relatively recent origin of the fast evolving microsporidian lineage. Therefore, the prokaryote-like ribosomal features, such as the absence of ITS2, may be derived rather than primitive characters.

A molecular phylogeny of the Aphidiinae (Hymenoptera: Braconidae)

Phylogenetic relationships within the Aphidiinae, and between this and other subfamilies of Braconidae (Hymenoptera), were investigated using sequence data from three genes: elongation factor-1alpha, cytochrome b, and the second expansion segment of the 28S ribosomal subunit. Variation in both protein-coding genes was characterized by a high level of homoplasy, but analysis of the expansion segment--robust over a range of alignment methods and parameters-resolved some of the older divergences. Parsimony analysis of the combined data suggests the following tribal relationships: (Ephedrini + (Praini + (Aphidiini + Trioxini))). In addition, the cyclostome subfamilies were found to form a clade separate from the Aphidiinae, but relationships between the Aphidiinae and the noncyclostome braconids could not be resolved. The inferred phylogeny also supported a secondary loss of internal pupation within the Praini and a polyphyletic origin of endoparasitism within the Braconidae.

Comparative variation of morphological and molecular evolution through geologic time: 28 S ribosomal RNA versus morphology in echinoids

The comparatively good fossil record of post-Palaeozoic echinoids allows rates of morphological change to be estimated over the past 260 million years and com pared with rates of molecular evolution. Parsimony analysis of morphological data, based predominantly on skeletal characteristics, and parsimony, distance and maximum likelihood analyses of molecular data, from the first 380 bases from the 5' end of the 28 S rRNA molecule, for 10 species of echinoid produce congruent phylogenies. The m olecular sequence chosen is dem onstrably far from saturation and sister groups have divergence times ranging from about 15 to 260 Ma. Parsimony analysis allows the great majority of molecular and morphological apomorphies to be placed in one of 18 independent geological time intervals, providing a direct measure of rates of evolution for periods in the geological past. Because most molecular fixed point m utations in our sequences cannot be polarized unambiguously by outgroup comparison (making the outgroup states effectively random), distance and parsimony analyses both tend spuriously to root the echinoid tree on the longest internal branch. A topology identical to that derived from morphological data is, however, obtained using Maximum Likelihood and also parsimony analysis where outgroup rooting is restricted to more conserved regions. This is taken as the correct topology for assessing rates of evolution. Overall, both morphological and molecular changes show a m oderately strong correlation with time elapsed, but a weaker correlation with one another. Statistically significant differences in evolutionary rate are found between some, but not all, pair-wise comparisons of sister lineages for both molecular and morphological data. The molecular clock rate for echinaceans is three times faster than that for cidaroids and irregular echinoids. Spearm an’s rank correlation test, which requires only relative m agnitude of changes to be known, suggests that morphological change has a slightly better correlation with time than does molecular change, averaged over all ten species. However, when just echinaceans are considered an extremely good correlation is found between the num ber of molecular changes and time elapsed, whereas morphological change remains poorly correlated. Thus, molecular rates approxim ate to a clocklike model within restricted echinoid elades, but vary significantly between clades. Averaging results over all echinoids produces a correlation that is no better than the correlation between morphological change and time elapsed.

Divergence towards a dead end? Cleavage of the divergent domains of ribosomal RNA in apoptosis

In several cases of apoptotic death the large ribosomal subunit 28S rRNA is specifically cleaved. The cleavages appear at specific sites within those domains of the rRNA molecule that have shown exceptional high divergence in evolution (D domains). The cleavages accompany rather than precede apoptosis, and there is a positive, but not complete, correlation between rRNA cleavage and internucleosomal DNA fragmentation. Most cell types studied so far show two alternative cleavage pathways that are mutually exclusive. Cleavage can either start in the D8 domain with secondary cuts within a subdomain of D2 (D2c), or in the D2 domain with subsequent excision of the D2c subdomain. The latter pathway is of particular interest since D2 (unlike D8) is normally inaccessible for RNase attack. That apoptosis specifically affects the ribosomal divergent domains suggests that these domains, which make up roughly 25% of total cellular RNA, might have evolved to serve functions related to apoptosis. Future studies will be directed to test the hypothesis that rRNA fragmentation may be part of an apoptotic program directed against the elimination of illegitimate (viral?) polynucleotides.

Secondary structures and features of the 18S, 5.8S and 26S ribosomal RNAs from the Apicomplexan parasite Toxoplasma gondii

The two major subunits of the ribosomal RNA (rRNA) of Toxoplasma gondii, 18S and 26S, as well as 5.8S, have been sequenced and folded according to known consensus and established secondary structures. Conserved and variable nucleotide (nt) regions were identified using multiple alignments with rRNA sequences of selected organisms. The 18S rRNA showed a well conserved core structure of 48 stems and a hypervariable V4 region identified four additional stems including a pseudoknot. The 18S rRNA contained an additional helix in the V2 region located between nt 204 to 258. We noted that T. gondii 18S does not have a true V6 region, but was organized as a motif of a simple stem. T. gondii 26S had a conserved core structure of 83 stems and its expansion segments, so-called divergent domains, demonstrated a high degree of similarity with secondary structures from rRNA of dinoflagellates and ciliates. For the T. gondii 26S sequence, we found two additional stems, D3d and D3e, composed of 140 nt having a higher deltaG value. These segments are absent from the prokaryotic rRNA structures, whereas the hypervariable V4 region of the small subunit is not as variable. The well preserved structures could indicate an additional function for the eukaryotic ribosome.

Sequence and characterisation of a ribosomal RNA operon from Agrobacterium vitis

One of the four ribosomal RNA operons (rrnA) from the Agrobacterium vitis vitopine strain S4 was sequenced, rrnA is most closely related to the rrn operons of Bradyrhizobium japonicum and Rhodobacter sphaeroides and carries an fMet-tRNA gene downstream of its 5S gene, as in the case of R. sphaeroides. The 16S rRNA sequence of S4 differs from the A. vitis K309 type strain sequence by only one nucleotide, in spite of the fact that S4 and K309 have very different Ti plasmids. The predicted secondary structure of the S4 23S rRNA shows several features that are specific for the alpha proteobacteria, and an unusual branched structure in the universal B8 stem. The 3' ends of the three other rrn copies of S4 were also cloned and sequenced. Sequence comparison delimits the 3' ends of the four repeats and defines two groups: rrnA/rrnB and rrnC/rrnD.

Simplicity-correlated size growth of the nuclear 28S ribosomal RNA D3 expansion segment in the crustacean order Isopoda

The expansion segments within the eukaryote nuclear 23S-like ribosomal RNA molecule are now well characterized in many diverse organisms. A different base compositional bias, a higher propensity for size variability, and an increased evolutionary rate distinguish these regions from the universally conserved "core" regions of the molecule. In addition, some expansion segments of higher eukaryotes exhibit significant sequence simplicity which is hypothesized to occur by slippage-mediated mutational processes. We describe the discovery of extreme size variation of the D3 expansion segment in the crustacean order Isopoda. Among 11 species D3 varies in size from 180 to 518 nucleotides but maintains a homologous secondary structure. The D3 size is significantly positively correlated to relative simplicity factor (RSF), indicating that growth is most likely by insertion of simple sequences. D3 size and RSF correlate approximately with a morphology-based phylogeny, and within oniscideans RSF increases as more recent divergences occur. The D3 of Armadillidium vulgare, with an RSF of 1.87, is the highest value recorded for any known expansion segment. Regions of high sequence simplicity in nuclear ribosomal RNA were previously only known from the higher vertebrate lineage. Here we demonstrate that this phenomenon occurs in a more extreme condition within a monophyletic invertebrate lineage. The extreme size changes identified could indicate that expansion segments are an extraneous element in the functioning ribosome.

Localization of nucleolin binding sites on human and mouse pre-ribosomal RNA

Nucleolin, a major RNA binding protein of the nucleolus is found associated mainly to the pre-ribosomal particles and is absent from the cytoplasmic mature ribosomes. The role of this protein in ribosome biogenesis remains largely unknown, and is likely to be reflected by its RNA binding properties. Nucleolin contains in its central domain four RNA recognition motifs (RRM, also called RBD for RNA binding domain) which are conserved among different species. RNA binding studies have revealed that nucleolin interacts specifically with a short stem loop structure called NRE (nucleolin recognition element). We show that nucleolin extracted from human, hamster and mouse cells interacts with the same specificity and affinity to a mouse 5'ETS (external transcribed spacer) RNA fragment which contains a NRE motif. A similar structure within the human 5'ETS is also efficiently recognized by mouse nucleolin. We identified putative NRE not only in the 5'ETS but also in the 3'ETS, ITS (internal transcribed spacer) and in the 18S and 28S RNA sequences. This is in agreement with in vivo cross-linking data and a previous immunocytological analysis of ribosomal transcription units. Interestingly, we found that all the NRE localized in the 28S region are within the variable domains. Despite considerable sequence divergence of these domains, several of the NRE have sequences perfectly conserved between these two species. This suggests that these nucleolin binding sites might be functionally important, in particular for ribosome biogenesis.

Structural constraints in expansion segments from a midge 26S rDNA

DNA sequences representing approximately 40% of the large-subunit rRNA gene from the lower dipteran Chironomus thummi were analyzed. Once aligned with their Drosophila counterparts, sequence and base content comparisons were carried out. Sequence identity was found to be high overall, except for six regions that displayed a local bias in nucleotide composition toward AT. These regions were identified as expansion segments D3, D4, D5, D6, D7a, and D12. Besides base sequence divergence, differences in length were observed between the respective variable domains of the two species, particularly for D7a. Prediction of secondary structure showed that the folding of the Chironomus expansion segments analyzed is in agreement with the general patterns proposed for eukaryotic LSU rRNA. The comparison with Drosophila revealed also that the Chironomus secondary structures of the variable domains are supported by multiple compensatory substitutions or even compensatory insertions. Chironomus D7a displayed an unusual structural feature with respect to the insect D7a models that have been inferred up to now. The structural constraint observed in the expansion segments of Diptera so distantly related as midges and Drosophila suggests that these regions contribute to some functional role. Concerning the D7a of insects so far analyzed, there can be, in addition to a conserved secondary structure, a nucleotide composition constraint that might be important for the process giving rise to the alpha and beta halves of the 26S rRNA.

Evolution of the rDNA spacer, ITS 2, in the ticks Ixodes scapularis and I. pacificus (Acari: Ixodidae)

Evolution of the rDNA spacer, ITS 2, is examined by comparing 17 DNA sequences of the ticks, Ixodes scapularis and I. pacificus. The distribution of fixed interspecific differences and the relative frequency of base changes vs. insertions/deletions (indels) matches the distribution and relative frequency for intraspecifically variable sites. This suggests that most intraspecific variation is not effectively selected against. The base composition of the ITS 2 transcript is G- and U-biased. But, 5-base regions enriched (> 80 per cent) for A or U occur more frequently than expected while G- and C-enriched regions occur less frequently than expected. Enriched sequences may be prone to replication slippage, accounting for the A/T bias in insertions. Slippage-mediated gains and losses of A/T-rich tandem repeats apparently account for most indels. Minimum-energy conformations of the two species' folded transcripts share major structural features. Structural inertia arises from intramolecular base pairing within stems that allows most mutations to be absorbed as new bulges off stems. Yet, there is evidence of selection to maintain the conformation. First, intraspecifically variable sites are concentrated at the ends of stems in loops and intersections, structures that do not contribute to intramolecular base pairing. Moreover, some indels that have become fixed in one species compensate for the presence of conformation-destabilizing indels. However, high rates of sequence evolution within stems and absence of compensatory base evolution contraindicates selective constraint. Degenerate dispersed and tandem copies of two subrepeats, each approximately 20 bases long, may account for much of the ITS 2 sequence. These are approximately inverses of each other and are, consequently, capable of significant intramolecular hydrogen bonding to produce folded transcripts of low energy. Evolution of the ITS 2 sequence may largely entail replication slippage-mediated gains and losses of these repeats or their composite subrepeats.

Fine mapping of 28S rRNA sites specifically cleaved in cells undergoing apoptosis

Bona fide apoptosis in rat and human leukemia cells, rat thymocytes, and bovine endothelial cells was accompanied by limited and specific cleavage of polysome-associated and monosome-associated 28S rRNA, with 18S rRNA being spared. Specific 28S rRNA cleavage was observed in all instances of apoptotic death accompanied by internucleosomal DNA fragmentation, with cleavage of 28S rRNA and of DNA being linked temporally. This indicates that 28S rRNA fragmentation may be as general a feature of apoptosis as internucleosomal DNA fragmentation and that concerted specific cleavage of intra- and extranuclear polynucleotides occurs in apoptosis. Apoptosis-associated cleavage sites were mapped to the 28S rRNA divergent domains D2, D6 (endothelial cells), and D8. The D2 cuts occurred in hairpin loop junctions considered to be buried in the intact ribosome, suggesting that this rRNA region becomes a target for RNase attack in apoptotic cells. D8 was cleaved in two exposed UU(U) sequences in bulge loops. Treatment with agents causing necrotic cell death or aging of cell lysates failed to produce any detectable limited D2 cleavage but did produce a more generalized cleavage in the D8 region. Of potential functional interest was the finding that the primary cuts in D2 exactly flanked a 0.3-kb hypervariable subdomain (D2c), allowing excision of the latter. The implication of hypervariable rRNA domains in apoptosis represents the first association of any functional process with these enigmatic parts of the ribosomes.

Phylogenetic relationship of the green alga Nanochlorum eukaryotum deduced from its chloroplast rRNA sequences

The marine green coccoidal alga Nanochlorum eukaryotum (N.e.) is of small size with an average diameter of 1.5 microns. It is characterized by primitive-appearing biochemical and morphological properties, which are considerably different from those of other green algae. Thus, it has been proposed that N.e. may be an early developed algal form. To prove this hypothesis, DNA of N.e. was isolated by a phenol extraction procedure, and the chloroplast DNA separated by preparative CsCl density-gradient centrifugation. The kinetic complexity of the nuclear and of the chloroplast DNA was evaluated by reassociation kinetics to 3 x 10(7) bp and 9 x 10(4) bp, respectively. Several chloroplast genes, including the rRNA genes, were cloned on distinct fragments. The order of the rRNA genes corresponds to the common prokaryotic pattern. The 16S rRNA gene comprises 1,548 bases and is separated from the 23S rRNA gene with its 2,920 bases by a short spacer of 460 bases, which also includes the tRNA(Ile) and tRNA(Ala) genes. The 5S rRNA gene has not been found; it must start further than 500 bases downstream from the 3'-end of the 23S rRNA gene. From the chloroplast rRNA sequences, we have deduced secondary structures of the 16S and 23S rRNAs, which are in agreement with standard models. The rRNA sequences were aligned with corresponding chloroplast sequences; phylogenetic relationships were calculated by several methods. From these calculations, we conclude that N.e. is most closely related to Chlorella vulgaris. Therefore, N.e. does not represent an early developed algal species; the primitive-appearing morphological and biochemical characteristics of N.e. must rather be explained by secondary losses.

Chemical secondary structure probing of two highly methylated regions in Xenopus laevis 28S ribosomal RNA

The large ribosomal subunit (LSU) RNA or 28S rRNA of vertebrates is characterized by two highly conserved and methylated regions towards the 3' end of the molecule that extend from domains IV to V of the molecule. In this report we describe the probing of the secondary structure of these two highly methylated regions in Xenopus laevis LSU RNA by chemical modification using the single-strand nucleotide specific probes; dimethyl sulphate (DMS) and 1-cyclo-hexyl-3-(2-morpholinoethyl-carbodiimide metho-p-toluene sulphonate (CMCT) followed by primer extension. The same regions in in vitro synthesized unmethylated X. laevis 28S rRNA were also probed for comparison. Our results in general tend to support the theoretically determined secondary structure model for the probed domains. From the results obtained, methylated cellular LSU RNA appears to be relatively more reactive than the in vitro transcript to the chemical probes. Accessibility to the probes was found to be similar at most sites for cellular and in vitro transcript LSU RNAs. This implies that structural destabilization due to 2'-O-methylations in cellular LSU RNA is not significant.

An rRNA variable region has an evolutionarily conserved essential role despite sequence divergence

Regions extremely variable in size and sequence occur at conserved locations in eukaryotic rRNAs. The functional importance of one such region was determined by gene reconstruction and replacement in Tetrahymena thermophila. Deletion of the D8 region of the large-subunit rRNA inactivates T. thermophila rRNA genes (rDNA): transformants containing only this type of rDNA are unable to grow. Replacement with an unrelated sequence of similar size or a variable region from a different position in the rRNA also inactivated the rDNA. Mutant rRNAs resulting from such constructs were present only in precursor forms, suggesting that these rRNAs are deficient in either processing or stabilization of the mature form. Replacement with D8 regions from three other organisms restored function, even though the sequences are very different. Thus, these D8 regions share an essential functional feature that is not reflected in their primary sequences. Similar tertiary structures may be the quality these sequences share that allows them to function interchangeably.

An rRNA variable region has an evolutionarily conserved essential role despite sequence divergence

Regions extremely variable in size and sequence occur at conserved locations in eukaryotic rRNAs. The functional importance of one such region was determined by gene reconstruction and replacement in Tetrahymena thermophila. Deletion of the D8 region of the large-subunit rRNA inactivates T. thermophila rRNA genes (rDNA): transformants containing only this type of rDNA are unable to grow. Replacement with an unrelated sequence of similar size or a variable region from a different position in the rRNA also inactivated the rDNA. Mutant rRNAs resulting from such constructs were present only in precursor forms, suggesting that these rRNAs are deficient in either processing or stabilization of the mature form. Replacement with D8 regions from three other organisms restored function, even though the sequences are very different. Thus, these D8 regions share an essential functional feature that is not reflected in their primary sequences. Similar tertiary structures may be the quality these sequences share that allows them to function interchangeably.

Molecular phylogeny of Drosophila based on ribosomal RNA sequences

Nucleotide sequences of 72 species of Drosophilidae were determined for divergent D1 and D2 domains (representing 200 and 341 nucleotides respectively in D. melanogaster) of large ribosomal RNA, using the rRNA direct sequencing method. Molecular phylogenetic trees were reconstructed using both distance and parsimony methods and the robustness of the nodes was evaluated by the bootstrap procedure. The trees obtained by these methods revealed four main lineages or clades which do not correspond to the taxonomical hierarchy. In our results, the genus Chymomyza is associated with the subgenus Scaptodrosophila of the genus Drosophila and their cluster constitutes the most ancient clade. The two other clades are constituted of groups belonging to the subgenus Sophophora of the genus Drosophila: the so-called Neotropical clade including the willistoni and saltans groups and the obscura-melanogaster clade itself split into three lineages: (1) obscura group + ananassae subgroup, (2) montium subgroup, and (3) melanogaster + Oriental subgroups. The fourth clade, the Drosophila one, contains three lineages. D. polychaeta, D. iri, and D. fraburu are branched together and constitute the most ancient lineage; the second lineage includes the annulimana, bromeliae, dreyfusi, melanica, mesophragmatica, repleta, robusta, and virilis groups. The third lineage is composed of the immigrans and the cardini, funebris, guaramunu, guarani, histrio, pallidipennis, quinaria, and tripunctata groups. The genera Samoaia, Scaptomyza, and Zaprionus are branched within the Drosophila clade. Although these four clades appear regularly in almost all tree calculations, additional sequencing will be necessary to determine their precise relationships.

Conserved secondary structures in the ITS2 of trematode pre-rRNA

Biological functions of transcribed spacer regions in eucaryotic pre-rRNAs remain elusive so far. Utilization of the comparative approach to analyse their secondary structure has been impeded by the extensive sequence divergence observed among most of the specimens available to date. However, we have taken advantage of the recent derivation of a set of largely similar sequences for different schistosome species to look for the presence of constrained secondary structures within the internally transcribed spacer 2 (ITS2). We show that several common features of secondary structure are shared by these species despite sequence variation, with some of them supported by compensatory changes, suggesting a significant role of ITS2 as an RNA domain during ribosome biogenesis.