Evolution of large-subunit rRNA structure. The diversification of divergent D3 domain among major phylogenetic groups

Comparative analysis of the mitogenomes of two Corydoras (Siluriformes, Loricarioidei) with nine known Corydoras, and a phylogenetic analysis of Loricarioidei

Corydoras is a speciose catfish genus from South America with widely investigated phylogenetic and evolutionary relationships. The complete mitogenomes of C. aeneus and C. paleatus were sequenced, assembled, and annotated using next-generation sequencing. The genome arrangements, gene contents, genome structures, base compositions, evolutionary features, codon usage, and tRNA structures of the two mitogenomes were compared and analyzed with nine published mitogenomes of Corydoras. Phylogenetic analysis was performed using concatenated nucleotide sequences with 13 protein-coding genes and two rRNAs with 44 mitogenomes of Siluriformes. These results provide information on the mitogenomes of eleven Corydoras species and evolutionary relationships within the suborder Loricarioidei, which may be applicable for further phylogenetic and taxonomic studies on Siluriformes and Loricarioidei.

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.

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.

Ribosomal cistrons in Bothrops neuwiedi (Serpentes) subspecies from Brazil

The ribosomal cistrons of six subspecies of Bothrops neuwiedi (Serpentes) were studied at both the cytogenetic and molecular levels. These subspecies populations occur in several Brazilian regions. The analysis of the nucleolar organizing region banding patterns showed variability in the chromosomal localization of rDNA cistrons. The rDNA clusters were found in two microchromosomes, in chromosome pair Number 6, or even in the one homologue of chromosome 6, and one microchromosome, according to the specimen examined but independent of the population from which it was selected. The organization of the rDNA repeat was studied by Southern blot hybridization using Xenopus laevis ribosomal DNA probes. The size of the repeat is 10.4 kilobases (kb), and the intergenic spacer (IGS), 2.6 kb long, is relatively small compared with the size of other vertebrate IGSs. Restriction mapping using the restriction enzymes EcoRI, BamHI, HindIII, and PvuII showed a highly conserved organization of the ribosomal repeats in the total genomic DNA of the subspecies studied.

Phylogenetic study of nine species of freshwater monogeneans using secondary structure and motif prediction from India

The present study was performed to identify and validate monogenean species from different piscine hosts using molecular tools. Nine species of freshwater monogeneans were collected from gills and skin of freshwater fishes at Hastinapur, Meerut, India. After microscopic examination, molecular analysis was performed utilizing 28S gene marker. Phylogenetic analysis indicated the validation and systematic position of these nine different monogeneans belongs to the Dactylogyridae and Gyrodactylidae families. The findings also confirm that the 28S rDNA sequence is highly conserved and may prove to be useful in taxonomic studies of parasitic platyhelminthes. Besides this, the study is also supplemented by molecular morphometrics that is based on 28S secondary structure homologies of nine monogenean species. The data indicate that 28S motifs i.e., ≤ 50bp in size can also be considered a promising tool for monogenean species identification and their validation.

Factors that affect large subunit ribosomal DNA amplicon sequencing studies of fungal communities: classification method, primer choice, and error

Nuclear large subunit ribosomal DNA is widely used in fungal phylogenetics and to an increasing extent also amplicon-based environmental sequencing. The relatively short reads produced by next-generation sequencing, however, makes primer choice and sequence error important variables for obtaining accurate taxonomic classifications. In this simulation study we tested the performance of three classification methods: 1) a similarity-based method (BLAST + Metagenomic Analyzer, MEGAN); 2) a composition-based method (Ribosomal Database Project naïve bayesian classifier, NBC); and, 3) a phylogeny-based method (Statistical Assignment Package, SAP). We also tested the effects of sequence length, primer choice, and sequence error on classification accuracy and perceived community composition. Using a leave-one-out cross validation approach, results for classifications to the genus rank were as follows: BLAST + MEGAN had the lowest error rate and was particularly robust to sequence error; SAP accuracy was highest when long LSU query sequences were classified; and, NBC runs significantly faster than the other tested methods. All methods performed poorly with the shortest 50-100 bp sequences. Increasing simulated sequence error reduced classification accuracy. Community shifts were detected due to sequence error and primer selection even though there was no change in the underlying community composition. Short read datasets from individual primers, as well as pooled datasets, appear to only approximate the true community composition. We hope this work informs investigators of some of the factors that affect the quality and interpretation of their environmental gene surveys.

