Dataset Information

Efficient algorithms for probing the RNA mutation landscape.

ABSTRACT: The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Although secondary structure predictions from thermodynamics-based algorithms are not as accurate as methods employing comparative genomics, the former methods are the only available tools to investigate novel RNAs, such as the many RNAs of unknown function recently reported by the ENCODE consortium. In this paper, we generalize the McCaskill partition function algorithm to sum over the grand canonical ensemble of all secondary structures of all mutants of the given sequence. Specifically, our new program, RNAmutants, simultaneously computes for each integer k the minimum free energy structure MFE(k) and the partition function Z(k) over all secondary structures of all k-point mutants, even allowing the user to specify certain positions required not to mutate and certain positions required to base-pair or remain unpaired. This technically important extension allows us to study the resilience of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. We have successfully applied RNAmutants to investigate deleterious mutations (mutations that radically modify the secondary structure) in the Hepatitis C virus cis-acting replication element and to evaluate the evolutionary pressure applied on different regions of the HIV trans-activation response element. In particular, we show qualitative agreement between published Hepatitis C and HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work also predicts other deleterious mutations, which could be verified experimentally. Finally, we provide evidence that the 3' UTR of the GB RNA virus C has been optimized to preserve evolutionarily conserved stem regions from a deleterious effect of pointwise mutations. We hope that there will be long-term potential applications of RNAmutants in de novo RNA design and drug design against RNA viruses. This work also suggests potential applications for large-scale exploration of the RNA sequence-structure network. Binary distributions are available at http://RNAmutants.csail.mit.edu/.

SUBMITTER: Waldispuhl J

PROVIDER: S-EPMC2475669 | biostudies-literature | 2008

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Efficient algorithms for probing the RNA mutation landscape.

Waldispühl Jérôme J Devadas Srinivas S Berger Bonnie B Clote Peter P

PLoS computational biology 20080808 8

The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. ...[more]

PMID: 18688270

Similar Datasets

Project description:In species with exalbuminous seeds such as crucifer oilseeds and legumes, the endosperm is eventually consumed and its space occupied by the embryo during seed development. However, the main constituent of the early developing seed is the liquid endosperm, and most of the carbon resources for the ensuing stages of seed development arrive at the embryo through the endosperm. In contrast to the extensive study of species with persistent endosperm, little is known about the global gene expression pattern in the endosperm of exalbuminous seed species. We took a multiparallel approach that combines ESTs, protein profiling and microarray analyses to look into the gene expression landscape in the endosperm of the oilseed crop Brassica napus. An EST collection of over 30,000 entries allowed us to detect close to 10,000 unisequences expressed in the endosperm. A protein profile analysis of more than 800 proteins corroborated several signature pathways uncovered by abundant ESTs. Using microarray analyses, we identified genes that are differentially or highly expressed across all developmental stages. These complementary analyses provided insight on several prominent metabolic pathways in the endosperm. We also discovered that LEC1 was highly expressed in the endosperm and that the regulatory cascade downstream of LEC1 operates in the endosperm. The endosperm EST collection and the microarray dataset provide a basic genomic resource for dissecting metabolic and developmental events important for oilseed improvement. Our findings on the featured metabolic processes and the LEC1 regulatory cascade offer new angles for investigation on the integration of endosperm gene expression with embryo development and storage product deposition in seed development. Keywords: seed development We defined the developmental phase of the endosperm based on the stages of the embedded embryos. Under our growth conditions, the globular-shape embryo stage is reached approximately six days after flowering (DAF), the heart-shape embryo stage is reached after eight days, and by day 14 the seed enters the cotyledon stage. The collected tissues allowed the construction of two cDNA libraries: one derived from endosperm of heart-shape embryo developing stage seeds and the other from endosperm of a mix of globular-shape embryo and cotyledon stage developing seeds. Heart-shaped-embryo endosperm is also ready for protein profiling analysis. For microarray analysis, endosperms were harvested from young seeds at the globular-shape embryo, heart-shape embryo and cotyledon stages, respectively. Two biological duplicates were set up, including four technological repeats in each duplicate. A dye-swap hybridization experiment was performed for each pair of target RNA comparisons. We performed four independent aRNA labelings for each tissue pair. For example, we used globular-shape embryo stage endosperm RNA coupled with Cy5 dye vs. heart-shape embryo stage endosperm RNA coupled with Cy3 dye, and the dye-swap experiment was with globular-shape embryo stage endosperm RNA with Cy3 dye vs. the heart-shape embryo stage endosperm RNA with Cy5 dye. Each tissue pair experiment was done with four repeats, using a total of 24 microarrays for the whole experiment.

Dataset Information

Efficient algorithms for probing the RNA mutation landscape.

Publications

Efficient algorithms for probing the RNA mutation landscape.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets