研究背景

纤毛虫是一类具有核双态性的单细胞真核生物，其生殖系小核（MIC）包含完整基因组，而体细胞大核（MAC）在发育过程中通过程序性DNA消除去除转座子等“垃圾DNA”。传统模型（如四膜虫、草履虫）认为生殖系来源的小RNA主要标记待清除序列，但在旋毛类纤毛虫（如Oxytricha）中，体细胞来源的sRNA反而起到保护体细胞DNA的作用。海洋游仆虫（Euplotes vannus）作为旋毛类的姐妹群，其基因组高度碎片化为数万个“纳米染色体”，但sRNA如何在此系统中调控DNA重组仍不清楚。

本研究通过PacBio HiFi长读长测序技术首次完成了Euplotes vannus高质量生殖系基因组组装（431 Mb），并发现约80%的生殖系DNA在体细胞发育中被消除。进一步研究揭示，一类体细胞来源的30nt小RNA在接合过程中大量积累，并通过精准靶向保护非转座子DNA不被清除。这些sRNA由Dicer-like酶（Dcl1）从长链双链RNA加工而来，其转录起始于保守的亚端粒染色体断裂序列（E-CBS，5′-TTGAA-3′）。结合Nanopore全长lncRNA测序与功能实验，本研究提出了一个全新的sRNA介导基因组重组模型（图8），不仅深化了我们对纤毛虫DNA消除机制的理解，也为真核生物中sRNA在遗传信息传递与基因组稳定性维持中的作用提供了新视角。

文章标题：Soma-derived 30-nt small RNAs are coupled with chromosome breakage and precisely target nontransposon DNA against elimination in Euplotes vannus

发表期刊：Science Advances（IF：12.5）

发表时间：2025年10月  

影响因子：12.5

发表单位：中国海洋大学高凤团队

贝纳基因参与了本研究中Nanopore全长lncRNA测序工作

主要研究结果

生殖系基因组特征与大规模DNA消除

研究成功组装了Euplotes vannus的431 Mb生殖系基因组（表1），并与85 Mb的体细胞基因组进行比较，发现高达80%的生殖系DNA在体细胞发育中被消除（图1B）。被清除的序列包括65.8%的转座子相关序列（TE GSS）和13.9%的非转座子序列（non-TE GSS）（图1A）。其中 ，Tec转座子（属于Tc1/mariner超家族）占据主导，约占生殖系基因组的31.3%。进一步分析显示，大多数被清除序列的边界含有5′-TA-3′二核苷酸重复（图1F），类似于转座子末端结构，揭示其可能起源于古代转座事件。此外，非转座子清除序列倾向于成簇分布（图1H），且多数位于蛋白质编码区内（图1G），突显了其精准切除的必要性。

表1 海洋游仆虫的基因组组装结果

图 1. 海洋游仆虫生殖系的基因组特征

接合特异性30nt小RNA的积累与起源  

通过对接合过程7个时间点的样品进行sRNA测序，研究发现一类30nt长度的小RNA在接合早期（8小时）开始积累，并在12小时达到峰值（占sRNA总量的60.9%）（图2A）。这些sRNA具有显著的5′-尿苷（U）偏好的特征（图2G），且31–33nt的sRNA在3′端常发生非模板依赖的腺苷化（图2H），表明存在转录后修饰。通过比对到基因组发现，超过85%的30nt sRNA源自体细胞基因组（图2B），且其覆盖度与基因组DNA拷贝数呈正相关，而与mRNA表达水平无关（图2C），表明其前体并非mRNA，而是由体细胞染色体独立转录产生的长链非编码RNA（lncRNA）。

图 2. 接合生殖特异性小RNA的发育动态

合成sRNA微注射证实其保护功能

为验证30nt sRNA的功能，研究团队设计了针对特定清除序列（GSS）的合成sRNA，并通过微注射技术导入接合细胞（图3A）。结果显示，双链sRNA或正义/反义链混合注射能导致目标GSS在子代细胞中被保留，而单链sRNA效果较弱（图3B）。此外，sRNA的长度灵活性也被证实：27–35 bp的sRNA均能诱导GSS保留，但25 bp的sRNA几乎无效（图3C）。这些实验不仅确认了sRNA的保护作用，还表明其机制对sRNA长度具有一定容忍度。

图 3. 合成小RNA诱导天然消除序列的保留

sRNA靶向精准且兼容杂合位点  