Fragmentation of the large subunit ribosomal RNA gene in oyster mitochondrial genomes

Background

Discontinuous genes have been observed in bacteria, archaea, and eukaryotic nuclei, mitochondria and chloroplasts. Gene discontinuity occurs in multiple forms: the two most frequent forms result from introns that are spliced out of the RNA and the resulting exons are spliced together to form a single transcript, and fragmented gene transcripts that are not covalently attached post-transcriptionally. Within the past few years, fragmented ribosomal RNA (rRNA) genes have been discovered in bilateral metazoan mitochondria, all within a group of related oysters.

Results

In this study, we have characterized this fragmentation with comparative analysis and experimentation. We present secondary structures, modeled using comparative sequence analysis of the discontinuous mitochondrial large subunit rRNA genes of the cupped oysters C. virginica, C. gigas, and C. hongkongensis. Comparative structure models for the large subunit rRNA in each of the three oyster species are generally similar to those for other bilateral metazoans. We also used RT-PCR and analyzed ESTs to determine if the two fragmented LSU rRNAs are spliced together. The two segments are transcribed separately, and not spliced together although they still form functional rRNAs and ribosomes.

Conclusions

Although many examples of discontinuous ribosomal genes have been documented in bacteria and archaea, as well as the nuclei, chloroplasts, and mitochondria of eukaryotes, oysters are some of the first characterized examples of fragmented bilateral animal mitochondrial rRNA genes. The secondary structures of the oyster LSU rRNA fragments have been predicted on the basis of previous comparative metazoan mitochondrial LSU rRNA structure models.

Non-monophyly of most supraspecific taxa of calcareous sponges (Porifera, Calcarea) revealed by increased taxon sampling and partitioned Bayesian analysis of ribosomal DNA

Calcareous sponges (Porifera, Calcarea) play an important role for our understanding of early metazoan evolution, since several molecular studies suggested their closer relationship to Eumetazoa than to the other two sponge 'classes,' Demospongiae and Hexactinellida. The division of Calcarea into the subtaxa Calcinea and Calcaronea is well established by now, but their internal relationships remain largely unresolved. Here, we estimate phylogenetic relationships within Calcarea in a Bayesian framework, using full-length 18S and partial 28S ribosomal DNA sequences. Both genes were analyzed separately and in combination and were further partitioned by stem and loop regions, the former being modelled to take non-independence of paired sites into account. By substantially increasing taxon sampling, we show that most of the traditionally recognized supraspecific taxa within Calcinea and Calcaronea are not monophyletic, challenging the existing classification system, while monophyly of Calcinea and Calcaronea is again highly supported.

Maintaining a conserved methylation in plant and insect U2 snRNA through compensatory mutation by nucleotide insertion

The extensive post-transcriptional modification of U2 snRNA is required for spliceosome assembly and pre-mRNA splicing in vertebrates. However, the rare modification of U2 snRNA in yeast implies a different mechanism for regulating spliceosome biogenesis in single-celled eukaryotes. To understand the evolutionary pattern of U2 snRNA methylation, we determined for the first time, the 2'-O-methylations of U2 snRNA in Oryza sativa, Arabidopsis thaliana and Drosophila melanogaster, and revealed two methylations which are conserved in a crucial region of U2 snRNA in plants. Interestingly, one of the methylations, U2-Cm29 is also methylated in D. melanogaster, but not in vertebrates. According to the methylation of U2-C29, computational analysis of databases identified three canonical box C/D snoRNAs, named OsmgU2-29, AtmgU2-29 and DmmgU2-28, as small methylation guides of U2 snRNA from O. sativa, A. thaliana and D. melanogaster, respectively. Although very divergent in their sequence, the three snoRNAs exhibit in common an 11 nucleotide-long sequence complementarity to corresponding U2 snRNA, implying a functional constraint on the modification during evolution. Interestingly, a nucleotide is found to be inserted both in U2 snRNA and DmmgU2-28 and maintains a perfect match of duplex specifying the methylation of C28 in Drosophila U2 snRNA. This is the first time a new model is being provided for compensatory mutations between a small guide RNA and its target by nucleotide insertion, instead of the known nucleotide substitution. In contrast to small Cajal body-specific RNAs (scaRNAs), the snoRNAs are similar to the reported singlet guide RNAs and are known to localize in nucleolus.

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.

