Efficient algorithms for reconstructing gene content by co-evolution

<p>Abstract</p> <p>Background</p> <p>In a previous study we demonstrated that co-evolutionary information can be utilized for improving the accuracy of ancestral gene content reconstruction. To this end, we defined a new computational problem, the Ancestral Co-Evolutionary (ACE) problem, and developed algorithms for solving it.</p> <p>Results</p> <p>In the current paper we generalize our previous study in various ways. First, we describe new efficient computational approaches for solving the ACE problem. The new approaches are based on reductions to classical methods such as linear programming relaxation, quadratic programming, and min-cut. Second, we report new computational hardness results related to the ACE, including practical cases where it can be solved in polynomial time.</p> <p>Third, we generalize the ACE problem and demonstrate how our approach can be used for inferring parts of the genomes of <it>non-ancestral</it> organisms. To this end, we describe a heuristic for finding the portion of the genome ('dominant set’) that can be used to reconstruct the rest of the genome with the lowest error rate. This heuristic utilizes both evolutionary information and co-evolutionary information.</p> <p>We implemented these algorithms on a large input of the ACE problem (95 unicellular organisms, 4,873 protein families, and 10, 576 of co-evolutionary relations), demonstrating that some of these algorithms can outperform the algorithm used in our previous study. In addition, we show that based on our approach a ’dominant set’ cab be used reconstruct a major fraction of a genome (up to 79%) with relatively low error-rate (<it>e.g.</it> 0.11). We find that the ’dominant set’ tends to include metabolic and regulatory genes, with high evolutionary rate, and low protein abundance and number of protein-protein interactions.</p> <p>Conclusions</p> <p>The <it>ACE</it> problem can be efficiently extended for inferring the genomes of organisms that exist today. In addition, it may be solved in polynomial time in many practical cases. Metabolic and regulatory genes were found to be the most important groups of genes necessary for reconstructing gene content of an organism based on other related genomes.</p

Chemomechanics of ionically conductive ceramics for electrical energy conversion and storage

Functional materials for energy conversion and storage exhibit strong coupling between electrochemistry and mechanics. For example, ceramics developed as electrodes for both solid oxide fuel cells and batteries exhibit cyclic volumetric expansion upon reversible ion transport. Such chemomechanical coupling is typically far from thermodynamic equilibrium, and thus is challenging to quantify experimentally and computationally. In situ measurements and atomistic simulations are under rapid development to explore how this coupling can be used to potentially improve both device performance and durability. Here, we review the commonalities of coupling between electrochemical and mechanical states in fuel cell and battery materials, illustrating with specific cases the progress in materials processing, in situ characterization, and computational modeling and simulation. We also highlight outstanding questions and opportunities in these applications – both to better understand the limiting mechanisms within the materials and to significantly advance the durability and predictability of device performance required for renewable energy conversion and storage.United States. Dept. of Energy (Basic Energy Sciences Division of Materials Sciences and Engineering, grant DE-SC0002633)United States. Dept. of Energy (Office of Science, Graduate Fellowship Program (DOE SCGF))United States. American Recovery and Reinvestment Act of 2009 (ORISE-ORAU, contract no. DE-AC05-06OR23100))United States. Dept. of Energy. Division of Materials Sciences and Engineering (MIT/DMSE Salapatas Fellowship)United States. Air Force Office of Scientific Research (Presidential Early Career Award in Science and Engineering (PECASE)

Small Polarons in Transition Metal Oxides

The formation of polarons is a pervasive phenomenon in transition metal oxide compounds, with a strong impact on the physical properties and functionalities of the hosting materials. In its original formulation the polaron problem considers a single charge carrier in a polar crystal interacting with its surrounding lattice. Depending on the spatial extension of the polaron quasiparticle, originating from the coupling between the excess charge and the phonon field, one speaks of small or large polarons. This chapter discusses the modeling of small polarons in real materials, with a particular focus on the archetypal polaron material TiO2. After an introductory part, surveying the fundamental theoretical and experimental aspects of the physics of polarons, the chapter examines how to model small polarons using first principles schemes in order to predict, understand and interpret a variety of polaron properties in bulk phases and surfaces. Following the spirit of this handbook, different types of computational procedures and prescriptions are presented with specific instructions on the setup required to model polaron effects.Comment: 36 pages, 12 figure

Reduced stability of mRNA secondary structure near the translation-initiation site in dsDNA viruses

<p>Abstract</p> <p>Background</p> <p>Recent studies have demonstrated a selection pressure for reduced mRNA secondary-structure stability near the start codon of coding sequences. This selection pressure can be observed in bacteria, archaea, and eukaryotes, and is likely caused by the requirement of efficient translation initiation in cellular organism.</p> <p>Results</p> <p>Here, we surveyed the complete genomes of 650 dsDNA virus strains for signals of reduced stability of mRNA secondary structure near the start codon. Our analysis included viruses infecting eukaryotic, prokaryotic, and archaeic hosts. We found that many viruses showed evidence for reduced mRNA secondary-structure stability near the start codon. The effect was most pronounced in viruses infecting prokaryotes, but was also observed in viruses infecting eukaryotes and archaea. The reduction in stability generally increased with increasing genomic GC content. For bacteriophage, the reduction was correlated with a corresponding reduction of stability in the phage hosts.</p> <p>Conclusions</p> <p>We conclude that reduced stability of the mRNA secondary structure near the start codon is a common feature for dsDNA viruses, likely driven by the same selective pressures that cause it in cellular organisms.</p

Synonymous Codon Ordering: A Subtle but Prevalent Strategy of Bacteria to Improve Translational Efficiency

Background: In yeast coding sequences, once a particular codon has been used, subsequent occurrence of the same amino acid tends to use codons sharing the same tRNA. Such a phenomenon of co-tRNA codons pairing bias (CTCPB) is also found in some other eukaryotes but it is not known whether it occurs in prokaryotes. Methodology/Principal Findings: In this study, we focused on a total of 773 bacterial genomes to investigate their synonymous codon pairing preferences. After calculating the actual frequencies of synonymous codon pairs and comparing them with their expected values, we detected an obvious pairing bias towards identical codon pairs. This seems consistent with the previously reported CTCPB phenomenon, since identical codons are certainly read by the same tRNA. However, among co-tRNA but non-identical codon pairs, only 22 were often found overrepresented, suggesting that many co-tRNA codons actually do not preferentially pair together in prokaryotes. Therefore, the previously reported co-tRNA codons pairing rule needs to be more rigorously defined. The affinity differences between a tRNA anticodon and its readable codons should be taken into account. Moreover, both within-gene-shuffling tests and phylogenetic analyses support the idea that translational selection played an important role in shaping the observed synonymous codon pairing pattern in prokaryotes. Conclusions: Overall, a high level of synonymous codon pairing bias was detected in 73 % investigated bacterial species

Ground, Proximal, and Satellite Remote Sensing of Soil Moisture

Soil moisture (SM) is a key hydrologic state variable that is of significant importance for numerous Earth and environmental science applications that directly impact the global environment and human society. Potential applications include, but are not limited to, forecasting of weather and climate variability; prediction and monitoring of drought conditions; management and allocation of water resources; agricultural plant production and alleviation of famine; prevention of natural disasters such as wild fires, landslides, floods, and dust storms; or monitoring of ecosystem response to climate change. Because of the importance and wide‐ranging applicability of highly variable spatial and temporal SM information that links the water, energy, and carbon cycles, significant efforts and resources have been devoted in recent years to advance SM measurement and monitoring capabilities from the point to the global scales. This review encompasses recent advances and the state‐of‐the‐art of ground, proximal, and novel SM remote sensing techniques at various spatial and temporal scales and identifies critical future research needs and directions to further advance and optimize technology, analysis and retrieval methods, and the application of SM information to improve the understanding of critical zone moisture dynamics. Despite the impressive progress over the last decade, there are still many opportunities and needs to, for example, improve SM retrieval from remotely sensed optical, thermal, and microwave data and opportunities for novel applications of SM information for water resources management, sustainable environmental development, and food security

Sequestration of Highly Expressed mRNAs in Cytoplasmic Granules, P-Bodies, and Stress Granules Enhances Cell Viability

Transcriptome analyses indicate that a core 10%–15% of the yeast genome is modulated by a variety of different stresses. However, not all the induced genes undergo translation, and null mutants of many induced genes do not show elevated sensitivity to the particular stress. Elucidation of the RNA lifecycle reveals accumulation of non-translating mRNAs in cytoplasmic granules, P-bodies, and stress granules for future regulation. P-bodies contain enzymes for mRNA degradation; under stress conditions mRNAs may be transferred to stress granules for storage and return to translation. Protein degradation by the ubiquitin-proteasome system is elevated by stress; and here we analyzed the steady state levels, decay, and subcellular localization of the mRNA of the gene encoding the F-box protein, UFO1, that is induced by stress. Using the MS2L mRNA reporter system UFO1 mRNA was observed in granules that colocalized with P-bodies and stress granules. These P-bodies stored diverse mRNAs. Granules of two mRNAs transported prior to translation, ASH1-MS2L and OXA1-MS2L, docked with P-bodies. HSP12 mRNA that gave rise to highly elevated protein levels was not observed in granules under these stress conditions. ecd3, pat1 double mutants that are defective in P-body formation were sensitive to mRNAs expressed ectopically from strong promoters. These highly expressed mRNAs showed elevated translation compared with wild-type cells, and the viability of the mutants was strongly reduced. ecd3, pat1 mutants also exhibited increased sensitivity to different stresses. Our interpretation is that sequestration of highly expressed mRNAs in P-bodies is essential for viability. Storage of mRNAs for future regulation may contribute to the discrepancy between the steady state levels of many stress-induced mRNAs and their proteins. Sorting of mRNAs for future translation or decay by individual cells could generate potentially different phenotypes in a genetically identical population and enhance its ability to withstand stress

Is There a Classical Nonsense-Mediated Decay Pathway in Trypanosomes?

In many eukaryotes, messenger RNAs with premature termination codons are destroyed by a process called “nonsense-mediated decay”, which requires the RNA helicase Upf1 and also, usually, an interacting factor, Upf2. Recognition of premature termination codons may rely on their distance from either a splice site or the polyadenylation site, and long 3′-untranslated regions can trigger mRNA decay. The protist Trypanosoma brucei relies heavily on mRNA degradation to determine mRNA levels, and 3′-untranslated regions play a major role in control of mRNA decay. We show here that trypanosomes have a homologue of Upf1, TbUPF1, which interacts with TbUPF2 and (in an RNA-dependent fashion) with poly(A) binding protein 1, PABP1. Introduction of a premature termination codon in either an endogenous gene or a reporter gene decreased mRNA abundance, as expected for nonsense-mediated decay, but a dependence of this effect on TbUPF1 could not be demonstrated, and depletion of TbUPF1 by over 95% had no effect on parasite growth or the mRNA transcriptome. Further investigations of the reporter mRNA revealed that increases in open reading frame length tended to increase mRNA abundance. In contrast, inhibition of translation, either using 5′-secondary structures or by lengthening the 5′-untranslated region, usually decreased reporter mRNA abundance. Meanwhile, changing the length of the 3′-untranslated region had no consistent effect on mRNA abundance. We suggest that in trypanosomes, translation per se may inhibit mRNA decay, and interactions with multiple RNA-binding proteins preclude degradation based on 3′-untranslated region length alone

Quantification of codon selection for comparative bacterial genomics

<p>Abstract</p> <p>Background</p> <p>Statistics measuring codon selection seek to compare genes by their sensitivity to selection for translational efficiency, but existing statistics lack a model for testing the significance of differences between genes. Here, we introduce a new statistic for measuring codon selection, the Adaptive Codon Enrichment (ACE).</p> <p>Results</p> <p>This statistic represents codon usage bias in terms of a probabilistic distribution, quantifying the extent that preferred codons are over-represented in the gene of interest relative to the mean and variance that would result from stochastic sampling of codons. Expected codon frequencies are derived from the observed codon usage frequencies of a broad set of genes, such that they are likely to reflect nonselective, genome wide influences on codon usage (<it>e.g</it>. mutational biases). The relative adaptiveness of synonymous codons is deduced from the frequency of codon usage in a pre-selected set of genes relative to the expected frequency. The ACE can predict both transcript abundance during rapid growth and the rate of synonymous substitutions, with accuracy comparable to or greater than existing metrics. We further examine how the composition of reference gene sets affects the accuracy of the statistic, and suggest methods for selecting appropriate reference sets for any genome, including bacteriophages. Finally, we demonstrate that the ACE may naturally be extended to quantify the genome-wide influence of codon selection in a manner that is sensitive to a large fraction of codons in the genome. This reveals substantial variation among genomes, correlated with the tRNA gene number, even among groups of bacteria where previously proposed whole-genome measures show little variation.</p> <p>Conclusions</p> <p>The statistical framework of the ACE allows rigorous comparison of the level of codon selection acting on genes, both within a genome and between genomes.</p