两个猫尾草(Phleum pratense)栽培种时间群体在关键冬季胁迫下的等位基因频率与基因表达差异

《Theoretical and Applied Genetics》：Allele frequency and gene expression differences under key winter stresses in temporal populations of two timothy cultivars

摘 要：猫尾草（Phleum pratense L.）是挪威北部最重要的多年生牧草物种，该区域因气候变化预计将面临多变的冬季天气条件。了解猫尾草栽培种如何应对气候变化，对于保障高纬度地区牧草生产至关重要。本研究以两个适应北部的猫尾草栽培种Engmo（地方品种）和Noreng（合成栽培种）的时间群体（种子世代）为对象，探究冰冻（freezing）和冰覆盖（ice encasement）胁迫下基因表达变化及SNP等位基因频率（allele frequency）偏移。总体而言，随时间推移冰冻耐受性（以LT50即致死50%群体的温度定义）下降，而冰覆盖耐受性（以LD50即致死50%群体的持续时间定义）上升。比较转录组分析鉴定到数个已知参与胁迫响应的基因——包括乙烯响应转录因子（ethylene-responsive transcription factors, ERFs）、脱水响应元件结合转录因子（dehydration-responsive element binding transcription factors, DREBs）、乙烯敏感性逆转蛋白1（reversion to ethylene sensitivity 1, RTE1）及脱落酸抑制因子1（abscisic acid repressor 1, ABR1）——在Noreng时间群体间于冰冻胁迫下呈差异表达。若干发生大幅等位频率变化的位点与Noreng中随时代呈现表达偏移的基因紧密相邻。两栽培种在冰覆盖胁迫下群体间极少基因表达差异，可能源于种子繁育过程中对此性状的选择压较弱。两栽培种群体遗传多样性均随时代逐渐下降。结果表明，植物激素介导的转录调控可能是高纬度地区适应多变冬季气候条件的关键机制之一。研究结果强调在种子繁育过程中监测遗传偏移以维持栽培种稳定性之重要性，并提示所鉴定的胁迫响应基因可作为培育气候适应性牧草作物的潜在靶点。

论文解读：两个猫尾草时间群体在关键冬季胁迫下的等位基因频率与基因表达差异

研究背景与立题依据

猫尾草（Phleum pratense L.）是挪威北部最重要的多年生牧草。气候变化导致高纬度地区冬季温度变率增大、秋季变暖延长及冬季偶发冰雪覆盖事件增加，影响牧草冷驯化（cold acclimation, CA）与越冬存活。越冬持久性主要由冰冻低温、长期冰覆盖、水涝及雪霉等复合冬季胁迫决定，其中冰冻胁迫每年发生而冰覆盖发生频率较低但呈上升趋势。挪威农民近年反映猫尾草栽培种Engmo（1953年发布的地方品种，源自北挪威Troms郡）和Noreng（2002年发布的合成栽培种，含10份Engmo克隆＋4份南部地方品种Grindstad克隆）越冬死亡率升高、持久性下降，可能与种子多代扩繁过程中的遗传漂变或选择导致的遗传偏移（genetic shift）、管理变更及气候暖化有关。目前缺乏对栽培种在多代种子扩繁（时间群体/temporal populations）过程中抗寒相关基因表达演化及基因组等位频率动态的直接证据。为此研究人员采用"复活研究"（resurrection study）策略，比较Engmo（种子批1988年早代、2003年中代、2020年现世代）和Noreng（1998年早代、2010年中代、2020年现世代）三个时间群体的存活表型、冰冻与冰覆盖胁迫下转录组差异及基于SLAF-seq（specific locus amplified fragment sequencing）的SNP等位频率变化，以揭示栽培种对变化冬季气候的适应/漂变机制。该论文发表于《Theoretical and Applied Genetics》。

主要关键技术方法概述

研究人员获取猫尾草栽培种Engmo三时间群体（1988早代、2003中代、2020现世代）与Noreng三时间群体（1998早代、2010中代、2020现世代）的种子批，进行：① 盆栽苗冷驯化后分别开展梯度低温冰冻处理（测定LT₅₀，取-18℃＜LT₅₀为T1、-24℃＞LT₅₀为T2及CA对照）与梯度时长冰覆盖处理（测定LD₅₀，取36天＜LD₅₀为T1、68天＞LD₅₀为T2及CA对照），记录存活率并以广义线性模型估算半致死温度/时间；② 对冻/冰胁迫下各处理植株提取总RNA进行Illumina双端转录组测序，比对猫尾草草图基因组，用DESeq2鉴定时间群体间差异表达基因（DEGs, differentially expressed genes），阈值|log₂FC|≥log₂(1.2)、BH校正p＜0.05，并行GO富集分析；③ 对各群体混池幼苗组织提取基因组DNA，采用RsaI/HaeIII酶切建SLAF-seq文库进行Illumina双端测序，比对参考基因组后用grenedalf估算1 kb滑动窗口核苷酸多样性（θπ）、Watterson's θ（θw）、Tajima's D及Hudson's F_st（群体间遗传分化系数），筛选F_st＞0.15且含≥5 SNP窗口，计算早—现世代绝对等位频率差|ΔAF|，关联高|ΔAF| SNP邻近（±10 kb）基因与DEGs。

