Phenotypic effects of targeted mutations in the small subunit rRNA gene of Tetrahymena thermophila.

Tetrahymena thermophila is an ideal organism with which to study functional aspects of the rRNAs in vivo since the somatic rRNA genes of T. thermophila can be totally replaced by cloned copies introduced via microinjection. In this study, we made small insertions into seven sites within the small subunit rRNA gene and observed their phenotypic effects on transformed cells. Two mutated genes coding for rRNA (rDNAs), both of which bear insertions in highly conserved sequences, failed to transform and are therefore believed to produce nonfunctional rRNAs. Three other altered rDNAs produce functional rRNAs that can substitute for most or all of the cellular rRNA. Two of these bear insertions in highly variable regions, and, surprisingly, the other has an insertion in a region that is well conserved for both sequence and secondary structure among eucaryotes. In addition, two other insertions appear to destabilize rRNAs that contain them. Our findings make predictions concerning the positions of some of these sites within the tertiary structure of the small ribosomal subunit and thus serve as an in vivo test of the existing tertiary structure models for the small subunit rRNA. Our results are in good agreement with expectations based on sequence comparison and in vitro work.     

Paramecium mitochondrial genes. I. Small subunit rRNA gene sequence and microevolution

The sequences of the small subunit mitochondrial rRNA genes from two divergent species of Paramecium (primaurelia and tetraurelia) were determined. The gene lies near the center of the linear mitochondrial genome, on the same strand as are all other currently identified genes. The sequences generally resemble their counterparts found in cytoplasmic, procaryotic, and other mitochondrial sources. The rDNA gene boundaries were located by nuclease S1 protection. Small subunit rDNA spans about 1680 nucleotides, including an extraneous 83-base pair sequence very near the 3' end which is unique to Paramecium mitochondria. This "insert" occurs at the apex of the highly variable in length penultimate helix, according to proposed models for small subunit rRNA secondary structure. A discontinuity occurs in isolated rRNA near the start of the insert, resulting in a stable 13 S RNA species and a small segment containing the remaining 3' portion of the gene. The overall rRNA gene sequence was 94% conserved between the two species, and the nucleotide differences consisted of 53% transitions, 37% transversions, and 9% insertions plus deletions. These substitutions were somewhat clustered, and the two most divergent regions coincided with the gene boundaries. The sequence was aligned with Escherichia coli 16 S rRNA for direct comparison of sequence and structure.     

The gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis contains a group I intron.

The nucleotide sequence of the gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis was determined. It revealed the presence of a group I intron with a length of 411 nucleotides. This is the third occurrence of such an intron discovered in a small subunit rRNA gene encoded by a eukaryotic nuclear genome. The other two occurrences are in Pneumocystis carinii, a fungus of uncertain taxonomic status, and Ankistrodesmus stipitatus, a green alga. The nucleotides of the conserved core structure of 101 group I intron sequences present in different genes and genome types were aligned and their evolutionary relatedness was examined. This revealed a cluster including all group I introns hitherto found in eukaryotic nuclear genes coding for small and large subunit rRNAs. A secondary structure model was designed for the area of the Ustilago maydis small ribosomal subunit RNA precursor where the intron is situated. It shows that the internal guide sequence pairing with the intron boundaries fits between two helices of the small subunit rRNA, and that minimal rearrangement of base pairs suffices to achieve the definitive secondary structure of the 18S rRNA upon splicing.     

The first complete monogenean ribosomal RNA gene operon: sequence and secondary structure of the Gyrodactylus salaris Malmberg, 1957, large subunit ribosomal RNA gene

The sequence of the Gyrodactylus salaris Malmberg, 1957, large subunit, or 28S, ribosomal RNA (rRNA) gene has been determined. This gene is the final portion of the Gyrodactylus rRNA gene operon to be sequenced and results in the first complete sequence of all rRNA genes and spacers from a monogenean. The nucleotide sequence was used to predict the secondary structure of the large subunit rRNA, and regions of conserved and variable sequence and structure were identified. The site where the 5' terminus of the 5.8S rRNA binds to a region within the large subunit rRNA was predicted and complements the anticipated interaction of the 3' terminus of the 5.8S with the 5' terminus of the large subunit rRNA. The large subunit gene may be useful in phylogenetic analysis of the Monogenea or Platyhelminthes and comparisons with other eukaryotes. The variable domains C and H may be most suitable for this purpose.     

