《International Journal of Molecular Sciences》：Engineering MSC Migration: Roles of Nanoparticles in Activating Migratory Pathways and Functions
Temuulen Batsaikhan,
Hyun Su Lee and
Young Joon Seo
编辑推荐：
这篇研究通过构建大黄鱼的高密度SNP遗传连锁图谱，系统解析了生长性状（体重、形态）与内脏白结节病抗性性状的遗传关系。研究发现，控制生长与抗性的数量性状位点(QTL)在很大程度上相互独立，仅存在有限的共定位区间，表明通过育种同步改良这两个重要性状是可行的，为该鱼种的平衡育种策略提供了重要的分子基础。
1. 引言
大黄鱼 (Larimichthys crocea) 是中国东南沿海重要的海水养殖鱼类，2022年全国年产量达25.7万吨。然而，由杀香鱼假单胞菌 (Pseudomonas plecoglossicida) 引起的内脏白结节病 (Visceral white nodules disease, VWND) 对该产业构成严重威胁，死亡率可超过80%，造成巨大经济损失。由于药物治疗效果有限，培育抗病品系成为一种更安全、可持续的解决方案。生长性能（包括产量性状和形态性状）和抗病力是核心育种性状，两者之间的遗传关系（拮抗或协同）是育种者关注的关键。传统育种方法如表型选择、家系选择等在世代周期长、表型可塑性高的大黄鱼中面临挑战。分子标记技术，特别是高密度遗传图谱和QTL (Quantitative Trait Locus) 定位，为解析这些复杂性状的遗传基础提供了更精确的工具。尽管已有关于大黄鱼遗传图谱和单一性状QTL定位的研究，但针对生长与VWND抗性间遗传关系的系统解析仍很缺乏。本研究旨在通过构建高密度SNP遗传图谱，对产量、形态和VWND抗性性状进行全面QTL定位，以阐明其遗传关系，为同步改良生长与抗病的育种策略提供分子依据。
2. 结果
2.1. 生长与VWND抗性性状的表型变异与正态性分析
对150尾大黄鱼F1个体的12个性状（9个生长相关，3个抗性相关）进行了表型测定。结果表明，大部分形态性状（如总长AL、体长BL等）及抗性性状中的存活时间(AT)近似服从正态分布。而产量性状(Wt)和两个病原载量性状（肝PPLL、脾PPSL）则呈现轻微偏态。所有性状均表现出良好的遗传变异，变异系数范围为8.92% (尾高TH) 至 87.50% (脾PPSL)，为后续分析奠定了基础。
2.2. 产量、形态与VWND抗性性状之间的表型相关性
相关性分析显示，所有形态性状之间以及形态性状与产量性状(Wt)之间均存在强正相关(r=0.73~0.99)。在抗性性状内部，存活时间(AT)与脾病原载量(PPSL)显著正相关(r=0.67)，而肝病原载量(PPLL)与两者相关性弱。值得注意的是，产量和所有形态性状与存活时间(AT)呈中等正相关(r=0.45~0.52)，与脾病原载量(PPSL)呈弱正相关，但与肝病原载量(PPLL)无显著相关。这表明生长表现优异的个体在攻毒后可能存活更久，脾脏病原载量也略高，但肝脏感染水平与生长表型无关。
2.3. 连锁分析与遗传图谱构建
本研究基于全基因组重测序和分箱标记，构建了包含24个连锁群(Linkage Groups, LGs)的高密度遗传图谱，对应大黄鱼的24条染色体。图谱共整合了32,429个分箱标记，总遗传距离为1417.17 cM，平均标记间隔仅为0.051 cM，表明图谱密度高、覆盖全面。
2.4. 生长相关性状的QTL定位
QTL定位在16个连锁群上共检测到53个显著的生长相关QTLs (LOD ≥ 3.0)，涉及产量性状(Wt)和总长(AL)、体长(BL)、体高(CH)等形态性状。Wt检测到12个QTLs，表型贡献率(PVE)在0.14%至5.83%之间。多个生长性状QTL在基因组特定区域（如LG6、LG16）存在共定位，表明这些区域对多个生长性状有协同调控作用。例如，LG6上的一个区间(20.4–28.4 Mb)包含了对Wt、AL、BL、CH都有显著效应的QTLs。LG16上的一个区间(15.84–16.23 Mb)甚至与所有9个被检测的生长性状QTL重叠，表明这是一个影响整体体型的基因组热点。
2.5. VWND抗性性状的QTL定位
针对三个抗性性状（存活时间AT，肝病原载量PPLL，脾病原载量PPSL）共鉴定出20个显著QTLs。其中，AT与5个QTL相关，PVE范围2.69%–8.93%；PPSL和PPLL分别与8个和7个QTL相关。研究发现，抗性性状QTL在某些基因组区域也存在共定位，例如LG9包含影响所有三个抗性性状的QTL区间，LG16和LG22也存在调控多个抗性性状的QTL簇，表明不同抗性成分可能受共同或连锁基因调控。
2.6. 生长与疾病抗性的QTL共定位
为了探究生长与抗性间的遗传关系，分析了两类性状QTL的共定位情况。结果显示，在全部53个生长相关QTL和20个抗性相关QTL中，仅发现两个基因组区间存在生长与抗性QTL的共定位。第一个位于LG16 (21.36–21.41 Mb)，被总长(AL)和存活时间(AT)共享。第二个位于LG22 (1.48–1.68 Mb)，被体长(BL)和存活时间(AT)共享。这两个区间对表型的贡献率均不高(0.66%–5.99%)。尤为关键的是，效应最大的生长QTL（主要在LG6、LG16）与效应最大的抗性QTL（主要在LG9、LG22）分别位于不同的连锁群。这表明，控制这两个重要经济性状的主要遗传位点在很大程度上是相互独立的。
3. 讨论
本研究构建的高密度遗传图谱实现了对大黄鱼生长和VWND抗性性状的精细QTL定位。研究发现，生长性状内部（产量与形态）遗传整合度高，而生长与抗性性状之间在基因组水平上主要表现为独立的遗传控制。有限的QTL共定位、候选基因功能（LG16上的轴突导向受体基因Unc5d和LG22上的钠离子通道基因SCN5A、DNA损伤检查点组分基因HUS1）均涉及基础细胞过程而非核心的生长或免疫通路，以及主要效应QTL的分离，共同支持了这一结论。这种遗传结构表明，生长性状与VWND抗性之间不存在强烈的遗传拮抗，这对于大黄鱼育种具有重要实践意义。育种者可以相对独立地针对高产性状和高抗性进行选择，有望在选育过程中同步改良这两个重要性状，克服传统育种中可能存在的“此消彼长”的瓶颈。本研究为开展旨在同时提高大黄鱼生长性能和抗病力的分子标记辅助育种提供了有价值的遗传信息和理论基础。
4. 材料与方法
本研究使用的实验鱼来自“岱衢”品系大黄鱼。通过人工授精建立了5个全同胞家系，从中选取一个家系的150尾8月龄健康个体进行后续研究。测定了1个产量性状（体重Wt）和8个形态性状。通过腹腔注射杀香鱼假单胞菌进行攻毒，测定3个VWND抗性相关性状：存活时间(AT)、肝(PPLL)和脾(PPSL)病原菌相对载量。对150尾个体进行了全基因组重测序，通过分箱标记过滤策略获得高质量SNP标记。利用Lepmap3软件，采用双亲假测交策略构建了高密度遗传连锁图谱。基于此图谱和F1代表型数据，使用MapQTL 6.0软件的区间作图法(Interval Mapping)对12个性状进行了QTL定位。统计分析使用R和SPSS软件完成。
5. 结论
本研究为大规模、系统地解析大黄鱼生长与VWND抗性的遗传结构奠定了基础。构建的高密度SNP遗传图谱实现了对两类性状QTL的精准定位。结果表明，生长与抗性性状主要由各自独立的QTLs控制，遗传相关性较弱，仅在两个基因组区间存在效应有限的共定位QTL。这种独立的遗传结构表明，在育种中同时提高大黄鱼的产量、改善体型和增强对VWND的抗性是可行的，不存在严重的遗传拮抗障碍。鉴定出的QTL热点区域和候选基因为未来的分子标记辅助选择和功能验证研究提供了明确靶点，有助于加速培育生长与抗病力兼备的优良大黄鱼新品种。
