Phylogenetic comparative analysis and the secondary structure of ribonuclease P RNA--a review

The most incisive a priori approach to inferring the higher order structure of large RNAs has proven to be the use of phylogenetic comparisons. This article provides guidelines to the method, using as an illustration the elucidation of the secondary structure of the catalytic RNA subunit of ribonuclease P (RNase P). The resultant structure is compared to the possibilities that are predicted thermodynamically for the RNase P RNA sequences of nine eubacteria.     

Circles: automating the comparative analysis of RNA secondary structure

SUMMARY: Circles is a program for inferring RNA secondary structure using maximum weight matching. The program can read in an alignment in FASTA, ClustalW, or NEXUS format, compute a maximum weight matching, and export one or more secondary structures in various file formats. AVAILABILITY: The program is available at no cost from http://taxonomy.zoology.gla.ac.uk/rod/circles/ and requires Windows 95/98/NT. CONTACT: r.page@bio.gla.ac.uk     

Phylogenetic origins of the plant mitochondrion based on a comparative analysis of 5S ribosomal RNA sequences

Summary The complete nucleotide sequences of 5S ribosomal RNAs fromRhodocyclus gelatinosa, Rhodobacter sphaeroides, andPseudomonas cepacia were determined. Comparisons of these 5S RNA sequences show that rather than being phylogenetically related to one another, the two photosynthetic bacterial 5S RNAs share more sequence and signature homology with the RNAs of two nonphotosynthetic strains.Rhodobacter sphaeroides is specifically related toParacoccus denitrificans andRc. gelatinosa is related toPs. cepacia.These results support earlier 16S ribosomal RNA studies and add two important groups to the 5S RNA data base. Unique 5S RNA structural features previously found inP. denitrificans are present also in the 5S RNA ofRb. sphaeroides; these provide the basis for subdivisional signatures. The immediate consequence of our obtaining these new sequences is that we are able to clarify the phylogenetic origins of the plant mitochondrion. In particular, we find a close phylogenetic relationship between the plant mitochondria and members of the alpha subdivision of the purple photosynthetic bacteria, namely,Rb. sphaeroides, P. denitrificans, andRhodospirillum rubrum.     

Phylogenetic comparative chemical footprint analysis of the interaction between ribonuclease P RNA and tRNA.

Ribonuclease P RNA is the catalytic moiety of the ribonucleoprotein enzyme that endonucleolytically cleaves precursor sequences from the 5' ends of pre-tRNAs. The bacterial RNase P RNA-tRNA complex was examined with a footprinting approach, utilizing chemical modification to determine RNase P RNA nucleotides that potentially contact tRNA. RNase P RNA was modified with dimethylsulfate or kethoxal in the presence or absence of tRNA, and sites of modification were detected by primer extension. Comparison of the results reveals RNase P bases that are protected from modification upon binding tRNA. Analyses were carried out with RNase P RNAs from three different bacteria: Escherichia coli, Chromatium vinosum and Bacillus subtilis. Discrete bases of these RNAs that lie within conserved, homologous portions of the secondary structures are similarly protected. One protection among all three RNAs was attributed to the precursor segment of pre-tRNA. Experiments using pre-tRNAs containing precursor segments of variable length demonstrate that a precursor segment of only 2-4 nucleotides is sufficient to confer this protection. Deletion of the 3'-terminal CCA sequence of tRNA correlates with loss of protection of a particular loop in the RNase P RNA secondary structure. Analysis of mutant tRNAs containing sequential 3'-terminal deletions suggests a relative orientation of the bound tRNA CCA to that loop.     

New insight into RNase P RNA structure from comparative analysis of the archaeal RNA