The nucleotide sequence of the 17S ribosomal RNA gene of Tetrahymena thermophila and the identification of point mutations resulting in resistance to the antibiotics paromomycin and hygromycin

We have determined the complete sequence of the nuclear gene encoding the small subunit (17 S) rRNA of the ciliated protozoan Tetrahymena thermophila. The gene encodes an RNA molecule which is 1753 nucleotides in length. The sequence of the Tetrahymena small subunit rRNA is homologous to those of other eukaryotes, and the predicted secondary structure for the molecule includes features which are characteristic of eukaryotic small subunit rRNAs. We have also determined the nature of two different mutations in the Tetrahymena 17 S gene which result in resistance to the aminoglycoside antibiotics paromomycin and hygromycin. In each case we have identified a single base change near the 3' end of the rRNA, within a region that is highly evolutionarily conserved in both sequence and secondary structure. Analysis of the effects of these mutations on rRNA structure, and of the impact of these drugs on translation, should help to elucidate the role of the small subunit ribosomal RNA in ribosome function.     

Structural-functional topography of human ribosomes based on the data from crosslinking with mRNA analogs--oligoribonucleotide derivatives

Reviewed are data on the position of template codons with respect to 18S rRNA and certain proteins on human ribosome obtained using a set of mRNA analogs, oligoribonucleotide derivatives carrying alkylating or photoactivatable group at different positions. A comparison of data on template position on the human and Escherichia coli ribosomes has revealed both the similarity in the structure of the mRNA-binding site of bacterial and mammalian ribosomes and the peculiarities of the functioning of mammalian (in particular, human) ribosomes. The similarity manifests itself in that the template codons at the A, P, and E sites of bacterial and human ribosomes are surrounded by similar nucleotides (occupying similar positions in the conserved regions of secondary structure) of small subunit rRNA. The template forms a loop whose foot is in proximity to the 530 stem-loop conserved region of rRNA. The specific features of mammalian ribosomes appear to be associated with their lower conformational mobility as compared with bacterial ribosomes, owing to which their interaction with the template involves a lesser number of molecular contacts.     

Structural and Functional Topography of Human Ribosomes Deduced from Crosslinking with mRNA Analogs,Oligoribonucleotide Derivatives

Reviewed are data on the position of template codons with respect to 18S rRNA and certain proteins on human ribosome obtained using a set of mRNA analogs, oligoribonucleotide derivatives carrying alkylating or photoactivatable groups at different positions. A comparison of data on the template position on the human and Escherichia coliribosomes has revealed both the similarity in the structure of the mRNA-binding site of bacterial and mammalian ribosomes and the peculiarities of the functioning of mammalian (in particular, human) ribosomes. The similarity manifests itself in that the template codons at the A-, P-, and E-sites of bacterial and human ribosomes are surrounded by similar nucleotides (occupying similar positions in the conserved regions of secondary structure) of small subunit rRNA. The template forms a loop whose foot is in proximity to the 530 stem–loop conserved region of rRNA. The specific features of mammalian ribosomes appear to be associated with their lower conformational mobility as compared with bacterial ribosomes, owing to which their interaction with the template involves a lesser number of molecular contacts.     

The ribosomal RNA genes of Drosophila mitochondrial DNA.

The nucleotide sequence of a segment of the mtDNA molecule of Drosophila yakuba which contains the A+T-rich region and the small and large rRNA genes separated by the tRNAval gene has been determined. The 5' end of the small rRNA gene was located by S1 protection analysis. In contrast to mammalian mtDNA, a tRNA gene was not found at the 5' end of the D. yakuba small rRNA gene. The small and large rRNA genes are 20.7% and 16.7% G+C and contain only 789 and 1326 nucleotides. The 5' regions of the small rRNA gene (371 nucleotides) and of the large rRNA gene (643 nucleotides) are extremely low in G+C (14.6% and 9.5%, respectively) and convincing sequence homologies between these regions and the corresponding regions of mouse mt-rRNA genes were found only for a few short segments. Nevertheless, the entire lengths of both of the D. yakuba mt-rRNA genes can be folded into secondary structures which are remarkably similar to secondary structures proposed for the rRNAs of mouse mtDNA. The replication origin-containing, A+T-rich region (1077 nucleotides; 92.8% A+T), which lies between the tRNAile gene and the small rRNA gene, lacks open reading frames greater than 123 nucleotides.     

