编辑推荐：
本综述系统地阐述了在耐缺氧物种——鲤鱼脑内，首次大规模注释长链非编码RNA（lncRNA）的研究。通过高通量测序，鉴定出6072个lncRNAs，并发现其中1321个在缺氧（anoxia）与复氧（reoxygenation）过程中发生差异表达。研究重点揭示了这些差异表达的lncRNAs与邻近基因在表达变化上存在显著正相关，并主要通过顺式（cis）作用发挥潜在调控功能。这项研究为理解鲤鱼在极端缺氧环境下的转录组适应性调节提供了新的分子视角，并丰富了鱼类基因组的功能注释。
引言：探索极端生存背后的RNA调控者
氧气缺失是脊椎动物生命面临的严峻威胁，但鲤鱼（Carassius carassius）却能在低温缺氧环境中生存数月，这得益于其一系列生理适应，如代谢抑制、转向无氧代谢并产生乙醇作为终产物。尽管生理适应机制已被广泛研究，但在缺氧与复氧过程中，表观遗传机制，特别是长链非编码RNA（lncRNA）如何调控转录组变化，仍不清楚。lncRNA是长度超过200个核苷酸、不编码蛋白质的RNA分子，它们能以顺式（cis，作用于邻近基因）或反式（trans，作用于远端基因）方式，在转录、转录后及表观遗传层面精密调控基因表达，参与发育、代谢、应激反应等多种过程。本研究旨在填补这一空白，系统注释鲤鱼脑中的lncRNAs，并探究其在缺氧应激转录组响应中的潜在作用。
材料与方法：多管齐下的鉴定与严谨的分析策略
研究使用了两批鲤鱼脑组织样本，分别在常氧（normoxia）、缺氧（anoxia，持续7天）和复氧（reoxygenation，缺氧后恢复24小时）条件下处理，每条件n=6。通过Illumina HiSeq平台进行链特异性转录组测序。原始 reads 经质量控制和修剪后，使用STAR比对到鲤鱼参考基因组。采用StringTie进行转录本组装，共获得145,264个转录本。为了高置信度地鉴定lncRNAs，研究同时使用了四种具有不同算法原理的工具：CPC2（基于比对）、CNCI（基于相邻核苷酸三联体）、CPAT（基于逻辑回归模型）和FEELnc（基于随机森林模型）。最终，被这四种工具共同预测为非编码的6072个转录本被认定为lncRNAs。利用FEELnc的分类器模块预测了这些lncRNAs的潜在相互作用伙伴（Interaction Partners, IPs），并分析了它们的基因组上下文（genic/intergenic）和亚型。差异表达分析使用DESeq2进行，筛选条件为调整后p值<0.05且|log2(倍数变化)| >0.38。研究还分析了差异表达的lncRNAs（DElncRNAs）与其预测的相互作用伙伴在表达变化上的相关性，并对显著变化的伙伴基因进行了基因本体（GO）富集分析。
结果
1. 长链非编码转录本的鉴定与差异表达
主成分分析（PCA）显示，不同处理组（常氧、缺氧、复氧）的样本能清晰区分。从组装的转录本中，四工具一致鉴定出6072个lncRNAs（）。在总共56,440个差异表达转录本（DETs）中，包含了1321个DElncRNAs。从常氧到缺氧（N to A）的转变引发了最强烈的转录组响应，DElncRNAs数量最多（1321个），其中上调和下调的数量相当。从常氧到复氧（N to R）的变化也很大，但从缺氧到复氧（A to R）的变化则最少，表明lncRNAs在进入缺氧状态时调整最为活跃。维恩图分析显示，不同比较组间的DElncRNAs存在部分重叠与特有集合，提示它们在不同生理阶段可能扮演特定角色。
2. 差异表达LncRNAs及其邻近基因的分类
DElncRNAs的长度分布与所有lncRNAs相似，中位长度为1672 bp。FEELnc为大多数lncRNAs预测了相互作用伙伴，其中约60%的lncRNAs位于基因间区（intergenic），40%位于基因区内（genic）（）。在基因间区，大多数lncRNAs与伙伴基因位于同一条链上（same strand），提示共转录调控。在基因区内，最常见的亚型是“嵌套型”（nested，lncRNA完全位于伙伴基因内部），且多数位于反义链。lncRNAs与其最佳预测伙伴之间的中位距离约为2 kb，许多作用发生在50 kb范围内，这强烈暗示这些相互作用主要是顺式（cis）调控。
3. 潜在相互作用伙伴的特性与差异表达
对DElncRNAs的预测相互作用伙伴进行表达聚类分析，揭示了在不同氧条件下具有明显表达模式的基因簇（）。在291个差异表达的相互作用伙伴（DEIPs）中，变化趋势与DElncRNAs类似，即N to A和N to R对比中变化最多。至关重要的是，DElncRNAs与其预测伙伴基因的log2(倍数变化)在三个对比（N to A, N to R, A to R）中均呈现显著正相关（Pearson相关系数分别为0.35, 0.57, 0.35）（）。这支持了DElncRNAs可能与其邻近基因共表达或对其进行调控的假设。对顶级变化的DElncRNAs-伙伴对的深入分析显示了三类关系：正相关（如STRG8312.1与跨膜蛋白TM275在缺氧时均上调）、无显著变化（如STRG.15579.1上调而其伙伴ZN729不变）以及负相关（如STRG.9408.4在缺氧时大幅上调而其伙伴SUFU下调）。
讨论与结论
本研究首次在耐缺氧的鲤鱼脑内进行了大规模的lncRNA注释，发现了6072个新型lncRNAs。与鲤科其他鱼类相比，本研究鉴定的lncRNAs数量较少但长度更长。在缺氧应激下，大量lncRNAs发生差异表达，尤其是在进入缺氧状态时，表明它们可能参与了转录组重编程以适应低氧环境。这些DElncRNAs在基因组上多位于基因间区，且与邻近的伙伴基因倾向于共表达，强烈提示其通过顺式作用机制调节附近基因的表达，从而在鲤鱼应对缺氧的分子适应中扮演重要角色。虽然具体的调控机制（如直接激活/抑制）有待未来功能实验验证，但观察到的表达相关性为理解lncRNAs在极端生理条件下的调控网络奠定了坚实基础。这项工作不仅扩展了鲤鱼基因组的功能注释，也为探索脊椎动物缺氧耐受的进化与分子机制提供了新的线索。