Forensic analysis of hallucinogenic fungi: a DNA-based approach

Hallucinogenic fungi synthesize two controlled substances, psilocin and psilocybin. Possession of the fungal species that contain these compounds is a criminal offence in North America. Some related species that are morphologically similar, do not contain the controlled substances. Therefore, unambiguous identification of fungi to the species level is critical in determining if a mushroom is illegal. We investigate a phylogenetic approach for the identification of species that contain the psychoactive compounds. We analyzed 35 North American specimens representing seven different genera of hallucinogenic and non-hallucinogenic mushrooms. We amplified and sequenced the internal transcribed spacer region of the rDNA (ITS-1) and a 5' portion of the nuclear large ribosomal subunit of rRNA (nLSU rRNA or 28S). ITS-1 locus sequence data was highly variable and produced a phylogenetic resolution that was not consistent with morphological identifications. In contrast, the nLSU rRNA data clustered isolates from the same species and separated hallucinogen containing and non-hallucinogen containing isolates into distinct clades. With this information, we propose an approach that combines the specificity of PCR detection and the resolving power of phylogenetic analysis to efficiently and unambiguously identify hallucinogenic fungal specimens for legal purposes.

The high diversity of snoRNAs in plants: identification and comparative study of 120 snoRNA genes from Oryza sativa

Using a powerful computer-assisted analysis strategy, a large-scale search of small nucleolar RNA (snoRNA) genes in the recently released draft sequence of the rice genome was carried out. This analysis identified 120 different box C/D snoRNA genes with a total of 346 gene variants, which were predicted to guide 135 2'-O-ribose methylation sites in rice rRNAs. Though not exhaustive, this analysis has revealed that rice has the highest number of known box C/D snoRNAs among eukaryotes. Interestingly, although many snoRNA genes are conserved between rice and Arabidopsis, almost half of the identified snoRNA genes are rice specific, which may highlight further the differences in rRNA methylation patterns between monocotyledons and dicotyledons. In addition to 76 singletons, 70 clusters involving 270 snoRNA genes were also found in rice. The large number of the novel snoRNA polycistrons found in the introns of rice protein-coding genes is in contrast to the one-snoRNA-per-intron organization of vertebrates and yeast, and of Arabidopsis in which only a few intronic snoRNA gene clusters were identified. Furthermore, due to a high degree of gene duplication, rice snoRNA genes are clearly redundant and exhibit great sequence variation among isoforms, allowing generation of new snoRNAs for selection. Thus, the large snoRNA gene family in plants can serve as an excellent model for a rapid and functional evolution.

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.

Cirripede phylogeny using a novel approach: molecular morphometrics

We present a new method using nucleic acid secondary structure to assess phylogenetic relationships among species. In this method, which we term "molecular morphometrics," the measurable structural parameters of the molecules (geometrical features, bond energies, base composition, etc.) are used as specific characters to construct a phylogenetic tree. This method relies both on traditional morphological comparison and on molecular sequence comparison. Applied to the phylogenetic analysis of Cirripedia, molecular morphometrics supports the most recent morphological analyses arguing for the monophyly of Cirripedia sensu stricto (Thoracica + Rhizocephala + Acrothoracica). As a proof, a classical multiple alignment was also performed, either using or not using the structural information to realign the sequence segments considered in the molecular morphometrics analysis. These methods yielded the same tree topology as the direct use of structural characters as a phylogenetic signal. By taking into account the secondary structure of nucleic acids, the new method allows investigators to use the regions in which multiple alignments are barely reliable because of a large number of insertions and deletions. It thus appears to be complementary to classical primary sequence analysis in phylogenetic studies.

Phylogenetic relationships in the mushroom genus Coprinus and dark-spored allies based on sequence data from the nuclear gene coding for the large ribosomal subunit RNA: divergent domains, outgroups, and monophyly

