Phylogenetic analysis of the structure of RNase MRP RNA in yeasts

RNase MRP is a ribonucleoprotein enzyme involved in processing precursor rRNA in eukaryotes. To facilitate our structure-function analysis of RNase MRP from Saccharomyces cerevisiae, we have determined the likely secondary structure of the RNA component by a phylogenetic approach in which we sequenced all or part of the RNase MRP RNAs from 17 additional species of the Saccharomycetaceae family. The structure deduced from these sequences contains the helices previously suggested to be common to the RNA subunit of RNase MRP and the related RNA subunit of RNase P, an enzyme cleaving tRNA precursors. However, outside this common region, the structure of RNase MRP RNA determined here differs from a previously proposed universal structure for RNase MRPs. Chemical and enzymatic structure probing analyses were consistent with our revised secondary structure. Comparison of all known RNase MRP RNA sequences revealed three regions with highly conserved nucleotides. Two of these regions are part of a helix implicated in RNA catalysis in RNase P, suggesting that RNase MRP may cleave rRNA using a similar catalytic mechanism.     

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.     

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.     

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 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 methods complement discriminant function analysis in ecomorphology

In ecomorphology, Discriminant Function Analysis (DFA) has been used as evidence for the presence of functional links between morphometric variables and ecological categories. Here we conduct simulations of characters containing phylogenetic signal to explore the performance of DFA under a variety of conditions. Characters were simulated using a phylogeny of extant antelope species from known habitats. Characters were modeled with no biomechanical relationship to the habitat category; the only sources of variation were body mass, phylogenetic signal, or random “noise.” DFA on the discriminability of habitat categories was performed using subsets of the simulated characters, and Phylogenetic Generalized Least Squares (PGLS) was performed for each character. Analyses were repeated with randomized habitat assignments. When simulated characters lacked phylogenetic signal and/or habitat assignments were random, <5.6% of DFAs and <8.26% of PGLS analyses were significant. When characters contained phylogenetic signal and actual habitats were used, 33.27 to 45.07% of DFAs and <13.09% of PGLS analyses were significant. False Discovery Rate (FDR) corrections for multiple PGLS analyses reduced the rate of significance to <4.64%. In all cases using actual habitats and characters with phylogenetic signal, correct classification rates of DFAs exceeded random chance. In simulations involving phylogenetic signal in both predictor variables and predicted categories, PGLS with FDR was rarely significant, while DFA often was. In short, DFA offered no indication that differences between categories might be explained by phylogenetic signal, while PGLS did. As such, PGLS provides a valuable tool for testing the functional hypotheses at the heart of ecomorphology. Am J Phys Anthropol 153:663–674, 2014. © 2013 Wiley Periodicals, Inc.     

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.     

RNA folding using quantum computers

The 3-dimensional fold of an RNA molecule is largely determined by patterns of intramolecular hydrogen bonds between bases. Predicting the base pairing network from the sequence, also referred to as RNA secondary structure prediction or RNA folding, is a nondeterministic polynomial-time (NP)-complete computational problem. The structure of the molecule is strongly predictive of its functions and biochemical properties, and therefore the ability to accurately predict the structure is a crucial tool for biochemists. Many methods have been proposed to efficiently sample possible secondary structure patterns. Classic approaches employ dynamic programming, and recent studies have explored approaches inspired by evolutionary and machine learning algorithms. This work demonstrates leveraging quantum computing hardware to predict the secondary structure of RNA. A Hamiltonian written in the form of a Binary Quadratic Model (BQM) is derived to drive the system toward maximizing the number of consecutive base pairs while jointly maximizing the average length of the stems. A Quantum Annealer (QA) is compared to a Replica Exchange Monte Carlo (REMC) algorithm programmed with the same objective function, with the QA being shown to be highly competitive at rapidly identifying low energy solutions. The method proposed in this study was compared to three algorithms from literature and, despite its simplicity, was found to be competitive on a test set containing known structures with pseudoknots.