IF50.5｜少数粗线期piRNA切割mRNA提升精子适应性

大多数粗线期piRNA无功能，仅约1%能切割靶mRNA，其中少数改变靶标丰度者维持雄性生育力。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

动物利用18–33核苷酸的PIWI相互作用RNA（piRNA）沉默生殖细胞中的转座子。胎盘哺乳动物还产生粗线期piRNA，这是一类源自长非编码RNA前体、在睾丸中高丰度表达的小RNA。尽管粗线期piRNA前体转录位点在哺乳动物中保守，但piRNA序列本身快速演化，甚至在人群中也是如此。此前研究提出了多种功能模型：类似miRNA的靶mRNA去稳定化、类似siRNA的切割、或激活部分互补mRNA的翻译，但也有观点认为它们只是降解产物。小鼠中六个主要粗线期piRNA簇的遗传敲除仅有两个导致不育，其余单个敲除未显示明显生育力表型，这引发了对其功能的争议。

科学问题

卡在哪里： 此前仅pi6和pi18两个簇的敲除导致雄性不育，而pi2、pi7、pi9、pi17单个敲除未显示生育力缺陷，这被解释为遗传冗余或功能缺失。同时，粗线期piRNA的作用机制（切割、miRNA样调控或翻译激活）尚无定论，且绝大多数piRNA序列缺乏保守性，其生物学意义不明。

本文要回答： 本文旨在确定粗线期piRNA是否具有生物学功能、其作用机制是什么，以及为何大多数piRNA序列不保守却能保留在哺乳动物进化中。

技术路线

作者主要用了：CRISPR/Cas9基因编辑、RNA-seq、5′单磷酸RNA测序、核糖体足迹测序、GRO-seq、体外切割实验。

作者首先构建了pi2、pi7、pi9、pi17的单敲除、双敲除和三敲除小鼠，通过连续交配评估生育力，并检测精子功能。随后对FACS分选的初级精母细胞进行RNA-seq，鉴定差异表达转录本，并通过5′单磷酸RNA测序检测piRNA切割产物。进一步利用核糖体足迹测序和GRO-seq分别评估翻译效率和转录速率，结合体外重组MIWI切割实验验证piRNA的直接切割活性。最后通过比较piRNA靶标与非靶标的保守性，提出进化模型。

图 1　Fig. 1。图1展示生育力和精子功能数据。a显示连续交配8个月后各基因型雄鼠的产仔数：对照中位数5窝，pi7−/−pi9−/−为3窝，pi9−/−pi17−/−仅0.5窝，pi2−/−pi9−/−pi17−/−几乎为零，Kruskal–Wallis检验P=7.3×10−7。b显示超活化精子比例：对照9.2%，pi9−/− 4.7%（P=0.004），pi17−/− 4.5%（P=0.01），表明单敲除即损害精子功能。c显示进行性运动精子比例，pi2−/−和pi9−/−显著降低。d显示交配后14.5天胚胎数，pi9−/−pi17−/−和pi2−/−pi9−/−pi17−/−雄鼠配偶中位数为0，对照为8。

核心亮点

亮点 1｜piRNA簇间存在冗余但非完全冗余

单敲除pi2、pi7、pi9或pi17不影响生育力，但双敲除pi7−/−pi9−/−和pi9−/−pi17−/−显著降低产仔数（中位数分别为3窝和0.5窝，对照为5窝），三敲除pi2−/−pi9−/−pi17−/−几乎不育。精子功能分析显示，pi2、pi9或pi17单敲除即导致超活化精子比例减半（对照9.2%，pi9−/− 4.7%，P=0.004），且pi2−/−和pi9−/−进行性运动精子减少。这表明这些piRNA簇对精子成熟有非冗余贡献。

图 2　Fig. 2。图2展示pi9和pi17靶标的鉴定。a左为C57BL/6与pi9−/−或pi17−/−初级精母细胞的转录组散点图，差异表达转录本（FDR<0.01）以颜色标出；右为pi9和pi17 piRNA直接靶标的列表，显示piRNA丰度、靶标在突变体中的增加倍数、靶位点位置、配对模式和结合能。b显示pi9和pi17靶标在pi9−/−、pi17−/−、pi9−/−pi17−/−和pi2−/−pi9−/−pi17−/−中的表达变化，表明靶标去抑制程度相似，且双敲除和三敲除导致更多非靶标转录本改变。