A case study to estimate the applicability of secondary structures of SSU‐rRNA gene in taxonomy and phylogenetic analyses of ciliates

Phylogenetic studies of ciliates are mainly based on the primary structure information of the nuclear genes. Some regions of the small subunit ribosomal RNA (SSU‐rRNA) gene have distinctive secondary structures, which have demonstrated value as phylogenetic/taxonomic characters. In the current work, we predict the secondary structures of four variable regions (V2, V4, V7 and V9) in the SSU‐rRNA gene of 45 urostylids. Structure comparisons indicate that the V4 region is the most effective in revealing interspecific relationships, while the V9 region appears suitable at the family level or higher. The V2 region also offers some taxonomic information, but is too conserved to reflect phylogenetic relationships at the family or lower level, at least for urostylids. The V7 region is the least informative. We constructed several phylogenetic trees, based on the primary sequence alignment and based on an improved alignment according to the secondary structures. The results suggest that including secondary structure information in phylogenetic analyses provides additional insights into phylogenetic relationships. Using urostylid ciliates as an example, we show that secondary structure information results in a better understanding of their relationships, for example generic relationships within the family Pseudokeronopsidae.     

Reverse splicing of the Tetrahymena IVS: evidence for multiple reaction sites in the 23S rRNA.

Group I introns in rRNA genes are clustered in highly conserved regions that include tRNA and mRNA binding sites. This pattern is consistent with insertion of group I introns by direct interaction with exposed regions of rRNA. Integration of the Tetrahymena group I intron (or intervening sequence, IVS) into large subunit rRNA via reverse splicing was investigated using E. coli 23S rRNA as a model substrate. The results show that sequences homologous to the splice junction in Tetrahymena are the preferred site of integration, but that many other sequences in the 23S rRNA provide secondary targets. Like the original splice junction, many new reaction sites are in regions of stable secondary structure. Reaction at the natural splice junction is observed in 50S subunits and to a lesser extent in 70S ribosomes. These results support the feasibility of intron transposition to new sites in rRNA genes via reverse splicing.     

The performance of several multiple-sequence alignment programs in relation to secondary-structure features for an rRNA sequence

The performances of five global multiple-sequence alignment programs (CLUSTAL W, Divide and Conquer, Malign, PileUp, and TreeAlign) were evaluated using part of the animal mitochondrial small subunit (12S) rRNA molecule. Conserved sequence motifs derived from an alignment based on secondary structural information were used to score how well each program aligned a data set of five vertebrate and five invertebrate taxa over a range of parameter values. All of the programs could align the motifs with reasonable accuracy for at least one set of parameter conditions, although if the whole sequence was considered, similarity to the structural alignment was only 25%-34%. Use of small gap costs generally gave more accurate results, although Malign and TreeAlign generated longer alignments when gap costs were low. The programs differed in the consistency of the alignments when gap cost was varied; CLUSTAL W, Divide and Conquer, and TreeAlign were the most accurate and robust, while PileUp performed poorly as gap cost values increased, and the accuracy of Malign fluctuated. Default settings for the programs did not give the best results, and attempting to select similar parameter values in different programs did not always result in more similar alignments. Poor alignment of even well-conserved motifs can occur if these are near sites with insertions or deletions. Since there is no a priori way to determine gap costs and because such costs can vary over the gene, alignment of rRNA sequences, particularly the less well conserved regions, should be treated carefully and aided by secondary structure and conserved motifs. Some motifs are single bases and so are often invisible to alignment programs. Our tests involved the most conserved regions of the 12S rRNA gene, and alignment of less well conserved regions will be more problematical. None of the alignments we examined produced a fully resolved phylogeny for the data set, indicating that this portion of 12S rRNA is insufficient for resolution of distant evolutionary relationships.     

DNA sequence,structure, and phylogenetic relationship of the small subunit rRNA coding region of mitochondrial DNA fromPodospora anserina