Phylogenetic relationships were investigated in the mushroom genus Coprinus based on sequence data from the nuclear encoded large-subunit rDNA gene. Forty-seven species of Coprinus and 19 additional species from the families Coprinaceae, Strophariaceae, Bolbitiaceae, Agaricaceae, Podaxaceae, and Montagneaceae were studied. A total of 1360 sites was sequenced across seven divergent domains and intervening sequences. A total of 302 phylogenetically informative characters was found. Ninety-eight percent of the average divergence between taxa was located within the divergent domains, with domains D2 and D8 being most divergent and domains D7 and D10 the least divergent. An empirical test of phylogenetic signal among divergent domains also showed that domains D2 and D3 had the lowest levels of homoplasy. Two equally most parsimonious trees were resolved using Wagner parsimony. A character-state weighted analysis produced 12 equally most parsimonious trees similar to those generated by Wagner parsimony. Phylogenetic analyses employing topological constraints suggest that none of the major taxonomic systems proposed for subgeneric classification is able to completely reflect phylogenetic relationships in Coprinus. A strict consensus integration of the two Wagner trees demonstrates the problematic nature of choosing outgroups within dark-spored mushrooms. The genus Coprinus is found to be polyphyletic and is separated into three distinct clades. Most Coprinus taxa belong to the first two clades, which together form a larger monophyletic group with Lacrymaria and Psathyrella in basal positions. A third clade contains members of Coprinus section Comati as well as the genus Leucocoprinus, Podaxis pistillaris, Montagnea arenaria, and Agaricus pocillator. This third clade is separated from the other species of Coprinus by members of the families Strophariaceae and Bolbitiaceae and the genus Panaeolus.

Evolutionarily conserved structural features in the ITS2 of mammalian pre-rRNAs and potential interactions with the snoRNA U8 detected by comparative analysis of new mouse sequences

Mechanisms of ITS2 excision from pre-rRNA remain largely elusive. In mammals, at least two endonucleolytic cleavages are involved, which result in the transient accumulation of precursors to 5.8S rRNA termed 8S and 12S RNAs. We have sequenced ITS2 in four new species of the Mus genus and investigated its secondary structure using thermodynamic prediction and comparative approach. Phylogenetic evidence supports an ITS2 folding organized in four domains of secondary structure extending from a preserved structural core. This folding is also largely conserved for the previously available mammalian ITS2 sequences, rat and human, despite their extensive sequence divergence relative to the Mus species. Conserved structural features include the structural core, containing the 3' end of 8S pre-rRNA within a single-stranded sequence, and a stem containing the 3' end of the 12S pre-rRNA species. A putative, phylogenetically preserved pseudoknot has been detected 1 nt downstream from the 12S 3' end. Two long complementarities have also been identified, in sequences conserved among vertebrates, between the pre-rRNA 32S and the snoRNA (small nucleolar RNA) U8 which is required for the excision of Xenopus ITS2. The first complementarity involves the 5.8S-ITS2 junction and 13 nt at the 5' end of U8, whereas the other one occurs between a mature 28S rRNA segment known to be required for ITS2 excision and positions 15-25 of snoRNA U8. These two potential interactions, in combination with ITS2 folding, could organize a functional pocket containing three cleavage sites and key elements for pre-rRNA processing, suggesting a chaperone role for the snoRNA U8.

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.

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

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

RAGA: RNA sequence alignment by genetic algorithm

We describe a new approach for accurately aligning two homologous RNA sequences when the secondary structure of one of them is known. To do so we developed two software packages, called RAGA and PRAGA, which use a genetic algorithm approach to optimize the alignments. RAGA is mainly an extension of SAGA, an earlier package for multiple protein sequence alignment. In PRAGA several genetic algorithms run in parallel and exchange individual solutions. This method allows us to optimize an objective function that describes the quality of a RNA pairwise alignment, taking into account both primary and secondary structure, including pseudoknots. We report results obtained using PRAGA on nine test cases of pairs of eukaryotic small subunit rRNA sequence (nuclear and mitochondrial).

Molecular cloning and cold-inducible gene expression of ferritin H subunit isoforms in rainbow trout cells

