基于混合分组分离分析揭示的加州猛水蚤（Tigriopus californicus）低氧耐受性自然变异遗传位点

《Frontiers in Ecology and Evolution》：Genetic loci associated with natural variation in hypoxia tolerance in the copepod Tigriopus californicus revealed by bulked segregant analysis

　　

猛水蚤属（Tigriopus californicus）是近期被发现丢失了编码低氧诱导因子（HIF）通路主调控因子关键基因的甲壳动物之一，而该通路被认为对动物氧稳态至关重要。该物种栖息于潮上带水洼，每日承受数小时低氧压力，常达缺氧状态。尽管缺乏典型HIF转录因子，T. californicus仍表现出水生动物中记录到的最高水平低氧耐受性之一。缺乏HIF组分的类群中低氧响应的遗传机制几乎完全未知。研究人员利用T. californicus中耐受性的地理变异，以耐受性较强的加利福尼亚种群与耐受性较弱的俄勒冈种群重组杂交后代为材料识别数量性状位点（QTL）。在利用临界氧张力（Pcrit）对近800只杂交猛水蚤进行个体水平耐受性量化后，研究人员执行Pool-seq与混合分组分离分析（BSA）以识别QTL。分析在5条染色体上揭示7个QTL，共含43,648个SNP。QTL内SNP对应721个基因，富集于应激响应、糖代谢、神经元信号与线粒体功能，并包含外骨骼与膜修饰相关基因。QTL中许多基因在T. californicus中已知于实验性低氧事件下呈差异表达，表明解释种群间Pcrit差异的遗传变异部分关联于急性转录组响应所涉位点。QTL中多数SNP位于内含子或基因上游1000 bp内，显示QTL偏好潜在调控影响变异。具已知转录响应的QTL基因较无急性低氧响应基因在这些区域拥有更多SNP。注释为高影响（含终止丢失或终止获得变异）的SNP影响十余个基因，包括一重要抗氧化调控基因。结果为缺乏关键氧感应调控基因物种中低氧耐受自然变异的遗传架构提供首窥，并凸显可能贡献极端环境应激韧性的多样生物学通路候选基因。

研究背景与立项依据

潮上带猛水蚤Tigriopus californicus栖居于北美西岸飞溅池，夏季夜间因共生藻类耗氧可经历长达6小时完全缺氧，演化出极强低氧耐受性，但种群间耐受性存在地理分化：南加州圣地亚哥（SD）种群优于北俄勒冈Strawberry Hill（SH）种群。该类群属于Harpacticoida目，已丢失HIF-α主调控基因及EGLN脯氨酰羟化酶，传统依赖HIF通路感知氧化的范式在此失效，无HIF背景下低氧耐受的遗传架构未知。研究人员此前证实SD种群在模拟夜间缺氧—复氧序列中产生涉及糖酵解、抗氧化、外骨骼修饰、色素代谢与糖衍生物保护的多网络转录重编程，但调控这些网络的遗传位点未明。鉴于SD与SH基因组分歧大且杂交多代可育，研究人员以SH×SD晚期重组杂交系为材料，采用混合分组分离分析（BSA）定位与临界氧张力（P_crit）相关的数量性状位点（QTL），以解析无HIF物种低氧耐受自然变异的遗传基础。

主要技术方法概括

研究人员建立SH雌×SD雄F1及延续至12代以上的多线重组杂交种群，自由混交约10代后形成外缘混合培养系；以分段断棍法（broken-stick segmented method）重复测定个体P_crit，筛选雌雄各前10%与后10%极端表型个体构建低P_crit（高耐受）与高P_crit（低耐受）两个DNA混合池（各78只）；对两池进行全基因组Illumina测序，比对至SD×SH杂合伪参考，用QTLseqr包以G’统计与ΔSNP指数识别QTL（窗口750 kb，BH校正q<0.01），SnpEff注释SNP影响；将QTL基因与既往SD种群缺氧—复氧转录组差异表达（DE）基因集交集，用quasipoisson GLM比较DE与非DE基因SNP分布，topGO做GO富集，PAML YN00估算d_N/d_S，并按低P_crit池中SD或SH等位基因频率偏移将QTL划分为SD偏向与SH偏向分别统计。

