Mccmb ’09 Ttcggtgatatccaggcggcgggcaatcatctt
Total Page:16
File Type:pdf, Size:1020Kb
TCAGAAAATGCGCTCCTGATGCACCCATACCGC TGCTTCCACGCGAGACTTGAGCTTCATTTTCTT CAGCATGTGCTTGACGTGCACTTTTACTGTGCT TTCGGTGATATCCAGGCGGCGGGCAATCATCTT GTTCGGCAAACCCTGGGCAATCAGCTTGAGAAT ATCGCGCTCGCGTGGGGTTAACTGGTTAACATC TCAGAAAATGCGCTCCTGATGCACCCATACCGC TGCTTCCACGCGAGACTTGAGCTTCATTTTCTT CAGCATGTGCTTGACGTGCACTTTTACTGTGCT MCCMB ’09 TTCGGTGATATCCAGGCGGCGGGCAATCATCTT POCEEDINGSR OF THE IN TERNATIONAL MSCOWCNFERENCEOO ONCMPUTATIONALO MLECULARBOLOGYOI July 20-23, 2009 Moscow, Russia Organizers ЕР И И Е Н И Ж Б Н И И О О И И Н Б Ф DepartmentofBioengineeringandBioinformatics О Т Р Е М Т Ь А Л ofM.V.LomonosovMoscowStateUniversity Т У И К К А И Ф 1930 BiologicalDepartment ofM.V.LomonosovMoscowStateUniversity У StateScientificCentreGosNIIGenetika Institutefor Information TrasnsmissionProblems,RAS The ScientificCouncil onBiophysicsRAS, EngelhardtInstituteofMolecular Biology Russian AcademyofSciences Sponsoredby Р И RussianFundofBasicResearch INRIA,France INRIA theFrenchNationalInstituteforResearch inComputerScienceandControl DepartmentofBioengineeringandBioinformatics ofM.V.LomonosovMoscowStateUniversity BiologicalDepartmentofM.V.LomonosovMoscowStateUniversity StateScientificCentreGosNIIGenetika InstituteforInformationTrasnsmissionProblems,Russian A cademyof S ciences TheScientificCouncilonBiophysics , RAS ussian cademyof ciences EngelhardtInstituteofMolecularBiology,Russian AcademyofSciences with financial support of RussianFundofBasicResearch INRIA,France (theFrenchNationalInstituteforResearchinComputerScienceandControl) POCEEDINGSR MCCMB ’09 Moscow, Russia July 20-23, 2009 4-TH MOSCOW CONFERENCE ON COMPUTATIONAL MOLECULAR BIOLOGY July 20–23, 2009 NEW METHOD TO IMPROVE ERROR PROBABILITY ESTIMATION APPLIED TO ILLUMINA SEQUENCING IRINA ABNIZOVA 1, TOM SKELLY 1, YUMI YAN 1, TONY COX 1 The new short read sequencing technique introduced new technological and computational challenges. It requires reconsideration of well-known error estimation algorithms, taking into account different sequencing platforms. 1 Wellcome Trust Sanger Institute, Hinxton, United Kingdom, [email protected] 1 4-TH MOSCOW CONFERENCE ON COMPUTATIONAL MOLECULAR BIOLOGY July 20–23, 2009 DETECTION OF GENES THAT UNDERWENT POSITIVE SELECTION IN DEEP-SEA ARCHAEBACTERIA OF PYROCOCCUS GENUS K.V. GUNBIN 1, D.A. AFONNIKOV 2, N.A.KOLCHANOV 2 Pressure is an environmental parameter of crucial importance for organisms. Archaeal species of the Pyrococcus genus live under both normal (~0,1MPa) and high pressures (>10MPa). To date, the genomes of three Pyrococcus species have been completely sequenced: P. furiosus bacteria live under normal pressure, whereas P. horikoshii and P. аbyssi are piezophilic (live in deep sea environment under high pressure at 14MPa and 20MPa, respectively). In this work we analyze the rate of nucleotide substitution in search for genes underwent positive selection in deep-sea species of Pyrococcus genus. A phylogenetic analysis was performed to determine the evolutionary relatedness of the piezophilic species of the Pyrococcus genus and T. kodekaraensis as outgroup. The analysis of phylogenetic tree demonstrates that piezophilic species have a common origin and the ancestor of piezophilic species emerged from archaebacteria phylogenetically close to the extant species of Pyrococcus genus inhabiting in normal pressure environments. Events of positive selection (PS) for adaptation of life under high pressure were searched for the set of 508 homologous genes which protein sequences are close homologs (amino acid sequence identity greater than 40%) and have no paralogs in genomes. We reconstructed genes and proteins of the most recent ancestor of piezophilic species of the Pyrococcus genus and the common ancestor of P. furiosus, P. horikoshii and P. аbyssi species. Reconstructed ancestral sequence of genes and proteins were compared with extant sequences using nonsynonymous to synonymous substitution rate ratio, radical to conservative amino acid replacement rate ratio, also amino acid dissimilarity measures. We use ArCOG functional classification of analyzed genes and demonstrated that positive selection events occurred in genes and proteins of ‘Coenzyme transport and metabolism’ and ‘Energy production and conversion’ functional groups (Table 1). The results suggest 1 Institute of Cytology and genetics SB RAS, [email protected] 2 Institute of Cytology and genetics SB RAS, Novosibirsk State University [email protected]; [email protected] 2 4-TH MOSCOW CONFERENCE ON COMPUTATIONAL MOLECULAR BIOLOGY July 20–23, 2009 that genes of these functional classes may be important for adaptation of piezophilic Pyrococcus species to deep-sea environment. Table 1. ArCOG group enrichment in the full set of analyzed genes and in genes with identified positive selection events. Last column represents the probablilty of difference in number of genes in full and PS sets observed by chance according to Monte Carlo shuffling test with 105 replicas. ArCOG groups with statistical significant difference (p<0.05) shown in bold. ArCOG group Number Number ppp-p---valuevalue of observing in full in PS by random chance dataset group Amino acid transport and metabolism 34 8 0.18225 Carbohydrate transport and metabolism 22 2 0.90471 Cell cycle control; cell division; chromosome partitioning 8 0 * Cell motility 7 2 0.32479 Cell wall/membrane/envelope biogenesis 13 1 0.90748 Coenzyme transport and metabolism 15 6 0.02416 Defense mechanisms 3 0 * Energy production and conversion 33 11 0.01072 Inorganic ion transport and metabolism 16 0 * Intracellular trafficking; secretion; and vesicular transport 6 0 * Lipid transport and metabolism 5 0 * Nucleotide transport and metabolism 24 5 0.36181 Posttranslational modification; protein turnover; chaperones 18 4 0.34395 Replication; recombination and repair 24 4 0.58199 Secondary metabolites biosynthesis; transport and catabolism 5 1 0.59636 Signal transduction mechanisms 3 1 0.42142 Transcription 28 4 0.70909 Translation; ribosomal structure and biogenesis 76 14 0.36997 Function unknown 87 6 0.99903 General function prediction only 75 14 0.34782 Not annotated 6 1 0.66633 Total 508 84 The work was supported by SB RAS integration proJect №109, Scientific School НШ-2447.2008.4, RAS program “Origin and evolution of Biosphere” and CRDF REC-008 grant. 3 4-TH MOSCOW CONFERENCE ON COMPUTATIONAL MOLECULAR BIOLOGY July 20–23, 2009 MATHEMATICAL MODELING OF THE MOLECULAR GENETIC SYSTEMS REGULATING A PLANT DEVELOPMENT ILYA AKBERDIN 1, FEDOR KAZANTSEV 1, STANISLAV FADEEV 2, IRINA GAINOVA 2, VITALY LIKHOSHVAI 1 Keywords: auxin metabolism, gene network, automatic generation, mathematical model, plant development Indole-3-acetic acid (IAA) is physiologically active in the form of the free acid, but can also be found in conJugated forms in plant tissues. IAA can be degraded and redundant pathways lead to its synthesis. Auxin participates in regulation of cell differentiation in development of embryo, leaves, vascular tissue, fruit, primary and lateral root and in controlling apical dominance and tropisms. The regulation of the IAA metabolism (synthesis, conJugation and degradations) is enough complex and may explain in some aspects how this simple substance is able to influence such diverse processes. Mathematical modeling of IAA metabolic gene network can help reveal the main factors governing this complex process. To reach this aim, we first reconstructed a gene network of auxin biosynthesis, conJugation degradation by annotating experimental data from 107 published papers into GeneNet computer system. This gene network after reduction was input into converter to generate the mathematical model of auxin metabolism. We have reconstructed the gene network and develop the mathematical model of auxin metabolism in arabidopsis shoots. The model allows to reproduce some phenomenological and molecular-genetic aspects of the auxin role in the plant development. The obtained results confirm adequacy of the developed model. In silico experiments testify to qualitatively rapid processes of the molecular genetic regulation of the systems homeostasis. The cumulative experimental data allowed starting construction of spatial distributed hierarchical model that describe both molecular genetic processes and processes on the level of cell- cell interactions simultaneously. So earlier we’ve developed the cellular automaton model that imitates morphodynamics of embryo development by means of regulation of signals produced by different embryonic cells is a first 1 The Institute of Cytology and Genetics SB RAS, Russian Federation, [email protected] 2 The Institute of Mathematics SB RAS, Russian Federation 4 4-TH MOSCOW CONFERENCE ON COMPUTATIONAL MOLECULAR BIOLOGY July 20–23, 2009 step in modelling the process of development in general and in modelling the gene network for morphogenesis in particular [1]. The next step in mathematical modeling application to studying of the plant development rules is integration of the spatial distributed hierarchical model with model of the intracellular auxin metabolism. Akberdin I.R., Ozonov E.A., Mironova V.V., Gorpinchenko D.N., Omelyanchuk N.A., Likhoshvai V.A., Kolchanov N.A. (2007). “A cellular automaton to model the development of shoot meristems of Arabidopsis thaliana”, Journal of Bioinformatics and Computational Biology Vol. 5, pp. 641-650. 5 4-TH MOSCOW CONFERENCE ON COMPUTATIONAL MOLECULAR BIOLOGY July 20–23, 2009 WATER-MEDIATED HYDROGEN BONDS ARE ESSENTIAL FOR LOOP STABILIZATION IN PROTEIN STRUCTURES EVGENIY AKSIANOV 1, SERGEI SPIRIN 1,2,