Dynamic scaling approach to study time series fluctuations

We propose a new approach for properly analyzing stochastic time series by mapping the dynamics of time series fluctuations onto a suitable nonequilibrium surface-growth problem. In this framework, the fluctuation sampling time interval plays the role of time variable, whereas the physical time is treated as the analog of spatial variable. In this way we found that the fluctuations of many real-world time series satisfy the analog of the Family-Viscek dynamic scaling ansatz. This finding permits to use the powerful tools of kinetic roughening theory to classify, model, and forecast the fluctuations of real-world time series.Comment: 25 pages, 7 figures, 1 tabl

Correlation between nucleotide composition and folding energy of coding sequences with special attention to wobble bases

Background: The secondary structure and complexity of mRNA influences its accessibility to regulatory molecules (proteins, micro-RNAs), its stability and its level of expression. The mobile elements of the RNA sequence, the wobble bases, are expected to regulate the formation of structures encompassing coding sequences. Results: The sequence/folding energy (FE) relationship was studied by statistical, bioinformatic methods in 90 CDS containing 26,370 codons. I found that the FE (dG) associated with coding sequences is significant and negative (407 kcal/1000 bases, mean +/- S.E.M.) indicating that these sequences are able to form structures. However, the FE has only a small free component, less than 10% of the total. The contribution of the 1st and 3rd codon bases to the FE is larger than the contribution of the 2nd (central) bases. It is possible to achieve a ~ 4-fold change in FE by altering the wobble bases in synonymous codons. The sequence/FE relationship can be described with a simple algorithm, and the total FE can be predicted solely from the sequence composition of the nucleic acid. The contributions of different synonymous codons to the FE are additive and one codon cannot replace another. The accumulated contributions of synonymous codons of an amino acid to the total folding energy of an mRNA is strongly correlated to the relative amount of that amino acid in the translated protein. Conclusion: Synonymous codons are not interchangable with regard to their role in determining the mRNA FE and the relative amounts of amino acids in the translated protein, even if they are indistinguishable in respect of amino acid coding.Comment: 14 pages including 6 figures and 1 tabl

Anomalous roughness with system size dependent local roughness exponent

We note that in a system far from equilibrium the interface roughening may depend on the system size which plays the role of control parameter. To detect the size effect on the interface roughness, we study the scaling properties of rough interfaces formed in paper combustion experiments. Using paper sheets of different width \lambda L, we found that the turbulent flame fronts display anomalous multi-scaling characterized by non universal global roughness exponent \alpha and the system size dependent spectrum of local roughness exponents,\xi_q, whereas the burning fronts possess conventional multi-affine scaling. The structure factor of turbulent flame fronts also exhibit unconventional scaling dependence on \lambda These results are expected to apply to a broad range of far from equilibrium systems, when the kinetic energy fluctuations exceed a certain critical value.Comment: 33 pages, 16 figure

Short ORF-Dependent Ribosome Shunting Operates in an RNA Picorna-Like Virus and a DNA Pararetrovirus that Cause Rice Tungro Disease

Rice tungro disease is caused by synergistic interaction of an RNA picorna-like virus Rice tungro spherical virus (RTSV) and a DNA pararetrovirus Rice tungro bacilliform virus (RTBV). It is spread by insects owing to an RTSV-encoded transmission factor. RTBV has evolved a ribosome shunt mechanism to initiate translation of its pregenomic RNA having a long and highly structured leader. We found that a long leader of RTSV genomic RNA remarkably resembles the RTBV leader: both contain several short ORFs (sORFs) and potentially fold into a large stem-loop structure with the first sORF terminating in front of the stem basal helix. Using translation assays in rice protoplasts and wheat germ extracts, we show that, like in RTBV, both initiation and proper termination of the first sORF translation in front of the stem are required for shunt-mediated translation of a reporter ORF placed downstream of the RTSV leader. The base pairing that forms the basal helix is required for shunting, but its sequence can be varied. Shunt efficiency in RTSV is lower than in RTBV. But in addition to shunting the RTSV leader sequence allows relatively efficient linear ribosome migration, which also contributes to translation initiation downstream of the leader. We conclude that RTSV and RTBV have developed a similar, sORF-dependent shunt mechanism possibly to adapt to the host translation system and/or coordinate their life cycles. Given that sORF-dependent shunting also operates in a pararetrovirus Cauliflower mosaic virus and likely in other pararetroviruses that possess a conserved shunt configuration in their leaders it is tempting to propose that RTSV may have acquired shunt cis-elements from RTBV during their co-existence

GC content around splice sites affects splicing through pre-mRNA secondary structures

