加州牙鲆（Paralichthys californicus）假白化症相关基因失调的发生时间与发育阶段特征研究

《Journal of Fish Biology》：Nature and timing of dysregulation associated with pseudoalbinism in California Halibut, Paralichthys californicus

　　

摘要专业翻译 比目鱼类的假白化症（pseudoalbinism），即有眼侧体表大部分区域缺乏黑色素细胞（melanophores）或黑色素（melanin）的现象，在育苗场中经常被观察到。假白化症的发病率可能非常高，在加州牙鲆（Paralichthys californicus）中有时超过90%的个体受到影响。假白化症在变态发育（metamorphosis）之前无法观察到，但其发病率受幼虫期饲养条件的影响。本项目采用定量实时聚合酶链反应（quantitative real time polymerase chain reaction, qRT-PCR）方法，在加州牙鲆的四个发育阶段（F、G、H和I期）测量了八个目标基因（sox10a、sox10b、gch2、cyp26b1、aldh1a2、aldh1a3、asip1和LOC109629521）的相对表达量。研究人员对具有高假白化症发病率（94%、97%）和中等发病率（42%、42%）的牙鲆批次进行了采样。作为参考，还采集了具有低发病率（1%、2%）的南方鲆（Southern Flounder，Paralichthys lethostigma）。在G期观察到基因表达的显著差异，尤其是gch2、asip1、cyp26b1和sox10b。综合所有阶段并纳入南方鲆的低发病率数据后发现，高发病率牙鲆组的gch2表达出现急剧下降（>10倍），而低发病率鲆中变化很小。中等发病率组的gch2表达呈双峰分布，表明最终色素表型是在牙鲆标准体长（standard length, SL）为7.2–7.6 mm的狭窄发育窗口期内被决定的。观察到的基因表达模式提出了一种假说，即加州牙鲆的假白化症发生在变态中期至晚期（G和H期）黑色素母细胞（melanoblasts）向成熟黑色素细胞转变的过程中，可能是由于黑皮质素-1受体（melanocortin-1 receptor, MC1R）受到拮抗以及正常无眼侧发育通路的异位表达所致。gch2的表达可能成为筛查早期假白化症的候选标志物，但有待正式验证。

论文解读文章

一、研究背景与科学问题

比目鱼类在育苗场生产过程中，变态发育后常出现皮肤色素异常的现象，其中假白化症（pseudoalbinism）——即有眼侧体表大部分区域缺乏黑色素细胞或黑色素——是最常见的异常色素类型，在牙鲆属（Paralichthys）中尤为突出。在加州牙鲆（Paralichthys californicus）中，假白化症的发病率可超过90%，严重影响苗种质量和养殖经济效益。正常状态下，比目鱼类有眼侧皮肤深色色素沉着，含有黑色素细胞、黄色素细胞（xanthophores）、红色素细胞（erythrophores）和虹彩细胞（iridophores），而无眼侧仅有虹彩细胞。假白化症在变态发育之前无法被观察，但其发病率与幼虫期的饲养条件密切相关，如饵料、养殖缸颜色、光周期、过度拥挤等因素均有报道。此外，研究也表明异常色素存在遗传贡献，遗传力估计值范围较大（h2 = 0.16–0.98），且已在多个物种中鉴定出与假白化症和黑化过度（hypermelanosis）相关的遗传标记或位点。然而，假白化症发生的分子机制和表型命运被决定的具体发育时间窗口仍不清楚，这构成了开展本研究的主要动因。

二、研究目的与总体发现

研究人员基于对色素细胞分化通路的已有认识，提出了四个先验假说来解释假白化症的可能发生机制：其一，神经嵴干细胞（neural crest stem cells）未能分化为色素细胞祖细胞（chromatophore progenitor cells）；其二，视黄酸（retinoic acid, RA）信号通路失调，阻止祖细胞分化为黑色素母细胞；其三，cyp26b1过度表达（双侧而非单侧），阻止两侧黑色素母细胞分化为成熟黑色素细胞；其四，成熟黑色素细胞内部的黑色素生成（melanogenesis）过程受到干扰。通过定量实时聚合酶链反应（qRT-PCR）技术测量八个目标基因在变态前后四个发育阶段的表达变化，研究结果不支持假说一和假说二，而支持假白化症发生在变态中期至晚期（G和H期）黑色素母细胞向成熟黑色素细胞转变的阶段，并提出了MC1R受体拮抗与无眼侧发育通路异位表达协同作用的机制假说。

三、主要关键技术方法

本研究采用定量实时聚合酶链反应（qRT-PCR）方法测量八个目标基因（sox10a、sox10b、gch2、cyp26b1、aldh1a2、aldh1a3、asip1和LOC109629521）的相对表达量，以延伸因子1α（ef1a）为内参基因，采用ΔΔC_t方法并经效率校正计算校准归一化相对量（calibrated normalized relative quantities, CNRQ）。样本来源方面，加州牙鲆来自Hubbs-SeaWorld研究所（加利福尼亚州圣地亚哥），南方鲆（Paralichthys lethostigma）来自德克萨斯大学海洋科学研究所（德克萨斯州阿兰萨斯港）。数据分析采用多元置换方差分析（PERMANOVA）、高斯线性模型、高斯位置-离散模型、层次聚类分析和基因共表达网络分析等方法。

四、研究结果

4.1 内参基因ef1a的验证