Cold acclimation of rainbow trout cells is considered to be mediated by alterations in the mRNAs and proteins present in cold-treated cells. A subtracted cDNA library from cold-treated rainbow trout RTG-2 cells was constructed and screened to isolate cDNA induced in the cold-treated cells in order to elucidate which genes are induced by cold acclimation. A set of cDNA clones encoding three members of ferritin H isoforms was isolated as cold-inducible genes. Northern blot analysis and nuclear run-on transcription assay showed that the transcription and accumulation of the ferritin H isoforms mRNA were enhanced by cold acclimation. Furthermore, the ferritin level in the trout cells increased on cold acclimation in response to a temperature shift from 22 degrees C to 4 degrees C. When the trout cells were subjected to 4 degrees C under the condition of a decreased ferritin H level obtained by the addition of an antisense oligonucleotide, cell growth was apparently inhibited. These findings indicate an association between the induction of ferritin H and cellular mechanisms during cold acclimation of trout cells.

Sequence variations in the large-subunit ribosomal RNA gene of Ammonia (Foraminifera, Protozoa) and their evolutionary implications

An unusually high divergence was observed in the ribosomal RNA genes of a free-living population of foraminifera belonging to the genus Ammonia. The sequences of a large-subunit (LSU) rDNA expansion segment D1 and flanking regions were obtained from 20 specimens named Ammonia sp. 1 and Ammonia sp. 2. The sequence divergence between the two species averages 14%. Within each species it ranges from 0.2% to 7.1% in Ammonia sp. 1 and from 0.7% to 2.3% in Ammonia sp. 2. We did not find two specimens having identical sequences. Moreover, in opposition to the generally accepted view, rDNA sequence variations were also found within a single individual. The variations among several rDNA copies in a single specimen of Ammonia may reach up to 4.9%. Most of the observed variations result from multiplication of CA or TA serial repeats occurring in two particularly variable regions. For single base changes, C-T transitions are most frequently observed. We discuss the evolution of expansion segments and their use for phylogenetic studies.

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.

Molecular phylogenetics of the four Schistosoma species groups determined with partial 28S ribosomal RNA gene sequences

Partial 28S ribosomal RNA (rRNA) gene sequences, including the variable domains D1, D2 and D3, were determined for representative species from the 4 Schistosoma species groups. On an alignment of 1345 bp from S. mansoni, S. haematobium, S. spindale and S. japonicum (with Heterobilharzia americana chosen as an outgroup), both maximum likelihood and maximum parsimony analyses provide a robust molecular phylogeny for the genus; ((((S. haematobium, S. spindale), S. mansoni), S. japonicum), H. americana). When analysed separately, both domain D1 and domain D2 yielded similarly informative data whereas D3 failed to resolve the phylogeny. These results confirm a phylogeny previously suggested by 18S rRNA gene sequences, corroborating the status of S. spindale as a sister taxon to S. haematobium, and demonstrate the utility of 28S rRNA gene sequence data for resolving phylogenies within the Schistosomatidae.

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.

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.

Extrachromosomal ribosomal DNA of Didymium iridis: sequence analysis of the large subunit ribosomal RNA gene and sub-telomeric region

The ribosomal DNA of the myxomycete Didymium iridis is organized as extrachromosomal linear molecules of about 20 kb, containing only one transcription unit of the ribosomal RNA genes. We have determined the sequence of the large subunit ribosomal RNA (LSU rRNA) gene as well as the sub-telomeric and telomeric regions. The LSU rRNA gene was found to encode a 3857 nucleotide-long LSU rRNA, interrupted by a transcribed spacer and two group I introns. A complete secondary structure model of D. iridis LSU rRNA has been constructed. The compact sub-telomeric region of D. iridis rDNA was found to contain several directly repeated sequence elements that include the simple telomere motif TTAGGG. Based on pairwise comparisons of LSU rRNA sequences, the time of divergence between the two myxomycete genera Didymium and Physarum was estimated.

Cell-free translation systems from different eukaryotes differ in their sensitivity to a template sugar-phosphate backbone

To study the role of a messenger sugar-phosphate backbone in the ribosomal decoding process, poly(U) and poly(dT) template activity in different eukaryotic systems has been compared. 80S ribosomes from Saccharomyces cerevisiae appeared to be able to translate poly(dT) both in the presence and in the absence of elongation factors, contrary to poly(U). However, ribosomes from higher eukaryotes (wheat germ, rabbit liver) are completely inefficient in poly(dT) translation. Moreover, rabbit liver ribosomes fail to bind effectively phenylalanyl-tRNA in the presence of poly(dT) although the polynucleotide seems to interact with the ribosomal decoding center. It is also of particular interest that hybrid ribosomes formed from the yeast and rabbit liver subunits can translate poly(dT) only when the large ribosomal subunit from yeast is used.