<p>Abstract</p> <p>Background</p> <p>Alternative splicing increases protein diversity by generating multiple transcript isoforms from a single gene through different combinations of exons or through different selections of splice sites. It has been reported that RNA secondary structures are involved in alternative splicing. Here we perform a genomic study of RNA secondary structures around splice sites in humans (<it>Homo sapiens</it>), mice (<it>Mus musculus</it>), fruit flies (<it>Drosophila melanogaster</it>), and nematodes (<it>Caenorhabditis elegans</it>) to further investigate this phenomenon.</p> <p>Results</p> <p>We observe that GC content around splice sites is closely associated with the splice site usage in multiple species. RNA secondary structure is the possible explanation, because the structural stability difference among alternative splice sites, constitutive splice sites, and skipped splice sites can be explained by the GC content difference. Alternative splice sites tend to be GC-enriched and exhibit more stable RNA secondary structures in all of the considered species. In humans and mice, splice sites of first exons and long exons tend to be GC-enriched and hence form more stable structures, indicating the special role of RNA secondary structures in promoter proximal splicing events and the splicing of long exons. In addition, GC-enriched exon-intron junctions tend to be overrepresented in tissue-specific alternative splice sites, indicating the functional consequence of the GC effect. Compared with regions far from splice sites and decoy splice sites, real splice sites are GC-enriched. We also found that the GC-content effect is much stronger than the nucleotide-order effect to form stable secondary structures.</p> <p>Conclusion</p> <p>All of these results indicate that GC content is related to splice site usage and it may mediate the splicing process through RNA secondary structures.</p

Identification and characterization of maize microRNAs involved in the very early stage of seed germination

<p>Abstract</p> <p>Background</p> <p>MicroRNAs (miRNAs) are a new class of endogenous small RNAs that play essential regulatory roles in plant growth, development and stress response. Extensive studies of miRNAs have been performed in model plants such as rice, <it>Arabidopsis thaliana </it>and other plants. However, the number of miRNAs discovered in maize is relatively low and little is known about miRNAs involved in the very early stage during seed germination.</p> <p>Results</p> <p>In this study, a small RNA library from maize seed 24 hours after imbibition was sequenced by the Solexa technology. A total of 11,338,273 reads were obtained. 1,047,447 total reads representing 431 unique sRNAs matched to known maize miRNAs. Further analysis confirmed the authenticity of 115 known miRNAs belonging to 24 miRNA families and the discovery of 167 novel miRNAs in maize. Both the known and the novel miRNAs were confirmed by sequencing of a second small RNA library constructed the same way as the one used in the first sequencing. We also found 10 miRNAs that had not been reported in maize, but had been reported in other plant species. All novel sequences had not been earlier described in other plant species. In addition, seven miRNA* sequences were also obtained. Putative targets for 106 novel miRNAs were successfully predicted. Our results indicated that miRNA-mediated gene expression regulation is present in maize imbibed seed.</p> <p>Conclusions</p> <p>This study led to the confirmation of the authenticity of 115 known miRNAs and the discovery of 167 novel miRNAs in maize. Identification of novel miRNAs resulted in significant enrichment of the repertoire of maize miRNAs and provided insights into miRNA regulation of genes expressed in imbibed seed.</p

Analysis of conserved microRNAs in floral tissues of sexual and apomictic Boechera species

<p>Abstract</p> <p>Background</p> <p>Apomixis or asexual seed formation represents a potentially important agronomic trait whose introduction into crop plants could be an effective way to fix and perpetuate a desirable genotype through successive seed generations. However, the gene regulatory pathways underlying apomixis remain unknown. In particular, the potential function of microRNAs, which are known to play crucial roles in many aspects of plant growth and development, remains to be determined with regards to the switch from sexual to apomictic reproduction.</p> <p>Results</p> <p>Using bioinformatics and microarray validation procedures, 51 miRNA families conserved among angiosperms were identified in <it>Boechera</it>. Microarray assay confirmed 15 of the miRNA families that were identified by bioinformatics techniques. 30 cDNA sequences representing 26 miRNAs could fold back into stable pre-miRNAs. 19 of these pre-miRNAs had miRNAs with <it>Boechera</it>-specific nucleotide substitutions (NSs). Analysis of the Gibbs free energy (ΔG) of these pre-miRNA stem-loops with NSs showed that the <it>Boechera</it>-specific miRNA NSs significantly (p ≤ 0.05) enhance the stability of stem-loops. Furthermore, six transcription factors, the Squamosa promoter binding protein like SPL6, SPL11 and SPL15, Myb domain protein 120 (MYB120), RELATED TO AP2.7 DNA binding (RAP2.7, TOE1 RAP2.7) and TCP family transcription factor 10 (TCP10) were found to be expressed in sexual or apomictic ovules. However, only SPL11 showed differential expression with significant (p ≤ 0.05) up-regulation at the megaspore mother cell (MMC) stage of ovule development in apomictic genotypes.</p> <p>Conclusions</p> <p>This study constitutes the first extensive insight into the conservation and expression of microRNAs in <it>Boechera </it>sexual and apomictic species. The miR156/157 target squamosa promoter binding protein-like 11 (SPL11) was found differentially expressed with significant (p ≤ 0.05) up-regulation at the MMC stage of ovule development in apomictic genotypes. The results also demonstrate that nucleotide changes in mature miRNAs significantly (p ≤ 0.05) enhance the thermodynamic stability of pre-miRNA stem-loops.</p