研究结果

分布与QTL关联P_crit

G’分析在染色体3、8、9、11、12上检出7个显著QTL，含43,648 SNP；93.09%为modifier影响（内含子、上下游、非编码区），7%为low/moderate/high影响，high影响28个SNP涉及26基因（含剪接位点、终止获得/丢失）。过滤至基因内或上游1 kb得721基因。QTL整体GO富集于呼吸爆发调控、碳水化合物衍生物代谢、葡萄糖刺激细胞响应、一氧化氮代谢调控、蛋白定位至线粒体等。

与转录组响应关联

721 QTL基因中101个为既往缺氧DE基因，重叠无显著富集（Fisher OR=1.18, p=0.129）；DE基因较非DE基因每kb SNP数更高且由内含子SNP驱动，modifier类SNP更多，但d_N/d_S与最大G’无差。DE基因GO偏负调转运体活性、质膜组织、内吞循环；非DE偏上皮形态发生负调、钙离子螯合、糖转运正调、一氧化氮代谢、渗透压应激。

低P_crit亲本等位频率偏向

染色体3/11/12的5个QTL在低P_crit池SD等位频率更高（SD偏向，498/721基因），8/9号染色体2个QTL SH等位更高（SH偏向）。SD偏向基因内含子SNP密度与G’更高，d_N/d_S略高。SD偏向GO富集中枢为碳水脂类代谢、质膜组织、激素刺激响应、神经发育调控、转录调控、精子发生；SH偏向为氧化供能、超分子纤维组织、膜蛋白定位、tRNA加工、肌肉发育。

QTL中低氧响应功能组

糖代谢路径获5基因（含海藻糖转运体Tret1、磷酸戊糖路径Sedohheptulokinase编码基因Shpk、丙酮酸路径Akr1a1）；线粒体相关23基因（含mt-ARS Hars/WARS2、mt-RP mRpL10/mRpS26、翻译因子Mtif2、ATP合酶亚基ATPsynCf6及BDKHA、线粒体HSP Trap1）；胆碱通路获slc44a2（ΔSNP最高，18错义）、Cha、ACHE1、ChT、slc44a1；抗氧化仅GSTT1及具stop-lost的Cyp6a14（SH等位末端多2密码子）；外骨骼组11基因含Laccase2（lcc2，SD偏向QTL，3错义，多内含子/上游SNP）及几丁质合成GALK1/GFPT1（SH偏向）；mesh（septate junction）获1剪接接受体high影响SNP。

讨论与结论翻译

研究人员此前证该潮间带猛水蚤极低氧耐受且缺HIF-α仍可快速转录重编程。以两表型分歧种群多代杂交系，研究人员定位与P_crit相关7个位点跨5染色体，涉及细胞信号、调控过程及应激通路（胆碱能通量、角质层重塑、线粒体功能）。典型低氧调控因子HIF-α非种群间P_crit变异所必需；全基因组广泛补偿演化足以在无此关键基因下产生耐受表型变异。SD偏向QTL主导自然种群P_crit变异，但8/9号SH偏向QTL表明SH种群独立演化出有益等位，最低P_crit杂交体同时携带双亲高耐受等位，符合无基因流异域种群多基因局部适应。SNP偏内含子/上游、DE基因内含子SNP密度高、QTL基因d_N/d_S不高于背景，支持调控区分歧而非编码松弛选择主导P_crit变异，SD种群调控演化权重更大。线粒体基因富集未必由核质兼容选择解释（杂交超22代且SD偏向占69%）。胆碱—乙酰胆碱通量、角质层Laccase2介导氧扩散屏障、Cyp6a14终止丢失等构成无HIF背景替代网络。该研究首次勾勒缺HIF-α物种低氧耐受自然变异遗传架构，指明应激韧性候选通路。