芦竹MADS-box基因家族全基因组鉴定及非生物胁迫表达谱与AdMADS33/AdMADS53激活AdNCED的候选调控研究

《Plants》：Genome-Wide Characterization and Expression Profiling of MADS-Box Family Genes Under Abiotic Stress in Arundo donax

　　

MADS-box转录因子（MADS-box TFs）在植物生长发育及非生物胁迫响应中发挥关键作用。然而，其在能源作物芦竹（Arundo donax，A. donax）中的抗非生物胁迫功能仍大部分未阐明。本研究在A. donax基因组中鉴定并系统分析了97个MADS-box基因（AdMADS）。系统发育分析（phylogenetic analysis）将这些基因分为I型（type I，44个成员）和II型（type II，53个成员）亚家族，并进一步细分为12个亚组。比较分析显示，I型和II型基因在编码蛋白理化性质及外显子-内含子结构方面存在明显差异。启动子分析鉴定出大量与光响应、植物激素和胁迫信号相关的顺式调控元件（cis-regulatory elements）。基因重复分析提示，全基因组重复（whole-genome duplication，WGD）/片段重复（segmental duplication）是已检测到的AdMADS基因家族扩张的主要贡献因素，其中II型基因主要参与这些重复事件。共线性分析（collinearity analysis）表明，A. donax、水稻（Oryza sativa，O. sativa）和谷子（Setaria italica，S. italica）之间的MADS-box基因高度保守，反映其在禾本科（Gramineae）内的进化保守性。通过分析公共RNA-seq数据，研究人员鉴定出45个AdMADS基因为水涝、干旱和盐胁迫下的胁迫响应候选基因，并通过RT-qPCR（reverse transcription-quantitative PCR）验证了候选基因。此外，研究人员发现AdMADS33和AdMADS53在烟草中可结合并激活脱落酸（abscisic acid，ABA）合成基因AdNCED的启动子，提示其可能通过促进AdNCED转录而参与耐盐性。本研究为未来功能基因挖掘提供了有价值资源，并为A. donax分子育种奠定了坚实基础。

研究背景：芦竹（Arundo donax L.）为多年生根茎禾草，具有快速生长、高生物量和在干旱、盐渍及重金属等边际环境适应性强等特点，是生物能源和植物修复候选作物；但其生物量组成和耐逆性等关键性状的分子机制仍大部分未阐明。MADS-box转录因子（MADS-box TFs，植物中调控发育与胁迫响应的一类转录因子）具有约60个氨基酸的保守MADS结构域，可结合CArG-box序列调控下游基因，在植物发育及非生物胁迫（abiotic stresses）响应中作用广泛；在拟南芥（Arabidopsis thaliana，A. thaliana）、水稻等中有耐盐、耐旱、耐冷等证据，但A. donax中MADS-box基因仍缺乏系统研究。因此，研究人员开展首个全基因组调查（genome-wide survey），以解析AdMADS基因家族特征、非生物胁迫响应及候选调控关系，为功能基因挖掘和分子育种提供基础。研究人员在A. donax基因组中鉴定97个AdMADS，分析其染色体定位、系统发育、基因结构、保守结构域与基序、顺式调控元件、基因重复及种间共线性；整合公共RNA-seq和RT-qPCR鉴定45个水涝、干旱和盐响应候选；进一步在异源系统中验证AdMADS33和AdMADS53结合并激活AdNCED启动子。结论提示AdMADS33/AdMADS53可能通过促进AdNCED转录参与耐盐，但尚需in planta遗传验证。其意义在于为AdMADS功能研究和耐逆A. donax分子育种提供资源与线索。

关键技术方法（≤250字）：研究人员采用BLAST v2.16.0与HMM v3.4检索A. donax基因组，并用Pfam确认MADS-box/SRF-TF与K-box结构域；以MUSCLE v5.1和IQ-TREE v2.4.0构建最大似然系统发育树；用CD-Search、MEME和PlantCARE分析保守结构域、基序及启动子顺式元件；用MCScanX、KaKs_Calculator 3.0和JCVI分析基因重复、Ka/Ks及种间共线性；公共RNA-seq数据来自NCBI SRA（水涝PRJNA1069878、干旱PRJEB10039、盐PRJNA497781），经HISAT2、featureCounts和DESeq2分析差异表达；以RT-qPCR验证候选基因，酵母单杂交（yeast one-hybrid，Y1H）和双荧光素酶（dual-luciferase，LUC）检测AdMADS33/AdMADS53与AdNCED启动子互作。植物材料源自中国湖北武汉单株（114.276°E，30.789°N）。

