The RNA-DNA world and the emergence of DNA-encoded heritable traits.

The RNA world hypothesis confers a central role to RNA molecules in information encoding and catalysis. Even though evidence in support of this hypothesis has accumulated from both experiments and computational modelling, the transition from an RNA world to a world where heritable genetic information is encoded in DNA remains an open question. Recent experiments show that both RNA and DNA templates can extend complementary primers using free RNA/DNA nucleotides, either non-enzymatically or in the presence of a replicase ribozyme. Guided by these experiments, we analyse protocellular evolution with an expanded set of reaction pathways made possible through the presence of DNA nucleotides. By encapsulating these reactions inside three different types of protocellular compartments, each subject to distinct modes of selection, we show how protocells containing DNA-encoded replicases in low copy numbers and replicases in high copy numbers can dominate the population. This is facilitated by a reaction that leads to auto-catalytic synthesis of replicase ribozymes from DNA templates encoding the replicase after the chance emergence of a replicase through non-enzymatic reactions. Our work unveils a pathway for the transition from an RNA world to a mixed RNA-DNA world characterized by Darwinian evolution, where DNA sequences encode heritable phenotypes.

The Effect of Environment on the Evolution and Proliferation of Protocells of Increasing Complexity.

The formation, growth, division and proliferation of protocells containing RNA strands is an important step in ensuring the viability of a mixed RNA-lipid world. Experiments and computer simulations indicate that RNA encapsulated inside protocells can favor the protocell, promoting its growth while protecting the system from being over-run by selfish RNA sequences. Recent work has also shown that the rolling-circle replication mechanism can be harnessed to ensure the rapid growth of RNA strands and the probabilistic emergence and proliferation of protocells with functionally diverse ribozymes. Despite these advances in our understanding of a primordial RNA-lipid world, key questions remain about the ideal environment for the formation of protocells and its role in regulating the proliferation of functionally complex protocells. The hot spring hypothesis suggests that mineral-rich regions near hot springs, subject to dry-wet cycles, provide an ideal environment for the origin of primitive protocells. We develop a computational model to study protocellular evolution in such environments that are distinguished by the occurrence of three distinct phases, a wet phase, followed by a gel phase, and subsequently by a dry phase. We determine the conditions under which protocells containing multiple types of ribozymes can evolve and proliferate in such regions. We find that diffusion in the gel phase can inhibit the proliferation of complex protocells with the extent of inhibition being most significant when a small fraction of protocells is eliminated during environmental cycling. Our work clarifies how the environment can shape the evolution and proliferation of complex protocells.

Network rewiring promotes cooperation in an aspirational learning model.

We analyze a cooperative decision-making model that is based on individual aspiration levels using the framework of a public goods game in static and dynamic networks. Sensitivity to differences in payoff and dynamic aspiration levels modulates individual satisfaction and affects subsequent behavior. The collective outcome of such strategy changes depends on the efficiency with which aspiration levels are updated. Below a threshold learning efficiency, cooperators dominate despite short-term fluctuations in strategy fractions. Categorizing players based on their satisfaction level and the resulting strategy reveal periodic cycling between the different categories. We explain the distinct dynamics in the two phases in terms of differences in the dominant cyclic transitions between different categories of cooperators and defectors. Allowing even a small fraction of nodes to restructure their connections can promote cooperation across almost the entire range of values of learning efficiency. Our work reinforces the usefulness of an internal criterion for strategy updates, together with network restructuring, in ensuring the dominance of altruistic strategies over long time-scales.

Evolution towards increasing complexity through functional diversification in a protocell model of the RNA world.

The encapsulation of genetic material inside compartments together with the creation and sustenance of functionally diverse internal components are likely to have been key steps in the formation of 'live', replicating protocells in an RNA world. Several experiments have shown that RNA encapsulated inside lipid vesicles can lead to vesicular growth and division through physical processes alone. Replication of RNA inside such vesicles can produce a large number of RNA strands. Yet, the impact of such replication processes on the emergence of the first ribozymes inside such protocells and on the subsequent evolution of the protocell population remains an open question. In this paper, we present a model for the evolution of protocells with functionally diverse ribozymes. Distinct ribozymes can be created with small probabilities during the error-prone RNA replication process via the rolling circle mechanism. We identify the conditions that can synergistically enhance the number of different ribozymes inside a protocell and allow functionally diverse protocells containing multiple ribozymes to dominate the population. Our work demonstrates the existence of an effective pathway towards increasing complexity of protocells that might have eventually led to the origin of life in an RNA world.

