《Biology》：Genome-Wide Identification and Characterization of the ABA Receptor PYL Gene Family in Broomcorn Millet (Panicum miliaceum L.)
编辑推荐：
脱落酸（abscisic acid，ABA）是协调植物生长、发育及逆境适应的重要植物激素。Pyrabactin抗性类似蛋白（Pyrabactin resistance-like，PYL）是ABA信号通路的核心受体，但在糜子（Panicum miliaceum
脱落酸（abscisic acid，ABA）是协调植物生长、发育及逆境适应的重要植物激素。Pyrabactin抗性类似蛋白（Pyrabactin resistance-like，PYL）是ABA信号通路的核心受体，但在糜子（Panicum miliaceum L.）中的家族组成及胁迫响应功能尚不清楚。研究人员在全基因组水平上鉴定了糜子PmPYL基因家族成员，并系统分析了其系统发育、基因结构、保守基序、染色体分布、共线性、选择压力、启动子顺式作用元件、预测调控网络、组织表达特征以及对ABA、盐和渗透胁迫的转录响应。共鉴定出15个PmPYL基因，其中8个位于A亚基因组，7个位于B亚基因组。系统发育分析将11个PmPYL蛋白归入Group II，4个归入Group III，Group I中无糜子成员。序列与结构分析显示，与ABA感知及蛋白磷酸酶2C（protein phosphatase 2C，PP2C）互作相关的受体特征高度保守。共线性分析鉴定出7对A/B同源基因对，所检测的复制基因对均呈现较低的Ka/Ks比值。启动子分析显示存在丰富的激素及胁迫响应元件，预测的蛋白互作网络支持保守的PYL–PP2C–蔗糖非发酵1相关蛋白激酶2（SNF1-related protein kinase 2，SnRK2）信号模块。同源基因对间的启动子及组织表达比较显示调控层面的分化。幼苗叶片qRT-PCR分析进一步显示不同处理和时间依赖的表达模式，其中PmPYL10、PmPYL8和PmPYL12在所检测条件下表现出显著的变化。这些发现揭示了糜子中保守的ABA受体框架及其调控多样性，并鉴定出可用于胁迫适应功能分析的候选PmPYL基因。
糜子ABA受体PYL基因家族研究解读
一、研究背景与科学问题
糜子（Panicum miliaceum L.）是欧亚干旱与半干旱地区栽培历史最悠久的谷类作物之一，作为一种C4、短生育期且耐旱的谷物，糜子适应于低降雨、土壤贫瘠和温度波动剧烈的边际土地。近年来，因其优异的胁迫耐受性、高水分利用效率以及富含膳食纤维、矿物质和生物活性物质的营养价值而重新受到关注。然而，与水稻和拟南芥等模式作物相比，糜子胁迫适应与激素信号的分子机制仍远未阐明。
脱落酸是协调植物生长、发育和逆境适应的核心植物激素，调控种子休眠与萌发、气孔关闭、根系发育、渗透调节以及对干旱和盐胁迫的响应。ABA核心信号通路由PYR/PYL/RCAR受体蛋白启动：ABA结合后，这些可溶性受体与A族蛋白磷酸酶2C互作并抑制其活性，从而解除PP2C对SnRK2激酶的抑制；活化的SnRK2进而磷酸化ABA响应转录因子、离子通道及其他下游靶标，重编程基因表达和生理响应。尽管该核心信号模块在陆生植物中高度保守，但不同PYL受体在表达模式、ABA敏感性、生化特性和生理功能上存在显著差异。拟南芥中的遗传学研究已证明PYR/PYL/RCAR受体间存在广泛的功能冗余，同时也揭示了受体特异的调控作用，如PYL8通过MYB77依赖的生长素响应转录促进侧根生长，PYL9则参与抗旱和ABA相关的叶片衰老。作物中的功能分化同样明显：水稻中过表达OsPYL5和OsPYL10可增强干旱或干旱与冷害复合耐受性，但ABA敏感性增强也可能带来生长代价；玉米中ZmPYL8、ZmPYL9和ZmPYL12有助于抗旱。谷子中SiPYL3的功能验证进一步表明，单个PYL受体在不同胁迫下可发挥相反的作用。
全基因组分析已在多种植物中鉴定出PYL基因家族，揭示了基因拷贝数、染色体组织、复制历史和胁迫响应表达方面的显著种间差异。基因复制在信号基因家族进化中尤为重要，因为复制后的受体可能保留保守的生化功能，而在空间或环境调控方面发生分化。这一过程对糜子可能尤为关键——糜子是异源四倍体物种，其基因组在两个亚基因组间保留了广泛的同源对应关系。然而，基因组复制是否促进了糜子PYL受体库的保留及随后的调控分化，此前尚不清楚。此外，PYL受体活性还受启动子结构、ABA结合后与PP2C的互作以及反馈调控机制等多层面的调节。尽管对模式作物和作物物种中PYL调控的理解日益深入，糜子PYL家族的组成、进化历史及调控特征仍 largely 未被探索，尤其是多倍化后PYL基因的保留方式、复制成员是否发生调控分化、以及哪些受体对ABA、盐和渗透胁迫响应最为强烈等问题均悬而未决。
二、研究内容与主要结论
针对上述问题，研究人员对糜子PYL基因家族进行了全基因组鉴定和系统表征，并利用时间序列qRT-PCR分析了PmPYL基因对50和100 μM ABA、200 mM NaCl及20% PEG-6000诱导渗透胁迫的转录响应。研究共鉴定出15个PmPYL基因，揭示了该家族在糜子中的组成、进化历史、调控特征和胁迫响应模式，为理解糜子ABA受体分化及筛选胁迫响应候选基因提供了框架。
三、关键技术方法
研究人员使用了以下主要技术方法：基于BLASTP和Hidden Markov Model（HMM，隐马尔可夫模型）双重策略的全基因组基因鉴定；利用MAFFT和IQ-TREE 2进行多序列比对和最大似然系统发育分析；采用MEME Suite进行保守基序鉴定；使用MCScanX进行共线性和基因复制分析；通过KaKs_Calculator 3.0计算非同义替换率（Ka）、同义替换率（Ks）及Ka/Ks比值评估选择压力；利用PlantCARE数据库进行启动子顺式作用元件预测；基于STRING数据库进行蛋白互作网络预测；采用qRT-PCR进行胁迫响应表达分析。植物材料为山西农业大学农业基因资源研究中心培育的糜子品种AJ 08，幼苗在人工气候箱中培养至三叶期后施加ABA、NaCl和PEG-6000处理。
四、研究结果
4.1 PmPYL基因家族的全基因组鉴定与特征分析
在AJ8参考基因组中，共鉴定出15个PmPYL基因，分布于A和B亚基因组的1、3、4、5和9号染色体上。A亚基因组含8个基因，B亚基因组含7个基因，整体分布均衡。大多数PmPYL基因以同源A/B染色体臂上的配对形式存在，表明基因组复制后PmPYL家族经历了较为协调的保留和分化。编码蛋白长度为196至223个氨基酸，预测分子量为21.26–24.02 kDa，理论等电点为5.08–8.40。
4.2 PmPYL蛋白的序列保守性与系统发育关系
基于糜子、拟南芥、水稻、谷子和高粱共59个PYL蛋白的最大似然系统发育树显示，11个PmPYL蛋白属于Group II，4个属于Group III，Group I中无糜子成员。7对A/B同源基因对均形成姐妹分支关系，与基于共线性的鉴定结果一致。多序列比对显示，15个PmPYL蛋白中13个含有SGLPA gate基序，PmPYL5和PmPYL7携带TGLPA变体；14个蛋白具有HRL latch基序。
4.3 PmPYL蛋白的二级结构与亚细胞定位预测
预测的二级结构中，α-螺旋占27.27%–34.29%，延伸链占25.23%–36.18%，无规卷曲占31.16%–43.20%。亚细胞定位预测显示，7个蛋白主要定位于细胞质，PmPYL4预测定位于细胞核和细胞质，PmPYL6、PmPYL8和PmPYL14预测定位于叶绿体，PmPYL9和PmPYL10预测为叶绿体和细胞质双重定位，PmPYL13和PmPYL15预测定位于细胞外空间和/或液泡。
4.4 PmPYL基因结构与保守基序组成
大多数PmPYL基因含有2至3个外显子，同一系统发育组的基因通常具有相似的外显子–内含子结构。MEME分析共鉴定出10个保守基序，大多数成员共享相似的基序组成和排列，其中Motif 1、3和5几乎存在于所有蛋白中。
4.5 PmPYL基因的种内与种间共线性分析
种内共线性分析鉴定出7对A/B同源基因对：PmPYL1–PmPYL4、PmPYL2–PmPYL3、PmPYL5–PmPYL7、PmPYL6–PmPYL8、PmPYL9–PmPYL10、PmPYL12–PmPYL14和PmPYL13–PmPYL15。PmPYL11未配对。种间比较显示糜子与水稻间的共线性关系多于糜子与拟南芥间的关系。
4.6 复制基因对的选择压力与分化时间
11个PmPYL共线性关系的Ka/Ks比值范围为0.0166–0.1563，均低于1.0，表明纯化选择在基因对进化保守中发挥了作用。7对A/B同源基因对的估计分化时间约为5.1–23.7百万年前（million years ago，MYA），4个额外非同源共线性关系的估计值约为68.9–104.9 MYA。
4.7 PmPYL基因启动子中的顺式调控元件
2 kb上游区域分析鉴定出激素响应、胁迫响应、防御相关、光响应和发育相关五大类顺式作用元件。激素响应元件包括ABA响应元件（ABRE）、生长素响应元件、赤霉素响应元件以及茉莉酸和水杨酸相关元件；胁迫相关元件包括MBS、MYB、MYC、LTR和ARE基序。
4.8 PmPYL相关蛋白互作网络预测
基于拟南芥参考网络的预测显示，PmPYL蛋白与A族PP2C（包括ABI1、ABI2、HAB1、PP2CA/HAB2、AHG1和HAI3）及SnRK2相关激酶存在潜在关联，支持糜子中保守的PYL–PP2C–SnRK2信号模块。
4.9 PmPYL基因的组织特异性表达模式
跨8个组织或发育阶段的表达分析显示，PmPYL2、PmPYL3、PmPYL5和PmPYL7在多个营养组织中广泛表达；PmPYL11在花序中优先表达；PmPYL10在根中表达相对较高；PmPYL1、PmPYL4、PmPYL8和PmPYL9在多数组织中转录丰度较低。
4.10 PmPYL基因的胁迫响应表达谱
在14个分析的PmPYL基因中，处理效应在13个基因中显著，时间效应和处理×时间互作在所有14个基因中均显著。PmPYL8在PEG-6000处理12小时时达到最大相对表达量61.45±10.50，对ABA和NaCl也表现出显著响应；PmPYL10在PEG-6000处理12小时达到10.06±1.25；PmPYL12在NaCl处理6小时达到8.41±2.85。PmPYL4的处理主效应不显著。
五、讨论与结论
讨论部分指出，糜子中15个PmPYL基因的鉴定及7对A/B同源基因对的保留表明，同源基因保留是当前PmPYL家族的重要特征。蛋白结构的强保守性与低Ka/Ks比值一致，反映了纯化选择对该受体家族的进化约束。然而，受体核心的保守并不意味着完全的功能冗余——同源基因对在启动子顺式元件组成和表达量方面存在变异，尽管这些观察在多重检验校正后未达到统计显著，不能确证调控亚功能化。分子定年结果需谨慎解读，因为估计值依赖于所假设的替换速率，且PmPYL11–PmPYL13对因Ks值过高而无法可靠定年。启动子预测元件数量不能解释qRT-PCR谱，PmPYL4含最多预测胁迫相关元件但其处理主效应不显著。组织表达差异和胁迫响应差异共同表明家族内存在功能分化潜力。由于实验设计缺乏时间匹配的未处理对照和0.1%乙醇溶剂对照，且培养托盘未随机分配，表达谱应视为探索性证据。预测的蛋白互作网络支持糜子中保守的PYL–PP2C–SnRK2信号模块，但属于计算推断而非实验验证。PYL受体活性还受磷酸化和泛素化等多层翻译后机制调控。
研究结论为：本研究在糜子中鉴定出15个PmPYL基因，包括8个A亚基因组和7个B亚基因组成员。序列、结构和共线性分析揭示了一个保守的受体框架，包含7对保留的A/B同源基因对，低Ka/Ks比值与纯化选择一致。家族成员在预测启动子组成和组织表达模式上存在差异。在幼苗叶片中，qRT-PCR分析揭示了ABA、NaCl和PEG-6000处理后相对于共享预处理基线的转录丰度变化，其中PmPYL8表现出最大的观察增幅。这些表达谱将PmPYL8、PmPYL10和PmPYL12鉴定为后续研究的候选基因，其在ABA信号和胁迫适应中的功能需实验验证。
