《Animals》：Advances in Poultry RNA-Omics Research: Technologies, RNA Information Layers, and Applications in Complex Traits
编辑推荐：
RNA组学技术已将家禽转录组研究拓展至组织水平基因丰度之外。这些技术现在能够解析完整转录本、RNA化学修饰、细胞来源和空间位置。本综述围绕一个核心问题：每种RNA组学技术究竟提供了哪些此前无法获取的特异性RNA信息？根据其首要测量靶标，研究人员将这些技术分为四
RNA组学技术已将家禽转录组研究拓展至组织水平基因丰度之外。这些技术现在能够解析完整转录本、RNA化学修饰、细胞来源和空间位置。本综述围绕一个核心问题：每种RNA组学技术究竟提供了哪些此前无法获取的特异性RNA信息？根据其首要测量靶标，研究人员将这些技术分为四类：用于分析基因表达和调控RNA网络的批量转录组及小RNA测序；用于解析转录本异构体的长读长和直接RNA测序；用于检测RNA修饰的表观转录组学方法；以及用于鉴定信号细胞来源和空间分布的单细胞及空间转录组学。解析RNA结构、RNA–蛋白质相互作用和翻译状态的方法（如结构探测、CLIP型定位和核糖体印迹）在家禽中仍处于起步阶段，导致输卵管、产脂肝脏和有核红细胞等组织中的翻译及RNA–蛋白质调控在很大程度上未被测量。针对每一类别，研究人员讨论了其直接测量靶标和主要输出、相较于既有方法所提供的额外信息、在家禽中的代表性应用及其固有局限性。在此框架基础上，研究人员探讨了如何根据具体研究问题选择和组合技术，并总结了这些技术在生产性能与产品品质、繁殖、健康与抗逆性研究中的应用。不同技术测量不同类型的RNA特征，不能简单互换。未来研究必须更加重视技术与科学问题的匹配、互补数据的整合，并夯实三个基础：家禽转录本和非编码RNA的注释、标准化分析流程以及功能验证体系。
论文主体部分总结如下：
3. RNA组学技术及其揭示的RNA信息层次
本文所讨论的技术提供互补的RNA信息，而非单向的简单到复杂的技术路径。短读长RNA测序可量化基因和外显子水平丰度，并识别未注释外显子和剪接位点，但片段化读段难以可靠重构全长转录本，组织水平信号不可避免受细胞组成影响。小RNA测序和全转录组测序将分析范围扩展至miRNA、lncRNA和circRNA，捕获常规基因表达分析常遗漏的调控RNA信息。长读长测序可解析完整或近完整转录本异构体、可变剪接和3'端使用。基于抗体富集的方法和天然RNA分析（如直接RNA测序）可用于绘制m6A等RNA修饰。单细胞和空间转录组学进一步确定表达信号的细胞来源和空间分布。RNA结构、RNA–蛋白质相互作用和翻译状态的信息目前仍处于家禽RNA组学研究前沿，应用有限。
3.1 批量转录组与调控RNA分析
短读长RNA测序仍是家禽RNA组学中应用最广泛的技术。通过修饰文库制备策略，分析靶标可扩展至各类非编码RNA：小RNA测序检测成熟miRNA，总RNA测序或circRNA富集方案用于鉴定circRNA。poly(A)富集主要捕获成熟mRNA和带poly(A)尾的lncRNA，但会遗漏大多数miRNA和许多circRNA；大小选择的小RNA文库可检测成熟miRNA但无法捕获其初级转录本。circRNA通常通过rRNA去除或RNase R外切酶富集环形分子后，基于反向剪接位点读段进行鉴定。基因水平分析聚合来自同一基因组位点的读段，较为稳定；转录本水平分析需将读段分配给不同异构体，高度依赖注释完整性、读段长度和异构体间序列差异。基于丰度估计可进一步进行差异表达、共表达模块和ceRNA网络分析。ceRNA模型最初在哺乳动物研究中提出，通常基于lncRNA、miRNA、circRNA和mRNA之间的序列关系和表达相关性构建。然而，标准差异表达分析仅证明RNA丰度在特定条件下变化，共表达分析反映基因间同步变化，二者均不能独立证明调控关系或直接物理相互作用。ceRNA关系基于结合位点预测和表达相关性，并非直接测量分子结合。批量组织的表达差异可能源于细胞类型比例变化，而非特定细胞类型内的转录调控。计算反卷积可部分校正，但不能完全消除不确定性。此外，lncRNA、circRNA和eRNA注释仍不稳定，缺乏标准化命名和跨版本映射，限制跨研究比较。
3.2 长读长与直接RNA测序
长读长测序将分析焦点从总基因表达转向特定转录本和异构体组成。PacBio Iso-Seq生成高精度共有读段，适合解析转录本结构；ONT cDNA测序可产生较长读段且通量和成本较优；ONT直接RNA测序通过纳米孔读取天然RNA分子，无需cDNA中间体，保留逆转录过程中可能丢失的分子信号。不同平台在通量、成本、单碱基准确度、定量可靠性和天然RNA信息之间需权衡。长读段不能保证覆盖每个RNA分子的全长；RNA降解、不完全逆转录和覆盖偏差均可导致截短转录本，“全长转录组”仍是近似描述。联合分析中，长读段用于定义转录本结构架构，高深度短读段用于丰度估计和剪接位点支持。ONT对19个鸡组织的测序鉴定出数万个异构体，其中许多不存在于现有注释中，异构体分类准确性高于短读段方法。鸭多组织PacBio Iso-Seq建立了全长参考转录组，揭示了数千个可变剪接事件和lncRNA。异构体分辨率可彻底改变性状相关表达变化的解释：基因水平“无变化”的基因可能发生异构体转换，影响编码能力或调控特征。但发现未注释异构体不等于证明其参与表型调控。直接RNA测序还可保留与RNA修饰相关的天然分子信号，但修饰判定依赖离子电流信号和计算模型，需独立方法验证。
3.3 表观转录组技术
家禽表观转录组研究主要集中于N6-甲基腺苷（m6A）。MeRIP-seq利用抗体富集含m6A的RNA片段，结果通常表现为修饰富集峰而非精确单碱基位点，每个峰可能跨越数十至数百个核苷酸。当多个异构体共享同一区域时，峰信号通常无法分配至特定转录本。MeRIP-seq必须与input RNA-seq文库联合解读，因为免疫沉淀信号受修饰水平和转录本丰度共同影响。高分辨率或单碱基方法可缩小候选位点范围，靶向实验可验证特定位点。天然RNA纳米孔测序通过离子电流偏移或碱基判定偏差推断RNA修饰，但结果依赖训练模型和算法，不同工具可能产生高度不一致的候选位点。抗体富集峰、高分辨率候选位点、独立验证位点和纳米孔模型推断信号应作为不同证据类型分别解读。目前家禽m6A研究涉及产蛋过程、胚胎性腺发育、排卵前卵泡选择、脂肪沉积差异等。许多研究测量METTL3和FTO等修饰酶的表达，但酶水平变化不能直接证明特定转录本修饰状态改变。目前最清晰的功能证据来自肌肉研究：METTL3促进circSIK2翻译并增强鸡成肌细胞增殖；METTL3依赖性GHR mRNA m6A修饰调节成肌细胞分化过程中的线粒体生物发生。该领域面临两个主要约束：峰水平MeRIP-seq分辨率不足以区分异构体或等位基因特异性修饰；m6A阅读蛋白对特定转录本命运的调控大多推断自哺乳动物研究，禽类RNA–蛋白质结合直接证据缺乏。
3.4 单细胞与空间转录组学
单细胞RNA测序（scRNA-seq）可记录单个细胞的表达谱；空间转录组学保留空间斑点、组织区域或单细胞的坐标信息。scRNA-seq和单核RNA测序（snRNA-seq）分析不同RNA库：液滴法scRNA-seq捕获细胞质RNA良好，但组织解离可能导致大而脆弱的细胞丢失；snRNA-seq分析核RNA，更适合冷冻组织及肌肉、脂肪等难以获得完整单细胞悬液的样本，通常产生更高比例的内含子和lncRNA读段。伪时间分析基于转录相似性排序细胞以映射潜在状态转变，但不是真实发育时间的直接记录。基于配体–受体数据库的细胞通讯分析仅提出潜在信号相互作用，不能证明实际通讯发生。标准液滴法scRNA-seq通常具有3'或5'端偏差，限制在基因水平分析。基于空间斑点的平台常在一个斑点内捕获多个细胞，需计算反卷积。细胞类型注释是家禽研究的另一挑战：目前常借用哺乳动物标记基因，但某些家禽细胞缺乏物种验证标记。禽类成熟红细胞有核且转录活跃，在液滴法scRNA-seq中产生大量血红蛋白相关读段，消耗测序通量并干扰聚类。荧光激活分选可最彻底去除红细胞，但会损失活性并诱导应激反应转录本；禽类特异性红细胞裂解速度较快但裂解不可预测；单核测序可绕过完整红细胞，适用于纤维化和脂质丰富组织，但仅读取核RNA。空间转录组平台分为两类：基于测序的空间条形码方法提供全转录组信息但分辨率受斑点限制；基于成像的原位杂交方法提供近单细胞分辨率但通常限于预先设计的基因面板。家禽细胞通讯推断严重依赖基于人和小鼠数据构建的配体–受体数据库，但禽类免疫配体和受体库与哺乳动物差异显著，直接应用易产生假阳性或假阴性推断。当前尚无经过整理的禽类配体–受体数据库。
3.5 RNA结构、RNA–蛋白质相互作用与翻译
转录组范围的化学探测方法（如DMS-seq和icSHAPE）通过测量核苷酸可及性研究RNA结构及其动态，但在家禽中尚无全转录组应用报道。RNA结合蛋白（RBP）对剪接、稳定性、定位和翻译调控不可或缺，但家禽缺乏系统性全转录组RNA–蛋白质相互作用图谱。使用CLIP或eCLIP直接定位家禽阅读蛋白结合位点的研究极为罕见，主要原因是缺乏能有效免疫沉淀靶蛋白的抗体。核糖体印迹（Ribo-seq）通过测序核糖体保护片段评估翻译状态和效率，在家禽中已有探索性应用，涉及鸡脾脏、生殖组织、脑组织、胚胎组织和原代成肌细胞。该技术要求严格控制核糖体状态保存、核酸酶消化、rRNA去除和物种特异性读段映射。核糖体印迹还可通过检测核糖体占用的短开放阅读框（sORF）发现功能性微肽。一项鸡研究结合RNA-seq和Ribo-seq鉴定了一个由lncRNA编码的74个氨基酸微肽，证明其调节成肌细胞增殖和分化。总体而言，这些方法在家禽中的系统数据仍匮乏，当前家禽RNA组学证据严重偏向表达丰度、异构体、修饰和细胞来源，结构与蛋白质结合层面有待拓展。
4. 家禽RNA组学技术的选择与组合
4.1 基于研究问题选择技术
RNA组学技术的选择取决于五个相互关联的因素：待解答的生物学问题、必须测量的RNA特征、可用的生物材料和采样条件、验证预期结论所需的方法，以及相对于所获信息的成本。技术适用性不取决于其新颖性，而取决于其测量读数是否直接回答研究问题，以及其固有局限性是否在实验设计和数据解读中得到充分处理。RNA丰度、转录本结构、化学修饰、细胞来源、空间位置、RNA–蛋白质相互作用和翻译状态均需要专门的测量方法，单一平台很少能同时包含所有信息。成本方面，批量短读长测序是最低层级；长读长测序和MeRIP-seq属中间层级；液滴法单细胞、单核平台和空间转录组学属最高层级。昂贵层级无法在个体水平性状推断所需样本量下应用。务实的设计是将昂贵层级保留给需要它的特定问题，用少量代表性动物，同时通过批量测序和反卷积维持整个队列的重复对比。
4.2 家禽中现有和潜在的技术组合
组合技术的根本目的是用不同方法填补互补信息空白。当前组合策略可分为三类：直接整合研究（在同一或配对样本上测量两种或多种RNA层）、跨数据集比较（独立生成的数据层相互解读）、以及技术可行但尚未实施的设计。MeRIP-seq标准分析需要配对的input RNA-seq文库以区分修饰富集变化与转录本丰度变化，但这应视为分析对照而非两个独立组学层的整合。ChickenGTEx项目汇集了基因型、基因表达、剪接和3'UTR使用数据，通过eQTL和sQTL分析将RNA特征与遗传变异关联。其组织面板、品种、年龄和生理状态来自贡献数据集而非设计采样框架，且效应具有背景依赖性，因此目录化关联应视为待检验的先验。短读段与长读段数据整合在方法学上已成熟：长读段定义转录本结构，更深短读段提供可靠定量和剪接位点支持，但在家禽中主要用于转录注释和资源构建。scRNA-seq与空间转录组学配对、单细胞与长读段异构体分析整合、CLIP/eCLIP与表达或修饰数据配对等设计在哺乳动物中已建立，但家禽实施证据极为有限。
4.3 整合在何种条件下加强证据
技术组合是否真正加强解读取决于其是否引入与问题相关的新信息。第一是生物学和样本对应性：不同测量应聚焦相同生物学对比，最好来自配对个体或匹配样本等分。第二是测量互补性：不同技术必须提供解决生物学难题所需的不同信息，如果两种测定最终仅提供基因水平丰度变化，结果可能只是数据量增加而非新RNA信息。第三是验证方法与结论一致：异构体结论需要异构体特异性检测或干扰；RNA修饰需要位点或区域水平验证；细胞类型结论需要细胞分辨定位或功能测定；RNA–蛋白质相互作用需要直接结合测定；生产或疾病性状结论必须与个体水平表型关联。当前家禽RNA组学文献中，独立生成的单层数据集多于来自同一或配对样本的多层数据。当务之急是围绕同一生物学问题设计互补测量，优先验证多组学整合提出的具体分子关系。
5. RNA组学在家禽复杂性状研究中的应用
5.1 生产与产品品质
肌肉生长、肌纤维类型、肉品质、脂肪沉积、饲料效率和蛋品质是家禽RNA组学中研究密集的生产性状。批量转录组和非编码RNA测序最初将研究从单个候选基因扩展至调控网络。circRNA、miRNA和mRNA的发育阶段联合分析已用于构建肌生成相关网络；lncRNA和circRNA谱可区分氧化型与糖酵解型肌纤维。脂肪沉积研究提出了一系列以候选ncRNA为中心的网络，涉及肝脏脂质沉积、脂肪肝、腹脂和肌内与腹脂差异。这些研究的常见输出是随性状或发育阶段共变的RNA集合，基于表达相关性和序列靶标预测组装，因此单个候选仍是提名而非已确立的控制点。肌肉研究显示了新层所增加的信息：表观转录组学报告转录本修饰状态；单核转录组学鉴定产生批量信号的细胞群；空间转录组学将木胸肌病起始定位到特定血管周围巨噬细胞生态位。证据在不同生产性状间分布不均，肌肉生长和脂肪沉积已有批量表达、ncRNA、修饰和细胞分辨数据，蛋成分、蛋壳质量、肉品质指标和饲料效率仍主要依赖少数品种的批量或ncRNA测序。从培养细胞验证的效应转化为活体动物生产性能的确认驱动因子仍是核心差距。
5.2 繁殖
家禽繁殖取决于生殖轴的发育和动态调控。开产日龄反映下丘脑–垂体–卵巢轴的发育轨迹，产蛋率和持续性衡量成熟轴的动态功能，卵泡选择取决于卵巢中重复的细胞命运决定，受精率和孵化率受卵母细胞质量和早期胚胎发育影响。目前家禽繁殖RNA组学严重依赖组织平均表达、调控RNA网络和m6A修饰。全转录组测序描绘了抱窝母鸡HPO轴和卵巢闭锁期间协调的ncRNA变化，但无法解析不同细胞群的贡献。m6A动力学与高、低产蛋鸡差异相关，特定m6A标记基因参与排卵前卵泡选择。长读段数据稀少，生殖轴的异构体动态基本未绘制，单细胞和空间证据有限。由于多数研究使用全卵巢、卵泡或输卵管样本，与产蛋性能相关的表达变化无法分配到颗粒细胞、膜细胞、生殖细胞、基质或驻留免疫细胞。卵泡组成随成熟显著变化，组织水平差异混合了组成变化和细胞内调控。miR-202-5p调节鹅颗粒细胞类固醇生成提供了坚实的细胞水平功能证据，但不直接代表个体鸟的体内产蛋率。
5.3 健康与抗逆性
解读家禽健康RNA组学研究时，必须区分感染或应激诱导的分子反应与抗性、耐受性和恢复力等个体性状。多数文献描述感染、环境应激或干预后的RNA变化。这类研究描述宿主如何反应，但前后测量很少能确定RNA特征是否决定个体先天抗病能力。单细胞转录组学将感染反应解析至特定细胞群，主要优势是区分细胞组成变化与细胞内调控。挑战、药物处理和免疫逃避研究捕获特定条件下的反应，而非明确的抗病性状。在相同挑战条件下对比抗病和易感个体能更直接评估RNA特征是否决定结局差异。马立克氏病抗性和易感鸡的比较将免疫细胞组成和状态与抗性表型联系起来；脾脏单细胞研究揭示了沙门氏菌感染的品种特异性细胞反应。抗性设计区别于反应设计之处在于随RNA测量一起记录的表型。耐受力要求在给定病原载量下解读损伤，恢复力需要通过恢复过程的重复测量。病原载量、存活、病变严重度和生产损失描述宿主结局的不同方面，并不总是一致变化。家禽健康RNA组学并不缺少感染后或应激后表达数据，缺少的是能分离上述健康性状的设计：遗传定义的抗性和易感群体、标准化配对挑战、精确结局表型、多个感染后时间点、控制细胞组成变化，以及基因型–环境互作下的纵向观察。
6. 家禽RNA组学的共同局限与优先方向
6.1 注释与参考资源
参考资源不足是家禽RNA组学的基础瓶颈。家禽转录本注释不完整，非编码RNA注释相对稀疏，异构体命名和跨版本映射缺乏稳定性。不同数据集中的“新”转录本无法可靠匹配。转录坐标和基因模型在不同基因组组装和注释版本间常发生变化，GC-rich区域和微染色体尤其易出现映射错误。长读段和lncRNA研究经常发现参考注释中缺失的转录本，但没有通用标识符和独立数据集中的重复检测，无法确定其是否为稳定生物实体。注释不能在不同家禽物种间无验证地转移。RefSeq和Ensembl对同一鸡组装的注释不同，基因模型、转录本集和标识符并非一一对应，定量结果不可互换。T2T无缺口基因组和禽类泛基因组对解决微染色体和GC-rich区域转录组映射问题尤为重要。
6.2 实验设计与平台偏差
样本处理、平台特性和研究设计共同影响RNA组学结果，但这些因素在家禽研究中有时未被充分报告。RNA完整性、保存方法和预处理决定下游检测保真度。scRNA-seq和snRNA-seq分析不同RNA区室，平台差异必须明确纳入方法选择和跨研究比较。组织解离可诱导应激反应基因表达；对于肌肉等难解离或纤维化组织，snRNA-seq通常是更优选择。MeRIP-seq结果对抗体批次和峰值调用参数敏感；长读段测序定量能力有限，异构体目录随过滤标准和分析流程波动；纳米孔修饰调用模型通常针对特定物种训练，应用前需严格验证。一些数据集生物学重复不足，表型记录粗略，以品种或发育阶段代替个体测量性状。技术重复、合并样本或来自单个动物的大量细胞有时被错误地视为独立生物学重复，低估个体间方差并人为膨胀统计显著性。已发表基准表明，每条件48个生物学重复的实验中，每组3个重复仅恢复20–40%的全集差异表达基因；当小倍数变化重要时至少需要6个重复，最好12个。家禽设计中常见每组3只鸟会遗漏许多真实差异。生物学重复是鸟本身，来自一只动物的细胞、技术重复和合并样本必须如此建模。禽类有核红细胞需要在样本制备中严格去除。禽类ZZ/ZW系统缺乏染色体范围剂量补偿，Z连锁基因在雄性（ZZ）中表达水平通常高于雌性（ZW）。混合性别队列中不将性别作为协变量会导致大量Z染色体基因出现假阳性。性别应记录并为每个样本建模，Z连锁和常染色体结果应分开报告。
6.3 计算分析与结果解读
ceRNA和调控网络预测、伪时间和细胞通讯分析、纳米孔修饰调用、长读段异构体鉴定和跨研究整合均依赖计算推断。其输出对工具选择、参数和参考版本敏感，不同方法常产生不一致的候选列表。可重复性取决于报告软件版本、参数、参考基因组和注释版本，以及原始数据和分析代码的可获取性。解读应区分直接读数与计算推断，充分报告分析参数和参考资源，并将未验证的关系标记为候选。计算输出不应承载超过其支持数据的生物学权重。
6.4 家禽RNA组学的优先方向
三个优先方向：参考与测量资源方面，建立完整家禽全长转录本注释、稳定跨版本非编码RNA和异构体标识符，构建物种和组织特异性细胞图谱及标准化鸡细胞类型本体，开发家禽特异性基准数据集并严格评估纳米孔修饰调用模型，系统建立RNA结构、RNA–蛋白质相互作用和翻译资源。整合与可重复性方面，从设计阶段就围绕同一生物学问题规划多层采样和测定，使用匹配样本等分、多组织谱或纵向采样，保持稳健生物学重复，不同数据层严格遵循统一参考基因组、注释版本和处理流程，完整记录软件版本和参数，公开原始数据和分析代码。机制与性状验证方面，下游实验必须靶向组学分析提出的具体RNA特征或分子关系，采用CRISPR干扰、异构体特异性干预、位点或区域水平修饰测定和RBP结合实验，在适当的细胞系、类器官或肠类器官中实施，并将细胞和分子读数与个体水平表型紧密关联。向分子育种的转化方面，eQTL、sQTL和3'UTR使用QTL目录可作为功能先验标记基因组区域，TWAS可将顺式调控基因表达与性状关联，贝叶斯或加权基因组预测模型可能提高预测准确性。但此类整合能否真正提升育种价值最终取决于家禽表达资源的规模和质量，以及在目标商业群体中的独立验证。调节信号的稳定性、基因型–环境互作和独立验证这三个条件目前均未满足。RNA特征与性状的相关性不是育种工具，只有遗传组分已被估计且效应在目标商业群体中验证过的变异或分子表型才能进入育种框架。家禽体内敲除模型极难生成，依赖原始生殖细胞操作，因此禽类类器官或肠类器官的正交功能验证不可或缺。最终目标不仅是阐明机制，而是将RNA特征作为分子表型，通过eQTL和sQTL整合入基因组选择系统，成为分子育种的转化工具。
7. 结论
家禽RNA组学已从分析组织平均丰度和调控RNA景观发展到解析转录本异构体、RNA修饰、细胞来源和空间分布。不同技术针对不同RNA特征且适合根本不同的生物学问题，因此多个数据层间的一致性本身不能证明因果关系。准确解读要求严格区分直接测量的RNA特征、计算推断的候选关系和实验验证的功能效应。未来研究必须针对具体生物学问题选择和组合技术，持续改进注释资源和分析流程，采用与结论相符的验证方法，在复杂RNA特征与精确定义的家禽性状之间建立可靠因果联系。对于育种应用，这意味着将RNA特征视为遗传组分已在目标群体中估计和验证的分子表型，而非等待重新解读的相关性。