ef1a的定量循环（Cq）值在四个发育阶段、两种饵料、两次产卵批次以及I期两种色素表型之间均保持稳定，牙鲆总体平均Cq值为15.4（变异系数CV = 2.1%），南方鲆为15.8（CV = 2.5%），验证了ef1a作为内参基因的可靠性。

4.2 最终色素表型比较（I期）

全局PERMANOVA分析显示，早期幼鱼（I期）正常与假白化个体之间的总体基因表达存在极显著差异（p = 0.0001）。单个基因分析表明，gch2、sox10b和LOC109629521三个基因的表达在两种表型之间存在显著差异。其中gch2差异最为突出，正常个体中位表达量是假白化个体的6.2倍（2.635 log₂ CNRQ单位），且两组分布几乎无重叠。sox10b在正常个体中表达量为假白化个体的1.9倍，而LOC109629521在假白化个体中表达量是正常个体的1.3倍。基因共表达网络分析显示，假白化个体中八个基因之间存在13个强相关性（|r| = 0.40–0.77），而正常个体中仅有6个，且正常个体中gch2与aldh1a3和aldh1a2的强相关性在假白化个体中完全消失。

4.3 变态前幼虫比较（F、G和H期）

F期：总体基因表达在发病率组之间无显著差异（p = 0.4137），但饵料效应显著（p = 0.002），仅sox10a和gch2两个基因的表达受饵料影响，差异较小。

G期：发病率对总体基因表达有显著影响（p = 0.0001），除aldh1a2外的所有基因平均表达在发病率组间均存在显著差异（adj. p < 0.021），其中gch2、asip1、cyp26b1和sox10b差异最为显著，高发病率组的表达量为中等组的1.5–2.7倍。gch2在高发病率组中的每个样本表达量均高于中等组的所有样本。cyp26b1、aldh1a3和asip1在高发病率组中表现出显著更大的离散度。高发病率组中除sox10b外的所有基因构成高度连通的共表达网络，cyp26b1、gch2和asip1之间呈正相关，而aldh1a3与这三个基因呈负相关。

H期：发病率对总体基因表达有显著影响（p = 0.0001），差异主要体现在均值而非方差。除sox10a和sox10b外，所有基因的表达在发病率组间均存在显著差异。gch2的差异尤为突出，中等发病率组的表达量是高发病率组的6.9倍。aldh1a3在中等组表达更高，而cyp26b1、aldh1a2、asip1和LOC109629521在高发病率组表达更高。

4.4 gch2表达的个体发育轨迹

三个发病率组（高、中、低）在最早采样的发育阶段（O_L = 78–88）具有相似的gch2表达水平，随后均上升至O_L = 85–87的最大值后迅速下降。低发病率组（南方鲆）在整个后续发育过程中呈现约5倍的缓慢稳定下降。高发病率组在O_L = 87–92期间出现>10倍（3.5 log₂ CNRQ单位）的急剧下降。中等发病率组的表达轨迹在高发病率组之后持续至O_L = 92，随后出现明显的双峰分布。层次聚类分析将O_L = 92–96的18个中等组幼虫分为两个主要聚类：7个幼虫具有较低gch2表达（与高发病率组水平相似），11个幼虫具有较高表达（与低发病率组相似）。低表达聚类样本的比例（39%）与中等组假白化症发病率（42%）高度吻合。

五、讨论与结论

研究结果表明，加州牙鲆假白化症并非源于神经嵴干细胞向色素细胞祖细胞分化失败（假说一）或RA信号通路阻止祖细胞向黑色素母细胞分化（假说二），因为高发病率组在F和G期的sox10a、sox10b和gch2表达并未降低，反而在G期更高。相反，基因表达模式提示假白化症发生在变态中期至晚期（G和H期）黑色素母细胞向成熟黑色素细胞转变的阶段。研究人员提出，假白化症是两个独立但协调的过程共同作用的结果：其一，asip1基因产物对黑皮质素-1受体（MC1R）的拮抗作用——高发病率组在G和H期asip1和LOC109629521表达显著升高，表明ASIP活性增强，抑制了MC1R信号通路对黑色素细胞成熟的促进作用；其二，正常无眼侧发育通路的异位表达——高发病率组cyp26b1表达升高且与aldh1a3呈强负相关，表明RA降解增强、RA合成减少，创造了不利于黑色素细胞分化的RA环境。cyp26b1、aldh1a3和asip1之间的强相关性提示RA降解、RA合成和MC1R拮抗可能构成一个协调的调控轴。

gch2表达的个体发育轨迹表明，与假白化症相关的基因失调最早在G期（O_L = 85–88，对应7.2–7.6 mm标准体长）可被检测到，而最终色素表型在O_L = 92–96（8.4–9.2 mm标准体长）时已被固定。这一发育窗口的精确界定为后续研究和干预措施的开发提供了重要基础。gch2表达可能成为筛查早期假白化症的候选标志物，但需进一步研究以确定其敏感性、特异性、预测价值和个体水平诊断准确性。研究人员进一步推测，多种环境胁迫因素可能通过激活应激轴（stress axis）汇聚于黑色素细胞分化过程，因为介导皮质醇合成和分泌的主要因子促肾上腺皮质激素（ACTH）在至少部分硬骨鱼类中可激活MC1R，异常色素可能是幼虫在关键发育窗口期经历应激的广义反应。未来研究建议采用转录组学方法、空间转录组学或显微解剖有眼侧和无眼侧组织的转录组分析，以及直接测量MC1R和RA信号通路活性，以进一步验证本研究的假说。