A detailed comparative analysis of archaeal RNase P RNA structure and a comparison of the resulting structural information with that of the bacterial RNA reveals that the archaeal RNase P RNAs are strikingly similar to those of Bacteria. The differences between the secondary structure models of archaeal and bacterial RNase P RNA have largely disappeared, and even variation in the sequence and structure of the RNAs are similar in extent and type. The structure of the cruciform (P7-11) has been reevaluated on the basis of a total of 321 bacterial and archaeal sequences, leading to a model for the structure of this region of the RNA that includes an extension to P11 that consistently organizes the cruciform and adjacent highly-conserved sequences.     

Phylogenetic comparative mutational analysis of the base-pairing between RNase P RNA and its substrate.

We have studied the base-pairing between the 3'-terminal CCA motif of a tRNA precursor and RNase P RNA by a phylogenetic mutational comparative approach. Thus, various derivatives of the Escherichia coli tRNA(Ser)Su1 precursor harboring all possible substitutions at either the first or the second C of the 3'-terminal CCA motif were generated. Cleavage site selection on these precursors was studied using mutant variants of M1 RNA, the catalytic subunit of E. coli RNase P, carrying changes at positions 292 or 293, which are involved in the interaction with the 3'-terminal CCA motif. From our data we conclude that these two C's in the substrate interact with the well-conserved G292 and G293 through canonical Watson-Crick base-pairing. Cleavage performed using reconstituted holoenzyme complexes suggests that this interaction also occurs in the presence of the C5 protein. Furthermore, we studied the interaction using various derivatives of RNase P RNAs from Mycoplasma hyopneumoniae and Mycobacterium tuberculosis. Our results suggest that the base-pairing between the 3'-terminal CCA motif and RNase P is present also in other bacterial RNase P-substrate complexes and is not limited to a particular bacterial species.     

Phylogenetic analysis and secondary structure of the Bacillus subtilis bacteriophage RNA required for DNA packaging

An unusual RNA molecule encoded by the Bacillus subtilis bacteriophage phi 29 is a structural component of the viral prohead and is required for the ATP-dependent packaging of DNA. Here we report a model of secondary structure for this prohead RNA developed from a phylogenetic analysis of the primary sequences of prohead RNAs of related phages. Twenty-nine phages related to phi 29 were found to produce prohead RNAs. These RNAs were analyzed by their ability to replace phi 29 RNA in in vitro phage assembly, by Northern blot hybridization with a probe complementary to phi 29 RNA, and by partial and complete sequence analyses. These analyses revealed four quite different sequences ranging in length from 161 to 174 residues. The secondary structure deduced from these sequences, in agreement with earlier observations, indicated that prohead RNA is organized into two domains. The larger 5'-domain (Domain I) is composed of 113-117 residues and contains four helices. Three of these helices appear to be organized into a central stem that is interrupted by two unpaired loops and the fourth helix and loop. The smaller 3'-domain (Domain II) is composed of 40-44 residues and consists of two helices. Domains I and II are separated by 8-13 unpaired residues. Nuclease cleavage occurs readily in this single-stranded joining region, and this cleavage allows the subsequent separation of the two RNA domains. The separated Domain I is fully active in DNA packaging in vitro. The functional significance and biological role of Domain II are unknown. The phylogenetic secondary structure model provides a basis for further analysis of the role of this RNA in bacteriophage morphogenesis.     

Phylogenetic relationships of ascomycetes and basidiomycetes based on comparative genomics analysis

The relationship between ascomycetes and basidiomycetes, the two main phyla of non-flagellated fungi, has rarely been investigated. In this study, we performed a comparative genomics analysis of genome sequences of 55 ascomycetes and 26 basidiomycetes species and detected 81 universal markers, 875 homologous genes and a conserved contig in the glucose-regulated protein gene. In dendrograms based on simple sequence repeat markers and homologous genes, ascomycetes and basidiomycetes formed distinct clusters, with each set of taxa having a high coefficient of relatedness. Ascomycetes and basidiomycetes also constituted distinct groups in a phylogenetic tree based on a conserved contig in the glucose-regulated protein gene. These results provide evidence that basidiomycetes may be derived from ascomycetes but are definitely genetically differentiated at the genomic level. The phylogenetic relationships of ascomycetes and basidiomycetes uncovered in this study provide new insights for future research related to fungal classification and evolution.     