研究结果

Freezing and ice encasement tolerance in temporal populations（时间群体的冰冻与冰覆盖耐受性）

两栽培种均表现为从早代至现世代LT₅₀升高（冰冻耐受性下降），LD₅₀延长（冰覆盖耐受性上升）；Engmo中代与现世代在-21℃存活率显著低于早代，Noreng各时代间冰冻存活率无显著差异，仅Noreng中代54天冰覆盖存活率高于早代。说明时间推移伴随抗寒力衰减与耐冰覆盖力微增，但两栽培种表型演化模式不同。

RNAseq read alignment and functional annotation（RNAseq比对与功能注释）

冰冻处理样品比对率＞85%，冰覆盖样品因细菌污染比对率随处理时长下降；同源性比对注释约92%–93%基因。数据质量满足下游分析。

Differential gene expression under freezing stress（冰冻胁迫下差异基因表达）

PCA显示处理与栽培种分离明显；Noreng时间群体间DEG数多于Engmo，且随时代间距增大而增多（早代vs中代＜早代vs现世代＞中代vs现世代），栽培种间DEG数也随时间增加。Noreng中早代vs中代/现世代对比富集激素介导信号通路（GO:0009755）与转录调节活性（GO:0011040），其中DREB1、开花位点T（Flowering Locus T, FT）、光响应基因LHY（Late Elongated Hypocotyl）、UV-B受体UVR8等在现中代表达下调，而乙烯响应转录因子（ERF, ethylene-responsive transcription factor）、RTE1、ABA抑制因子ABR1（abscisic acid repressor 1）、冷响应蛋白激酶1（CRPK1, cold-responsive protein kinase 1）、NAC6及MADS-box 51上调，锌指LSD1（负调控细胞死亡）下调、NAC6（正调控程序性细胞死亡/programmed cell death, PCD）上调。表明Noreng冰冻胁迫转录程序发生方向性偏移，涉及昼夜节律–光温耦合、乙烯/ABA激素信号及PCD调控。

Differential gene expression under ice encasement（冰覆盖胁迫下差异基因表达）

时间群体间DEG未呈现方向性偏移模式，无显著GO富集，且重叠栽培种间DEG极少。推断冰覆盖胁迫下缺乏强选择压致随机表达漂变，未见适应性转录演化信号。

Allele frequency changes and genetic diversity estimates（等位频率变化与遗传多样性估算）

SNP-based PCA按栽培种背景、时代顺序分离；合成栽培种Noreng的θπ（核苷酸多样性）高于地方品种Engmo（p＜0.05）。两栽培种θπ、θw、Tajima's D均随时代下降。早—中、中—现、早—现世代间基因组F_st均约0.026–0.031。筛选得Engmo 303个、Noreng 268个高分化（F_st≥0.15）窗口；Noreng中14个与高F_st窗口重叠的DEGs（含皮质细胞界定蛋白编码基因下调＞2倍）及81个邻近大|ΔAF| SNP的胁迫响应基因（含DREB1B、ABR1、bHLH148-like）被鉴定，提示等位频率偏移可能与部分基因表达演化相关联。

讨论与结论翻译

Noreng时间群体间冰冻胁迫下DEG增多且具方向性表达模式（早代＞中代≥现世代或反向），部分DEG与栽培种间DEG高度重叠，说明栽培种间冰冻转录分歧主要来自Noreng演化偏移。Noreng中激素信号（乙烯/ABA途径ERF、RTE1、ABR1）、CRPK1–DREB1调控模块下调及PCD相关LSD1/NAC6表达重编程暗示植物激素介导转录调控是高纬度多变冬季适应关键。Engmo虽SNP等位频率变化明显却未见冰冻DEG方向性偏移，其LT₅₀下降或与不同纬度种子扩繁环境（北vs南）影响冷驯化相关。冰覆盖胁迫下无表达偏移符合弱田间选择压预期。基因组θπ下降反映纯化选择或奠基者效应致遗传侵蚀。综上，研究表明猫尾草栽培种在种子扩繁中出现方向性及发散性等位频率变化，Noreng经历冰冻胁迫转录程序演化，植物激素介导转录重调可能是适应变化冬季气候的机制之一；复活研究与混池SLAF-seq联用可有效解析作物对气候变化的适应演化，需在种子繁育中监控遗传偏移以保持栽培种稳定性，所鉴定胁迫响应基因为培育气候韧性牧草提供候选靶标。