How strategy environment and wealth shape altruistic behaviour: cooperation rules affecting wealth distribution in dynamic networks.

Societies rely on individual contributions to sustain public goods that benefit the entire community. Several mechanisms, that specify how individuals change their decisions based on past experiences, have been proposed to explain how altruists are not outcompeted by selfish counterparts. A key aspect of such strategy updates involves a comparison of an individual's latest payoff with that of a random neighbour. In reality, both the economic and social milieu often shapes cooperative behaviour. We propose a new decision heuristic, where the propensity of an individual to cooperate depends on the local strategy environment in which she is embedded as well as her wealth relative to that of her neighbours. Our decision-making model allows cooperation to be sustained and also explains the results of recent experiments on social dilemmas in dynamic networks. Final cooperation levels depend only on the extent to which the strategy environment influences altruistic behaviour but are largely unaffected by network restructuring. However, the extent of wealth inequality in the community is affected by a subtle interplay between the environmental influence on a person's decision to contribute and the likelihood of reshaping social ties, with wealth-inequality levels rising with increasing likelihood of network restructuring in some situations.

Emergence of ribozyme and tRNA-like structures from mineral-rich muddy pools on prebiotic earth.

The RNA world hypothesis, although a viable one regarding the origin of life on earth, has so far failed to provide a compelling explanation for the synthesis of RNA enzymes from free nucleotides via abiotic processes. To tackle this long-standing problem, we develop a realistic model for the onset of the RNA world, using experimentally determined rates for polymerization reactions. We start with minimal assumptions about the initial state that only requires the presence of short oligomers or just free nucleotides and consider the effects of environmental cycling by dividing a day into a dry, semi-wet and wet phases that are distinguished by the nature of reactions they support. Long polymers, with maximum lengths sometimes exceeding 100 nucleotides, spontaneously emerge due to a combination of non-enzymatic, non-templated polymer extension and template-directed primer extension processes. The former helps in increasing the lengths of RNA strands, whereas the later helps in producing complementary copies of the strands. Strands also undergo hydrolysis in a structure-dependent manner that favour breaking of bonds connecting unpaired nucleotides. We identify the most favourable conditions needed for the emergence of ribozyme and tRNA-like structures and double stranded RNA molecules, classify all RNA strands on the basis of their secondary structures and determine their abundance in the population. Our results indicate that under suitable environmental conditions, non-enzymatic processes would have been sufficient to lead to the emergence of a variety of ribozyme-like molecules with complex secondary structures and potential catalytic functions.

RiboD: a comprehensive database for prokaryotic riboswitches.

SUMMARY: Riboswitches are cis-regulatory non-coding genomic segments that control the expression of downstream genes by undergoing conformational change upon ligand binding. We present a comprehensive database of prokaryotic riboswitches that allows the user to search for riboswitches using multiple criteria, extract information about riboswitch location and gene/operon it regulates. RiboD provides a very useful resource that can be utilized for the better understanding of riboswitch-based gene regulation in bacteria and archaea. AVAILABILITY AND IMPLEMENTATION: RiboD can be freely accessed on the web at http://ribod.iiserkol.ac.in/.

The evolution of antibiotic production rate in a spatial model of bacterial competition.

We consider competition between antibiotic producing bacteria, non-producers (or cheaters), and sensitive cells in a two-dimensional lattice model. Previous work has shown that these three cell types can survive in spatial models due to the presence of spatial patterns, whereas coexistence is not possible in a well-mixed system. We extend this to consider the evolution of the antibiotic production rate, assuming that the cost of antibiotic production leads to a reduction in growth rate of the producers. We find that coexistence occurs for an intermediate range of antibiotic production rate. If production rate is too high or too low, only sensitive cells survive. When evolution of production rate is allowed, a mixture of cell types arises in which there is a dominant producer strain that produces sufficient to limit the growth of sensitive cells and which is able to withstand the presence of cheaters in its own species. The mixture includes a range of low-rate producers and non-producers, none of which could survive without the presence of the dominant producer strain. We also consider the case of evolution of antibiotic resistance within the sensitive species. In order for the resistant cells to survive, they must grow faster than both the non-producers and the producers. However, if the resistant cells grow too rapidly, the producing species is eliminated, after which the resistance mutation is no longer useful, and sensitive cells take over the system. We show that there is a range of growth rates of the resistant cells where the two species coexist, and where the production mechanism is maintained as a polymorphism in the producing species and the resistance mechanism is maintained as a polymorphism in the sensitive species.

Bribery games on interdependent complex networks.

