通过端粒长度分布数学建模推断端粒动力学参数

端粒是线性染色体末端的重复序列，对基因组完整性不可或缺：它随每次细胞分裂而缩短，又由逆转录酶端粒酶在特定细胞中加以延伸^[1]。在酿酒酵母中，端粒延伸受 Rap1p 负调控，这套反馈机制能辨别末端结合的 Rap1p 分子数目，从而把端粒维持在一个较窄的长度范围内^[2]；裂殖酵母中 Taz1p 的缺失同样引起端粒急剧变长^[4]。动力学实验还表明，延伸速率随端粒变长而递减，而降解速率近似恒定^[6]；全基因组缺失筛选则鉴定出150余个影响端粒长度的基因^[3]。然而端粒长度本质上是异质的，理应用整条分布而非单一均值来刻画；如今长读长测序已能高分辨率地实测这一分布，但如何由分布形态反推控制缩短与延伸的生化参数，一直缺乏定量的方法框架。

新工作由此切入：作者提出假设，端粒长度分布本身编码着与其生化机制相关的信息。他们建立数学模型，拟合长读长测序所得的酿酒酵母端粒长度分布，证明模型中的端粒长度可表示为一个遍历马尔可夫链，并给出稳态分布的显式表达式；进而证明生化参数可由分布唯一推断，并在已知真实参数的合成数据上完成验证。将这一参数推断流程用于野生型和端粒调控因子突变体的实测分布，得到了生物学上合理的参数值；再对66株天然酵母菌株拟合，揭示出端粒长度调控参数的多样性；此外，个别染色体末端有时对应于不同的调控参数。

相比以往以遗传学手段刻画反馈机制、以平均长度描述表型的研究路径^[2][3][6]，这项研究把分布形状转化为可推断的生化学量，使不同菌株之间、乃至不同染色体末端之间的调控差异可以被定量比较。末端参数的差异也呼应了端粒酶缺失小鼠中的观察：2号染色体端粒显著短于全基因组平均，染色体特异性端粒长度在一定程度上决定了染色体不稳定的发生^[5]。在群体层面，据报道天然酵母端粒长度分布多样性的长读长测序刻画此前已有预印本（Garrido 等）发表^[7]，并据称无需特殊 DNA 制备即可获得这类分布数据^[8]。

局限同样清楚：模型的构建与验证目前限于酿酒酵母，能否外推到端粒维持异常表型与小鼠并不吻合的人类^[1]仍是未知数；参数的“生物学合理性”目前只是间接印证，缺少独立的生化实测对照；个别染色体末端为何适用不同参数，其分子基础也不明了。而围绕这一模型的数学理论仍在发展，据报道相关项目（涉及长时行为分析，J. Olayé 与 M. Tomašević）以“进行中”状态列于法国数学年会页面^[9]。