本研究由北京生命科学研究所（National Institute of Biological Sciences, Beijing）的杜立林（Li-Lin Du）课题组主导完成，第一作者为胡雯（Wen Hu）。研究团队还包括蒋兆迪（Zhao-Di Jiang）、索芳（Fang Suo）、郑金鑫（Jin-Xin Zheng）以及何万中（Wan-Zhong He）。该研究成果于2017年6月20日发表于国际知名期刊《eLife》，题为《A large gene family in fission yeast encodes spore killers that subvert Mendel’s law》（一个庞大的粟酒裂殖酵母基因家族编码颠覆孟德尔定律的孢子杀手）。本文属于基于原始研究的研究论文（类型A），以下将从研究背景、工作流程、主要结果、研究结论及价值与亮点等方面进行详细报告。
孟德尔的均等分离定律（Mendel’s law of equal segregation）认为，等位基因在遗传给后代时应享有平等的机会。然而，被称为减数分裂驱动元件（meiotic drivers）的自私遗传元件（selfish genetic elements）能通过扭曲分离比来使自身获得遗传优势。在真菌中，此类元件表现为孢子杀手（spore killers），特异性地破坏未继承该基因的孢子。然而，此前尚不明确到底有多少真菌物种携带孢子杀手，以及它们作用机制的多样性程度。
粟酒裂殖酵母（*Schizosaccharomyces pombe*）作为模式生物，其在杂交实验中常表现出严重的内种杂交不育（intraspecific hybrid sterility），孢子存活率往往低于5%。尽管此前研究发现，从康普茶中分离的“kambucha”菌株与实验室菌株杂交时，至少有三种孢子杀手导致了杂交不育，但具体的杀手基因身份始终未被揭开。因此，本研究的根本目标在于从分子水平揭示导致内种杂交不育的孢子杀手元件，并深入阐明其作用机制、演化特征与生物学意义。
1. 孢子杀手基因的初步定位 研究团队首先选取了一株从西班牙葡萄中分离的天然菌株CBS5557，与包含2.23 Mb染色体倒位的实验室菌株进行杂交。尽管引入相同的倒位使孢子存活率从14%提高到22%，仍有大量不育现象存在。为探究剩余的不育原因，作者利用二代测序（NGS）辅助的集群分离分析（bulk segregant analysis）技术，对存活后代的全基因组SNP等位基因频率进行检测。结果显示，实验室菌株来源的III号染色体上的SNP等位基因频率显著偏低，形成了左臂和右臂两个宽的传播扭曲（transmission distortion）峰，表明至少有两个孢子杀手基因位于CBS5557的III号染色体上。
2. 杀手基因的精细定位与鉴定 为了精确定位这些杀手基因，研究团队对含有部分CBS5557染色体片段的回交菌株进行了多轮回交与NGS辅助集群分离分析。通过检测回交一代（backcrossed-1）和回交二代（backcrossed-2）菌株的分离模式，他们将候选区域缩小到特定的CBS5557来源的基因组片段。随后，他们利用Illumina和PacBio测序数据对差异最大的序列进行检查，发现两个基因，分别为*cw9*和*cw27*，均属于此前功能未知的*wtf*（with Tf LTRs）基因家族。*cw9*与参考基因组中的*wtf9*基因序列一致性低至42%；而*cw27*则在参考基因组中对应的位置因长末端重复序列（LTR）介导的重组而完全缺失。
3. 基因功能验证与遗传学分析 为验证*cw9*和*cw27*的孢子杀手功能，研究者在CBS5557菌株背景下构建了杂合缺失或纯合缺失突变体。通过在CBS5557自身的h+/h-二倍体中进行四分子分析（tetrad analysis），他们发现纯合缺失对孢子存活无影响，而杂合缺失则导致孢子活力显著下降，且死亡的孢子正是那些继承缺失等位基因的孢子。当两个杀手基因同时处于杂合状态时，双缺失孢子死亡尤为严重，表明继承一种杀手只能对抗其自身的杀伤效应，而无法对抗另一种杀手。
为进一步验证这些杀手基因在不同基因组环境下的功能独立性，研究者构建了携带*cw9*或*cw27*的整合型质粒，将其整合到实验室菌株基因组的不同位置（his3、lys1、leu1等位点）。杂交实验显示，只要一方亲本携带杀手基因，未继承该基因的孢子便大幅死亡，且杀手的效应不受插入位点影响。当双亲分别携带不同的杀手基因时，继承任一杀手基因的孢子均会经历严重的活力丧失，进一步证实二者之间不存在交叉抗性。
4. 超微结构观察与截短分析 为了从形态学层面观察杀伤过程，作者对同步进行减数分裂与孢子形成的二倍体细胞进行了电子显微镜（electron microscopy）分析。在转入氮源缺乏培养基18小时后，含*cw27*的细胞形成的子囊中，两个孢子显示出细胞质电子密度低、细胞壁电子透明以及非球形形态等异常，表明未受保护的孢子无法顺利完成成熟过程。
为解析杀手基因的结构-功能关系，作者进一步对*cw9*和*cw27*进行了一系列截短突变分析（ta、tb、tc、td、te形式）。结果显示：位于上游的LTR对杀伤活性和保护活性均非必需；然而，若缺失紧邻保守上游序列（conserved_up sequence）上游的约90 bp序列（tb形式），携带该基因的孢子自身也遭受一定的活力丧失，表明保护活性受损；而缺失C末端最后10个氨基酸的截短形式（td形式）完全丧失了杀伤活性，但保留了完整的保护活性。这种“活性分离”突变表明，一个基因可通过不同的表达产物分别实现杀伤与保护功能，有力支持了随后由Nuckolls等人提出的毒药-解毒药双重转录模型。
5. 比较基因组学与系统发育分析 最后，研究者对参考基因组中的25个*wtf*基因和CBS5557中的29个同源*wtf*基因进行了详细的比较分析。他们发现在多个位点上存在基因拷贝数的差异（如从单拷贝变为双拷贝或三拷贝），并有三个仅在CBS5557中存在的singleton *wtf*基因（cw26~*cw28*）在实验室菌株中因LTR介导的重组而丢失。利用最大似然法构建的系统发育树显示，许多基因存在异位基因转换（ectopic gene conversion）的迹象，而同一基因簇内的*wtf*基因却保持高度分化，提示多样化选择（diversifying selection）阻止了序列同质化。
本研究取得了一系列层次递进且互有关联的结果。首先，通过集群分离分析确定了CBS5557菌株的III号染色体上存在至少两个孢子杀手区域，并成功将杀手基因锁定为*wtf*家族的*cw9*和*cw27*。缺失分析、回交分析以及不同遗传背景下的杂交实验均充分证实，这两个基因以独立遗传、独立功能和独立于基因组位置的方式发挥作用。电子显微镜的结果为杀手导致孢子成熟障碍提供了直观的超微结构证据。截短实验结果不但揭示杀手基因的C末端对其杀伤活性至关重要，更通过td活性分离突变体首次证明保护活性与杀伤活性可以在同一编码序列中被突变分离。最后，比较基因组学与系统发育学数据全面展示了*wtf*基因家族在不同菌株间的快速分化、基因丢失以及非等位基因转换的演化图景。
该研究的核心结论是：粟酒裂殖酵母的*wtf*基因家族成员编码了能颠覆孟德尔分离定律的孢子杀手，其中*cw9*和*cw27*以单基因、自足的方式行事，通过产生不同的蛋白行使毒药（poison）和解毒药（antidote）的双重功能。本研究的科学价值在于首次从分子水平揭示了裂殖酵母内种杂交不育的关键遗传元件，为理解减数分裂驱动元件的演化机制提供了全新的理想模型。此外，该研究发现*wtf*基因家族的高度多样性和快速演化直接导致了不同天然菌株间广泛的生殖隔离（reproductive isolation），对物种形成和基因组演化研究具有深远的启发性。同时，由于该家族的庞大及其仅存于S. pombe的特性，也为研究物种特异性基因家族的扩张和功能分化提供了一个极其珍贵的案例。
本研究的亮点主要体现在以下几个方面：其一，发现新型杀手元件，揭示了此前功能未明的*wtf*基因家族的生物学本质——作为强大的减数分裂驱动者；其二，解析“杀手”与“保护”双重活性，通过巧妙的截短实验获得活性分离突变体，奠定后续“毒药-解毒药”双重转录模型的基础，这属于机制研究上的重要突破；其三，展现基因家族的演化动力学，综合运用PacBio测序和系统发育分析，清晰展示了*wtf*基因通过基因复制、异位基因转换和LTR介导重组等方式实现快速的拷贝数和序列演化，以形成不同菌株间互不兼容的杀手体系；其四，方法学严谨性，结合了集群分离分析、回交策略、四分子分析、电子显微镜以及多种基因操作技术，构成了从遗传定位到分子机制再到进化分析的完整证据链。这些成果不仅强有力地回应该领域长期以来的关键科学问题，也为更广泛的遗传冲突与物种形成研究范式带来了诸多启发。