进化引导的酵母互补分析揭示人类PSPH变体的功能差异

《FEBS Open Bio》：Evolution-guided yeast complementation reveals functional differences in human PSPH variants

　　

解析人类遗传变异如何影响保守代谢酶，对于理解其进化和临床意义至关重要。在此，研究人员将时间分层人类基因组的序列分析与定量酿酒酵母（Saccharomyces cerevisiae）互补分析相结合，在缺乏SER2（酵母中L-丝氨酸生物合成最后一步所需的基因）的菌株中，检验人类磷酸丝氨酸磷酸酶（human phosphoserine phosphatase, PSPH）变体之间的功能差异。古代狩猎采集者与现今人类之间的群体基因组比较鉴定出两个PSPH外显子在核苷酸多样性方面存在较高差异，指导选择两个古代基因组优先变体（R27S和Q83H）进行功能测试。为了将其效应置于功能背景下，研究人员分别表达每个变体，并在多种环境条件下与现代PSPH等位基因及两个疾病相关等位基因（D32N和A35T）进行比较互补。人类PSPH增强了SER2缺失突变体的生长，并揭示了等位基因之间可重复的定量差异。现代等位基因通常赋予最强的互补能力，而古代基因组变体支持可测量但更具条件依赖性的拯救，疾病相关等位基因则显示出最弱的互补能力。这些差异在不同条件下大体一致，而特定环境扰动揭示了效应大小的背景依赖性变化。总之，研究结果建立了一个可扩展的框架，将进化基因组学与实验功能分析联系起来，以识别和评估具有可测量体内效应的人类代谢酶变体。

L-丝氨酸（L-serine）支持蛋白质合成、氮和氨基酸代谢以及磷脂酰丝氨酸、磷脂酰乙醇胺和鞘脂等必需膜脂质的产生。人类磷酸丝氨酸磷酸酶（human phosphoserine phosphatase, PSPH）基因编码磷酸丝氨酸磷酸酶，其变异与L-丝氨酸缺乏症及神经发育表型相关。已有功能研究表明，PSPH替换如D32N和M52T可产生从轻微或背景依赖性变化到酶活性、蛋白稳定性、相互作用或细胞表型显著降低等一系列效应。Sacchi et al.最近报道，在阿尔茨海默病海马样本中检测到的R27S/D32G单倍型可降低催化效率，但杂合条件下不显著降低L-丝氨酸产生，同时改变亚细胞定位。相反，PSPH过表达与癌症增殖和代谢表型相关，表明该通路对PSPH剂量、活性和细胞背景敏感。鉴于PSPH的核心代谢作用及与人类疾病的关联，理解遗传变异如何塑造其功能仍是一个重要挑战。虽然全基因组群体方法和大型测序项目已发现大量变异，古DNA数据集也增加了时间维度，但很少有古代或古老等位基因被功能表征。因此，将古代和现代基因组的进化信号与可扩展功能分析相结合，对于解释人类遗传变异的生物学后果至关重要。

在这项研究中，研究人员将时间分层人类基因组的序列分析与定量酿酒酵母互补分析相结合，在缺乏SER2（酵母中L-丝氨酸生物合成最后一步所需基因）的菌株中，检验人类PSPH变体之间的功能差异。古代与现今人类基因组的群体比较鉴定出两个PSPH外显子具有升高的核苷酸多样性差异，指导选择R27S和Q83H两个古代基因组优先变体，并与现代PSPH等位基因及两个疾病相关等位基因D32N和A35T进行比较。结果表明，人类PSPH能够增强SER2缺失突变体的生长，并表现出可重复的等位基因特异性定量差异：现代等位基因通常互补最强，古代变体支持可测量但更具条件依赖性的拯救，疾病相关等位基因最弱。该研究建立了一个可扩展框架，将进化基因组学与实验功能分析相结合，用于识别和评估具有可测量体内效应的人类代谢酶变体，论文发表在《FEBS Open Bio》。

主要技术方法包括：分析8个古代人类基因组（MA-1、Anzick-1、Kennewick、Loschbour、Saqqaq、Bichon、Kotias和La Bra?a-1）、高覆盖度尼安德特人和丹尼索瓦人基因组，以及1000 Genomes Project Phase 3的2504个现代人类基因组；使用Massive Tajima软件在非重叠100-bp窗口中估计核苷酸多样性π，并计算基因组间差异Δπ；利用Ensembl VEP注释候选编码变异，并用ChimeraX进行人PSPH与酵母Ser2的结构比对。在酿酒酵母BY4741中通过PCR介导的同源重组敲除SER2，将人PSPH等位基因克隆到pCM189载体并进行定点突变，转化酵母后，先通过稀释点板实验定性验证互补，再在96孔板中测量液体生长动力学。在35种环境条件下测试生长，计算归一化拯救系数R，并使用单因素ANOVA和Tukey HSD检验评估统计差异。

