[论文解读] Adaptive evolution of transcription factor binding sites

Johannes Berg, Stana Willmann+1 more|arXiv (Cornell University)|Jan 29, 2003

Evolution and Genetic Dynamics参考文献 34被引用 4

一句话总结

本文提出了一种生物物理种群遗传学模型，用于研究通过点突变实现转录因子结合位点的适应性进化，整合了基于DNA-蛋白质结合能量学的适应度景观。结果表明，对功能性结合位点的自然选择导致特定的核苷酸频率相关性，并识别出两个关键约束条件——短基序长度和长启动子区域——这些条件使真核生物调控网络能够实现快速适应性进化。

ABSTRACT

The regulation of a gene depends on the binding of transcription factors to specific sites located in the regulatory region of the gene. The generation of these binding sites and of cooperativity between them are essential building blocks in the evolution of complex regulatory networks. We study a theoretical model for the sequence evolution of binding sites by point mutations. The approach is based on biophysical models for the binding of transcription factors to DNA. Hence we derive empirically grounded fitness landscapes, which enter a population genetics model including mutations, genetic drift, and selection. We show that the selection for factor binding generically leads to specific correlations between nucleotide frequencies at different positions of a binding site. We demonstrate the possibility of rapid adaptive evolution generating a new binding site for a given transcription factor by point mutations. The evolutionary time required is estimated in terms of the neutral (background) mutation rate, the selection coefficient, and the effective population size. The efficiency of binding site formation is seen to depend on two joint conditions: the binding site motif must be short enough and the promoter region must be long enough. These constraints on promoter architecture are indeed seen in eukaryotic systems. Furthermore, we analyse the adaptive evolution of genetic switches and of signal integration through binding cooperativity between different sites. Experimental tests of this picture involving the statistics of polymorphisms and phylogenies of sites are discussed.

研究动机与目标

• 理解调控DNA中转录因子结合位点适应性进化的分子机制。
• 量化在选择、突变和遗传漂变作用下，新功能结合位点出现所需的进化时间。
• 通过选择与突变熵的相互作用，解释观察到的结合位序列模糊性。
• 研究多个位点之间的协同结合如何实现信号整合和调控网络中遗传开关的形成。
• 识别促进复杂调控系统高效进化所需的启动子结构约束条件。

提出的方法

• 基于转录因子与DNA序列之间的结合能，构建一个具有生物学物理基础的适应度景观。
• 采用主方程方法对有限种群中的突变动力学进行序列进化建模。
• 应用包含点突变、遗传漂变和选择的种群遗传学框架，其中选择系数由结合亲和力推导得出。
• 利用稳态近似和突变算符的本征函数展开，估算位点形成的中性等待时间。
• 推导出平均适应性等待时间 T₀ = Q_stat^(L₁)(r_s+1)/(L₁λ)，以量化功能位点出现的时间。
• 分析基序长度、启动子长度和选择强度对调控网络可演化性联合影响。

实验结果

研究问题

• RQ1 在选择作用下，通过点突变形成新功能转录因子结合位点的进化时间尺度是什么？
• RQ2 由于对功能性结合的自然选择，结合位点不同位置之间的核苷酸频率相关性是如何产生的？
• RQ3 为了实现调控位点的高效适应性进化，启动子长度和基序长度需要满足哪些结构约束？
• RQ4 多个位点之间的结合协同性如何促进复杂调控逻辑（如信号整合）的进化？
• RQ5 突变熵在功能结合位点观察到的序列模糊性中起到了多大程度的作用？

主要发现

• 通过点突变形成新功能转录因子结合位点的时间尺度为 T_s ~ 1/(sμN)，显著短于中性时间尺度 T_0 ~ 1/μ。
• 由于对功能性结合的自然选择，即使在无直接相互作用的情况下，结合位点不同位置之间的核苷酸频率也出现特定相关性。
• 结合位点形成的效率关键取决于两个联合约束条件：基序必须足够短，启动子区域必须足够长，以提供足够的序列空间使功能位点得以出现。
• 该模型预测，调控网络可通过过去的选择积累冗余的调控相互作用，从而增强对扰动的鲁棒性。
• 观察到的结合位序列模糊性可由对最优结合的选择与亚最优序列的高突变熵之间的平衡来解释。
• 结合位点适应性进化的时标强烈依赖于有效种群大小 N、突变率 μ 和选择系数 s，其中更大的 N 和 s 会加速位点形成。