在针对Dcl1基因中一个长GSS（2904 bp）的实验中，作者设计了两条重叠5 bp的30nt sRNA，覆盖GSS的3′端区域（图4A）。注射后，子代细胞中出现了11 bp或36 bp的GSS部分保留（图4C），证实sRNA能精确指导切除边界。值得注意的是，即便sRNA与目标序列存在单碱基不匹配（如杂合位点），其保护作用依然有效，说明sRNA与模板的配对并不要求完全互补，这一特性可能有助于维持遗传多样性并在演化中发挥作用。

图 4. 小RNA引导的精准DNA切割与Dicer样核糖核酸酶的关键功能

Dcl1酶是sRNA成熟的关键且为生存所必需  

系统发育分析识别出五个Dcr/Dcl同源蛋白，其中Dcl1（仅含两个RNase III结构域）在接合早期高表达（图4B, E），与sRNA积累时间吻合。通过sRNA诱导的Dcl1基因突变体研究发现，突变细胞无法产生成熟的30nt sRNA（图4D），且在接合过程中大量死亡（≥75%）（图4F），表明Dcl1在sRNA生物发生与基因组重组中不可或缺。

sRNA前体为双向转录lncRNA，起始于E-CBS位点  

免疫荧光染色与lncRNA测序显示，接合早期在体细胞大核中积累大量双链RNA（图5），其反义链转录本在12小时覆盖了绝大多数纳米染色体（图6A）。研究进一步采用Nanopore全长lncRNA测序技术，去除了占多数的rRNA干扰，并能够直接获取来自正义与反义链的完整单分子lncRNA序列。分析结果显示，这些作为sRNA前体的lncRNA，转录起始位点（TSS）均位于距离端粒约17–20 bp的E-CBS（5′-TTGAA-3′） motif处（图6B, E）。这些转录本保留内含子（图6A），且与mRNA无关，支持其作为sRNA前体。研究提出，RdRP可能参与将单链转录本转化为双链RNA，供Dcl1切割生成sRNA。

图 5. 野生型与Dcl1突变体细胞在接合过程中双链RNA的免疫荧光染色

图 6. 接合早期发育性双链RNA的特征

DNA清除相关基因的表达谱分析  

通过7个时间点的mRNA-seq分析，识别出大量在接合过程中差异表达的基因（图7A），包括sRNA生物发生相关基因（如RPB2、Spt4mB、RdRP、Piwi）、染色质修饰因子（如PRC2组件Ezl1、去甲基酶Jmj1）以及DNA切割与修复蛋白（如Ku70/80、SUMO通路成员）。此外，28个MAC编码的转座酶同源物在体细胞发育中高表达（图7C），表明驯化转座酶在DNA清除中发挥作用。

图 7. 海洋游仆虫接合过程中可能与DNA消除相关基因的表达谱

研究总结

本研究通过整合高质量基因组组装、全转录组测序和功能实验，系统揭示了海洋游仆虫中体细胞来源30nt小RNA在程序性DNA清除中的保护作用。这些sRNA并非标记清除对象，而是通过精准靶向保护非转座子DNA，确保其在染色体断裂与重组过程中得以保留。进一步利用Nanopore全长lncRNA测序技术，成功解析了sRNA前体的精确转录边界，阐明sRNA的前体为从E-CBS位点起始的双向转录本，经Dcl1加工后装载至Piwi蛋白，指导发育中新大核的基因组重塑。该机制不仅解决了“如何从数万纳米染色体快速产生sRNA”的难题，也展现了sRNA在真核生物基因组稳定性维持与演化适应中的深远意义。

图 8. 体细胞来源的30nt小RNA在海洋游仆虫基因组重排过程中作用的模型

参考文献：  

Lyu L, Ding B, Fu J, et al. Soma-derived 30-nt small RNAs are coupled with chromosome breakage and precisely target nontransposon DNA against elimination in *Euplotes vannus*. *Sci Adv*. 2025;11(41):eadx3690. doi:10.1126/sciadv.adx3690

全长lncRNA测序——全球首发

全长lncRNA测序是指无需对去除rRNA后的RNA进行打断，直接获取逆转录cDNA的5’到3’全长序列。不依赖于Poly(A)尾富集，同时对含有Poly(A)尾和不含有Poly(A)尾的lncRNA转录本进行鉴定和定量，并系统全面的进行可变剪切分析。

技术优势

1. 无需拼接，直接获得全长lncRNA序列

2. 对lncRNA转录本进行定量

3. 鉴定lncRNA的可变剪切

4. 捕获mRNA、tRNA、小RNA等其它类型的RNA