亮点 2｜pi9和pi17 piRNA切割少数mRNA

在pi9−/−和pi17−/−初级精母细胞中，分别仅有7个和16个转录本显著改变（FDR<0.01），其中17个mRNA上调1.4–4.4倍。14个上调mRNA被鉴定为pi9或pi17 piRNA的直接切割靶标，其3′切割产物在对照中可检测，在相应突变体中减少≥8倍。体外重组MIWI实验证实这些piRNA可切割靶RNA。pi9−/−pi17−/−双敲除导致额外324个转录本改变，但仅Asb1同时受pi9和pi17 piRNA调控，表明不同piRNA簇很少共享靶标。

图 3　Fig. 3。图3展示piRNA切割对靶标丰度的影响有限。a左为具有可检测切割产物的pi9和pi17靶标在pi9−/−pi17−/−中的丰度变化，与对照piRNA靶标相比显著增加（KS检验P=0.002）；右为不考虑切割产物时，两组无显著差异（P=0.21）。b比较靶标丰度变化≤1.25倍和>1.25倍的切割效率、piRNA丰度和RNA PolII密度，显示低效切割和高转录速率是丰度不变的原因。c显示pachytene piRNA靶基因的PolII密度比非靶基因高约3倍。d显示有切割靶标的piRNA与无靶标piRNA的保守性评分无显著差异。

亮点 3｜大多数piRNA切割不影响靶标丰度

通过5′单磷酸RNA测序，鉴定出约112个pi9和pi17 piRNA的切割靶标，但其中位数靶标丰度仅增加4.9%。比较切割效率、piRNA丰度和靶标转录速率发现，低效切割和高转录速率共同解释了为何大多数切割事件不改变稳态mRNA水平。靶基因的RNA PolII密度（GRO-seq）比非靶基因高约3倍（中位数370 ppm vs 123 ppm，P=7×10−7），且切割效率低的靶标其piRNA丰度也较低。

图 4　Fig. 4。图4展示piRNA不通过miRNA或翻译调控机制作用。a显示具有种子互补性的潜在miRNA样靶标在pi9−/−pi17−/−中的丰度变化与对照无差异。b显示假定翻译激活靶标在pi9−/−pi17−/−中的翻译效率变化与对照无差异（KS检验P值均>0.05）。c为模型图，提出大多数piRNA无功能，少数切割靶标且改变丰度的piRNA维持精子适应性，从而确保整个piRNA库的保留。

亮点 4｜piRNA不通过miRNA或翻译调控机制作用

对具有种子互补性的潜在miRNA样靶标进行分析，未发现其在pi9−/−pi17−/−中上调。核糖体足迹测序显示，pi9或pi17 piRNA的假定翻译激活靶标与对照靶标在翻译效率变化上无显著差异（初级精母细胞P=0.64，次级精母细胞P=0.09，圆形精子细胞P=0.37）。此前报道的翻译激活靶标Tbpl1、Cnot4和Spesp1的核糖体占位或mRNA丰度均无变化。因此，粗线期piRNA既不激活也不抑制翻译。

亮点 5｜功能性piRNA快速演化且多源自转座子

在14个pi9或pi17 piRNA–靶标对中，仅5个在胎盘哺乳动物间保守，其余9个仅存在于鼠亚科。超过一半的改变靶标丰度的piRNA源自转座子序列，其中单个URR1B转座子衍生的piRNA调控4个mRNA。重复序列来源的piRNA在切割靶标但不改变丰度的类别中富集2.5倍。这表明piRNA介导的基因调控在进化中快速周转，转座子来源的piRNA在新靶标获取中起重要作用。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： CRISPR/Cas9基因编辑、RNA-seq、5′单磷酸RNA测序、核糖体足迹测序、GRO-seq、体外切割实验
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 研究piRNA、精子发生或小RNA调控机制的生信与实验人员。

带走一句

做piRNA功能分析时，不要默认所有piRNA都有靶标或功能；应结合切割产物、转录速率和piRNA丰度来评估实际调控效应，并关注少数关键piRNA的进化保守性。