《Genes》：Pan-Genome-Wide Identification, Evolutionary Analysis, and Stress-Responsive Expression Profiling of the Maize DUF668 Gene Family
编辑推荐：
背景/目的：含有未知功能域668（Domain of Unknown Function 668, DUF668）的植物蛋白已被证实参与生长、发育及非生物胁迫响应，但该家族在玉米中的群体水平多样性和潜在功能仍知之甚少。因此，研究人员在玉米泛基因组范围内对该家族进
背景/目的：含有未知功能域668（Domain of Unknown Function 668, DUF668）的植物蛋白已被证实参与生长、发育及非生物胁迫响应，但该家族在玉米中的群体水平多样性和潜在功能仍知之甚少。因此，研究人员在玉米泛基因组范围内对该家族进行了表征，并检测了其在非生物胁迫下的表达情况。方法：研究人员对从26个玉米基因组中鉴定的DUF668基因进行了存在-缺失变异（Presence-Absence Variation, PAV）、系统发育关系、结构变异、功能注释及预测功能关联评估。利用公共RNA-seq数据及实时荧光定量PCR（quantitative real-time PCR, qRT-PCR）技术，分析了B73基因组驻留基因在涝渍、温度、渗透、盐度和紫外-B（UV-B）胁迫下的表达谱。结果：研究人员共鉴定出51个非冗余位点，包括8个核心位点、5个近核心位点、9个非必需位点和29个私有位点。对69个玉米、拟南芥和水稻蛋白的系统发育分析将其分为两个亚家族。大多数位点表现出纯化选择的证据，而ZmDUF668-16具有较高的非同义替换率与同义替换率比值（Ka/Ks）；在ZmDUF668-4和ZmDUF668-15中检测到大的插入或缺失。B73基因组驻留的19个位点表现出独特的发育、时序及胁迫响应表达模式。ZmDUF668-1也位于宽泛的2号染色体L-ABA/root-ABA1区间内。结论：玉米DUF668家族在泛基因组、进化和转录水平上均表现出显著的变异。这些发现为优先选择候选基因进行玉米发育及非生物胁迫响应的功能验证奠定了基础。
研究背景与科学问题
植物蛋白中广泛存在一类被称为“未知功能域”（Domain of Unknown Function, DUF）的保守序列。尽管其名称暗示功能未知，但越来越多的证据表明，许多含DUF的蛋白在植物生长发育及逆境适应中扮演关键角色。其中，DUF668家族是一个植物特异性的高度保守蛋白家族，先前的研究已在多个物种中揭示了其潜在重要性。然而，对于重要粮食作物玉米而言，对该家族的认识主要局限于单一参考基因组（B73）的分析，这极大地限制了对该家族遗传多样性、进化动态及功能潜力的全面理解。传统的单参考基因组分析方法无法有效捕捉不同种质间的存在-缺失变异（PAV）和结构变异（SV），而这些变异正是群体遗传多样性的重要来源。因此，为了突破这一瓶颈，深入挖掘玉米DUF668基因家族的群体遗传学特征及其在非生物胁迫响应中的潜在功能，研究人员利用玉米泛基因组资源开展了此项研究。
研究方法概述
研究人员首先利用隐马尔可夫模型（HMM）从26个高质量玉米基因组组装中鉴定出所有DUF668基因，并根据泛基因分配方案将其归类为非冗余位点。随后，采用最大似然法（Maximum Likelihood）构建了包含玉米、拟南芥和水稻DUF668蛋白的系统发育树。通过计算非同义替换率与同义替换率的比值（Ka/Ks）评估进化压力，并利用ANNOVAR软件注释与基因区域重叠的结构变异（SV）。基因表达分析整合了公共RNA-seq数据集（来源于NCBI GEO和SRA数据库，如GSE50191等）以及自主进行的实时荧光定量PCR（qRT-PCR）实验，所用植物材料为玉米自交系B73。
主要研究结果
3.1 泛基因组范围的ZmDUF668鉴定
研究人员在26个玉米基因组中鉴定出51个非冗余的ZmDUF668泛基因组位点，每个品系的位点数在13至21个不等。根据在不同品系中的出现频率，这些位点被划分为8个核心位点、5个近核心位点、9个非必需位点和29个私有位点，揭示了该基因家族在群体水平上显著的PAV。
3.2 玉米DUF668基因的泛基因组分析
蛋白长度比较显示，部分成员如ZmDUF668-10在不同品系间存在显著的长度变异。系统发育分析将来自三个物种的69个蛋白分为两个亚家族，且每个亚家族均包含玉米中不同PAV类别的成员。
3.3 ZmDUF668经历不同的选择压力
Ka/Ks比值分析表明，大部分ZmDUF668基因受到纯化选择（Ka/Ks < 1），功能趋于保守。然而，ZmDUF668-16是一个明显的例外，其在多个品系中表现出升高的Ka/Ks比值，暗示其经历了放松的纯化选择或正向选择，序列分化程度更高。
3.4 结构变异与玉米ZmDUF668基因表达的关联
通过比较分析，研究人员发现两个结构变异（SV）与基因表达差异显著相关。在ZmDUF668-4的启动子区发现一个25.2 kb的插入，其表达量在携带该SV的品系（Il14H）中高于不携带的品系（B73）。而在ZmDUF668-15的3‘侧翼区发现一个7.8 kb的插入，其表达趋势则相反。尽管编码区和保守基序在多数品系中保持保守，但这些调控区的SV与表达变化相偶联。
3.5 B73驻留ZmDUF668位点与玉米ABA相关QTL的对应关系
通过与已知的脱落酸（ABA）相关数量性状位点（QTL）进行比较，发现ZmDUF668-1定位于2号染色体上一个与叶片ABA浓度和根系性状相关的宽泛区间（L-ABA/root-ABA1）内，使其成为一个值得优先研究的候选基因。
3.6 ZmDUF668蛋白的结构特征、基因结构与预测功能关联
三维结构预测显示，代表性蛋白具有保守的α螺旋折叠结构，但ZmDUF668-19因发生大规模序列截短而呈现简化的构象。基因本体论（GO）富集分析表明B73驻留的DUF668基因显著富集于“生长正调控”。STRING网络分析揭示了一个包含31个蛋白和245条预测关联的功能网络，其中ZmDUF668-10具有较高的连接度。基因结构分析发现，B73驻留的19个基因中存在两种截然不同的架构：10个为多外显子基因，9个为无内含子基因。
3.7 启动子顺式作用元件分析
对B73和Il14H两个品系的启动子分析显示，两者预测的顺式作用元件总数相似，但在光响应、激素响应和胁迫响应等类别上的分布存在差异。
3.8 B73驻留ZmDUF668基因在不同条件下的表达模式
公共RNA-seq数据分析显示，B73驻留的19个基因表现出组织特异性表达模式，例如ZmDUF668-27/28/29在花粉中高表达，ZmDUF668-4/12在花丝中高表达。在非生物胁迫下，这些基因呈现出多样的时间动态响应，如ZmDUF668-5/6/12/16在涝渍早期（2-4小时）上调，而ZmDUF668-2/3/8则在后期（8-10小时）上调。在温度、干旱、盐和UV-B胁迫下，不同成员也展现出特异的表达变化。
3.9 通过qRT-PCR揭示的温度响应表达
针对12个B73驻留基因的qRT-PCR验证进一步确认了其温度响应性。ZmDUF668-4和-16在低温（4°C和10°C）下显著上调，被认为是冷响应候选基因；而ZmDUF668-8、-12和-14在高温（37°C和42°C）下显著上调，被认为是热响应候选基因。
3.10 涝渍期间的时间表达模式
对15个基因在涝渍时间进程中的qRT-PCR分析显示，ZmDUF668-4/7/15/28在2小时达到峰值，而ZmDUF668-5/6/10/11/16在4小时达到峰值，ZmDUF668-2/3/8/14则在8小时达到峰值，揭示了基因特异性的时序表达模式。
3.11 多种非生物胁迫下的qRT-PCR分析
qRT-PCR分析证实，ZmDUF668-4、-11、-12和-28在UV-B辐射下表达增加，其中ZmDUF668-4增幅最大。ZmDUF668-3和-5在冷胁迫下上调，ZmDUF668-1和-28在盐胁迫下上调。
3.12 水分亏缺胁迫下的表达模式
在水分子缺胁迫下，ZmDUF668-8和-17在严重缺水24小时后表现出最大的表达增量，被认为是晚期严重缺水响应候选基因。而ZmDUF668-1、-2、-4、-6、-11和-19则在轻度缺水24小时后显著上调。
讨论与结论
讨论部分指出，利用泛基因组策略克服了单参考基因组分析的局限性，揭示了玉米DUF668家族在群体水平上惊人的PAV多样性。虽然大部分基因受纯化选择约束，但少数基因（如ZmDUF668-16）和特定结构变异（如ZmDUF668-4和-15的插入）可能驱动了功能和表达的分化。结合B73驻留基因的组织特异性表达、对不同非生物胁迫（特别是水分子缺、盐害和UV-B）的动态响应，以及与ABA相关QTL区间的位置对应关系，这些结果为后续的功能验证提供了丰富的候选基因资源。研究结论明确指出，通过对26个高质量玉米基因组的分析，研究人员整合了51个非冗余的ZmDUF668泛基因组位点。综合的生物信息学分析显示核心成员保留了保守的架构，而结构变异则与基因结构和表达模式的差异相关联。公共RNA-seq数据和qRT-PCR分析共同证明了所选B73驻留基因在多种非生物胁迫下的响应表达。这些结果为深入研究玉米发育和非生物胁迫响应机制提供了重要的候选基因基础，并优先推荐ZmDUF668-1作为位于2号染色体L-ABA/root-ABA1区间内的位置候选基因进行后续功能验证。