研究结果：2.1 鉴定、性质及染色体定位：通过BLAST+HMM从A. donax基因组鉴定97个AdMADS，分布于30条染色体，3号（11）、1号（9）、7号（8）最多，21、23、27、31、34、35号无分布；Pfam分为44个I型（仅MADS-box）和53个II型（MADS-box+K-box）；氨基酸65–633，分子量15.5–71.2 kDa，pI 4.68–9.85，不稳定指数32.9–77.66，均为亲水；I型较II型氨基酸多、分子量大、pI低、不稳定指数略低、GRAVY较不负。2.2 系统发育关系：用A. donax、O. sativa、A. thaliana的I/II型蛋白构建系统发育树；I型分Mα（25）、Mβ（8）、Mγ（11）；II型分ALG12（6）、ANR1-like（5）、B（7）、BS（2）、C/D（4）、E（4）、FUL-like（12）、SOC1-like（6）、SVP-like（7）；A. donax与O. sativa分布高度相似，Mα和FUL-like成员最多，符合近缘关系。2.3 基因结构及保守结构域与基序：基于系统发育框架分析motif、结构域和基因结构；motif 1–4为核心，分别出现于92.8%、89.7%、77.3%、76.3%；motif 5–7主要见于II型；I型无motif 5/6，仅11.6%含motif 7；所有蛋白含MADS-box/SRF-TF结构域，II型另有K-box；I型77.3%含0或1个内含子，II型均含5个以上外显子，II型内含子平均5898.9 bp、I型1989.3 bp。2.4 启动子顺式调控元件：用PlantCARE预测，鉴定50类顺式元件，平均26.1个/基因；AdMADS20最多（64），AdMADS79（58）、AdMADS85（52）次之，AdMADS87和AdMADS88无预测元件；光响应元件最丰富，G-box（247）见于82.5%启动子，Box 4（108）见于67.0%；激素相关MeJA元件（CGTCA/TGACG）346个、ABA响应ABRE 329个；胁迫相关厌氧诱导ARE/GC-motif 242个、干旱MBS 95个、低温LTR 57个，提示受外部环境和内部生理信号组合调控，但仍需实验验证。2.5 种内与种间共线性：用MCScanX等鉴定97个共线基因对，来自72个AdMADS，均预测源于WGD/片段重复；AdMADS28同源基因最多（6），AdMADS26和AdMADS84各5；重复对均限于同一型，提示I/II型分化早于A. donax的WGD；II型参与重复比例79.2%高于I型67.4%；Ka/Ks均<1.0，显示强纯化选择；重复时间28.86–0.35百万年，与第二次WGD一致；与O. sativa有81对（涉及68个AdMADS和30个OsMADS），与S. italica有79对（涉及74个AdMADS和38个SiMADS）；多数OsMADS（83.3%）和SiMADS（65.8%）在A. donax有2个以上直系同源，符合异源多倍体特征。2.6 非生物胁迫相关AdMADS鉴定：用公共RNA-seq分析水涝、干旱和盐胁迫表达，平均比对率85.93%；水涝30 d叶18、根19个差异表达，45 d叶15、根18个；叶中8个基因一致调控，根中8个基因一致调控；干旱10% PEG根无变化、仅叶AdMADS26下调，20% PEG根AdMADS37上调、叶4个基因上调；盐256.67 mM叶7个、419.23 mM叶16个差异表达，6个基因在两种盐浓度下共同响应；共45个AdMADS至少在一种胁迫下差异表达；RT-qPCR验证出胁迫特异性表达，RNA-seq与RT-qPCR的fold-change线性回归为y=0.416x+1.025，R²=0.667，呈显著正相关，支持RNA-seq可靠性。2.7 AdMADS33和AdMADS53在异源系统激活AdNCED启动子：200/400 mM NaCl处理A. donax幼苗，存活率分别为72.2%和25.0%；检测6个盐响应AdMADS和4个ABA合成基因，AdMADS33、AdMADS47、AdMADS53上调，AdNCED在盐胁迫下显著上调，其他ABA合成基因变化不明显；AdNCED启动子含2个CArG box；Y1H显示AdMADS33和AdMADS53可结合该启动子；双LUC在N. benthamiana中显示AdMADS33或AdMADS53增强AdNCED启动子驱动LUC活性，突变CArG后激活消失；结论为AdMADS33和AdMADS53可在异源系统结合并激活AdNCED启动子，但未建立in planta调控，需遗传证据。

讨论部分总结：讨论指出，A. donax的MADS-box基因在亚组分布上具有独特性；I/II型在蛋白理化性质和基因结构上差异显著，II型内含子更多，提示不同进化轨迹及I型可能的广泛内含子丢失。与水稻、谷子比较显示I型分布保守，II型中AGL6和MIKC*亚组缺失；这些亚组与花发育、花粉成熟相关，而A. donax具有不育和配子发生异常，但原文强调不育是多因素复杂现象，不能仅归因于MADS-box基因丢失。胁迫响应方面，水涝响应AdMADS最多（35个），干旱6个、盐17个；水涝响应基因多数在叶和根均差异表达，也有组织特异模式，提示转录网络复杂；45个胁迫响应基因中II型占91.1%（41/45），与其他物种趋势一致，但差异表达本身不能证明功能，需要遗传验证。盐响应候选方面，AdMADS33和AdMADS53可结合并激活AdNCED启动子，提示AdMADS-AdNCED-ABA模块可能参与盐响应；但现有证据基于异源报告系统，未建立in planta调控和耐盐表型，需过表达、敲除或沉默等直接遗传证据，并需测定ABA含量。

研究结论翻译：作为有前景的能源作物，A. donax需要深入的分子表征以支持其遗传改良。本研究首次对A. donax的MADS-box基因家族进行了全面的全基因组鉴定和表征；研究人员系统分析了其系统发育关系、染色体分布、基因结构、保守蛋白基序和结构域以及共线性关系。此外，通过整合RNA-seq和RT-qPCR分析，研究人员鉴定出45个响应水涝、干旱和盐胁迫的AdMADS基因。其中，AdMADS33和AdMADS53在异源系统中被证明可结合并激活AdNCED启动子，提示其可能通过促进AdNCED转录而参与耐盐性。然而，仍需直接的in planta遗传验证以确认其生理功能。研究结果为进一步功能研究奠定了基础，并为耐逆A. donax品种的分子育种提供了有价值的遗传资源。未来工作应聚焦于关键候选基因的功能表征、其调控网络的阐明以及通过转基因试验验证其作用。