Bribe demands present a social conflict scenario where decisions have wide-ranging economic and ethical consequences. Nevertheless, such incidents occur daily in many countries across the globe. Harassment bribery constitute a significant sub-set of such bribery incidents where a government official demands a bribe for providing a service to a citizen legally entitled to it. We employ an evolutionary game-theoretic framework to analyse the evolution of corrupt and honest strategies in structured populations characterized by an interdependent complex network. The effects of changing network topology, average number of links and asymmetry in size of the citizen and officer population on the proliferation of incidents of bribery are explored. A complex network topology is found to be beneficial for the dominance of corrupt strategies over a larger region of phase space when compared with the outcome for a regular network, for equal citizen and officer population sizes. However, the extent of the advantage depends critically on the network degree and topology. A different trend is observed when there is a difference between the citizen and officer population sizes. Under those circumstances, increasing randomness of the underlying citizen network can be beneficial to the fixation of honest officers up to a certain value of the network degree. Our analysis reveals how the interplay between network topology, connectivity and strategy update rules can affect population level outcomes in such asymmetric games.

Phylogenomic and comparative analysis of the distribution and regulatory patterns of TPP riboswitches in fungi.

Riboswitches are metabolite or ion sensing cis-regulatory elements that regulate the expression of the associated genes involved in biosynthesis or transport of the corresponding metabolite. Among the nearly 40 different classes of riboswitches discovered in bacteria so far, only the TPP riboswitch has also been found in algae, plants, and in fungi where their presence has been experimentally validated in a few instances. We analyzed all the available complete fungal and related genomes and identified TPP riboswitch-based regulation systems in 138 fungi and 15 oomycetes. We find that TPP riboswitches are most abundant in Ascomycota and Basidiomycota where they regulate TPP biosynthesis and/or transporter genes. Many of these transporter genes were found to contain conserved domains consistent with nucleoside, urea and amino acid transporter gene families. The genomic location of TPP riboswitches when correlated with the intron structure of the regulated genes enabled prediction of the precise regulation mechanism employed by each riboswitch. Our comprehensive analysis of TPP riboswitches in fungi provides insights about the phylogenomic distribution, regulatory patterns and functioning mechanisms of TPP riboswitches across diverse fungal species and provides a useful resource that will enhance the understanding of RNA-based gene regulation in eukaryotes.

Comparative genomics and phylogenomic analyses of lysine riboswitch distributions in bacteria.

Riboswitches are cis-regulatory elements that regulate the expression of genes involved in biosynthesis or transport of a ligand that binds to them. Among the nearly 40 classes of riboswitches discovered so far, three are known to regulate the concentration of biologically encoded amino acids glycine, lysine, and glutamine. While some comparative genomics studies of riboswitches focusing on their gross distribution across different bacterial taxa have been carried out recently, systematic functional annotation and analysis of lysine riboswitches and the genes they regulate are still lacking. We analyzed 2785 complete bacterial genome sequences to systematically identify 468 lysine riboswitches (not counting hits from multiple strains of the same species) and obtain a detailed phylogenomic map of gene-specific lysine riboswitch distribution across diverse prokaryotic phyla. We find that lysine riboswitches are most abundant in Firmicutes and Gammaproteobacteria where they are found upstream to both biosynthesis and/or transporter genes. They are relatively rare in all other prokaryotic phyla where if present they are primarily found upstream to operons containing many lysine biosynthesis genes. The genome-wide study of the genetic organisation of the lysine riboswitches show considerable variation both within and across different Firmicute orders. Correlating the location of a riboswitch with its genomic context and its phylogenetic relationship with other evolutionarily related riboswitch carrying species, enables identification and annotation of many lysine biosynthesis, transporter and catabolic genes. It also reveals previously unknown patterns of lysine riboswitch distribution and gene/operon regulation and allows us to draw inferences about the possible point of origin of lysine riboswitches. Additionally, evidence of horizontal transfer of riboswitches was found between Firmicutes and Actinobacteria. Our analysis provides a useful resource that will lead to a better understanding of the evolution of these regulatory elements and prove to be beneficial for exploiting riboswitches for developing targeted therapies.

Bribery games on inter-dependent regular networks.