Phylogenetic comparative assembly

Background  

Recent high throughput sequencing technologies are capable of generating a huge amount of data for bacterial genome sequencing projects. Although current sequence assemblers successfully merge the overlapping reads, often several contigs remain which cannot be assembled any further. It is still costly and time consuming to close all the gaps in order to acquire the whole genomic sequence.     

Phylogenetic comparative assembly

Background

Biclustering is an important analysis procedure to understand the biological mechanisms from microarray gene expression data. Several algorithms have been proposed to identify biclusters, but very little effort was made to compare the performance of different algorithms on real datasets and combine the resultant biclusters into one unified ranking.

Results

In this paper we propose differential co-expression framework and a differential co-expression scoring function to objectively quantify quality or goodness of a bicluster of genes based on the observation that genes in a bicluster are co-expressed in the conditions belonged to the bicluster and not co-expressed in the other conditions. Furthermore, we propose a scoring function to stratify biclusters into three types of co-expression. We used the proposed scoring functions to understand the performance and behavior of the four well established biclustering algorithms on six real datasets from different domains by combining their output into one unified ranking.

Conclusions

Differential co-expression framework is useful to provide quantitative and objective assessment of the goodness of biclusters of co-expressed genes and performance of biclustering algorithms in identifying co-expression biclusters. It also helps to combine the biclusters output by different algorithms into one unified ranking i.e. meta-biclustering.     

Phylogenetic analysis of tmRNA secondary structure.

The bacterial tmRNA acts with dual tRNA-like and mRNA-like character to tag incomplete translation products for degradation. Comparative analysis of 17 tmRNA genes (including eight new sequences) has allowed us to deduce conserved features of the tmRNA secondary structure. Except in a segment that includes the first codon of the tag reading frame, tmRNA is highly structured, with four pseudoknots and a total of 11 conserved base pairing regions. The previously identified tRNA minihelix structure is connected by a long base paired region to a large structured domain composed of a pseudoknot, followed by the tag reading frame and a string of three rather similar pseudoknots. The conservation of numerous structural elements among diverse eubacterial species indicates that these elements have important function beyond simply forming an endonuclease-resistant link between the reading frame and the tRNA-like domain.     

The accuracy of ribosomal RNA comparative structure models

The determination of the 16S and 23S rRNA secondary structure models was initiated shortly after the first complete 16S and 23S rRNA sequences were determined in the late 1970s. The structures that are common to all 16S rRNAs and all 23S rRNAs were determined using comparative methods from the analysis of thousands of rRNA sequences. Twenty-plus years later, the 16S and 23S rRNA comparative structure models have been evaluated against the recently determined high-resolution crystal structures of the 30S and 50S ribosomal subunits. Nearly all of the predicted covariation-based base pairs, including the regular base pairs and helices, and the irregular base pairs and tertiary interactions, were present in the 30S and 50S crystal structures.     

Macintosh sequence analysis software

The analysis of information in nucleotide and amino acid sequence data from an investigator’s own laboratory, or from the ever-growing worldwide databases, is critically dependent on well planned and written software. Although the most powerful packages previously have been confined to workstations, there has been a dramatic increase over the last few years in the sophistication of the programs available for personal computers, as the speed and power of these have increased. A wide choice of software is available for the Macintosh, including the LaserGene suite of programs from DNAStar. This review assesses the strengths and weaknesses of LaserGene and concludes that it provides a useful and comprehensive range of sequence analysis tools.     

Phylogenetic comparative analysis of microsporogenesis in angiosperms with a focus on monocots