Summary DNA sequence analysis and the localization of the 5 and 3 termini by S1 mapping have shown that the mitochondrial (mt) small subunit rRNA coding region fromPodospora anserina is 1980 bp in length. The analogous coding region for mt rRNA is 1962 bp in maize, 1686 bp inSaccharomyces cerevisiae, and 956 bp in mammals, whereas its counterpart inEscherichia coli is 1542 bp. TheP. anserina mt 16S-like rRNA is 400 bases longer than that fromE. coli, but can be folded into a similar secondary structure. The additional bases appear to be clustered at specific locations, including extensions at the 5 and 3 termini. Comparison with secondary structure diagrams of 16S-like RNAs from several organisms allowed us to specify highly conserved and variable regions of this gene. Phylogenetic tree construction indicated that this gene is grouped with other mitochondrial genes, but most closely, as expected, with the fungal mitochondrial genes.     

How slippage-derived sequences are incorporated into rRNA variable-region secondary structure: implications for phylogeny reconstruction

We analyzed the type and frequency of mutational changes in hypervariable rRNA regions, using the highly length-variable region V4 of the small subunit rRNA locus of tiger beetles (Cicindelidae) as an example. Phylogenetic analysis of indels in closely related species showed that (1) most indels are single nucleotides (usually A or T and sometimes G) or di-nucleotides of A and T. These occur at numerous foci, and they exhibit a strong bias for duplication of 5' single and di-nucleotide motifs but not 3' motifs. (2) Insertions/deletions in stem-forming regions affected paired and unpaired bases with about equal frequency but they did not disrupt the secondary structure. (3) Recurring mutations involving short repeats of the same bases caused parallel evolution of similar sequence motifs in the rRNA of different lineages. The observed types of change are consistent with the propostion that slippage is the main mutational mechanism. Slippage-derived sequences tend to be self-complementary, and therefore the stem-loop structure could be self-organizing as a consequence of the underlying mutational mechanism. Thus, the secondary structure in the cicindelid V4 region may be conserved due to the dynamics of the mutational mechanism rather than to functional constraints. These processes may also have a tendency to produce similar primary sequences irrespective of phylogenetic associations. The findings have implications for sequence alignment in phylogenetic analysis and should caution against the use of secondary structure to improve the determination of positional homology in hypervariable regions.     

Sequence and secondary structure variation in the Gyrodactylus (Platyhelminthes: Monogenea) ribosomal RNA gene array

Nucleotide sequences were determined for the rRNA internal transcribed spacers 1 and 2 (ITS1 and 2) and the 5' terminus of the large subunit rRNA in selected Gyrodactylus species. Examination of primary sequence variation and secondary structure models in ITS2 and variable region V4 of the small subunit rRNA revealed that structure was largely conserved despite significant variation in sequence. ITS1 sequences were highly variable, and models of structure were unreliable but, despite this, show some resemblance to structures predicted in Digenea. ITS2 models demonstrated binding of the 3' end of 5.8S rRNA to the 5' end of the large subunit rRNA and enabled the termini of these genes to be defined with greater confidence than previously. The structure model shown here may prove useful in future phylogenetic analyses.     

Cytochrome oxidase subunit III gene in Neurospora crassa mitochondria. Location and sequence

We have located and sequenced the gene for cytochrome oxidase subunit III (CoIII) in Neurospora crassa mitochondria. The CoIII gene is located downstream from the small rRNA gene within a cluster of tRNA genes and is coded by the same strand as the tRNA and the rRNA genes. Like the tRNA and the rRNA genes, the CoIII gene is also flanked by the GC-rich palindromic DNA sequences which are highly conserved in N. crassa mitochondria. The CoIII coding sequence predicts a protein 269 amino acids long including 8 tryptophan residues. All 8 tryptophan residues are coded for by UGA. This supports our previous conclusion based on the anticodon sequence of N. crassa mitochondrial tryptophan tRNA and provides evidence for the notion that use of UGA as a codon for tryptophan rather than chain termination may be a feature common to most mitochondrial protein synthesis systems. The close correspondence between the amino acid composition of N. crassa CoIII and that of the protein predicted by the CoIII gene sequence suggests that unlike in mammalian mitochondria, AUA is a codon for isoleucine and not for methionine in N. crassa mitochondria. The N. crassa CoIII sequence shows strong homologies to the corresponding yeast and human proteins (53 and 47%, respectively). The overall hydrophobic character of the protein is consistent with suggestions that most of CoIII is embedded in the mitochondrial inner membrane.     

