《Biology》：Comparative Transcriptome Profiling Reveals Developmental Regulatory Networks in Rainbow Trout Eggs and Hatched Larvae
编辑推荐：
摘要：虹鳟（Oncorhynchus mykiss）是全球重要的冷水养殖物种，但从胚胎向自由生活仔鱼转变的分子事件仍未充分解析。本研究采用真核有参转录组测序，比较了 O. mykiss 眼卵阶段与刚孵化仔鱼的全基因组基因表达谱。共鉴定出618个差异表达基因（d
摘要：虹鳟（Oncorhynchus mykiss）是全球重要的冷水养殖物种，但从胚胎向自由生活仔鱼转变的分子事件仍未充分解析。本研究采用真核有参转录组测序，比较了 O. mykiss 眼卵阶段与刚孵化仔鱼的全基因组基因表达谱。共鉴定出618个差异表达基因（differentially expressed genes, DEGs），其中297个上调、321个下调。基因本体（Gene Ontology, GO）富集分析显示，DEGs显著富集于免疫系统过程、刺激响应和生物调节。在免疫系统过程类别中，仔鱼有4个DEGs上调、11个下调，表明免疫相关基因表达重塑较单纯激活更复杂。京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）富集分析仅鉴定出糖酵解/糖异生（Glycolysis/Gluconeogenesis, map00010）显著富集（padj < 0.05）。氮代谢（Nitrogen metabolism, map00910）和色氨酸代谢（tryptophan metabolism, map00380）显示富集趋势，但经多重检验校正后无统计学显著性。由于这些通路包含的DEGs数量少且调控方向混合，未推断通路水平的整体上调或下调。这些结果提示，在从卵黄依赖合成代谢向仔鱼分解利用转变过程中可能存在代谢重塑。转录本组装进一步鉴定出134,770个潜在新异构体和2902个基因间区转录本，归类为推定新转录本，丰富了 O. mykiss 基因组注释。趋势聚类表明，免疫和应激响应基因协同调控，同时氨基酸和能量代谢基因下调，揭示了该发育过渡中的程序性时序调控。这些发现为理解虹鳟早期发育提供了基础数据，并为健康导向的水产养殖管理提供了候选分子靶标。
论文解读文章
鱼类胚胎发育是决定个体存活和种群补充的关键生命阶段，其有序推进依赖合子基因组激活、胚层分化、组织形态发生和器官原基形成等高度保守且时序耦合的细胞与分子事件。模式物种斑马鱼（Danio rerio）研究已表明，母源因子与合子基因组之间的精确时序转换是正常胚胎发生的前提，关键节点失调可导致发育停滞或表型异常；代谢微环境和小分子代谢物也紧密调控孵化与仔鱼生长。然而，虹鳟（Oncorhynchus mykiss）作为具有重要商业价值的鲑科鱼类，其早期发育的精细分子知识仍有限。近年来，高通量测序推动转录组学成为解析非模式鱼类胚胎发育分子基础的核心手段，可从基因表达动态、表观修饰和信号通路调控等层面系统识别关键基因及共调控网络。虹鳟是全球广泛养殖的冷水鱼，经济价值高、养殖产量大，且具有较高质量参考基因组；其胚胎在低温下孵化期较长。研究中，在16 °C恒温孵化条件下，孵化约发生在受精后24天（days post-fertilization, dpf），约384度日。从胚胎（卵）向刚孵化仔鱼的转变具有重要生物学意义，涉及能量代谢从卵黄源营养分解向主动摄食和消化外源食物的根本转变，并伴随消化、感觉和运动系统功能成熟。现有鲑科转录组研究多集中于幼鱼和成鱼的应激耐受、养殖性能及病原响应，对早期胚胎和仔鱼发育全过程转录动态关注有限；尤其是卵与刚孵化仔鱼直接比较、免疫激活与代谢重塑协调调控尚未在虹鳟中系统表征。因此，研究人员提出两个假设：孵化伴随内源免疫相关基因表达重塑，反映从母源保护向自主防御转变；该转变涉及从卵黄依赖合成储存向分解利用的代谢转变。研究目标是鉴定阶段特异DEGs、富集通路和新转录本，为虹鳟早期发育提供转录组框架。
研究人员采用真核有参转录组测序比较 O. mykiss 眼卵与刚孵化仔鱼的全基因组表达谱。共鉴定618个DEGs，297个上调、321个下调。GO富集显示DEGs显著富集于免疫系统过程、刺激响应和生物调节；免疫系统过程类别中4个上调、11个下调。KEGG富集仅糖酵解/糖异生显著，氮代谢和色氨酸代谢呈趋势但校正后不显著。转录组装鉴定134,770个潜在新异构体和2902个基因间区转录本为推定新转录本，丰富了基因组注释。趋势聚类显示免疫和应激响应基因协同调控，同时氨基酸和能量代谢基因下调，揭示程序性时序调控。该研究为理解虹鳟早期发育提供基础数据，并为健康导向的水产养殖管理提供候选分子靶标。
关键方法概括（不超过250字）：样本来自中国水产科学研究院东海水产研究所养殖基地。研究人员对 O. mykiss 进行三倍体诱导，16 °C恒温孵化，取眼卵（24 dpf）和新孵化仔鱼（31 dpf，约7 dph）样本，每组3个生物学重复。采用真核有参转录组测序，NovaSeq X Plus平台PE150。fastp质控，HISAT2比对 O. mykiss 参考基因组，StringTie组装并用Class code预测新转录本，DIAMOND进行NR、Swiss-Prot、GO、KEGG、EggNOG和Pfam注释；RSEM计算TPM（Transcripts Per Million，每百万转录本数）定量，DESeq2鉴定DEGs，GOATOOLS进行GO富集，KOBAS进行KEGG富集，并进行GSEA（Gene Set Enrichment Analysis，基因集富集分析）及层次聚类和趋势聚类。
研究结果如下。
3.1 测序数据质量控制与参考基因组比对：通过质控和比对分析，所有样本Q30均高于92%，测序错误率低；比对到 O. mykiss 参考基因组的总体比对率为94.97%–95.57%，唯一比对 reads 占87.25%–87.66%，多比对 reads 占7.48%–8.31%，表明数据质量高、外源污染少，参考基因组选择可靠。
3.2 转录本组装与新转录本预测：通过StringTie组装并与参考注释比较，多数组装转录本含1–10个外显子，多数基因有1–3个转录本异构体。新转录本预测鉴定134,770个潜在新异构体（class code“j”），其中2902个为基因间区转录本（class code“u”），定义为推定新转录本；另有134,606个与已知转录本完全匹配，“o”类1170个。结果表明 O. mykiss 发育过程中存在既往未注释的转录事件。
3.3 功能注释概览：研究人员对全部44,665个在所有样本中可检测到转录信号的基因进行NR、SwissProt、GO、KEGG、EggNOG和Pfam注释；其中25,955个基因满足组内平均TPM ≥ 1，定义为活跃表达基因，用于后续样本关系分析。
3.4 表达分布与样本关系：通过箱线图、密度曲线、Pearson相关、PCA（Principal Component Analysis，主成分分析）和维恩分析，发现组内Pearson相关系数均高于0.95，重复性好；组间相关略低但仍较高；PCA显示两组分离，转录组谱可区分；维恩图显示共有和特有基因数，支持组间差异。
3.5 差异表达基因鉴定：以B组仔鱼对A组眼卵进行比较，使用DESeq2，阈值|log
2(fold change)| ≥ 1且p
adj < 0.05，共鉴定618个DEGs，其中297个上调、321个下调；火山图和热图显示代表性表达差异。
3.6 DEGs的GO功能富集分析：使用GOATOOLS进行GO富集（p
adj < 0.05）。生物学过程（Biological Process, BP）中，DEGs显著富集于免疫系统过程、生物调节、稳态过程、生殖过程、细胞过程、发育过程、多细胞生物过程、定位和刺激响应；细胞组分（Cellular Component, CC）中富集膜界细胞器和细胞骨架；分子功能（Molecular Function, MF）中蛋白结合和信号受体结合突出。免疫系统过程中4个上调、11个下调，提示免疫防御、应激感知和细胞内组织发生重编程。
3.7 DEGs的KEGG通路富集分析：在218个映射通路中，仅糖酵解/糖异生（map00010）显著富集（p
adj < 0.05）。氮代谢（map00910）和色氨酸代谢（map00380）显示富集趋势，但多重检验校正后不显著；三者均属代谢类别。通路内DEGs呈上调和下调混合模式，且每个通路DEGs数量有限，因此未推断通路整体上调或下调。结果提示卵到仔鱼转变中能量和氮代谢可能重塑。
基因集富集分析：为补充基于阈值的DEG富集并在全转录组水平检测协调变化，研究人员以GO术语和KEGG通路为预定义基因集进行GSEA，以B组对A组的信号噪声指标排序，1000次置换，FDR < 0.05。GO基因集中，A组眼卵显著富集的基因集主要关联合成和代谢过程，如糖异生、脂蛋白代谢、己糖生物合成、血红蛋白复合物形成、胶原三聚体组装和活性氧代谢；B组仔鱼显著富集的主要关联视觉系统发育和神经功能，如眼晶状体结构成分、晶状体发育、视觉感知、突触后膜特化、谷氨酸能突触传递和神经系统发育。KEGG通路中，A组有9个代谢通路显著富集，包括糖酵解/糖异生、淀粉和蔗糖代谢、丙酮酸代谢、磷酸戊糖途径和卟啉代谢；B组无KEGG通路达到FDR < 0.05，仅黑色素生成、心肌肾上腺素信号和Wnt信号等呈弱富集趋势。
3.8 DEGs聚类分析：通过趋势聚类分析，618个DEGs被分为10个表达亚簇。亚簇2包含266个基因，从眼卵到仔鱼表达增加，代表早期虹鳟发育中程序性时序表达模式；EggNOG功能分类显示该亚簇基因的直系同源类别分布。聚类结果为识别发育阶段特异核心调控网络提供线索。
讨论部分总结：本研究通过转录组测序系统比较 O. mykiss 卵和刚孵化仔鱼。DEGs显著富集于免疫系统过程、刺激响应和生物调节，其中免疫系统过程4个上调、11个下调，提示免疫相关基因协同重塑而非简单单向激活。胚胎主要依赖母源抗体和卵膜内非特异免疫因子保护；孵化后仔鱼进入外部水环境，面临病原微生物、温度波动和渗透挑战，因此免疫重塑对存活关键。11个下调、4个上调提示免疫系统建立涉及精细稳态调节，而非单纯激活。KEGG中糖酵解/糖异生显著，但GSEA显示该通路负NES，即眼卵A组整体转录活性更高；DEGs混合调控，不能推断仔鱼中简单激活。氮代谢和色氨酸代谢趋势不显著，且DEGs少，应谨慎解释。鱼卵富含卵黄，储存蛋白质、脂质和糖原，胚胎期合成代谢占优；孵化后外源摄食尚未完全建立，内部卵黄储备继续消耗，逐渐转向分解利用。本研究中7 dph仔鱼仍有可见残留卵黄囊，且未提供外源饲料，首次外源摄食预计在约10–14 dph，因此7 dph代谢富集反映内源卵黄依赖阶段。代谢重塑为后续外源摄食建立基础。新转录本方面，134,770个潜在新异构体和2902个基因间区转录本提示 O. mykiss 基因组可能存在大量既往未注释转录单元，可能包括长非编码RNA，但未进行功能验证。能量代谢与转录调控关联方面，持续上调和下调趋势子集主导，持续上调基因主要富集免疫和应激响应，持续下调基因主要集中于氨基酸和能量代谢，反映资源分配策略和发育重编程。DEG比例约1.4%，与鲑科及其他硬骨鱼类似发育转变报道的1%–5%范围一致。局限性包括：样本为三倍体虹鳟，外推到二倍体需谨慎；仅两个时间点，不能区分孵化直接触发的短暂响应与程序性发育变化；未进行基因敲除或药理抑制等功能验证；新转录本未实验验证；全胚胎/仔鱼样本不能反映器官特异表达；实验室条件结果对自然种群适用性需验证；全卵和全仔鱼组织组成不同，部分差异可能反映新细胞类型出现而非发育重编程本身。
研究结论翻译：总之，本研究通过RNA-seq分析系统揭示了 O. mykiss 从卵向刚孵化仔鱼发育转变伴随的转录组重编程特征。差异表达基因主要关联两个核心模块：免疫相关基因表达重塑和代谢通路变化。KEGG富集分析仅鉴定出糖酵解/糖异生显著富集，而氮代谢和色氨酸代谢显示富集趋势，但多重检验校正后无统计学显著性。免疫系统过程DEGs中4个上调、11个下调，表明复杂重塑而非简单激活。由于这些代谢通路仅含少量DEGs且呈混合调控，未推断通路水平整体方向。这些发现提示该发育转变可能涉及从合成储存向分解利用的转变，并伴随免疫相关基因复杂调节。推定新转录本的发现丰富了 O. mykiss 基因组注释，为未来发育生物学研究提供候选分子资源。尽管因果关系需实验验证，本研究为理解硬骨鱼早期存活策略的分子基础提供了有用线索。未来研究应整合时间采样、组织特异转录组和功能验证，以进一步解析调控 O. mykiss 早期发育的关键基因网络。本研究对现有鲑科发育知识有三方面贡献：建立虹鳟眼卵和刚孵化仔鱼阶段的全基因组转录组基线资源；系统表征驱动孵化期免疫重塑和代谢重编程的关键DEGs和功能通路，揭示两个核心生物过程间的协调模式；提供完整DEG列表作为未来鲑科比较和功能发育研究的参考数据。
