《Discover Plants》：Genome-wide and transcriptomic characterization of the GRAS transcription factor family in cowpea (Vigna unguiculata)
编辑推荐：
GRAS（GAI-RGA-SCR，植物特有转录因子）家族以参与植物发育与胁迫信号转导著称，但在豇豆（Vigna unguiculata）中尚未被系统表征。研究人员通过整合BLASTp与HMMER（基于PF03514.17隐马尔可夫模型）在豇豆参考蛋白组中鉴定并
GRAS（GAI-RGA-SCR，植物特有转录因子）家族以参与植物发育与胁迫信号转导著称，但在豇豆（Vigna unguiculata）中尚未被系统表征。研究人员通过整合BLASTp与HMMER（基于PF03514.17隐马尔可夫模型）在豇豆参考蛋白组中鉴定并解析该家族（VuGRAS），结合RNA-Seq与qPCR探究其在非生物（根区脱水，root-zone desiccation, RZD）与生物（创伤辅助CABMV或CPSMV接种）胁迫下的表达动态。共获得57个VuGRAS基因；基因结构显示多数无内含子，蛋白预测提示核定位、偏酸性及负GRAVY值（亲水特性）。该家族不均分布于11条染色体，扩张由分散、串联及近端重复驱动。系统发育将其归入16个亚族，各含保守基序。启动子区富集与胁迫响应TF相关的顺式作用元件（cis-regulatory elements, CREs），尤其DOF因子。GO与KEGG注释关联激素信号等过程。RZD下，VuGRAS在RZD25′普遍诱导、RZD150′普遍抑制，呈时间依赖性；VuGRAS6.1与VuGRAS8.1跨阶段持续诱导，为脱水响应优先候选。创伤病毒接种初期（≤16 h）对VuGRAS表达影响有限，暗示早期参与弱或具迟发性调控角色。该研究深化了对VuGRAS在胁迫适应中作用的理解，为豇豆遗传改良的功能研究奠基。
研究背景与立项逻辑
GRAS（GAI-RGA-SCR，植物特有转录因子家族）是植物发育与胁迫信号网络的关键调控节点，其C端含LHRI、VHIID、LHRII、PFYRE、SAW等保守区域。尽管在拟南芥、大豆、木豆、紫花苜蓿等物种中已有全基因组表征，但豇豆（Vigna unguiculata）作为耐瘠薄、固氮能力强且广泛分布于半干旱区的重要豆类作物，其GRAS家族长期处于空白状态。豇豆生产受根区干旱等非生物胁迫，以及CABMV（cowpea aphid-borne mosaic virus，豇豆蚜传花叶病毒）和CPSMV（cowpea severe mosaic virus，豇豆重型花叶病毒）等生物胁迫双重制约，而转录因子层面的调控机制不清。研究人员据此开展首项豇豆VuGRAS系统鉴定与胁迫表达研究，旨在绘制该家族进化、结构与功能响应图谱。
主要技术方法概览
研究人员以Phytozome v1.2豇豆蛋白组为对象，采用BLASTp与HMMER（PF03514.17）双策略鉴定，经Batch CD-Search校验GRAS结构域；利用MCScanX解析染色质定位与重复类型，IQ-TREE构建与菜豆、月季的最大似然系统树，MEME与CD-Search刻画基序/结构域；启动子1 kb区用MEME/TomTom比对JASPAR_CORE_2022_plant推断CREs；功能注释接PANNZER2与BlastKOALA。表达层面复用Cowpea Genomics Consortium（CpG）RNA-Seq：根区脱水用cv. Pingo de Ouro水培根样（RZD25′/RZD150′），病毒胁迫用抗性基因型IT85F-2687（CABMV）与BR-14 Mulato（CPSMV）叶片创伤辅助接种60 min/16 h样本，de novo组装后经RSEM定量、edgeR差异分析，并以qPCR（REST 2009相对表达）验证。
研究结果浓缩
3.1 鉴定与特征：双策略获58非冗余序列，剔除1个结构域不完整者，定57个VuGRAS（VuGRAS1–57）。多数无内含子；蛋白391–833 aa，MW 44.6–90.5 kDa，91.2%偏酸（pI<7），GRAVY全为负，82.4%预测核定位。表明家族具精简结构、亲水性与核内转录调控倾向。
3.2 染色体分布与扩张：57基因不均布11染色体，Vu05（11）与Vu11（10）最密。扩张机制以分散重复为主（77%），次为串联（16%）与近端（7%）；串联对Ka/Ks<1，受纯化选择，且多归LISCL亚族。
3.2.1 系统发育与基序：联合菜豆、月季共166蛋白建树，分16亚族（SCL3、DELLA、RAM1、Rc2等），PAT1（11）与LISCL（8）最庞大，Rc2首次在豆科外近缘出现。MEME获10基序，映射至LHRI–VHIID–LHRII–PFYRE–SAW；LISCL/PAT1保全套基序，NSP2缺SAW相关motif 4，SCR缺LHRII/SAW相关motif 4/8，显示亚族级结构分化。
3.3 启动子CREs：1 kb启动子分析（严格e-value<10?2）检出DOF、TCP、NAC、bHLH类结合位点，DOF（MA1274.1）出现125次最频，提示DOF可能为VuGRAS转录调控重要上游因子；3基因未见检出元件，归因为阈值与近端区域局限。
3.4 GO与KEGG：15个GO term锚定36基因，涵盖转录调控（32基因）与共生真菌响应（27基因）；VuGRAS17注释最宽（8 term）。KEGG仅map04075植物激素信号转导，专属DELLA亚族（VuGRAS4/33/39），指向激素通路专门化。
3.5 RNA-Seq胁迫表达：三类处理共84转录本。RZD25′ 22差异（偏诱导），RZD150′ 19差异（偏抑制）；WaCABMV60′仅4抑制，16 h及CPSMV各时点无显著变化。VuGRAS6.1与VuGRAS8.1跨RZD两时点持续诱导，为脱水核心候选。
3.6 qPCR验证：根脱水下VuGRAS6.1、8.1、14.1、22.1、28.1、54.1在RZD25′诱导；VuGRAS6.1/8.1等在RZD150′仍诱导，与RNA-Seq一致。病毒处理中已检靶点均无显著波动，印证早期病毒响应中GRAS参与有限。
讨论与结论翻译
讨论指出，57个VuGRAS规模与菜豆、苜蓿相近而少于大豆，反映豆科谱系特异扩增；无内含子结构利于胁迫下快速转录。核定位与亲水剖面符合GRAS共性。分散/串联重复与纯化选择维持家族保守性。Rc2跨科恢复提示古老分支留存。CRE中DOF富集将VuGRAS置于更广的胁迫转录网络。表达上，RZD早期诱导—晚期抑制的时间反转，说明VuGRAS在急性脱水中有阶段性重编程功能；病毒创伤初期不显效，不代表后期无用，可能介入系统抗性迟发层。
研究结论部分译写：本研究完成豇豆GRAS家族首次全基因组与转录组学刻画，明确57个VuGRAS成员的进化分类、结构特征与染色体动力学；揭示其在根区脱水下具明显时序表达重排，并筛选出VuGRAS6.1与VuGRAS8.1两类优先功能验证靶标；同时表明创伤辅助病毒侵染早期对该家族扰动较弱。成果为理解豇豆胁迫适应的转录调控基础提供框架，并支撑后续基因功能验证与抗逆育种设计。
