Evolutionary trends in the prokaryotic community and prokaryotic community-phage systems

The Evolutionary Constructor software has been used for computer simulation of the life and evolution of communities of unicellular haploid organisms (prokaryotic cells). Opposite trends of the community evolution (simplification and complication of the genome) have been studied. It has been demonstrated that species with reduced genomes tend to replace genetically and metabolically rich species under highly favorable environmental conditions. Under unfavorable conditions, the opposite tendency is observed. It has also been shown that introduction of phages capable of killing the cells into the system may radically change the current evolutionary trend.     

Theories of biological evolution from the viewpoint of the modern systemic biology

Theories of biological evolution advanced in the last 200 years are reviewed from the viewpoint of advances of modern genetics. The theory of gene networks as a key direction of systemic biology is a link connecting different evolutionary theories.     

Prediction of Bacterial and Archaeal Allergenicity with AllPred Program

Nowadays, allergic disorders have become one of the most important social problems in the world. This can be related to the advent of new allergenic agents in the environment, as well as an increasing density of human contact with known allergens, including various proteins. Thus, the development of computer programs designed for the prediction of allergenic properties of proteins becomes one of the urgent tasks of modern bioinformatics. Previously we developed a web accessible Allpred Program (http://www-bionet.sscc.ru/ psd/cgi-bin/programs/Allpred/allpred.cgi) that allows users to assess the allergenicity of proteins by taking into account the characteristics of their spatial structure. In this paper, using AllPred, we predicted the allergenicity of proteins from 462 archaea and bacteria species for which a complete genome was available. The segregation of considered proteins on archaea and bacteria has shown that allergens are predicted more often among archaea than among bacteria. The division of these proteins into groups according to their intracellular localization has revealed that the majority of allergenic proteins were among the secreted proteins. The application of methods for predicting the level of gene expression of microorganisms based on DNA sequence analysis showed a statistically significant relationship between the expression level of the proteins and their allergenicity. This analysis has revealed that potentially allergenic proteins were more common among highly expressed proteins. Sorting microorganisms into the pathogenic and nonpathogenic groups has shown that pathogens can potentially be more allergenic because of a statistically significant greater number of allergens predicted among their proteins.     

Bacteriophages affect evolution of bacterial communities in spatially distributed habitats: a simulation study

Background

Bacteriophages are known to be one of the driving forces of bacterial evolution. Besides promoting horizontal transfer of genes between cells, they may induce directional selection of cells (for instance, according to more or less resistance to phage infection). Switching between lysogenic and lytic pathways results in various types of (co)evolution in host-phage systems. Spatial (more generally, ecological) organization of the living environment is another factor affecting evolution. In this study, we have simulated and analyzed a series of computer models of microbial communities evolving in spatially distributed environments under the pressure of phage infection.

Results

We modeled evolving microbial communities living in spatially distributed flowing environments. Non-specific nutrient supplied in the only spatial direction, resulting in its non-uniform distribution in environment. We varied the time and the location of initial phage infestation of cells as well as switched chemotaxis on and off. Simulations were performed with the Haploid evolutionary constructor software (http://evol-constructor.bionet.nsc.ru/).

Conclusion

Simulations have shown that the spatial location of initial phage invasion may lead to different evolutionary scenarios. Phage infection decreases the speciation rate by more than one order as far as intensified selection blocks the origin of novel viable populations/species, which could carve out potential ecological niches. The dependence of speciation rate on the invasion node location varied on the time of invasion. Speciation rate was found to be lower when the phage invaded fully formed community of sedentary cells (at middle and late times) at the species-rich regions. This is especially noticeable in the case of late-time invasion.

Our simulation study has shown that phage infection affects evolution of microbial community slowing down speciation and stabilizing the system as a whole. This influencing varied in its efficiency depending on spatially-ecological factors as well as community state at the moment of phage invasion.

     

Haploid evolutionary constructor: new features and further challenges

In this paper we consider the recent advances in methodology for modeling of prokaryotic communities evolution and new features of the software package "Haploid evolutionary constructor" (http://evol-constructor.bionet.nsc.ru). We show the principles of building complex computer models in our software tool. These models describe several levels of biological organization: genetic, metabolic, population, ecological. New features of the haploid evolutionary constructor include the modeling of gene networks and phage infections.     

Relationship between a Gene Network Graph and Qualitative Modes of Its Functioning

Theoretical investigation of properties of assumed gene networks constructed from elementary units of two types, genetic elements and control links, was carried out. A test was formulated for a subclass of such networks with cyclic structure called S(n,k) networks allowing calculation-free prediction of the network limit properties (the presence/absence and number of stationary and/or cyclic functioning modes) from a graph of the network structure. The data obtained can be useful for constructing gene networks with predefined properties.     

Generalized Chemokinetic Method for Gene Network Simulation

Development of methods for mathematical simulation of biological systems and building specific simulations is an important trend of bioinformatics. Here we describe the method of generalized chemokinetic simulation generating flexible and adequate simulations of various biological systems. Adequate simulations of complex nonlinear gene networks—control system of cholesterol by synthesis in the cell and erythrocyte differentiation and maturation—are given as examples. The simulations were expressed in terms of unit processes—biochemical reactions. Optimal sets of parameters were determined and the systems were numerically simulated under various conditions. The simulations allow us to study the possible functional conditions of these gene networks, calculate the consequences of mutations, and define optimal strategies for their correction including therapeutic ones. A graphical user interface for these simulations is available at http://wwwmgs.bionet.nsc.ru/systems/MGL/GeneNet/     

NotI flanking sequences: a tool for gene discovery and verification of the human genome

A set of 22 551 unique human NotI flanking sequences (16.2 Mb) was generated. More than 40% of the set had regions with significant similarity to known proteins and expressed sequences. The data demonstrate that regions flanking NotI sites are less likely to form nucleosomes efficiently and resemble promoter regions. The draft human genome sequence contained 55.7% of the NotI flanking sequences, Celera’s database contained matches to 57.2% of the clones and all public databases (including non-human and previously sequenced NotI flanks) matched 89.2% of the NotI flanking sequences (identity ≥90% over at least 50 bp, data from December 2001). The data suggest that the shotgun sequencing approach used to generate the draft human genome sequence resulted in a bias against cloning and sequencing of NotI flanks. A rough estimation (based primarily on chromosomes 21 and 22) is that the human genome contains 15 000–20 000 NotI sites, of which 6000–9000 are unmethylated in any particular cell. The results of the study suggest that the existing tools for computational determination of CpG islands fail to identify a significant fraction of functional CpG islands, and unmethylated DNA stretches with a high frequency of CpG dinucleotides can be found even in regions with low CG content.