This paper presents the first broad overview of three main features of microsporogenesis (male meiosis) in angiosperms: cytokinesis (cell division), intersporal wall formation, and tetrad form. A phylogenetic comparative approach was used to test for correlated evolution among these characters and to make hypotheses about evolutionary trends in microsporogenesis. The link between features of microsporogenesis and pollen aperture type was examined. We show that the pathway associated with successive cytokinesis (cytoplasm is partitioned after each meiotic division) is restricted to wall formation mediated by centrifugally developing cell plates, and tetragonal (or decussate, T-shaped, linear) tetrads. Conversely, much more flexibility is observed when cytokinesis is simultaneous (two meiotic divisions completed before cytoplasmic partitioning). We suggest that the ancestral type of microsporogenesis for angiosperms, and perhaps for all seed plants, associated simultaneous cytokinesis with centripetal wall formation, resulting in a large diversity in tetrad forms, ranging from regular tetrahedral to tetragonal tetrads, including rhomboidal tetrads. From this ancestral pathway, switches toward successive cytokinesis occurred among basal angiosperms and monocots, generally associated with a switch toward centrifugal intersporal wall formation, whereas eudicots evolved toward an almost exclusive production of regular tetrahedral tetrads. No straightforward link is found between the type of microsporogenesis and pollen aperture type.     

A comprehensive comparison of comparative RNA structure prediction approaches

Background  

An increasing number of researchers have released novel RNA structure analysis and prediction algorithms for comparative approaches to structure prediction. Yet, independent benchmarking of these algorithms is rarely performed as is now common practice for protein-folding, gene-finding and multiple-sequence-alignment algorithms.     

'DNA Strider': a 'C' program for the fast analysis of DNA and protein sequences on the Apple Macintosh family of computers

DNA Strider is a new integrated DNA and Protein sequence analysis program written with the C language for the Macintosh Plus, SE and II computers. It has been designed as an easy to learn and use program as well as a fast and efficient tool for the day-to-day sequence analysis work. The program consists of a multi-window sequence editor and of various DNA and Protein analysis functions. The editor may use 4 different types of sequences (DNA, degenerate DNA, RNA and one-letter coded protein) and can handle simultaneously 6 sequences of any type up to 32.5 kB each. Negative numbering of the bases is allowed for DNA sequences. All classical restriction and translation analysis functions are present and can be performed in any order on any open sequence or part of a sequence. The main feature of the program is that the same analysis function can be repeated several times on different sequences, thus generating multiple windows on the screen. Many graphic capabilities have been incorporated such as graphic restriction map, hydrophobicity profile and the CAI plot- codon adaptation index according to Sharp and Li. The restriction sites search uses a newly designed fast hexamer look-ahead algorithm. Typical runtime for the search of all sites with a library of 130 restriction endonucleases is 1 second per 10,000 bases. The circular graphic restriction map of the pBR322 plasmid can be therefore computed from its sequence and displayed on the Macintosh Plus screen within 2 seconds and its multiline restriction map obtained in a scrolling window within 5 seconds.     

A comparative analysis of the triloops in all high-resolution RNA structures reveals sequence structure relationships

Despite an increasing number of experimentally determined RNA structures, the gap between the number of structures and that of RNA families is still growing. To overcome this limitation, efficient and reliable RNA modeling methodologies must be developed. In order to reach this goal, here, we show how triloop sequence-structure relationships have been inferred through a systematic analysis of all triloops found in available high-resolution structures. The structural annotation of all triloops allowed us to define discrete states of the triloop's conformational space, and therefore an explicit sequence-to-structure relation. The sequence-structure relationships inferred from this explicit relation are presented in a convenient modeling table that provides a limited set of possible three-dimensional structures given any triloop sequence. The table is indexed by the two nucleotides that form the triloop's flanking base pair, since they are shown to provide the most information about the triloop three-dimensional structures. We also report the observations in the X-ray crystallographic structures of important conformational variations, which we believe might be the result of RNA dynamic.