《Molecular Ecology》：Ontogenetic Sequence of Differential Gene Expression in Predator-Induced Daphnia pulex
编辑推荐：
许多大型溞（Daphnia）物种在捕食者存在时会发育出诱导型形态防御。这些性状具有表型可塑性，即在无捕食者时其产生受到抑制。尽管先前研究已鉴定出多种参与防御发育的神经体液因子和候选基因，但特定基因和因子的参与仍不确定，且关于发育过程中基因表达变化时序的信息极为
许多大型溞（Daphnia）物种在捕食者存在时会发育出诱导型形态防御。这些性状具有表型可塑性，即在无捕食者时其产生受到抑制。尽管先前研究已鉴定出多种参与防御发育的神经体液因子和候选基因，但特定基因和因子的参与仍不确定，且关于发育过程中基因表达变化时序的信息极为有限。因此，本研究对大型溞（D. pulex）幼体发育阶段防御发育进行了非候选基因依赖的表达分析。研究人员分析了微甲壳动物对捕食者释放的天然混合化合物以及近期在实验室鉴定的纯化信息素（kairomone）的转录组响应。所得数据表明，与大型溞全局基因表达模式相关的主要因素是蜕皮周期。信息素处理对全局基因表达的影响仅在特定阶段显现。不同阶段通常表现出独特的信息素诱导基因表达模式。然而，对于自然释放与化学合成的信息素，在各阶段的表达谱具有相似性。多种具有调控、结构及解毒功能的基因在信息素处理下差异表达。最一致的响应发现于编码胰岛素样肽（insulin-like peptide）的ilp-3表达水平。基因敲低实验提示该激素参与防御产物的形成。许多响应信息素处理的基因尚无预测功能，强调了对大型溞基因功能进行研究的必要性。
大型溞（Daphnia）作为淡水生态系统中的关键浮游动物，在面对捕食压力时展现出多样化的表型可塑性防御机制，包括形态结构的改变如头棘（neckteeth）、头盔（helmets）和脊状突起（crests）等。其中，由幽蚊（Chaoborus）幼虫诱导产生的头棘是研究最为深入的系统之一。尽管先前研究通过候选基因方法鉴定了部分神经体液因子和相关基因，但防御发育的完整分子机制，特别是基因表达变化的时空动态，仍远未阐明。此外，早期研究多基于混合样本的RNA测序，难以捕捉个体间变异和精细的发育阶段特征。
为系统解析大型溞防御发育的分子基础，研究人员设计了一项精细的时序转录组研究。研究采用无偏倚的方法，聚焦于从信息素感知到效应衰退的关键发育窗口，选取四个时间点：诱导后18小时的预释放胚胎（e18）、释放后30小时的一龄幼体（jI30）、42小时的一龄幼体（jI42）以及56小时的二龄幼体（jII56）。为区分信息素本身与其他捕食者伴生化合物的效应，研究设置了对照、天然信息素混合物和化学合成纯信息素（N-亚油酸-L-谷氨酰胺）三组处理。每个处理-阶段组合设置六个生物学重复，以单个母体的同胞 clutch 为单位，确保样本间发育同步性。
研究材料来源于加拿大安大略省分离的R9克隆系，该克隆属于北美或泛北极大型溞。实验前通过七代严格单性生殖培养以减少个体间变异。信息素敏感期的确定基于早期形态学分期标准。样本采集后使用Tag-Seq技术进行数字基因表达谱分析，该技术较传统RNA-Seq在差异表达分析中具有更高的统计效力。
研究发现，样本聚类的主要驱动因素是发育阶段而非信息素处理，蜕皮周期是解释样本间变异的最主要因素。jI42阶段呈现异质性，可细分为两个亚群（jI42.1和jI42.2），分别对应蜕皮前后的不同发育时相。仅jII56阶段在聚类中显示出处理间的明显分离，合成信息素处理样本与部分天然信息素处理样本聚为一类，区别于对照组。
差异基因表达分析采用两阶段检验策略，以最大化检测对两种信息素均有响应的基因。结果共鉴定出313个在至少一个阶段有强差异表达证据的基因。这些差异表达基因（DEGs）中，绝大多数为上调表达，且在不同阶段间呈现一致的表达变化方向。值得注意的是，后蜕皮阶段（jI30和jII56）的差异表达基因数量显著多于前蜕皮阶段（e18和jI42.2），提示信息素响应存在蜕皮周期依赖性。
在跨阶段持续差异表达的基因中，ilp-3编码的胰岛素样肽3（ILP-3）最为突出。该基因在胚胎阶段响应最强，随发育逐渐减弱，至二龄幼体转为下调。系统发育分析表明，ILP-3属于节肢动物胰岛素样生长因子（aIGF）家族，大型溞基因组中存在ilp-1、ilp-3和ilp-4三个旁系同源基因，源于两次基因复制事件。研究人员通过双链RNA干扰（dsRNAi）技术对ilp-3进行功能验证，结果显示基因敲低导致一龄幼体头棘表达量降低64.3%，二龄幼体降低32.3%，证实ILP-3参与防御结构的形成调控。同时，细胞因子信号抑制因子2（SOCS2）同源基因和尿苷二磷酸糖基转移酶基因ugt209B1也呈现类似的早期响应、逐渐衰退的表达模式。
另一类具有晚期表达特征的差异基因包括：血小板衍生生长因子/血管内皮生长因子（PDGF/VEGF）相关生长因子、几丁质结合蛋白以及细胞色素P450基因cyp4AP7。尤为特殊的是一个长末端重复序列（LTR）转座元件（KRLE），该元件在信息素诱导的幼体中高度表达，含有gag多蛋白和逆转录胃蛋白酶（retropepsin）两个开放阅读框，但缺乏pol基因，属于枝角类特有的Ty3/Gypsy型元件。这是首例报道的水生甲壳动物中差异表达的LTR元件。
胚胎期特异性响应的基因包括：编码latrophilin家族蛋白的cirl基因和编码胆碱转运体的chT1基因，两者均与神经体液调节相关。此外，含BTB结构域和螺旋-转角-螺旋（HTH）DNA结合结构域的转录调控因子XM_046591911.1~GN，以及小核仁蛋白IMP3编码基因imp-3也在该阶段差异表达。
基因集富集和过表达分析揭示，含几丁质结合结构域的蛋白、蛋白酶（样）结构域蛋白及昆虫表皮蛋白家族是最显著富集的功能类别。这些蛋白的表达呈现蜕皮周期依赖性波动，部分上调、部分下调，反映了表皮结构的深度重塑而非简单的组分线性增加。
研究讨论部分指出，这是首次对化学合成信息素进行转录组响应分析的研究，证实了脂肪酸-L-谷氨酰胺共轭物是触发大型溞抗捕食反应的活性成分。与先前基于混合样本的研究相比，单个clutch的分析策略虽然牺牲了部分测序深度，但能够评估个体间变异并识别发育阶段的精细结构。研究强调了将精细发育视角纳入环境胁迫遗传响应研究的重要性。
研究结论部分翻译如下：大型溞基因表达模式在发育过程中高度动态，许多参与信息素转录组响应的基因亦遵循蜕皮周期，导致不同蜕皮时相差异表达基因数量的偏倚。因此，在信息素响应基因中发现大量编码表皮组分的基因。这类基因的差异化表达具有阶段依赖性，表明表皮经历深度重塑而非其组分的线性增加。ILP-3的持续差异上调是幼体阶段头棘产生的最强相关因子，为胰岛素样信号参与抗捕食响应提供了证据。尽管敲低结果提示ILP-3可能参与头棘形成的调控，但该基因在属内的保守性表明其在大型溞中负责信息素的系统性响应，并具有其他生理功能。SOCS2与ILP-3信号的关系仍有待研究。胚胎期差异表达的短列表富含神经体液因子，包括与乙酰胆碱生物合成相关的chT1和与多巴胺信号相关的cirl。 unexpectedly发现的LTR元件KRLE是首个报道的水生甲壳动物中差异表达的LTR元件，其近期丢失pol开放阅读框的特征提示其可能通过劫持相关元件的逆转录 machinery 而保持移动性。葡萄糖-6-磷酸酶上调和Methuselah样黏附G蛋白偶联受体（aGPCR）下调可能与生命史变化相关。许多决定头棘产生的细胞决定因子在未被诱导的大型溞中已存在，因此信息素化学感受器等通路组分可能组成型表达而不在信息素暴露后改变表达水平。Tag-Seq方法的局限性在于需要增加测序深度和样本量以扩展差异表达基因列表，而组织特异性变化需单细胞转录组技术解决。随着基因表达数据集的积累，特别是本次研究所提供的数据，将有助于更好地理解大型溞的基因功能与相互作用。