We examine a scenario of social conflict that is manifest during an interaction between government servants providing a service and citizens who are legally entitled to the service, using evolutionary game-theory in structured populations characterized by an inter-dependent network. Bribe-demands by government servants during such transactions, called harassment bribes, constitute a widespread form of corruption in many countries. We investigate the effect of varying bribe demand made by corrupt officials and the cost of complaining incurred by harassed citizens, on the proliferation of corrupt strategies in the population. We also examine how the connectivity of the various constituent networks affects the spread of corrupt officials in the population. We find that incidents of bribery can be considerably reduced in a network-structured populations compared to mixed populations. Interestingly, we also find that an optimal range for the connectivity of nodes in the citizen's network (signifying the degree of influence a citizen has in affecting the strategy of other citizens in the network) as well as the interaction network aids in the fixation of honest officers. Our results reveal the important role of network structure and connectivity in asymmetric games.

Deciphering a survival strategy during the interspecific competition between Bacillus cereus MSM-S1 and Pseudomonas sp. MSM-M1.

Interspecific competition in bacteria governs colony growth dynamics and pattern formation. Here, we demonstrate an interesting phenomenon of interspecific competition between Bacillus cereus MSM-S1 and Pseudomonas sp. MSM-M1, where secretion of an inhibitor by Pseudomonas sp. is used as a strategy for survival. Although B. cereus grows faster than Pseudomonas sp., in the presence of Pseudomonas sp. the population of B. cereus reduces significantly, whereas Pseudomonas sp. do not show any marked alteration in their population growth. Appearance of a zone of inhibition between growing colonies of two species on nutrient agar prevents the expanding front of the MSM-S1 colony from accessing and depleting nutrients in the region occupied by MSM-M1, thereby aiding the survival of the slower growing MSM-M1 colonies. To support our experimental results, we present simulations, based on a chemotactic model of colony growth dynamics. We demonstrate that the chemical(s) secreted by Pseudomonas sp. is responsible for the observed inhibition of growth and spatial pattern of the B. cereus MSM-S1 colony. Our experimental results are in excellent agreement with the numerical results and confirm that secreted inhibitors enable Pseudomonas sp. to survive and coexist in the presence of faster growing B. cereus, in a common niche.

Finite population analysis of the effect of horizontal gene transfer on the origin of an universal and optimal genetic code.

The origin of a universal and optimal genetic code remains a compelling mystery in molecular biology and marks an essential step in the origin of DNA and protein based life. We examine a collective evolution model of genetic code origin that allows for unconstrained horizontal transfer of genetic elements within a finite population of sequences each of which is associated with a genetic code selected from a pool of primordial codes. We find that when horizontal transfer of genetic elements is incorporated in this more realistic model of code-sequence coevolution in a finite population, it can increase the likelihood of emergence of a more optimal code eventually leading to its universality through fixation in the population. The establishment of such an optimal code depends on the probability of HGT events. Only when the probability of HGT events is above a critical threshold, we find that the ten amino acid code having a structure that is most consistent with the standard genetic code (SGC) often gets fixed in the population with the highest probability. We examine how the threshold is determined by factors like the population size, length of the sequences and selection coefficient. Our simulation results reveal the conditions under which sharing of coding innovations through horizontal transfer of genetic elements may have facilitated the emergence of a universal code having a structure similar to that of the SGC.

Riboswitch Scanner: an efficient pHMM-based web-server to detect riboswitches in genomic sequences.

UNLABELLED: Riboswitches are non-coding RNA located in the 5' untranslated regions where they bind a target metabolite used to specify the riboswitch class and control the expression of associated genes. Accurate identification of riboswitches is the first step towards understanding their regulatory and functional roles in the cell. In this article, we describe a new web application named Riboswitch Scanner which provides an automated pipeline for pHMM-based detection of riboswitches in partial as well as complete genomic sequences rapidly, with high sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: Riboswitch Scanner can be freely accessed on the web at http://service.iiserkol.ac.in/∼riboscan/ CONTACT: mukherjee.sumit89@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An Efficient Minimum Free Energy Structure-Based Search Method for Riboswitch Identification Based on Inverse RNA Folding.

