《PLOS Genetics》：Estimating the distribution of fitness effects of loss of heterozygosity (LOH) events using an engineered library of Saccharomyces cerevisiae
编辑推荐：
本研究通过构建一个携带随机诱导杂合性缺失（Loss of Heterozygosity, LOH）事件的条形码标记酵母二倍体文库，系统地探讨了LOH事件在七种不同环境（包括化学、温度胁迫和活体宿主模型）中对适应度影响的分布（Distribution of Fitness Effects, DFE）。研究发现，LOH事件的DFE主要表现为中性到有害，且较长的LOH片段通常与更低的适应度相关。然而，在新型或胁迫条件下，有益的LOH事件虽然罕见但可带来显著的适应优势。一个关键发现是普遍存在的拮抗性多效性（antagonistic pleiotropy），超过75%的菌株表现出适应度权衡，即在一个环境获得优势的代价是在其他环境中表现不佳。这项工作揭示了LOH在进化中的双重角色：通常是一种遗传负担，但也为快速、环境特异性的适应和生态位特化提供了强大动力。
实验设计与LOH文库构建
为了系统地研究全基因组范围内LOH事件对适应度影响的分布，本研究开发了一种高通量方法，在二倍体酿酒酵母中生成随机的LOH事件，以便易于测量其适应度。研究的主要目标是：构建一个多样化、带有条形码标记的二倍体酿酒酵母菌株文库，每个菌株携带不同的LOH事件；通过竞争性适应度实验和条形码测序（Bar-Seq）定量测定这些LOH事件相关的DFE；通过测量跨环境的适应度来确定LOH事件的多效性效应。
为了构建这样的文库，研究采用了源自高效基因组编辑的方法，利用稀有限制性内切酶I-SceI在BY4741 × SK1杂交株中引入DNA双链断裂。DSB后续以同源染色体为模板进行修复，从而在基因组中产生随机的LOH事件。所得菌株组成了一个文库，其中包含具有广泛LOH事件的二倍体酵母，这些事件在长度和基因组位置上各不相同，并且每个菌株都标记有一个可用于追踪基因型和进行混合适应度实验的条形码序列。通过对文库中195个FUdR抗性菌株进行全基因组测序，最终确定了150个二倍体LOH菌株和7个二倍体对照菌株用于后续分析。
构建的LOH文库菌株的特性
所构建的LOH菌株文库显示出了LOH事件的多样性，在长度、位置和亲本来源上均存在差异。绝大多数（82.6%）的LOH事件为SK1型。在单个菌株中，SK1型LOH的累积长度也占主导地位，其中74.0%的菌株没有BY4741型LOH。观察到的LOH片段覆盖了大部分染色体区域。末端LOH（T-LOH）少于中间LOH（I-LOH），比例约为38%对62%。T-LOH倾向于产生更长的LOH片段，其中位数长度为200,241 bp，而I-LOH的中位数长度为10,478 bp。每个染色体上的LOH事件数量与染色体长度强相关，表明DSB在整个基因组中是随机诱导的。
不同环境中诱导LOH事件的适应度效应分布
通过Bar-seq在多种环境条件下定量测定了LOH文库菌株和对照菌株的适应度，以确定其DFE。测试的环境包括六种不同的液体生长环境，代表了酵母实验进化中使用的各种胁迫源和程度，以及使用大蜡螟（Galleria mellonella）的无脊椎动物宿主模型。具体包括：标准合成完全培养基（SC）、30°C的酵母提取物蛋白胨葡萄糖培养基（YPD-30C）、酸性合成培养基（pH 2.6）、含4.5 mM过氧化氢（H2O2）的YPD、含2 mg/ml咖啡因的YPD以及大蜡螟体内（20°C）环境。
总体适应度图谱显示，LOH的影响在测试环境中各不相同，从中性到负面。在SC和咖啡因环境中，平均适应度变化接近中性。在YPD-30C、酸性、YPD-37C、大蜡螟-20°C和H2O2条件下，观察到更明显的平均适应度损害。除了YPD-30C和SC，其他五种处理都显示出双峰分布模式，存在一个接近中性的峰和一个低适应度峰。
通过数值最大似然估计进一步推断单个LOH片段的底层分布模型。结果显示，咖啡因数据集为负偏的威布尔分布，而SC为轻尾广义正态分布。其他处理的右尾（极端值）分析表明，YPD-30C和YPD-37C环境中适应度增益的右尾存在截断（即处于威布尔域），最大选择系数（s）小于0.020。而其他五种处理（酸性、H2O2、SC、咖啡因、大蜡螟-20°C）则表现出指数右尾，观察到的最大s值在0.061到0.193之间。
一个普遍趋势是，菌株中LOH的累积长度增加通常与更负面的适应度效应相关。在七个环境中的三个以及跨处理的平均适应度中，都观察到了显著的负相关。类似但较弱的相关性也出现在LOH事件数量上，但仅与T-LOH的数量显著相关，这表明是同合子区域本身对适应度产生了负面影响。累积LOH长度与全基因组杂合性之间呈负相关。定量性状位点（QTL）作图表明，所研究环境中的适应度大多是多基因控制的。在SC培养基中观察到一个相反的趋势，即较高的LOH事件数量和累积长度增加了菌株的适应度。
LOH的多效性效应
尽管许多LOH事件赋予中性或负面适应度效应，但特定的LOH菌株在特定环境中表现出显著的适应度增益，展示了它们驱动适应性进化的潜力。在某些处理中，适应度最高的菌株生长速率在H2O2中增加了高达19.3%，在SC中增加18.7%，在咖啡因中增加14.4%，在酸性环境中增加8.65%。在这些环境中表现最佳的菌株，其基因组中具有较长的LOH片段，这与YPD中表现最佳但LOH片段最少的菌株形成了对比。
然而，这种适应度增益是环境特异性的，并且以其他环境中表现更差为代价，表明存在拮抗性多效性。例如，菌株NFS-2A在咖啡因环境中表现出显著的持续适应度增益，但在酸性环境中表现明显更差；菌株P1-4C在酸性环境中表现出增益，但在咖啡因环境中表现更差。绝大多数菌株（75.7%）表现出这种权衡，即在某些环境中具有正适应度，而在其他环境中具有负适应度。一些菌株（23.5%）在所有测试环境中都表现出持续的负适应度，显示了同合子区域带来的潜在负担。只有一个菌株（P1-5A）在所有处理中均具有正适应度变化。此外，研究观察到，跨越更大基因组区域的LOH事件更可能导致显著的适应度权衡。分析发现，菌株的适应度权衡指标与菌株的累积LOH长度之间存在统计学上显著的正相关。使用其他多效性指标（如跨环境适应度的标准差或范围）也得到了一致的正相关结果。这表明，尽管长LOH片段在大多数环境中可能是有害的，但其产生的变异可导致在特定环境中具有较高的相对适应度。
讨论与结论
LOH是克隆二倍体生物中普遍存在的现象，但随机LOH事件产生的表型变异量尚不清楚。本研究通过一种遗传工程方法，成功探究了随机LOH在多样化LOH文库中的适应度图谱。结果显示，与点突变等其他变异类型类似，大多数LOH事件是中性的或有害的，其DFE偏向负面。然而，与环境相关，相当一部分LOH事件可能是有益的。尽管LOH的DFE可能是双峰的，但没有像自发突变那样存在大量高度有害的等位基因峰。LOH较低的危害性DFE与它重组了已经历自然选择的等位基因这一事实一致，这些等位基因的遗传负荷预期较低。
研究的主要发现支持LOH是快速、环境特异性适应的强大来源。在测试的每个环境中，都有一部分LOH菌株表现出显著的适应度增益，其中一些菌株在所有重复中均显示出稳健的正选择系数，增益高达19.2%，这些菌株很可能在几代之内就在种群中占据主导地位。此外，一个显著成果是拮抗性多效性的普遍存在，绝大多数菌株在至少一种环境中具有适应性。当基因型专门适应于特定环境而以牺牲普遍性能力为代价时，可以加速局部适应或生态位分化，从而导致谱系多样化和物种形成。
这项系统的全基因组评估揭示了LOH在产生适应度变异方面极其多方面的作用。研究表明，LOH的适应度效应并非单一，而是遗传背景、片段大小和环境之间复杂的相互作用。我们将其最常见的作用描述为一种遗传负担的来源，其中同合性增加带来了显著的适应度成本。然而，我们也证明，其多效性效应可带来环境特异性的适应度增益，这可以克服其固有成本，并通过牺牲普遍稳健性来换取生态位特异性优势，从而成为特化的驱动力。
