In silico evolution of the Drosophila gap gene regulatory sequence under elevated mutational pressure

BMC Evolutionary Biology

Research

Joshua S. Schiffman¹  Maria G. Samsonova²  Aleksandra A. Chertkova²  Konstantin N. Kozlov²  Sergey V. Nuzhdin³  Vitaly V. Gursky⁴ 

[1] Molecular and Computational Biology, University of Southern California, 90089, Los Angeles, CA, USA;Systems Biology and Bioinformatics Laboratory, Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya, 29, 195251, St. Petersburg, Russia;Systems Biology and Bioinformatics Laboratory, Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya, 29, 195251, St. Petersburg, Russia;Molecular and Computational Biology, University of Southern California, 90089, Los Angeles, CA, USA;Systems Biology and Bioinformatics Laboratory, Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya, 29, 195251, St. Petersburg, Russia;Theoretical Department, Ioffe Institute, Polytechnicheskaya, 26, 194021, St. Petersburg, Russia;

关键词: Evolution;    Regulatory sequence;    Gap gene regulatory network;    Thermodynamic expression model;    Drosophila development;   

DOI  :  10.1186/s12862-016-0866-y

来源: Springer

【 摘 要 】

BackgroundCis-regulatory sequences are often composed of many low-affinity transcription factor binding sites (TFBSs). Determining the evolutionary and functional importance of regulatory sequence composition is impeded without a detailed knowledge of the genotype-phenotype map.ResultsWe simulate the evolution of regulatory sequences involved in Drosophila melanogaster embryo segmentation during early development. Natural selection evaluates gene expression dynamics produced by a computational model of the developmental network. We observe a dramatic decrease in the total number of transcription factor binding sites through the course of evolution. Despite a decrease in average sequence binding energies through time, the regulatory sequences tend towards organisations containing increased high affinity transcription factor binding sites. Additionally, the binding energies of separate sequence segments demonstrate ubiquitous mutual correlations through time. Fewer than 10% of initial TFBSs are maintained throughout the entire simulation, deemed ‘core’ sites. These sites have increased functional importance as assessed under wild-type conditions and their binding energy distributions are highly conserved. Furthermore, TFBSs within close proximity of core sites exhibit increased longevity, reflecting functional regulatory interactions with core sites.ConclusionIn response to elevated mutational pressure, evolution tends to sample regulatory sequence organisations with fewer, albeit on average, stronger functional transcription factor binding sites. These organisations are also shaped by the regulatory interactions among core binding sites with sites in their local vicinity.

【 授权许可】

CC BY   
© The Author(s) 2017