Riboswitches are RNA genetic control elements that were originally discovered in bacteria and provide a unique mechanism of gene regulation. They work without the participation of proteins and are believed to represent ancient regulatory systems in the evolutionary timescale. One of the biggest challenges in riboswitch research is to find additional eukaryotic riboswitches since more than 20 riboswitch classes have been found in prokaryotes but only one class has been found in eukaryotes. Moreover, this single known class of eukaryotic riboswitch, namely the TPP riboswitch class, has been found in bacteria, archaea, fungi and plants but not in animals. The few examples of eukaryotic riboswitches were identified using sequence-based bioinformatics search methods such as a combination of BLAST and pattern matching techniques that incorporate base-pairing considerations. None of these approaches perform energy minimization structure predictions. There is a clear motivation to develop new bioinformatics methods, aside of the ongoing advances in covariance models, that will sample the sequence search space more flexibly using structural guidance while retaining the computational efficiency of sequence-based methods. We present a new energy minimization approach that transforms structure-based search into a sequence-based search, thereby enabling the utilization of well established sequence-based search utilities such as BLAST and FASTA. The transformation to sequence space is obtained by using an extended inverse RNA folding problem solver with sequence and structure constraints, available within RNAfbinv. Examples in applying the new method are presented for the purine and preQ1 riboswitches. The method is described in detail along with its findings in prokaryotes. Potential uses in finding novel eukaryotic riboswitches and optimizing pre-designed synthetic riboswitches based on ligand simulations are discussed. The method components are freely available for use.

Bribe and Punishment: An Evolutionary Game-Theoretic Analysis of Bribery.

Harassment bribes, paid by citizens to corrupt officers for services the former are legally entitled to, constitute one of the most widespread forms of corruption in many countries. Nation states have adopted different policies to address this form of corruption. While some countries make both the bribe giver and the bribe taker equally liable for the crime, others impose a larger penalty on corrupt officers. We examine the consequences of asymmetric and symmetric penalties by developing deterministic and stochastic evolutionary game-theoretic models of bribery. We find that the asymmetric penalty scheme can lead to a reduction in incidents of bribery. However, the extent of reduction depends on how the players update their strategies over time. If the interacting members change their strategies with a probability proportional to the payoff of the alternative strategy option, the reduction in incidents of bribery is less pronounced. Our results indicate that changing from a symmetric to an asymmetric penalty scheme may not suffice in achieving significant reductions in incidents of harassment bribery.

Pathways of Genetic Code Evolution in Ancient and Modern Organisms.

There have been two distinct phases of evolution of the genetic code: an ancient phase--prior to the divergence of the three domains of life, during which the standard genetic code was established--and a modern phase, in which many alternative codes have arisen in specific groups of genomes that differ only slightly from the standard code. Here we discuss the factors that are most important in these two phases, and we argue that these are substantially different. In the modern phase, changes are driven by chance events such as tRNA gene deletions and codon disappearance events. Selection acts as a barrier to prevent changes in the code. In contrast, in the ancient phase, selection for increased diversity of amino acids in the code can be a driving force for addition of new amino acids. The pathway of code evolution is constrained by avoiding disruption of genes that are already encoded by earlier versions of the code. The current arrangement of the standard code suggests that it evolved from a four-column code in which Gly, Ala, Asp, and Val were the earliest encoded amino acids.

Two perspectives on the origin of the standard genetic code.

The origin of a genetic code made it possible to create ordered sequences of amino acids. In this article we provide two perspectives on code origin by carrying out simulations of code-sequence coevolution in finite populations with the aim of examining how the standard genetic code may have evolved from more primitive code(s) encoding a small number of amino acids. We determine the efficacy of the physico-chemical hypothesis of code origin in the absence and presence of horizontal gene transfer (HGT) by allowing a diverse collection of code-sequence sets to compete with each other. We find that in the absence of horizontal gene transfer, natural selection between competing codes distinguished by differences in the degree of physico-chemical optimization is unable to explain the structure of the standard genetic code. However, for certain probabilities of the horizontal transfer events, a universal code emerges having a structure that is consistent with the standard genetic code.

Revisiting the physico-chemical hypothesis of code origin: an analysis based on code-sequence coevolution in a finite population.

The origin of the genetic code marked a major transition from a plausible RNA world to the world of DNA and proteins and is an important milestone in our understanding of the origin of life. We examine the efficacy of the physico-chemical hypothesis of code origin by carrying out simulations of code-sequence coevolution in finite populations in stages, leading first to the emergence of ten amino acid code(s) and subsequently to 14 amino acid code(s). We explore two different scenarios of primordial code evolution. In one scenario, competition occurs between populations of equilibrated code-sequence sets while in another scenario; new codes compete with existing codes as they are gradually introduced into the population with a finite probability. In either case, we find that natural selection between competing codes distinguished by differences in the degree of physico-chemical optimization is unable to explain the structure of the standard genetic code. The code whose structure is most consistent with the standard genetic code is often not among the codes that have a high fixation probability. However, we find that the composition of the code population affects the code fixation probability. A physico-chemically optimized code gets fixed with a significantly higher probability if it competes against a set of randomly generated codes. Our results suggest that physico-chemical optimization may not be the sole driving force in ensuring the emergence of the standard genetic code.