白纹佛蝗线粒体全基因组序列

通过长PCR扩增线粒体全基因组进行保守引物步移法结合克隆测序技术,对白纹佛蝗mtDNA 全序列进行了测定和分析.结果表明:白纹佛蝗线粒体基因组全长15 657 bp,包含13 个蛋白编码基因、22个tRNA 基因和2 个rRNA 基因以及1个非编码的控制区域,它们的长度分别是11 202 bp,1 486 bp,2 156 bp 和 728 bp.37个基因的位置与飞蝗的一致,有9对基因间存在41 bp重叠,重叠碱基数在 1～8 bp之间;基因间隔序列共计21处 126 bp,间隔长度从 1～20 bp不等,最大的基因间隔是20 bp,是在tRNALys 和 ATP8 基因之间.还对lrRNA和srRNA二级结构进行了预测,同时也对tRNA反密码子臂的碱基对类型以及不同链上蛋白编码基因的A/T,C/G组成偏向性进行了详细的讨论.     

Characterization of mitochondrial ribosomal RNA genes in gadiformes: sequence variations,secondary structural features,and phylogenetic implications

Secondary structure features of mitochondrial ribosomal RNAs (mt-rRNAs) of bony fishes were investigated by a DNA sequence alignment approach. The small subunit (SSU) and large subunit (LSU) mt-rRNA genes were found to contain several additional variable regions compared to their mammalian counterparts. Fish mt-LSU rRNA genes were found to be longer than the mammalians due to increased length of some of the variable regions. The 5' and 3' ends of Atlantic cod mt-rRNAs were precisely mapped. The 3' ends of mt-SSU rRNAs were found to be homogenous and mono-adenylated, whereas that of the mt-LSU rRNAs were heterogenous and oligo-adenylated. The 5' ends of mt-SSU rRNAs appeared to be heterogenous, corresponding to the presumed first and second positions of the gene. Sequences of the central domain and the D-domain of the mt-SSU and mt-LSU rRNA genes, respectively, were determined and characterized for 11 gadiform species (representing the families Gadidae, Lotidae, Ranicipitidae, Merlucciidae, Phycidae, and Macrouridae) and one Lophiidae species. Detailed secondary structure models of the RNA regions are presented for the Atlantic cod (Gadus morhua) and Roundnose grenadier (Coryphaeonides rupestris). Saturation plots revealed that DNA nucleotide positions corresponding to unpaired RNA regions become saturated with transitions at sequence divergence levels about 0.15. Phylogenetic analyses revealed some aspects of gadiform relationships. Gadidae was identified as the most derived of the gadiform families. Lotidae was found to be the family closest related to Gadidae, and Ranicipitidae was also recognized as a derived gadiform taxon.     

Multisequence comparisons in protein coding genes. Search for functional constraints

A very powerful method for detecting functional constraints operative in biological macromolecules is presented. This method entails performing a base permanence analysis of protein coding genes at each codon position simultaneously in different species. It calculates the degree of permanence of subregions of the gene by dividing it into segments, c codons long, counting how many sites remain unchanged in each segment among all species compared. By comparing the base permanence among several sequences with the expectations based on a stochastic evolutionary process, gene regions showing different degrees of conservation can be selected. This means that wherever the permanence deviates significantly from the expected value generated by the simulation, the corresponding regions are considered "constrained" or "hypervariable". The constrained regions are of two types: alpha and beta. The alpha regions result from constraints at the amino acid level, whereas the beta regions are those probably involved in "control" processing. The method has been applied to mitochondrial genes coding for subunit 6 of the ATPase and subunit 1 of the cytochrome oxidase in four mammalian species: human, rat, mouse, and cow. In the two mitochondrial genes a few regions that are highly conserved in all codon positions have been identified. Among these regions a sequence, common to both genes, that is complementary to a strongly conserved region of 12S rRNA has been found. This method can also be of great help in studying molecular evolution mechanisms.