-
降香黄檀（Dalbergia odorifera）隶属于豆科（Fabaceae）黄檀属（Dalbergia），俗称海南黄花梨，为中国特有的珍稀常绿半落叶乔木，亦是国家二级重点保护野生植物[1]。该树种环境适应性强、抗风性能优异，具备良好的水土保持与防风固沙生态功能，生态应用价值突出。其心材呈深红褐色，质地坚硬致密、纹理独特、自带芳香，兼具极佳的耐腐与耐湿特性，是制作高档家具、工艺摆件与乐器的顶级名贵木材，经济价值极高。同时，降香黄檀也是传统名贵中药材，具有理气止痛、降压降脂等药用功效，在中医药领域应用广泛[2-3]。目前，降香黄檀全基因组测序与基础注释工作已完成，为其基因家族挖掘与功能解析提供了数据支撑。但相较于完善的基因组信息，调控降香黄檀生长发育、木材形成及环境适应性的核心分子机制仍不清晰，关键功能基因亟待挖掘，严重限制了其分子生物学机制研究与遗传改良进程。因此，系统开展降香黄檀关键转录因子家族的鉴定与功能解析，对揭示其生长发育调控机理具有重要的理论与应用价值。
MYB转录因子是植物中数量庞大、功能多样的核心转录因子家族之一，广泛参与植物生长发育、次生代谢调控、激素信号传导及生物与非生物胁迫应答等多项生命进程。MYB家族蛋白的典型特征为含有高度保守的MYB结构域，由51～52个保守氨基酸残基构成螺旋−转角−螺旋（helix-turn-helix）空间结构，可特异性结合下游靶基因启动子的顺式作用元件，从而精准调控靶基因的转录表达[4]。依据N端保守MYB结构域的数量与排布特征，植物MYB家族可划分为1R-MYB（MYB-related）、2R-MYB（R2R3-MYB）、3R-MYB（R1R2R3-MYB）与4R-MYB四个亚类[5]，各亚类成员呈现显著的功能分化特征。其中，1R-MYB主要参与植物生物钟调控、光周期响应、花青素合成及根毛发育等生理过程[6]；2R-MYB为植物MYB家族最主要的亚类，广泛调控次生代谢产物合成、细胞分化、激素信号通路及逆境应答过程，是植物生长发育调控的核心成员[7-8]；3R-MYB则主要参与细胞分裂调控与非生物胁迫适应性响应[9]。
近年来，MYB基因家族的生物学功能已在多种模式植物与木本植物中得到广泛验证，其功能多样性与调控保守性已形成普遍共识。在模式植物拟南芥（Arabidopsis thaliana）中，AtMYB75正向调控花青素生物合成，AtMYB21与AtMYB24特异性参与花药发育调控[10-11]。在木本植物与经济作物中，MYB转录因子是木材形成、器官发育、逆境耐受及次生代谢的关键调控因子。毛果杨（Populus trichocarpa）中PtrMYB3与PtrMYB20可通过激活纤维素与木质素合成通路，主导次生细胞壁发育与木材形成进程[12]；水稻（Oryza sativa）OsMYB110能够响应磷素营养信号，正向促进茎秆伸长与增粗生长[13]。此外，MYB家族成员可通过多条信号通路参与植物环境适应与代谢调控：月季（Rosa hybrida）RhMYB108依托乙烯与茉莉酸信号通路调控花瓣衰老进程[14]；苹果（Malus domestica）MdMYB88与MdMYB124可通过激活CBF信号通路增强植株低温胁迫耐受性[15]；水稻OsMYB1则负调控花青素的生物合成，参与植株次生代谢稳态调控[16]。综上可见，MYB转录因子在植物生长发育、代谢调控及逆境响应中发挥不可替代的作用，但目前关于降香黄檀MYB基因家族的系统鉴定、进化特征及生物学功能的研究仍较为匮乏。
基于此，本研究依托降香黄檀全基因组数据，系统完成DoMYB基因家族的成员鉴定、系统进化分析、基因结构与保守基序分析、染色体定位、共线性分析及启动子顺式作用元件预测，并结合不同组织的表达模式解析其功能分化特征。同时，基于前期研究基础，筛选核心2R-MYB亚类基因DoMYB48，进一步探究其组织表达特异性与亚细胞定位特征。本研究旨在系统阐明降香黄檀MYB基因家族的结构特征与表达规律，为深入解析其生长发育的分子调控机制、挖掘优良功能基因提供理论基础与候选基因资源。
-
3年生降香黄檀植株样本采自海南省海南大学儋州校区试验基地（19°30'40" N，109°29'25" E），取样时间为2021年11月；7年生降香黄檀植株样本采自海南省定安县人工培育苗圃（19°38'56" N，110°14'29" E），取样时间为2020年7月。分别采集3年生与7年生降香黄檀植株的韧皮部、边材、茎干过渡区及叶片组织，用于转录组测序（RNA-seq）分析，其中，茎干过渡区定义为植株茎干心材与边材之间的过渡生长区域。本氏烟草（Nicotiana benthamiana）作为亚细胞定位试验材料，常规培养条件设置为：温度24℃、光照周期16 h光照/8 h黑暗。
-
本研究所用降香黄檀基因组数据为本课题组前期组装注释的未公开高质量基因组，基因组总大小约 1.2 Gb，Contig N50 为 58.21 Mb，共注释得到56 679 个蛋白编码基因，基因组组装与注释完整性良好。
从 Pfam 数据库（
https://pfam.xfam.org/ ）获取 MYB 保守结构域对应的隐马尔可夫模型 HMM 文件 PF00249 [17]，借助 TBtools 软件中 Simple HMM Search 工具对降香黄檀全基因组蛋白序列进行比对检索，筛选标准设定 E-value＜1×10−5。同时，从拟南芥专用数据库 TAIR（http://www.arabidopsis.org/ ）下载拟南芥全部 MYB 家族蛋白序列作为查询集 [18]，采用 BLASTP 程序比对降香黄檀蛋白库，阈值同样设为 E-value＜1×10−5。合并 HMM 检索与 BLASTP 比对两套筛选结果，剔除重复冗余序列；再分别利用 Pfam 在线工具与 NCBI 保守结构域批量检索平台 Batch CD-search（
https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi?SEQUENCE ）校验候选蛋白的 MYB 结构域完整度，删除结构域残缺的序列，最终获得可靠的 DoMYB 家族成员 [19]。采用 Expasy ProtParam 工具（
https://web.expasy.org/protparam/ ）批量预测 DoMYB 蛋白的分子量、理论等电点、亲疏水性、不稳定系数等基础理化特征 [20]；通过 WoLF PSORT 在线平台（https://wolfpsort.hgc.jp/ ）完成所有 MYB 蛋白的亚细胞定位生物信息学预测，物种分类选择植物（Plant）模式开展分析 [21]。 -
根据基因组注释文件，利用TBtools软件将降香黄檀DoMYB基因家族成员在染色体上的物理位置进行可视化展示，并对基因结构进行分析。使用NCBI Batch CD-search进行保守结构域的预测；运用MEME Suite（
https://meme-suite.org/meme/ ）在线分析工具对保守基序进行预测，基序识别的最大数量被设定为20。使用TBtools进行保守结构域、保守基序及基因结构的可视化分析 [22]。 -
从TAIR数据库、Uniprot数据库（
https://www.uniprot.org/ ）和Phytozome数据库（https://phytozome-next.jgi.doe.gov/ ）下载拟南芥、大豆（Glycine max）、水稻、毛果杨的家族的数据。利用ClustalW算法对拟南芥与降香黄檀的序列进行比对，使用MEGA12（v12.0）软件基于邻接法（Neighbor-Joining，NJ）构建系统发育树，bootstrap重复次数为1 000次，其余参数为默认值。降香黄檀种内系统进化树划分为23个亚组（A1～A23），与拟南芥MYB基因构建的系统发育树划分为25个亚组（C1～C25）。使用iTOL（https://itol.embl.de/ ）进行美化和可视化。在生成的系统发育树中，拟南芥的基因用方框标记，降香黄檀的基因用星形标记。此外，不同的MYB亚家族成员用颜色清晰区分：1R-MYB蛋白以棕色标示，2R-MYB蛋白以黑色标示，3R-MYB蛋白以红色标示，拟南芥中非典型的4R-MYB蛋白则以蓝色标示。使用TBtools中的MCScanX插件对降香黄檀种内的MYB基因进行共线性分析，同时计算非同义替换率（Ka）、同义替换率（Ks）以及二者的比值（Ka/Ks），其中Ka/Ks>1表明正选择，Ka/Ks=1表明中性进化，Ka/Ks<1表明纯化选择；使用TBtools中的MCScanX插件对降香黄檀（Do）、拟南芥（At）、水稻（Os）、大豆（Gm）、毛果杨（Ptr）的基因组内的MYB基因进行共线性分析，绘制物种间共线性图谱并可视化展示共线性基因对。
-
本研究采用TBtools程序，对所有转录起始位点上游2 000 bp范围内的启动子片段进行批量提取，再将得到的序列导入PlantCARE（https://bioinformatics.psb.ugent.be/webtools/plantcare/html）在线网站进行顺式调控元件的识别与注释，对预测结果进行整理筛选与精简归类，使用Excel进行绘制。
-
利用降香黄檀3年生树木的韧皮部（B01～B05）、边材（S01～S05）、过渡区（T01～T05）及叶片（L01～L05）四个组织的转录组表达数据（PRJCA056155）的TPM值，其中每个组织的样本取自5株已形成心材的独立植株，使用Z-score标准化数据进行热图绘制。采用RNA提取试剂盒（CW0581S），从降香黄檀7年生植株的韧皮部、边材、过渡区和叶片组织中提取RNA。利用1 μg RNA，按照PrimeScript TM RT reagent Kit with gDNA Eraser试剂盒配置反转录体系，反转录为cDNA；以TB Green Premix Ex Taq Ⅱ试剂盒配置qPCR体系，设计DoMYB48基因特异性引物，以Actin为内参基因进行RT-qPCR反应 [23]。实验设置3次生物学重复，每次重复包含3次技术重复。其中RT-qPCR反应体系如下：引物为0.4 μL、TB Green Premix Ex TaqII为10 μL、cDNA为2 μL、无菌水为7.2 μL；反应程序使用两步法PCR扩增程序：95℃预变性30 s，95℃条件下5 s，60℃条件下34 s，循环数为40。使用2-∆∆Ct法计算相对表达量，使用Graphpad Prism制作柱状图，采用单因素方差（One-way ANOVA）结合Bonferroni校正，确定各组织间表达差异的显著性，使用星号标注显著水平，其中****表示P≤0.000 1。
表 1 RT-qPCR引物序列
Table 1. Primer sequences used for RT-qPCR
基因名称
Gene基因ID
Gene ID引物F（5′−3′）
Forward primer引物R（5′−3′）
Reverse primerActin Daodo09aG0171500 CAATGAATTGCGTGTTGCT ATACCAGTTGTGCGACCACTT DoMYB48 Daodo02aG0332400 CCCTTCACCACCACCAGTTTC TTACCATCACCGTTATCATTAGCC -
全基因组关联分析所用的303份降香黄檀植株群体栽植于海南大学儋州校区，连续测定了12个月份的株高数据。基于全基因组重测序数据，使用GEMMA软件结合亲缘关系与群体结构，采用混合线性模型（Mixed Linear Model，MLM）对303份样本的株高表型数据进行关联分析，显著性阈值采用Bonferroni校正法确定，设置P value<1×10−5的SNP（Single Nucleotide Polymorphism，SNP）为显著的SNP。
-
选取降香黄檀叶片cDNA为模板，利用同源重组方法将扩增片段构建至pCAMBIA1300载体中。将构建好的pCAMBIA1300-DoMYB48质粒转入GV3101菌株中，在LB液体培养基中进行振荡培养，直至OD600达到0.5～0.6，离心收集菌体，使用由10 mmol·L−1 MgCl2、10 mmol·L−1 MES、100 μmol·L−1 乙酰丁香酮（Acetosyringone，AS）组成的侵染液重悬菌体，室温暗培养2～3 h。将重悬液注射到烟草叶片中，黑暗培养12 h后光照培养24 h，在激光共聚焦显微镜（LSM900）下对荧光信号进行亚细胞定位情况判断。
-
本研究使用Pfam数据库中MYB基因的HMM模型在降香黄檀基因组中进行筛选，经Pfam及NCBI结构域验证并去除冗余后，共筛选鉴定到180个MYB家族成员（DoMYB1～DoMYB180）。根据R结构域数量分为1R-MYB（38个）、2R-MYB（137个）、3R-MYB（5个）3种类型，未检测到4R-MYB类型成员。染色体定位显示180个DoMYB基因不均匀地分布于全部10条染色体上，不同染色体上的基因数量介于14至29之间，其中2号、3号、8号及9号染色体上存在明显的基因簇（图1）。
图 1 DoMYB基因家族在降香黄檀染色体上的基因分布
Figure 1. Gene distribution of the DoMYB gene family on D. odorifera chromosomes
理化性质分析表明，DoMYB蛋白相对分子质量范围为8.17～183.08 kDa；氨基酸数目差异显著，最少的有72个氨基酸，最多的有1 683个；等电点范围为4.27～9.92；不稳定系数介于30.86～88.91之间，仅有6个蛋白的不稳定系数小于40，多数DoMYB蛋白为不稳定蛋白；脂肪系数范围为50.00～91.25；除DoMYB130外，其他蛋白的平均亲水性小于0，表明该家族成员多为不稳定亲水性蛋白。3R-MYB亚类的平均氨基酸数目和分子质量均高于1R-MYB和2R-MYB，其余理化性质在各亚类间无显著差异。亚细胞定位预测显示，全部成员表现出细胞核特异性定位的特点。
-
以180条DoMYB蛋白序列为基础构建系统发育树，根据系统发育关系将DoMYB基因划分为23个亚组（A1～A23）。其中，1R-MYB的38个成员集中分布于A10、A11、A13、A14和A21亚组中，其余分散在A1、A7、A9、A18、A22和A23亚组中；3R-MYB主要分布在A22亚组中，2R-MYB类型的137个成员分布最为广泛，遍布多个亚组。
保守结构域分析显示，所有的DoMYB基因都包含典型的MYB结构域（PLN03091、PLN03212、SANT），其中PLN03091占比最高（97个，53.9%），其次为PLN03212结构域（39个，21.7%）与SANT结构域（18个，10.0%）。其中DoMYB9、DoMYB105及DoMYB122不仅含有典型的MYB结构域，还包含COG5114、RSC8等其他类型的结构域（图2-A）。MEME分析鉴定出20个不同的保守基序（Motif1～Motif20），每个蛋白包含1到10个不等的基序。Motif3（91.7%）、Motif1（85.0%）和Motif5（77.8%）出现频率最高，其余的基序出现频率为3.9%～26.1%。同一亚组内的基因通常具有相似的基序组成，如A5亚组均含Motif1、Motif2、Motif3、Motif5和Motif12，A11亚组中均含Motif8、Motif10、Motif18。此外，亚组中存在不同分支具有不同Motif数目的情况，A22亚组10个基因均含Motif1、Motif2、Motif3、Motif5、Motif12，其中分支1的5个基因（DoMYB8、DoMYB29、DoMYB121、DoMYB141、DoMYB158）相较于同亚组分支2（DoMYB9、DoMYB42、DoMYB48、DoMYB122、DoMYB151）多出Motif8（图2-B）。基因结构解析结果表明，外显子数量范围为1～12个，内含子数量范围为0～11个，其中120个基因含3个外显子，29个基因含有2个外显子，仅有9个基因含有1个外显子，其余的22个基因的外显子数量范围为4～12个。同亚组中亲缘关系紧密的基因，其基因结构模式较为一致，而不同亚组间UTR、CDS和内含子在数量和位置上存在差异（图2-C）。
图 2 DoMYB基因家族保守结构域、保守基序、基因结构分析
Figure 2. Analysis of conserved domains, conserved motifs and gene structure of the DoMYB gene family
-
基于180个降香黄檀MYB蛋白与168个拟南芥的MYB蛋白序列共同构建系统发育树，被划分为25个亚组（C1～C25）（图3）。1R-MYB成员主要分布在C1～C5亚组，3R-MYB蛋白主要集中于C23亚组，2R-MYB蛋白分布于整个进化树范围内。各亚组成员数量差异显著，C5亚组最多（60个），C13亚组（37个）和C12亚组（35个）次之。根据拟南芥中MYB转录因子的已知功能，可推测降香黄檀同源基因功能。C9亚组中的DoMYB112、DoMYB113、DoMYB148、DoMYB97、DoMYB90、DoMYB27等与拟南芥基因AtMYB123聚类于同一分支中，这些基因可能在根系发育与形态建成发挥作用[24]；C20亚组中的DoMYB48与AtMYB104聚类于同一分支，属于GAMYB亚类成员，可能在植物的种子萌发、茎的伸长及花药分化过程中发挥重要的作用[25]。
图 3 降香黄檀和拟南芥MYB基因家族系统进化树分析
Figure 3. Phylogenetic analysis of the MYB gene families in D. odorifera and A. thaliana
-
对DoMYB基因进行种内共线性分析，共鉴定出56对共线性关系的基因对。主要分布在2号、3号及5号染色体中，其中2R/2R类型基因对46对，1R/1R类型基因对4对，1R/2R类型基因对4对，3R/3R类型基因对2对。计算了56对重复基因对的Ka、Ks以及Ka/Ks比值以评估选择压力，结果显示所有可分析基因对的Ka/Ks比值均未超过1（0.106～0.465），表明这些基因对受纯化选择作用（图4、表2）。
图 4 DoMYB基因共线性分析
Figure 4. Intra-species collinearity analysis of DoMYB genes
表 2 DoMYB基因家族种内共线性Ka/Ks计算
Table 2. Intraspecific co-linearity Ka/Ks calculation of the DoMYB gene family
基因1（Gene 1） 基因2（Gene 2） Ka Ks Ka/Ks DoMYB1 DoMYB39 0.149 0.763 0.195 DoMYB2 DoMYB38 0.119 0.689 0.172 DoMYB3 DoMYB37 0.465 1.063 0.438 DoMYB9 DoMYB48 0.560 1.351 0.415 DoMYB2 DoMYB75 0.342 3.228 0.106 DoMYB17 DoMYB110 0.303 0.887 0.342 DoMYB15 DoMYB115 0.186 0.527 0.353 DoMYB14 DoMYB118 0.109 0.367 0.297 DoMYB21 DoMYB127 0.186 1.052 0.177 DoMYB6 DoMYB124 0.135 0.463 0.291 DoMYB8 DoMYB121 0.148 0.428 0.344 DoMYB32 DoMYB62 0.176 0.510 0.345 DoMYB50 DoMYB80 0.171 0.550 0.311 DoMYB51 DoMYB82 0.215 0.536 0.400 DoMYB30 DoMYB79 0.078 0.563 0.139 DoMYB24 DoMYB88 0.129 0.792 0.163 DoMYB26 DoMYB89 0.079 0.575 0.137 DoMYB27 DoMYB90 0.160 0.616 0.260 DoMYB28 DoMYB91 0.069 0.464 0.149 DoMYB35 DoMYB83 0.318 0.683 0.465 DoMYB44 DoMYB104 0.129 0.504 0.256 DoMYB45 DoMYB103 0.133 0.578 0.230 DoMYB46 DoMYB102 0.115 0.500 0.229 DoMYB40 DoMYB153 0.496 2.227 0.223 DoMYB53 DoMYB77 0.189 0.877 0.216 DoMYB55 DoMYB98 0.173 0.821 0.211 DoMYB56 DoMYB99 0.209 0.683 0.305 DoMYB66 DoMYB165 0.091 0.563 0.161 DoMYB67 DoMYB166 0.166 0.823 0.202 DoMYB57 DoMYB173 0.163 0.758 0.215 DoMYB59 DoMYB174 0.321 0.838 0.383 DoMYB60 DoMYB175 0.118 0.445 0.265 DoMYB61 DoMYB176 0.284 0.821 0.346 DoMYB63 DoMYB177 0.171 0.596 0.287 DoMYB64 DoMYB178 0.163 0.553 0.295 DoMYB65 DoMYB179 0.155 0.711 0.218 DoMYB53 DoMYB172 0.380 3.518 0.108 DoMYB53 DoMYB173 0.300 2.451 0.122 DoMYB70 DoMYB137 0.12 0.486 0.253 DoMYB71 DoMYB136 0.199 0.459 0.433 DoMYB78 DoMYB149 0.170 0.772 0.220 DoMYB84 DoMYB139 0.111 0.707 0.157 DoMYB92 DoMYB150 0.149 0.518 0.288 DoMYB93 DoMYB152 0.164 0.619 0.266 DoMYB94 DoMYB154 0.153 0.627 0.245 DoMYB108 DoMYB119 0.101 0.470 0.214 DoMYB107 DoMYB120 0.094 0.589 0.160 DoMYB123 DoMYB134 0.328 1.594 0.206 DoMYB143 DoMYB160 0.190 0.483 0.393 DoMYB146 DoMYB162 0.187 1.181 0.158 DoMYB140 DoMYB159 0.089 0.461 0.194 DoMYB141 DoMYB158 0.157 0.483 0.326 进一步将降香黄檀与拟南芥、大豆、水稻和毛果杨4个物种进行种间共线性分析，降香黄檀与4个物种均存在共线性基因对，与拟南芥之间鉴定出201对同源基因对，与大豆之间鉴定出699对同源基因对，与水稻之间鉴定出92对同源基因对，与毛果杨之间鉴定到386对同源基因对。进一步分析表明，DoMYB133与其他物种的20个基因存在同源共线性，该基因与拟南芥中5个基因同源，与大豆中8个基因同源，与毛果杨中6个基因同源，与水稻中1个基因同源，是具有最多共线性对的基因。此外，DoMYB19、DoMYB53、DoMYB98等47个基因在拟南芥、大豆、水稻和毛果杨之间都存在共线性关系，表明这些基因在MYB基因家族的长期进化过程中发挥了关键的调控作用（图5）。
图 5 降香黄檀与拟南芥、大豆、水稻、毛果杨MYB家族基因的种间共线性关系
Figure 5. Interspecific syntenic relationships of MYB family genes among D. odorifera, A. thaliana, soybean, rice, and P. trichocarpa
-
为探究降香黄檀DoMYB基因转录调控潜力，对其上游启动子区域内的顺式作用元件进行预测，共筛选到41种类型各异的顺式调控功能元件，对应的元件总数为4 541个（图6）。将这些功能序列根据其生物学功能归属，划分为4大功能类别：光响应元件（13种）、生长发育相关元件（12种）、植物激素响应元件（10种）以及生物及非生物胁迫响应元件（6种）。光响应元件的占比与数量高于其他类别，179个基因含有光响应元件，其中Box 4和G-box出现频率较高，表明DoMYB基因家族普遍受到光照信号的调控。胁迫响应元件以厌氧诱导元件（ARE）分布最广（136个基因），耐旱元件（MBS）和低温响应元件（LTR）分别分布于约1/3和1/4的成员中。生长发育相关元件如分生组织表达元件（CAT-box）分布于63个基因中，昼夜节律控制元件（Circadian）分布于42个基因中。此外，植物激素响应元件中，脱落酸响应元件（ABRE）分布于136个基因中，茉莉酸甲酯响应元件（TGACG-motif、CGTCA-motif）分布于125个基因中，响应赤霉素元件（TATC-box、GARE-motif、P-box）分布于90个基因中。上述结果表明，DoMYB基因家族的启动子区域富含各类顺式作用元件，表明其成员在多个生物学过程中可能发挥重要的转录调控功能。
图 6 降香黄檀MYB基因家族成员启动子顺式作用元件分析
Figure 6. Analysis of the cis-acting elements of the promoter regions of MYB gene family members in D. odorifera
-
利用降香黄檀3年生植株的韧皮部（B01～B05）、边材（S01～S05）、过渡区（T01～T05）及叶片（L01～L05）组织的转录组测序数据，对DoMYB在各组织样本中的表达模式进行统计解析，综合其相对组织特异性及表达量数据，不同成员呈现显著的组织差异表达（图7）。在叶片中的高表达基因数量最多，占总基因数的52.7%，其中DoMYB56和DoMYB155在叶片中的表达水平显著高于其他组织；韧皮部高表达基因占总基因数的15.6%，包括DoMYB83、DoMYB140、DoMYB159等，其中DoMYB159在韧皮部中表现出明显的组织特异性；边材高表达基因占总基因数的15.6%，包括DoMYB30、DoMYB65、DoMYB79、DoMYB169等；过渡区高表达基因占总基因数的7.8%，包括DoMYB6、DoMYB84、DoMYB89、DoMYB138、DoMYB179等。DoMYB基因在不同组织中发挥重要且独特的作用，可能参与组织特异性的基因调控网络。降香黄檀种内发育树中的A22亚组中包含2R-MYB成员和3R-MYB成员，其两个分支在保守基序和表达模式上存在显著差异。分支1中DoMYB8、DoMYB29、DoMYB121在边材和过渡区中高表达，DoMYB141和fgvbDoMYB158在叶片中高表达；而分支2基因整体主要在叶片中高表达，在其他组织中表达受到抑制，其中DoMYB48在L01与L02样本中高表达。
图 7 DoMYB基因家族表达模式分析
Figure 7. Analysis of the expression pattern of the DoMYB gene family
前期课题组通过全基因组关联分析，筛选到与株高显著关联的DoMYB48基因（P=4.151×10−6）。为进一步研究DoMYB48在不同组织样本中的表达模式，采用RT-qPCR检测该基因在7年生降香黄檀韧皮部、边材、过渡区以及叶片中的表达量，结果显示，DoMYB48在7年生降香黄檀过渡区组织中表达量最高，与叶片组织中的表达量相比呈现出583倍的显著差异（图8）。
图 8 降香黄檀DoMYB48基因在7年生降香黄檀中组织差异表达
Figure 8. Tissue-specific expression of the DoMYB48 gene in 7-year-old D. odorifera
-
为了进一步探究DoMYB48基因功能，本研究通过构建重组pCAMBIA1300-eGFP-DoMYB48载体，转化农杆菌后瞬时转化至烟草叶片中，在激光共聚焦显微镜下观察DoMYB48融合蛋白的荧光信号，并以pCAMBIA1300-eGFP空载作为对照。注射空载菌液的烟草叶片中，绿色荧光在细胞核和细胞膜上分布；注射pCAMBIA1300-eGFP-DoMYB48载体菌液的烟草叶片中，绿色荧光分布在细胞核中，与细胞核标记蛋白的红色荧光重叠，表明DoMYB48蛋白定位于细胞核，符合转录因子的调控功能特征（图9）。
图 9 DoMYB48亚细胞定位分析
Figure 9. Subcellular localization analysis of DoMYB48
-
MYB转录因子广泛参与各类生长发育进程、次生代谢、逆境响应等关键生物学过程[26]。本研究在降香黄檀中共鉴定到180个DoMYB成员，分为1R-MYB（38个）、2R-MYB（137个）和3R-MYB（5个）三种类型。与拟南芥（198个）[27]、大豆（244个）[28]、毛果杨（192个）[29]等物种数量相当，反映了MYB家族在植物进化中的保守性；与小麦（Triticum aestivum）（719个）[30]、棉花（Gossypium hirsutum）（812个）[31]等相比，MYB基因家族的数量差异反映在进化过程中全基因组多倍体化事件驱动的MYB家族扩张变化规律。DoMYB基因在10条染色体上呈不均匀分布，部分基因呈簇状分布，这种分布与串联重复事件相关，是基因家族扩张的重要机制[32]。系统发育分析将180个DoMYB成员划分为23个亚组，与拟南芥MYB家族的分类系统具有良好的对应关系。基因组内共鉴定出56对共线性基因对，Ka/Ks比值均小于1，重复基因经历了强烈的纯化选择后得以保留[33]。种间共线性分析显示，降香黄檀与大豆的共线性对远多于与水稻的共线性对，这与豆科内部的近缘关系一致。47个DoMYB基因在拟南芥、大豆、水稻和毛果杨中均存在同源基因，可能具有潜在的跨物种功能保守性，在MYB家族进化中具有重要作用。
理化性质分析表明，家族成员多为不稳定亲水性蛋白，其氨基酸长度差异显著，该家族可能参与从基础代谢到发育等广泛与复杂的调控[34]。所有成员均预测定位于细胞核，这与MYB转录因子在核内调控基因表达的功能一致。对保守基序进行分析显示，20个保守基序在不同基因间呈现不均匀分布，具有不同保守基序的基因可能参与特定的生理途径或胁迫响应。A22亚组内出现基序分化：分支1含有Motif8，而分支2缺失该基序，基序的丢失可能导致DNA结合特异性或蛋白互作能力的改变，进而驱动功能分化[35]。基因结构分析显示，多数基因的结构为3个外显子和2个内含子的模式，这与油菜（Brassica napus）中的研究一致[36]。启动子区域顺式作用调控元件分析显示，光响应元件最为丰富，表明MYB基因在光信号响应及光形态建成中扮演重要角色[37]。
组织表达分析显示，大多数DoMYB基因在叶片中的表达量较高，这与启动子区域富含光响应元件结果相符，也与叶片作为光合作用主要场所的生理功能相一致。在叶片中高表达的DoMYB155位于C9分支，与AtMYB123聚类，可能在花青素合成与种子发育中发挥关键作用[24]。在韧皮部中高表达的DoMYB159位于C13分支，与HAG亚家族的AtMYB28和AtMYB29聚类，其在拟南芥中参与硫代葡萄糖苷的生物合成调控，表明DoMYB159可能在韧皮部防御反应或次生代谢物合成中发挥重要作用[38]。在边材中高表达的DoMYB30、DoMYB65、DoMYB79和DoMYB169位于C15分支，与AtMYB26、AtMYB46、AtMYB83、AtMYB103等次生细胞壁及木质素生物合成重要基因聚类，可能参与边材的次生细胞壁形成和木质化过程[39-40]。在过渡区中高表达的DoMYB179与AtMYB46、AtMYB83、AtMYB103等次生细胞壁相关基因聚类于C15分支，可能参与边材向心材过渡区中细胞壁组分的调控与代谢物积累过程[41]。不同DoMYB基因分别调控相应组织中的特异基因表达网络，为理解MYB转录因子在植物生长发育中的复杂调控机制提供了重要线索。
基于课题组前期研究，筛选到与株高显著相关的2R-MYB成员DoMYB48，该基因与拟南芥C20亚组中AtMYB104聚类，属于与赤霉素调控相关的GAMYB分支，参与花药发育及GA信号传导，调节茎伸长及植株开花[42-44]。组织表达分析显示，DoMYB48在3年生降香黄檀植株中叶片表达量最高，可能主要参与营养生长阶段的叶片发育及茎伸长过程；在7年生植株中过渡区表达量最高，表明其功能可能随发育阶段转换，从早期调控叶片发育转向参与心材形成的次生代谢调控，这种转变可能与降香黄檀心材形成的树龄相关，7年生植株进入心材发育阶段，过渡区成为代谢活跃区[45]。在其他木本植物中，MYB转录因子随着发育阶段转换的现象已有报道；在杉木中，基因在幼龄茎部表达较高，但随着茎的成熟表达量逐渐下降[46]。DoMYB48基因的表达模式随着植株的发育过程产生了变化。
综上所述，在降香黄檀基因组中共鉴定到180个DoMYB基因家族成员，在保守基序、基因结构及理化性质方面表现出明显的特异性与多样性。启动子区域内光响应元件最为丰富。组织表达分析显示DoMYB基因呈现组织特异表达模式，其中DoMYB48在降香黄檀的不同发育阶段具有不同的差异表达模式。本研究为解析MYB转录因子在降香黄檀生长发育中调控机制及后续功能验证提供了理论基础。
降香黄檀DoMYB基因家族的鉴定及表达模式分析
DOI: 10.15886/j.cnki.rdswxb.20260089
CSTR: 32425.14.j.cnki.rdswxb.20260089
Genome-wide identification and expression pattern analysis of DoMYB gene family in Dalbergia odorifera
-
摘要: 降香黄檀（Dalbergia odorifera）兼具极高的社会经济价值、生态功能与药用研究潜力。MYB转录因子是植物重要的调控因子，可通过多样化调控模式广泛参与植物生长发育进程。为揭示降香黄檀生长发育的分子调控网络，本研究基于其全基因组数据对DoMYB基因家族开展系统的鉴定、结构特征及表达模式分析。结果显示，共鉴定获得180个降香黄檀MYB家族基因成员，依据R结构域数量可划分为1R-MYB（38个）、2R-MYB（137个）和3R-MYB（5个）三类；该家族基因不均匀分布于降香黄檀10条染色体上，且部分基因呈簇状分布特征。系统进化分析将180个DoMYB基因归为23个亚组，同一亚组成员的基因结构与保守基序高度保守，进化特征相似。启动子顺式作用元件预测结果表明，DoMYB家族基因广泛参与光响应、激素信号传导以及生物与非生物胁迫应答过程。组织表达模式分析发现，DoMYB各成员在降香黄檀韧皮部、边材、茎干过渡区及叶片中存在显著的组织特异性表达特征。其中，关键基因DoMYB48在7年生降香黄檀植株的茎干过渡区表达量最高，且该基因定位于细胞核。综上，降香黄檀DoMYB基因家族在基因结构与表达模式上具有丰富的多样性，研究结果为阐明降香黄檀生长发育的分子调控机制提供了关键的候选基因资源与理论依据。Abstract: Dalbergia odorifera possesses significant socioeconomic value, ecological functions, and potential for medicinal research, making it an important subject of study. MYB transcription factors are important regulatory elements in plants and are extensively involved in plant growth and development through diverse regulatory mechanisms. To elucidate the molecular regulatory network governing the growth and development of D. odorifera, a systematic analysis of the DoMYB gene family, including its identification, structural characteristics and expression patterns, was conducted by using whole-genome data. The results showed that 180 members of the D. odorifera MYB gene family were identified and could be classified into three categories based on the number of R domains: 1R-MYB (38), 2R-MYB (137), and 3R-MYB (5). These genes are unevenly distributed across the 10 chromosomes of D. odorifera, with some occurring in clusters. Phylogenetic analysis grouped 180 DoMYB genes into 23 subgroups; genes within the same subgroup share highly conserved gene structures and conserved motifs and exhibit similar evolutionary characteristics. Predictions of promoter cis-acting elements indicate that DoMYB family genes are extensively involved in light response, hormone signaling, and responses to biotic and abiotic stresses. Analysis of tissue expression patterns revealed that members of the DoMYB family exhibit significant tissue-specific expression in the phloem, sapwood, stem transition zone, and leaves of D. odorifera.. Among them, the key gene DoMYB48 showed the highest expression level in the stem transition zone of 7-year-old D. odorifera plants and was localized to the cell nucleus. In summary the DoMYB gene family in D. odorifera exhibits rich diversity in both gene structure and expression patterns. All these findings provide key candidate gene resources and a theoretical basis for elucidating the molecular regulatory mechanisms underlying the growth and development of D. odorifera.
-
图 2 DoMYB基因家族保守结构域、保守基序、基因结构分析
Fig. 2 Analysis of conserved domains, conserved motifs and gene structure of the DoMYB gene family
图 3 降香黄檀和拟南芥MYB基因家族系统进化树分析
注：图中At代表拟南芥AtMYB蛋白，Do代表降香黄檀DoMYB蛋白
Fig. 3 Phylogenetic analysis of the MYB gene families in D. odorifera and A. thaliana
Note: In the figure, At refers to A. thaliana AtMYB protein, and Do refers to D. odorifera DoMYB protein.
图 5 降香黄檀与拟南芥、大豆、水稻、毛果杨MYB家族基因的种间共线性关系
注：灰色线条代表共线区块，红色及蓝色线条代表共线MYB基因对。
Fig. 5 Interspecific syntenic relationships of MYB family genes among D. odorifera, A. thaliana, soybean, rice, and P. trichocarpa
Note: Gray lines represent collinear blocks, red and blue lines represent collinear MYB gene pairs.
图 6 降香黄檀MYB基因家族成员启动子顺式作用元件分析
Fig. 6 Analysis of the cis-acting elements of the promoter regions of MYB gene family members in D. odorifera
图 7 DoMYB基因家族表达模式分析
Fig. 7 Analysis of the expression pattern of the DoMYB gene family
图 8 降香黄檀DoMYB48基因在7年生降香黄檀中组织差异表达
注：采用单因素方差分析结合Bonferroni校正检验各组织间表达差异的显著性，****P<0.0001表示与过渡区组织相比的显著性差异
Fig. 8 Tissue-specific expression of the DoMYB48 gene in 7-year-old D. odorifera
Note: One-way ANOVA analysis of variance combined with Bonferroni correction was used to test the significance of expression differences among the various tissues. ****P < 0.0001 indicated a significant difference compared to the transitional zone tissues.
图 9 DoMYB48亚细胞定位分析
注：GFP：绿色荧光；mCherry：红色荧光细胞核marker；Bright：明场；Merge：叠加图像。
Fig. 9 Subcellular localization analysis of DoMYB48
Note: GFP: Green Fluorescent protein; mCherry: Cell nucleus marker - red fluorescence; Bright: Bright; Merge: Image showing the superimposition.
表 1 RT-qPCR引物序列
Table 1 Primer sequences used for RT-qPCR
基因名称
Gene基因ID
Gene ID引物F（5′−3′）
Forward primer引物R（5′−3′）
Reverse primerActin Daodo09aG0171500 CAATGAATTGCGTGTTGCT ATACCAGTTGTGCGACCACTT DoMYB48 Daodo02aG0332400 CCCTTCACCACCACCAGTTTC TTACCATCACCGTTATCATTAGCC 
下载: 导出CSV
表 2 DoMYB基因家族种内共线性Ka/Ks计算
Table 2 Intraspecific co-linearity Ka/Ks calculation of the DoMYB gene family
基因1（Gene 1） 基因2（Gene 2） Ka Ks Ka/Ks DoMYB1 DoMYB39 0.149 0.763 0.195 DoMYB2 DoMYB38 0.119 0.689 0.172 DoMYB3 DoMYB37 0.465 1.063 0.438 DoMYB9 DoMYB48 0.560 1.351 0.415 DoMYB2 DoMYB75 0.342 3.228 0.106 DoMYB17 DoMYB110 0.303 0.887 0.342 DoMYB15 DoMYB115 0.186 0.527 0.353 DoMYB14 DoMYB118 0.109 0.367 0.297 DoMYB21 DoMYB127 0.186 1.052 0.177 DoMYB6 DoMYB124 0.135 0.463 0.291 DoMYB8 DoMYB121 0.148 0.428 0.344 DoMYB32 DoMYB62 0.176 0.510 0.345 DoMYB50 DoMYB80 0.171 0.550 0.311 DoMYB51 DoMYB82 0.215 0.536 0.400 DoMYB30 DoMYB79 0.078 0.563 0.139 DoMYB24 DoMYB88 0.129 0.792 0.163 DoMYB26 DoMYB89 0.079 0.575 0.137 DoMYB27 DoMYB90 0.160 0.616 0.260 DoMYB28 DoMYB91 0.069 0.464 0.149 DoMYB35 DoMYB83 0.318 0.683 0.465 DoMYB44 DoMYB104 0.129 0.504 0.256 DoMYB45 DoMYB103 0.133 0.578 0.230 DoMYB46 DoMYB102 0.115 0.500 0.229 DoMYB40 DoMYB153 0.496 2.227 0.223 DoMYB53 DoMYB77 0.189 0.877 0.216 DoMYB55 DoMYB98 0.173 0.821 0.211 DoMYB56 DoMYB99 0.209 0.683 0.305 DoMYB66 DoMYB165 0.091 0.563 0.161 DoMYB67 DoMYB166 0.166 0.823 0.202 DoMYB57 DoMYB173 0.163 0.758 0.215 DoMYB59 DoMYB174 0.321 0.838 0.383 DoMYB60 DoMYB175 0.118 0.445 0.265 DoMYB61 DoMYB176 0.284 0.821 0.346 DoMYB63 DoMYB177 0.171 0.596 0.287 DoMYB64 DoMYB178 0.163 0.553 0.295 DoMYB65 DoMYB179 0.155 0.711 0.218 DoMYB53 DoMYB172 0.380 3.518 0.108 DoMYB53 DoMYB173 0.300 2.451 0.122 DoMYB70 DoMYB137 0.12 0.486 0.253 DoMYB71 DoMYB136 0.199 0.459 0.433 DoMYB78 DoMYB149 0.170 0.772 0.220 DoMYB84 DoMYB139 0.111 0.707 0.157 DoMYB92 DoMYB150 0.149 0.518 0.288 DoMYB93 DoMYB152 0.164 0.619 0.266 DoMYB94 DoMYB154 0.153 0.627 0.245 DoMYB108 DoMYB119 0.101 0.470 0.214 DoMYB107 DoMYB120 0.094 0.589 0.160 DoMYB123 DoMYB134 0.328 1.594 0.206 DoMYB143 DoMYB160 0.190 0.483 0.393 DoMYB146 DoMYB162 0.187 1.181 0.158 DoMYB140 DoMYB159 0.089 0.461 0.194 DoMYB141 DoMYB158 0.157 0.483 0.326
下载: 导出CSV
-
[1] 麦有专, 杨玉坤, 陈彧, 等. 热带珍贵树种降香黄檀研究进展[J]. 热带林业, 2018, 46(4): 25−28. https://doi.org/10.3969/j.issn.1672-0938.2018.04.007 doi: 10.3969/j.issn.1672-0938.2018.04.007 [2] 杨素丹, 王楠, 刘红梅, 等. 降香黄檀的开发价值及园林应用[J]. 湖北农业科学, 2019, 58(6): 78−81. https://doi.org/10.14088/j.cnki.issn0439-8114.2019.06.017 doi: 10.14088/j.cnki.issn0439-8114.2019.06.017 [3] 刘衡, 施福军, 黎韦水, 等. 降香黄檀与东京黄檀木材解剖特征及抽提物化学成分的可区别性[J]. 林业科学, 2021, 57(2): 103−114. https://doi.org/10.11707/j.1001-7488.20210211 doi: 10.11707/j.1001-7488.20210211 [4] 解雪华, 秦晨, 管崧淳, 等. 七叶树R2R3-MYB转录因子基因家族鉴定及表达分析[J]. 中草药, 2025, 56(20): 7512−7525. https://doi.org/10.7501/j.issn.0253-2670.2025.20.021 doi: 10.7501/j.issn.0253-2670.2025.20.021 [5] Wu X Y, Xia M, Su P, et al. MYB transcription factors in plants: a comprehensive review of their discovery, structure, classification, functional diversity and regulatory mechanism [J]. International Journal of Biological Macromolecules, 2024, 282: 136652. https://doi.org/10.1016/j.ijbiomac.2024.136652 doi: 10.1016/j.ijbiomac.2024.136652 [6] He Z R, Qin X F, Jia T, et al. Genome-wide identification of 1R-MYB transcription factors family and functional characterization of TrMYB130 under drought stresses in Trifolium repens (L. ) [J]. Gene, 2025, 943: 149247. https://doi.org/10.1016/j.gene.2025.149247 doi: 10.1016/j.gene.2025.149247 [7] Li S N, Wang W Y, Gao J L, et al. MYB75 phosphorylation by MPK4 is required for light-induced anthocyanin accumulation in Arabidopsis [J]. The Plant Cell, 2016, 28(11): 2866−2883. https://doi.org/10.1105/tpc.16.00130 doi: 10.1105/tpc.16.00130 [8] Zhang P, Wang R L, Ju Q, et al. The R2R3-MYB transcription factor MYB49 regulates cadmium accumulation [J]. Plant Physiology, 2019, 180(1): 529−542. https://doi.org/10.1104/pp.18.01380 doi: 10.1104/pp.18.01380 [9] Yang W B, Cortijo S, Korsbo N, et al. Molecular mechanism of cytokinin-activated cell division in Arabidopsis [J]. Science, 2021, 371(6536): 1350−1355. https://doi.org/10.1126/science.abe2305 doi: 10.1126/science.abe2305 [10] Teng S, Keurentjes J, Bentsink L, et al. Sucrose-specific induction of anthocyanin biosynthesis in Arabidopsis requires the MYB75/PAP1 gene [J]. Plant Physiology, 2005, 139(4): 1840−1852. https://doi.org/10.1104/pp.105.066688 doi: 10.1104/pp.105.066688 [11] Song S S, Qi T C, Huang H, et al. The Jasmonate-ZIM domain proteins interact with the R2R3-MYB transcription factors MYB21 and MYB24 to affect Jasmonate-regulated stamen development in Arabidopsis [J]. The Plant Cell, 2011, 23(3): 1000−1013. https://doi.org/10.1105/tpc.111.083089 doi: 10.1105/tpc.111.083089 [12] McCarthy R L, Zhong R Q, Fowler S, et al. The poplar MYB transcription factors, PtrMYB3 and PtrMYB20, are involved in the regulation of secondary wall biosynthesis [J]. Plant and Cell Physiology, 2010, 51(6): 1084−1090. https://doi.org/10.1093/pcp/pcq064 doi: 10.1093/pcp/pcq064 [13] Wang T T, Jin Y, Deng L X, et al. The transcription factor MYB110 regulates plant height, lodging resistance, and grain yield in rice [J]. The Plant Cell, 2024, 36(2): 298−323. https://doi.org/10.1093/plcell/koad268 doi: 10.1093/plcell/koad268 [14] Zhang S, Zhao Q C, Zeng D X, et al. RhMYB108, an R2R3-MYB transcription factor, is involved in ethylene- and JA-induced petal senescence in rose plants [J]. Horticulture Research, 2019, 6(1): 131. https://doi.org/10.1038/s41438-019-0221-8 doi: 10.1038/s41438-019-0221-8 [15] Xie Y P, Chen P X, Yan Y, et al. An atypical R2R3 MYB transcription factor increases cold hardiness by CBF-dependent and CBF-independent pathways in apple [J]. New Phytologist, 2018, 218(1): 201−218. https://doi.org/10.1111/nph.14952 doi: 10.1111/nph.14952 [16] Yin M M, Wei C H, Du H M, et al. Comprehensive analysis of R2R3-MYB transcription factors reveals OsMYB1 as a key regulator of anthocyanin biosynthesis in rice [J]. Plant Science, 2025, 360: 112732. https://doi.org/10.1016/j.plantsci.2025.112732 doi: 10.1016/j.plantsci.2025.112732 [17] Paysan-Lafosse T, Andreeva A, Blum M, et al. The Pfam protein families database: embracing AI/ML [J]. Nucleic Acids Research, 2025, 53(D1): D523−D534. https://doi.org/10.1093/nar/gkae997 doi: 10.1093/nar/gkae997 [18] Reiser L, Proia A, Bakker E, et al. Recent major changes to TAIR: updates to the database, website, and Arabidopsis genome [J]. Genetics, 2026, 232(4): iyaf248. https://doi.org/10.1093/genetics/iyaf248 doi: 10.1093/genetics/iyaf248 [19] Wang J Y, Chitsaz F, Derbyshire M K, et al. The conserved domain database in 2023 [J]. Nucleic Acids Research, 2023, 51(D1): D384−D388. https://doi.org/10.1093/nar/gkac1096 doi: 10.1093/nar/gkac1096 [20] Duvaud S, Gabella C, Lisacek F, et al. Expasy, the Swiss Bioinformatics Resource Portal, as designed by its users [J]. Nucleic Acids Research, 2021, 49(W1): W216−W227. https://doi.org/10.1093/nar/gkab225 doi: 10.1093/nar/gkab225 [21] 刘艳丽, 周媛, 曹丹, 等. 基于茶蛋白质组学数据分析植物亚细胞定位预测软件的应用[J]. 植物科学学报, 2020, 38(5): 671−677. https://doi.org/10.11913/PSJ.2095-0837.2020.50671 doi: 10.11913/PSJ.2095-0837.2020.50671 [22] Chen C J, Wu Y, Li J W, et al. TBtools-II: a “one for all, all for one” bioinformatics platform for biological big-data mining [J]. Molecular Plant, 2023, 16(11): 1733−1742. https://doi.org/10.1016/j.molp.2023.09.010 doi: 10.1016/j.molp.2023.09.010 [23] Zhao W X, Meng X X, Xu J H, et al. Integrated mRNA and small RNA sequencing reveals microRNAs associated with xylem development in Dalbergia odorifera [J]. Frontiers in Genetics, 2022, 13: 883422. https://doi.org/10.3389/fgene.2022.883422 doi: 10.3389/fgene.2022.883422 [24] 刘忠丽. 拟南芥AtMYB123与AtKOR1基因突变体的筛选鉴定及其功能研究[D]. 杭州: 浙江大学, 2011. [25] 王梦琦, 解振强, 孙欣, 等. 葡萄miR159及其靶基因VvGAMYB在花发育过程中的作用分析[J]. 园艺学报, 2017, 44(6): 1061−1072. https://doi.org/10.16420/j.issn.0513-353x.2016-0816 doi: 10.16420/j.issn.0513-353x.2016-0816 [26] Zhang D W, Zhou H P, Zhang Y, et al. Diverse roles of MYB transcription factors in plants [J]. Journal of Integrative Plant Biology, 2025, 67(3): 539−562. https://doi.org/10.1111/jipb.13869 doi: 10.1111/jipb.13869 [27] Chen Y H, Yang X Y, He K, et al. The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the Rice MYB family [J]. Plant Molecular Biology, 2006, 60(1): 107−124. https://doi.org/10.1007/s11103-005-2910-y doi: 10.1007/s11103-005-2910-y [28] Du H, Yang S S, Liang Z, et al. Genome-wide analysis of the MYB transcription factor superfamily in soybean [J]. BMC Plant Biology, 2012, 12(1): 106. https://doi.org/10.1186/1471-2229-12-106 doi: 10.1186/1471-2229-12-106 [29] Wilkins O, Nahal H, Foong J, et al. Expansion and diversification of the Populus R2R3-MYB family of transcription factors [J]. Plant Physiology, 2009, 149(2): 981−993. https://doi.org/10.1104/pp.108.132795 doi: 10.1104/pp.108.132795 [30] Sukumaran S, Lethin J, Liu X, et al. Genome-wide analysis of MYB transcription factors in the wheat genome and their roles in salt stress response [J]. Cells, 2023, 12(10): 1431. https://doi.org/10.3390/cells12101431 doi: 10.3390/cells12101431 [31] 朱英洁, 赵航, 包颖. 棉花MYB转录因子的研究进展[J]. 曲阜师范大学学报(自然科学版), 2021, 47(4): 105−111. https://doi.org/10.3969/j.issn.1001-5337.2021.4.105 doi: 10.3969/j.issn.1001-5337.2021.4.105 [32] Hanada K, Zou C, Lehti-shiu M D, et al. Importance of lineage-specific expansion of plant tandem duplicates in the adaptive response to environmental stimuli [J]. Plant Physiology, 2008, 148(2): 993−1003. https://doi.org/10.1104/pp.108.122457 doi: 10.1104/pp.108.122457 [33] Radwan N S, Lamlom S F, Emwas A H, et al. Genome-wide identification, characterization, and evolutionary analysis of the HSP70 gene family in rice (Oryza sativa L. ) [J]. BMC Plant Biology, 2025, 25(1): 1671. https://doi.org/10.1186/s12870-025-07574-8 doi: 10.1186/s12870-025-07574-8 [34] 陈艳丽, 朱桐杉, 徐国鑫, 等. 柴胡R2R3-MYB转录因子家族鉴定及其盐胁迫下的表达分析[J]. 中草药, 2026, 57(7): 2704−2717. https://doi.org/10.7501/j.issn.0253-2670.2026.07.022 doi: 10.7501/j.issn.0253-2670.2026.07.022 [35] Bai X T, Xu J M, Shao X M, et al. A novel gene coding γ-aminobutyric acid transporter may improve the tolerance of Populus euphratica to adverse environments [J]. Frontiers in Plant Science, 2019, 10: 1083. https://doi.org/10.3389/fpls.2019.01083 doi: 10.3389/fpls.2019.01083 [36] Li J, Lin K Y, Zhang S, et al. Genome-wide analysis of Myeloblastosis-related genes in Brassica napus L. and positive modulation of osmotic tolerance by BnMRD107 [J]. Frontiers in Plant Science, 2021, 12: 678202. https://doi.org/10.3389/fpls.2021.678202 doi: 10.3389/fpls.2021.678202 [37] 阙秋霞, 赖恭梯, 潘若, 等. 刺葡萄MYB基因家族鉴定及其在不同光质下的表达模式分析[J]. 农业生物技术学报, 2023, 31(2): 298−310. https://doi.org/10.3969/j.issn.1674-7968.2023.02.008 doi: 10.3969/j.issn.1674-7968.2023.02.008 [38] Gigolashvili T, Engqvist M, Yatusevich R, et al. HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana [J]. New Phytologist, 2008, 177(3): 627−642. https://doi.org/10.1111/j.1469-8137.2007.02295.x doi: 10.1111/j.1469-8137.2007.02295.x [39] Ko J H, Jeon H W, Kim W C, et al. The MYB46/MYB83-mediated transcriptional regulatory programme is a gatekeeper of secondary wall biosynthesis [J]. Annals of Botany, 2014, 114(6): 1099−1107. https://doi.org/10.1093/aob/mcu126 doi: 10.1093/aob/mcu126 [40] Öhman D, Demedts B, Kumar M, et al. MYB103 is required for FERULATE-5-HYDROXYLASE expression and syringyl lignin biosynthesis in Arabidopsis stems [J]. The Plant Journal, 2013, 73(1): 63−76. https://doi.org/10.1111/tpj.12018 doi: 10.1111/tpj.12018 [41] Shi J W, Zhang Q X, Yan X, et al. A conservative pathway for coordination of cell wall biosynthesis and cell cycle progression in plants [J]. The Plant Journal, 2021, 106(3): 630−648. https://doi.org/10.1111/tpj.15187 doi: 10.1111/tpj.15187 [42] Yang X T, Wang K, Bu Y N, et al. The transcription factor TaGAMYB modulates tapetum and pollen development of TGMS wheat YanZhan 4110S via the gibberellin signaling [J]. Plant Science, 2022, 324: 111447. https://doi.org/10.1016/j.plantsci.2022.111447 doi: 10.1016/j.plantsci.2022.111447 [43] Aya K, Ueguchi-Tanaka M, Kondo M, et al. Gibberellin modulates anther development in rice via the transcriptional regulation of GAMYB [J]. The Plant Cell, 2009, 21(5): 1453−1472. https://doi.org/10.1105/tpc.108.062935 doi: 10.1105/tpc.108.062935 [44] Gocal G F W, Sheldon C C, Gubler F, et al. GAMYB-like genes, flowering, and gibberellin signaling in Arabidopsis [J]. Plant Physiology, 2001, 127(4): 1682−1693. https://doi.org/10.1104/pp.010442 doi: 10.1104/pp.010442 [45] 陈水莲, 梁远楠, 周莹, 等. 降香黄檀心材与树龄关系研究[J]. 南方农业, 2016, 10(4): 43−45. https://doi.org/10.19415/j.cnki.1673-890x.2016.04.016 doi: 10.19415/j.cnki.1673-890x.2016.04.016 [46] Zhuang H B, Chong S L, Priyanka B, et al. Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir) [J]. BMC Plant Biology, 2021, 21(1): 581. https://doi.org/10.1186/s12870-021-03322-w doi: 10.1186/s12870-021-03322-w -
点击查看大图
计量
- 文章访问数: 6
- HTML全文浏览量: 3
- 被引次数: 0