A Universal Trend of Reduced mRNA Stability near the Translation-Initiation Site in Prokaryotes and Eukaryotes

Recent studies have suggested that the thermodynamic stability of mRNA secondary structure near the start codon can regulate translation efficiency in Escherichia coli, and that translation is more efficient the less stable the secondary structure. We survey the complete genomes of 340 species for signals of reduced mRNA secondary structure near the start codon. Our analysis includes bacteria, archaea, fungi, plants, insects, fishes, birds, and mammals. We find that nearly all species show evidence for reduced mRNA stability near the start codon. The reduction in stability generally increases with increasing genomic GC content. In prokaryotes, the reduction also increases with decreasing optimal growth temperature. Within genomes, there is variation in the stability among genes, and this variation correlates with gene GC content, codon bias, and gene expression level. For birds and mammals, however, we do not find a genome-wide trend of reduced mRNA stability near the start codon. Yet the most GC rich genes in these organisms do show such a signal. We conclude that reduced stability of the mRNA secondary structure near the start codon is a universal feature of all cellular life. We suggest that the origin of this reduction is selection for efficient recognition of the start codon by initiator-tRNA

The Proteomic Code: a molecular recognition code for proteins

<p>Abstract</p> <p>Background</p> <p>The Proteomic Code is a set of rules by which information in genetic material is transferred into the physico-chemical properties of amino acids. It determines how individual amino acids interact with each other during folding and in specific protein-protein interactions. The Proteomic Code is part of the redundant Genetic Code.</p> <p>Review</p> <p>The 25-year-old history of this concept is reviewed from the first independent suggestions by Biro and Mekler, through the works of Blalock, Root-Bernstein, Siemion, Miller and others, followed by the discovery of a Common Periodic Table of Codons and Nucleic Acids in 2003 and culminating in the recent conceptualization of partial complementary coding of interacting amino acids as well as the theory of the nucleic acid-assisted protein folding.</p> <p>Methods and conclusions</p> <p>A novel cloning method for the design and production of specific, high-affinity-reacting proteins (SHARP) is presented. This method is based on the concept of proteomic codes and is suitable for large-scale, industrial production of specifically interacting peptides.</p

The surprising negative correlation of gene length and optimal codon use - disentangling translational selection from GC-biased gene conversion in yeast

<p>Abstract</p> <p>Background</p> <p>Surprisingly, in several multi-cellular eukaryotes optimal codon use correlates negatively with gene length. This contrasts with the expectation under selection for translational accuracy. While suggested explanations focus on variation in strength and efficiency of translational selection, it has rarely been noticed that the negative correlation is reported only in organisms whose optimal codons are biased towards codons that end with G or C (-GC). This raises the question whether forces that affect base composition - such as GC-biased gene conversion - contribute to the negative correlation between optimal codon use and gene length.</p> <p>Results</p> <p>Yeast is a good organism to study this as equal numbers of optimal codons end in -GC and -AT and one may hence compare frequencies of optimal GC- with optimal AT-ending codons to disentangle the forces. Results of this study demonstrate in yeast frequencies of GC-ending (optimal AND non-optimal) codons decrease with gene length and increase with recombination. A decrease of GC-ending codons along genes contributes to the negative correlation with gene length. Correlations with recombination and gene expression differentiate between GC-ending and optimal codons, and also substitution patterns support effects of GC-biased gene conversion.</p> <p>Conclusion</p> <p>While the general effect of GC-biased gene conversion is well known, the negative correlation of optimal codon use with gene length has not been considered in this context before. Initiation of gene conversion events in promoter regions and the presence of a gene conversion gradient most likely explain the observed decrease of GC-ending codons with gene length and gene position.</p