全基因组核苷酸多样性分析优先筛选出古代和现代人类数据集中的PSPH候选变异。通过比较8个古代人类基因组与1000 Genomes Phase 3现代人类数据集的π，研究人员在PSPH基因内发现两个编码外显子（外显子4和外显子5）的Δπ值显著升高，峰值分别为3.94×10^?3和3.55×10^?3，位列全基因组窗口前1.42%和1.86%。这些区域包含14个变异，其中外显子4有2个非同义变异，外显子5有6个。研究优先考虑至少两个古代基因组支持且产生非同义替换的变异，最终选择R27S、D32G和Q83H。R27S和D32G也存在于尼安德特人和丹尼索瓦人基因组中。功能预测显示Q83H和R27S为中性，而D32G预测有害，与残基32与疾病关联的证据一致。这些结果确定了用于功能测试的候选等位基因。

酵母互补揭示人类PSPH变体间等位基因特异性功能差异。研究人员开发了定量异源互补分析，通过测量生长动态而非终点拯救来检测分级功能差异。稀释点板实验确认人类PSPH可恢复ser2Δ突变体在需要内源丝氨酸生物合成条件下的生长，表明磷酸丝氨酸磷酸酶功能在真核生物中保守。在1011个酿酒酵母基因组面板中，SER2的编码序列核苷酸多样性高于基因组中位数，且是主要丝氨酸生物合成基因中最高的，支持其作为内源参考。随后在单倍体ser2Δ背景中比较现代PSPH、R27S、Q83H、D32N和A35T的互补能力。每个构建体使用5个独立转化菌落作为生物学重复，每个重复测量6次。单因素ANOVA显示构建体间差异显著（P<10^?23），Tukey HSD检验将酵母SER2、现代PSPH、古代变体和疾病变体分为不同互补组。酵母SER2拯救最强，与野生型无差异；现代PSPH也提供强拯救，但显著低于酵母SER2；疾病相关等位基因最弱，且彼此无差异；古代变体居中，Q83H与现代等位基因有显著差异，R27S效应更温和，两者之间也有轻微差异。补充外源L-丝氨酸（375 mg·L^?1）可恢复生长并减少等位基因间差异，表明观察到的表型主要反映丝氨酸生物合成限制。

环境谱分析揭示PSPH等位基因效应的背景依赖性调节。研究人员在35种环境条件下测量互补能力，包括替代碳源、离子和金属胁迫、乙醇暴露、DNA复制胁迫、细胞骨架扰动等。5个生物学重复间高度一致（归一化RMSD为6.8%），野生型携带空质粒与SER2互补菌株的生长率高度相关（r=0.985），表明系统可靠。使用归一化拯救系数R（R=0表示无拯救，R=1表示完全拯救）量化互补能力，跨条件的累积分布显示一致的互补层级：现代等位基因显著强于疾病相关变体和Q83H，但与R27S无显著差异；疾病相关变体最弱且彼此无差异，但显著弱于两个古代变体；两个古代变体之间也有轻微差异。少数环境条件改变了效应大小：10 μm苯菌灵、0.12 mm CaCl₂和3%乙醇加大现代等位基因与部分变体的差距；0.5 mm CuSO₄、2 mm乙醇胺和棉子糖则减轻某些变体的缺陷。这些结果表明，虽然等位基因特异性趋势总体可重复，环境扰动可调节效应幅度并揭示基因型-环境互作。

讨论部分强调，该研究建立了人类PSPH对酵母ser2Δ突变体的功能互补，作为跨物种体内分析系统，并比二元互补读数更灵敏。现代参考等位基因通常比测试变体支持更强互补，疾病相关等位基因差异最明显，古代变体R27S和Q83H显示可测量但更条件依赖的拯救。特别值得注意的是，这两个古代变体被序列工具预测为中性，但酵母分析捕捉到标准计算预测无法检测的细微体内效应。进化分析用于优先排序，但古代与现代数据集大小和采样结构差异大，多样性模式可能反映选择、人口统计、漂变或采样历史。R27S、D32G和Q83H在1000 Genomes Phase 3中未观察到多态性，但在gnomAD等资源中存在于现代人群，尤其是非洲和非洲裔美国人以及部分东亚人群，提示其可能为古老的站立变异。结果与先前生化研究一致，即PSPH变体可通过多种机制影响酶功能。该研究的局限性包括未测试D32G和M52T，且生长拯救不能区分催化效率、蛋白稳定性或相互作用机制，也不能反映人体组织特异性调控。研究结论翻译如下：该研究将古代基因组分析、临床变异优先级排序和实验功能生物学整合在一个工作流程中。古代与现代基因组比较有助于优先筛选那些在临床注释、现代频率或计算机预测中不明显的变异，而酵母互补可快速测试这些变异是否在体内改变PSPH依赖性生长。随着古代基因组数据集和临床变异目录不断扩展，整合进化基因组学与可扩展功能分析，对于识别那些仅靠现代数据集可能无法显现其生物学意义的变异至关重要。