Evolution of compensatory substitutions through G.U intermediate state in Drosophila rRNA

It has often been suggested that the frequently observed Watson-Crick base-pair compensatory substitutions in RNA helical structures occur mainly through a slightly deleterious G.U intermediate state. We have scored base substitutions in a set of 82 related Drosophila species for the D1 and D2 variable domains of the large rRNA subunit. In all locations where a G-C in equilibrium with A-U compensatory base change occurred, a G.U pair has been observed in one or several species. As this dominant process implies two transitions, their rate was far higher in paired regions (92%) than in unpaired regions (47%). The other types of compensation were rarer and no intermediate states were observed. Most of the G.U base pairs observed in a species are not slightly deleterious. The rate of evolution of compensatory substitution is close to that predicted by a simple model of compensatory substitution through slightly deleterious or slightly advantageous G.U pairs, although some exceptions are presented.

Secondary structure constraints on the evolution of Drosophila 28 S ribosomal RNA expansion segments

Eukaryotic ribosomal RNA genes contain rapidly evolving regions of unknown function termed expansion segments. We present the comparative analysis of the primary and secondary structure of two expansion segments from the large subunit rRNA gene of ten species of Drosophila and the tsetse fly species Glossina morsitans morsitans. At the primary sequence level, most of the differences observed in the sequences obtained are single base substitutions. This is in marked contrast with observations in vertebrate species in which the insertion or deletion of repetitive motifs, probably generated by a DNA-slippage mechanism, is a major factor in the evolution of these regions. The secondary structure of the two regions, supported by multiple compensatory base changes, is highly conserved between the species examined and supports the existence of a general folding pattern for all eukaryotes. Intriguingly, the evolutionary rate of expansion segments is very slow relative to other genic and non-genic regions of the Drosophila genome. These results suggest that the evolution of expansion segments in the rDNA multigene family is a balance between the homogenization of new mutations by unequal crossing over and a combination of selection against some such mutations per se and selection for subsequent compensatory mutations, in order to maintain a particular RNA secondary structure.

rRNA genes from the lower chordate Herdmania momus: structural similarity with higher eukaryotes

Ascidians, primitive chordates that have retained features of the likely progenitors to all vertebrates, are a useful model to study the evolutionary relationship of chordates to other animals. We have selected the well characterized ribosomal RNA (rRNA) genes to investigate this relationship, and we describe here the cloning and characterization of an entire ribosomal DNA (rDNA) tandem repeat unit from a lower chordate, the ascidian Herdmania momus. rDNA copy number and considerable sequence differences were observed between two H. momus populations. Comparison of rDNA primary sequence and rRNA secondary structures from H. momus with those from other well characterized organisms, demonstrated that the ascidians are more closely related to other chordates than invertebrates. The rDNA tandem repeat makes up a larger percentage (7%) of the genome of this animal than in other higher eukaryotes. The total length of the spacer and transcribed region in H. momus rDNA is small compared to most higher eukaryotes, being less than 8 kb, and the intergenic spacer region consists of smaller internal repeats. Comparative analysis of rDNA sequences has allowed the construction of secondary structures for the 18S, 5.8S and 26S rRNAs.

A proposal for the secondary structure of a variable area of eukaryotic small ribosomal subunit RNA involving the existence of a pseudoknot

Eukaryotic small ribosomal subunit RNAs contain an area of variable structure, V4, which comprises about 250 nucleotides in most species, whereas the corresponding area in bacterial small ribosomal subunit RNAs consists of about 64 nucleotides folded into a single hairpin. There is no consensus on the secondary structure of area V4 in eukaryotes, about 10 different models having been proposed. The prediction of a model on a comparative basis poses special problems because, due to the variability of the area in length as well as sequence, a dependable alignment is very difficult to achieve. A new model was derived by systematic examination of all combinations of helices that have been hitherto proposed, plus some new ones. The following properties of the helices were examined: transposability to all presently known sequences, presence of compensating substitutions, and thermodynamic stability. A model was selected by ranking all possible combinations of transposable helices according to the number of compensating substitutions scored. The optimal model comprises a pseudoknot and four hairpin structures. Certain species contain additional hairpins inserted between these structural elements, while in others the structure is partially or entirely deleted.