-
冬瓜（Benincasa hispida（Thunb.）Cogn.）是葫芦科（Cucurbitaceae）冬瓜属（Benincasa Savi）一年生蔓性草本植物[1]。我国冬瓜种植面积常年维持在33万 hm2以上[2]。近几年，由于复种指数的提高和冬瓜品种的单一化趋势，冬瓜土传病害的发生越来越严重和频繁，其中冬瓜枯萎病发病较为严重[3]，该病发病轻时减产20%～30%，重时减产在60%以上[4]。因此，如何防治枯萎病已经成为保证冬瓜产业健康可持续发展的关键。
冬瓜枯萎病的发生涉及病原菌、寄主和环境条件三者的相互作用。枯萎病致病菌一旦在土壤成功繁殖几乎不可能根除[5]，因此防治难度极大。嫁接是防治冬瓜枯萎病的有效措施之之一。本课题组前期研究发现‘海砧1号’南瓜（Cucurbita spp.）具高抗枯萎病特性[6]，用其做砧木的嫁接冬瓜也具有高抗枯萎病特性，病情指数为0[7]，且产量[8]、氮吸收效率[9]比传统黑籽南瓜砧木嫁接的冬瓜要佳。课题组前期分别从化感作用和生理生化的角度研究了该砧木嫁接提高冬瓜枯萎病抗性的原因[6,10 − 11]，但枯萎病抗性提高的分子机制目前还不清楚。
植物对病原物侵染的抗性是一系列基因活化和表达的结果。在寄主-病原菌互作中，通过识别、信号传递和诱导植物防卫反应基因的表达，从而使植物抗病性得到表现。转录组测序技术是挖掘植物抗逆基因，揭示植物分子抗性机制的有效手段之一。为此，本研究以嫁接冬瓜作为研究对象，自接冬瓜与自接南瓜为对照，通过转录组测序技术分析其幼苗根部和叶片中基因表达的变化，挖掘差异抗病相关基因，为解析南瓜嫁接提高冬瓜枯萎病抗性的分子机制，同时为冬瓜新品种的选育提供理论基础。
-
冬瓜品种选用海南主栽品种‘铁柱2号’，南瓜品种选用‘海砧1号’，由海南农业科学院蔬菜研究所提供。
育苗地点设在海南省农业科学院澄迈基地。2023年5月24日浸种冬瓜种子用于培养自接冬瓜的砧木，6月2日浸种冬瓜种子用于培养自接冬瓜和嫁接冬瓜的接穗，6月6日浸种南瓜种子用于培养自接南瓜和嫁接冬瓜的砧木，6月9日浸种南瓜种子用于培养自接南瓜的接穗，浸种时间均为12 h，砧木种子露白长度达到种子本身1/2播于60孔穴盘，接穗种子露白长度达到种子本身1/2播于托盘中。待接穗长至两片子叶平展，砧木第一片真叶成梭状时采用插接法嫁接与自接[12]，得到嫁接冬瓜、自接冬瓜和自接南瓜3种幼苗，分别用GW、SW和SP表示。3种幼苗培养至两叶一心期备用。
-
实验菌种由广西农业科学院蔬菜研究所提供。该菌种为尖孢镰刀菌冬瓜专化型（Fusarium oxysporum Schl. F. sp. benincasae）生理小种1号（后用FOB表示），从冬瓜‘铁柱2号’中分离获得。
PDA与PDB培养基的制备参照[13 − 14]。FOB于PDA培养基中培养7 d后，挑取菌丝转至PDB培养基，置25 ℃、150 r·min−1转速的恒温振荡器培养3 d，脱脂棉过滤得到孢子悬浮液，测定得到孢子浓度为2.75×107 CFU·g−1，按1∶4的比例稀释孢子悬浮液得稀释液，备用。
-
土壤采自澄迈县瑞溪镇罗浮村，基本理化性质为：pH 5.94、有机质5.1 g·kg−1、碱解氮28.1 mg·kg−1、速效磷103.3 mg·kg−1与速效钾129.6 mg·kg−1。拣出实验土壤杂物，过0.18 mm筛子，灭菌3次，每次2 h[15]，风干后得灭菌土壤，备用。
将孢子稀释液和水分别按土壤∶液体= 6 g∶1 ml混合拌匀，得接菌土壤。测土壤FOB浓度，为8.0×105 CFU·g−1和0 CFU·g−1（用F和0表示）。装接菌土壤，每钵120 g，备用。
-
实验设3类幼苗，分别为GW、SW与SP，2种浓度接菌土壤，分别为F和0。选长势一致的3种幼苗，清洗根系后，定植在装有接菌土壤的方形塑料钵中，每钵1株。6种接菌处理的幼苗分别用GWF、GW0、SWF、SW0、SPF、SP0表示。GWF、SWF和SPF处理幼苗各25株，GW0、SW0和SP0处理幼苗各10株。
定植后将幼苗置光照培养箱培养，培养条件：光照周期12 h/12 h，温度26 ℃/21 ℃（昼/夜），光强为
2000 Lux，培养前2 d无需浇水，后每2 d补充植株所蒸发的水分，实验过程中不施肥。 -
课题组前期通过不同浓度病菌的接菌实验发现，接菌土壤FOB浓度为2.8 × 106 CFU·g−1时，自接冬瓜幼苗第6 d开始出现枯萎的症状，接菌土壤FOB浓度为8.0 × 105 CFU·g−1时，第10 d开始出现枯萎的症状，而嫁接冬瓜均未出现枯萎的症状[11]，故本研究选择病菌入侵阶段（第5 d）开展幼苗根系与叶片的转录组分析。
培养5 d后对GWF、GW0、SWF和SW0幼苗叶片（分别用GWFL、GWL、SWFL和SWL表示）和GWF、GW0、SPF和SP0幼苗根部（分别用GWFR、GWR、SPFR和SPR表示）进行采样，样本均设4个重复，每组4个生物学重复命名为 1、2、3、4。先用自来水冲洗第3～4片叶片的正反面（根部用自来水冲洗根部表面），再用去离子水再次冲洗、立即用滤纸吸干残留的水分，用灭菌后的剪刀剪下叶片或根部后置于锡箔纸中，去除叶片主脉后剪碎（根部首尾各1 cm，将中间根部剪碎），置于2 mL冻存管中，随后立即于液氮中速冻，后转到−80 ℃超低温冰箱中保存，随后立即将样品交由上海美吉生物医药科技有限公司完成RNA-seq测序工作。
-
使用fastp v 0.19.5软件（
https://github.com/OpenGene/fastp ）对原始数据进行过滤，过滤掉带接头、首尾端低质量、碱基N含量过高以及舍弃去adapter及质量修剪后长度小于30 bp的序列，从而得到高质量的测序数据（clean data）以保证后续分析的顺利进行。使用HiSat2v2.1.0软件（
http://ccb.jhu.edu/software/hisat2/index.shtml ）将质控后的clean data （reads），与参考基因组比对获得用于后续转录本组装、表达量计算等的mapped data （reads），叶片参考基因来源: Benincasa_hispida （https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_009727055.1/ ），根系参考基因来源：Cucurbita_moschata （https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_002738365.1/ ），同时对该次转录组测序的比对结果进行质量评估，主要包括测序饱和度、基因覆盖度、Reads在参考基因组不同区域分布以及Reads在不同染色体分布分析。 -
获得基因的Read Counts数后，使用DESeq2 v 1.10.1软件（
http://biocon-ductor.org/pac-kages/stats/bioc/DESeq2 ）进行样本间基因的表达差异分析，将|log2FC| ≥ 1且FDR （False Discovery Rate）<0.05的基因定义为差异表达基因（differentially expressed genes, DEGs）。 -
利用GO数据库，将基因按照生物学过程（Biological Process, BP）、细胞组分（Cellular Component, CC）和分子功能（Molecular Function, MF）进行分类，对差异表达基因的主要功能进行GO注释。
利用KEGG数据库，将基因按照pathway通路或行使的功能进行分类，对差异表达基因进行KEGG注释。
利用植物抗性基因数据库（PRGdb;
http://prgdb.org/prgdb4/ ），将差异表达基因的蛋白序列提交到DRAGO 3工具中预测R基因的类别及结构域，进一步筛选抗病基因。此外，GO功能富集和KEGG功能富集分析分别通过软件Goatools v 0.6.5 （
https://git-hub.com/tanghaibao/GOatool ）和Python scipy （https://scipy.org/install/ ）进行，使用 4种多重检验方法 （BH, BY, holm 和 bonferroni）对p值进行了校正，当p值≤0.05时，认为GO与KEGG功能存在显著富集情况。 -
选取叶片4个、根部2个差异表达基因进行实时荧光定量PCR，验证转录组结果的准确性。
以F-box protein基因作为叶片的内参基因[16]，以Actin基因作为根部的内参基因[17]进行qRT-PCR反应。实验设3个重复，使用2−∆∆ct法计算样本基因相对表达量，使用Origin 2024软件进行绘图。利用NCBI中的Pick Primers设计引物（表1），使用Blast进行引物特异性检验，将设计好的引物交于擎科生物科技有限公司海南分公司合成。
表 1 叶片和根部基因qRT-PCR引物
Table 1. qRT-PCR primers for leaf and root genes
部位
Tissue基因
Gene上游引物序列
Forward primer sequence（5'→3'）下游引物序列
Reverse primer sequence（5'→3'）叶片
LeafLOC120076313 CGCTTGTAGCACAGCGAAAG TCGATTATCTTGGGCGGGAC LOC120067870 TGCAGACGGCATAGACAGTG TGCTATTCGGCCGTCTCTTC LOC120086723 TCACTGCTGATCCCCGAAAC GAAACTCAATGCTCGCCACC LOC120075868 GAAACGCATCTGGAAGTGCC ACGTAGGGTCCTGCTCATCA F-box protein CAAACTCCCCAATGCTCAGT CAAACAAGACCACCAGATGAAC 根部
RootLOC111431282 AACACCAAGCCAACTCGGAA AGAGGTTGCCGAAGGAGTTC LOC111434154 AGATGACCCGGAAGCCTTTG GTTGACTGGGTCAGCAAGGT Actin AGCCATCTCTCATCGGTAT CATGGTTGAACCACCACTG -
16个叶片样品测序共获得原始序列700 019 546条，原始数据经过滤后获得高质量的reads数共649 983 334条，共获得103.28 Gb Clean Data，各样品Clean Data均达到6.04 Gb以上。Q20碱基百分比含量的平均值为97.73％，Q30碱基百分比含量的平均值为95.88％，GC的含量的平均值为45.89％（表2）。16个根部样品测序共获得原始序列796 045 970条，原始数据经过滤后获得高质量的reads数共788 820 340条，共获得118.44 Gb Clean Data，各样品Clean Data均达到6.29 Gb以上。Q20碱基百分比含量的平均值为97.90％，Q30碱基百分比含量的平均值为93.75％，GC的含量的平均值为45.93％（表3），说明文库构建质量较高，可满足下一步的数据分析。
表 2 叶片测序数据统计表
Table 2. Statistical table of leaf sequencing data
处理
Treatment生物学重复
Biological repeat高质量测序读段
Clean reads高质量碱基数
Clean bases错误率
Error rate（%）Q20（%）
Q20（%）Q30（%）
Q30（%）GC含量
GC content（%）SWFL 1 43 179 818 6 458 447 311 0.0253 97.30 95.21 46.33 2 42 651 922 6 376 573 009 0.0245 97.70 95.82 45.16 3 42 875 394 6 406 161 037 0.0249 97.50 95.51 45.24 4 45 009 320 6 752 776 974 0.0263 96.83 94.48 45.98 GWFL 1 43 320 074 6 459 906 009 0.0241 97.92 96.16 45.05 2 43 053 652 6 427 466 383 0.0241 97.92 96.15 44.86 3 42 666 922 6 322 273 436 0.0241 97.91 96.17 45.65 4 42 812 064 6 373 731 835 0.0241 97.91 96.18 45.34 SWL 1 42 549 842 6 340 332 373 0.0241 97.91 96.16 45.45 2 43 429 956 6 458 514 075 0.0241 97.91 96.19 45.85 3 47 369 184 7 007 856 484 0.0239 97.98 96.35 45.87 4 43 195 524 6 390 480 048 0.0247 97.58 95.76 46.07 GWL 1 49 023 592 7 288 104 998 0.0243 97.83 96.02 46.81 2 40 544 988 6 041 735 660 0.0244 97.8 95.94 46.68 3 40 588 474 6 060 309 165 0.0243 97.84 96.02 47.10 4 41 032 682 6 117 046 724 0.0244 97.80 95.92 46.78 表 3 根部测序数据统计表
Table 3. Statistical table of root sequencing data
处理
Treatment生物学重复
Biological repeat高质量测序读段
Clean reads高质量碱基数
Clean bases错误率
Error rate（%）Q20（%）
Q20（%）Q30（%）
Q30（%）GC含量
GC content（%）SPFR 1 55 068 966 8 274 790 336 0.0253 97.95 93.88 46.00 2 41 816 552 6 288 523 427 0.0250 98.05 94.16 45.84 3 56 613 890 8 512 784 710 0.0252 98.00 93.97 45.69 4 53 390 314 8 004 266 521 0.0254 97.91 93.79 45.94 GWFR 1 49 987 502 7 494 885 674 0.0250 98.03 94.1 45.93 2 45 607 438 6 842 619 711 0.0251 97.99 94.01 46.14 3 42 957 918 6 456 667 868 0.0253 97.92 93.84 46.02 4 51 699 286 7 766 345 815 0.0256 97.83 93.57 45.95 SPR 1 47 556 512 7 125 387 830 0.0250 98.03 94.09 46.07 2 48 154 996 7 207 863 128 0.0260 97.64 93.17 45.96 3 53 055 568 7 969 322 838 0.0256 97.83 93.57 45.94 4 52 423 320 7 867 678 851 0.0255 97.86 93.70 45.87 GWR 1 50 296 086 7 558 923 882 0.0256 97.81 93.51 45.69 2 46 241 656 6 952 236 167 0.0257 97.76 93.42 45.92 3 45 286 002 6 803 395 798 0.0256 97.82 93.55 45.98 4 48 664 334 7 309 752 018 0.0254 97.90 93.71 45.93 将叶片组织样品比对到冬瓜参考基因组后，利用HiSat2软件分别将各样品的Clean Reads与指定的参考基因组进行序列比对，比对率从96.65%到97.64%不等。该比对率较高，可以保证后续分析的准确性。与参考基因组相比，最终鉴定得到检测到表达基因共22 808个，其中已知基因有21 121个，新基因1 687个。将根部组织样品比对到南瓜参考基因组后，利用HiSat2软件分别将各样品的Clean Reads与指定的参考基因组进行序列比对，比对率从72.43%到86.66%不等。该比对率较高，可以保证后续分析的准确性。与参考基因组相比，最终鉴定得到检测到表达基因共29 436个，其中已知基因有28 566个，新基因870个。
-
为了检验样本的合理性，分别对叶片与根部各样品间所有基因表达量进行相关性分析。如图1所示，叶片生物学重复样品间Pearson相关系数的平方值（R2）均在0.8以上（图1−A），根部生物学重复样品间Pearson相关系数的平方值（R2）均在0.9以上（图1−B），说明二者样品之间生物学重复性良好，可以满足后续分析工作。
图 1 叶片（A）和根部（B）样品相关性热图
Figure 1. Sample correlation heat map of leaf（A）and root（B）
-
将|log2FC| ≥ 1且FDR < 0.05的基因定义为差异表达基因。在叶片当中，GWFL vs SWFL组共鉴定出4 310个DEGs，其中包含1 956个上调DEGs和2 354个下调DEGs；GWFL vs GWL组共鉴定出6 808个DEGs，其中包含3 130个上调DEGs和3 678个下调DEGs；SWFL vs SWL组共鉴定出3 759个DEGs，其中包含2 290个上调DEGs和1 469个下调DEGs。在根部当中，GWFR vs SPFR组共鉴定出56个DEGs，其中包含49个上调DEGs和7个下调DEGs；GWFR vs GWR组共鉴定出803个DEGs，其中包含403个上调DEGs和401个下调DEGs；SPFR vs SPR组共鉴定出242个DEGs，其中包含115个上调DEGs和127个下调DEGs。
为进一步筛选关键抗枯萎病基因，本研究以接菌后在嫁接冬瓜中上调，且在自接冬瓜或南瓜中下调的差异表达基因定义为潜在抗病基因[18]，并对其进行进一步分析。对叶片GWFL vs GWL组的3 130个上调DEGs、GWFL vs SWFL组的1 956个上调DEGs与SWFL vs SWL组的
1469 个下调DEGs进行Venn分析，筛选出284个叶片潜在抗病基因（图2−A）。对根部GWFR vs SPFR组的49个上调DEGs、GWFR vs GWR组的402个上调DEGs与SPFR vs SPR组的127个下调DEGs进行Venn分析，筛选出2个根部潜在抗病基因（图2−B）。
图 2 差异表达基因维恩图（A：叶片，B：根部）
Figure 2. Venn map of DEGs in leaf（A）and root（B）
-
对所筛选出来的284叶片关键抗病基因进行GO功能分析，在生物过程中，叶片的关键抗病基因主要富集在细胞过程、代谢过程与生物调节；在细胞组分中，叶片的关键抗病基因主要富集在细胞组分、膜组分和细胞器；在分子功能中，叶片的关键抗病基因主要富集在催化活性与结合（图3−A）。
图 3 潜在抗病基因的GO功能分类（A：叶片，B：根部）
Figure 3. GO functional classification of potential disease resistance genes in leaf（A）and root（B）
同时对所筛选出来的2个根部关键抗病基因进行GO功能分析，在生物过程中，根部的关键抗病基因主要富集在代谢过程、细胞过程和应激反应；在细胞组分中，根部的关键抗病基因主要富集在含蛋白质复合物、膜组分、细胞器组分、细胞膜和细胞组分；在分子功能中，根部的关键抗病基因主要富集在催化活性与结合（图3−B）。
-
叶片关键抗病基因在丙氨酸、天门冬氨酸和谷氨酸代谢（Alanine, aspartate and glutamate metabolism）、苯丙烷生物合成（Phenylpropanoid biosynthesis）、角质、亚蜡和蜡的生物合成（Cutin, suberine and wax biosynthesis）与氮代谢（Nitrogen metabolism）等通路显著富集（图4−A）。
图 4 潜在抗病基因的KEGG富集分析（A：叶片，B：根部）
Figure 4. KEGG enrichment analysis of potential disease resistance genes in leaf（A）and root（B）
根部的关键抗病基因在芪类、二芳基庚烷类和姜酚生物合成（Stilbenoid, diarylheptanoid and gingerol biosynthesis）、光合作用天线蛋白（photosynthesis-antenna proteins）、类黄酮生物合成（Flavonoid biosynthesis）与苯丙烷生物合成（Phenylpropanoid biosynthesis）等通路显著富集（图4−B）。
叶片中丙氨酸、天门冬氨酸和谷氨酸代谢（Alanine, aspartate and glutamate metabolism）通路涉及到3个基因，分别是LOC120076313、LOC120089431与LOC120083554。苯丙烷生物合成（Phenylpropanoid biosynthesis）涉及到5个基因，分别是LOC120067870、LOC120073918、LOC120081818、LOC120085310与LOC120091171。角质、亚蜡和蜡生物合成（Cutin, suberine and wax biosynthesis）涉及到2个基因，分别是LOC120086723与LOC120071639。氮代谢（Nitrogen metabolism）也涉及到2个基因，分别为LOC120074685与LOC120075868（表4）。
表 4 叶片KEGG富集分析详情表
Table 4. Detailed table of leaf KEGG enrichment analysis
代谢通路描述
Pathway description通路编号
Pathway ID基因数量
Gene number基因
Gene基因功能描述
Gene description丙氨酸、天门冬氨酸和谷氨酸代谢
Alanine, aspartate and
glutamate metabolismmap00250 3 LOC120076313 Asparagine synthetase [glutamine-hydrolyzing]1-like LOC120089431 Asparagine synthetase [glutamine-hydrolyzing] 1-like, transcript variant X1 LOC120083554 Probable aldehyde dehydrogenase, transcript variant X1 苯丙烷生物合成
Phenylpropanoid biosynthesismap00940 5 LOC120067870 Peroxidase 2-like LOC120073918 Berberine bridge enzyme-like 15 LOC120081818 Peroxidase 11, transcript variant X1 LOC120085310 4-coumarate--CoA ligase 2-like LOC120091171 Caffeoylshikimate esterase-like, transcript variant X1 角质、亚蜡和蜡生物合成
Cutin, suberine and wax biosynthesismap00073 2 LOC120086723 Omega-hydroxypalmitate O-feruloyl transferase-like LOC120071639 Cytochrome P450 86B1-like 氮代谢
Nitrogen metabolismmap00910 2 LOC120074685 Uncharacterized LOC120074685 LOC120075868 Beta carbonic anhydrase 5, chloroplastic-like, transcript variant X1 根部中芪类、二芳基庚烷类和姜酚生物合成（Stilbenoid, diarylheptanoid and gingerol biosynthesis）、光合作用天线蛋白（photosynthesis-antenna proteins）、类黄酮生物合成（Flavonoid biosynthesis）与苯丙烷生物合成（Phenylpropanoid biosynthesis）等通路均涉及到1个基因LOC111431282，光合作用天线蛋白（photosynthesis-antenna proteins）涉及到1个基因为LOC111434154（表5）。因此推测LOC111431282基因可能在抵御尖孢镰刀菌侵染发挥重要的作用。
表 5 根部KEGG富集分析详情表
Table 5. Detailed table of root KEGG enrichment analysis
代谢通路描述
Pathway description通路编号
Pathway ID基因数量
Gene number基因
Gene基因功能描述
Gene description芪类、二芳基庚烷类和姜酚生物合成 Stilbenoid, diarylheptan-oid and gingerol biosynthesis map00945 1 LOC111431282 Vinorine synthase-like 光合作用天线蛋白Photosynthesis-antenna proteins map00196 1 LOC111434154 Chlorophyll a-b binding protein of LHCII type 1-like 类黄酮生物合成 Flavonoid biosynthesis map00941 1 LOC111431282 Vinorine synthase-like 苯丙烷生物合成Phenylpropanoid biosynt-hesis map00940 1 LOC111431282 Vinorine synthase-like 从KEGG富集分析可知，苯丙烷生物合成是冬瓜叶片与南瓜根部共有的富集代谢通路。因此对苯丙烷生物合成进行代谢通路分析（ko00904）。该通路以苯丙氨酸化学分子为底物，经过苯丙氨酸解氨酶（Phenylalanine Ammonia-Lyase, PAL）催化作用变成肉桂酸（Cinnamic acid），而后经过3条不同的代谢通路形成对-羟基苯基木质素（p-hydroxyphenyl lignin, H木质素）、愈创木基木质素（Guaiacyl lignin, G木质素）和紫丁香基木质素（Syringyl lignin, S木质素）3种不同的木质素（图5）。
图 5 苯丙烷生物合成通路图（A：叶片，B：根部）
Figure 5. Map of phenylpropanoid biosynthesis pathway in leaf（A）and root（B）
-
分别对嫁接冬瓜叶片284个及根部2个潜在抗病基因进行蛋白序列提取，利用PRGdb网站的DRAGO 3工具对蛋白结构域进行预测。在叶片284个潜在抗病基因中，共预测到16与R基因结构域相关的基因，根据其不同的结构域主要分为7类，分别为CK、CL、KIN、LECRK、N、RLP、TRAM（表6），而根部2个潜在抗病基因没有预测到与R基因相关的结构域。
表 6 叶片潜在抗病基因蛋白结构域预测
Table 6. Prediction of protein domain of potential disease resistance genes in leaves
基因 Gene 分类 Class 结构域 Domain 基因功能描述 Description LOC120087570 CK CC,TM,Kinase Contains a coiled-coil domain（CC）,a transmembrane
domain（TM）and a kinase domainLOC120090355 CL CC,LRR,TM Contains a coiled-coil domain,a leucine-rich repeat
domain（LRR）and a kinase domain（Kinase）LOC120076759 KIN TM,Kinase Contains a transmembrane domain and a kinase domain LOC120083806 KIN TM,Kinase LOC120082723 KIN TM,Kinase LOC120088593 KIN TM,Kinase LOC120089817 KIN TM,Kinase LOC120092573 KIN TM,Kinase LOC120072525 KIN TM,Kinase LOC120090305 LECRK LECM,TM,Kinase Lectin-receptor kinase, contains a lectin-like motifs（LECM）
instead of LRR,a transmembrane domain and a kinase domainLOC120090267 LECRK LECM,TM,Kinase LOC120090948 N NBS,TM Contains a nucleotide-binding site domain（NBS）domain and a
transmembrane domainLOC120085520 N NBS,TM LOC120072345 RLP LRR,TM Receptor like proteins, contains a leucine-rich repeat domain and
a kinase domainLOC120085476 RLP LRR,TM LOC120085949 TRAM TM Contains a transmembrane domain only -
对叶片所有转录组序列进行TFs分析，共鉴定出2 828个差异表达的TFs，这些TFs属于48个基因家族，主要包括MYB_related、bHLH、MYB、HB-other、ERF、bZIP、MIKC、B3、WRKY、NAC与M_type等（图6−A）。对根部所有转录组序列进行TFs分析，共鉴定出5 426个差异表达的TFs，这些TFs也属于48个基因家族，主要包括MYB_related、bHLH、MYB、ERF、NAC、bZIP、HB-other、ARF、MIKC、B3、WRKY、M_type、FAR1、HD-ZIP、C3H、GRAS与Dof等（图6−B）。
图 6 尖孢镰刀菌侵染下响应的TFs（A：叶片，B：根部）
Figure 6. TFs in response to Fusarium oxysporum infection in leaf（A）and root（B）
叶片与根部共有的基因家族有MYB_related、bHLH、MYB、ERF、NAC、bZIP、HB-other、MIKC、B3、WRKY与M_type。这些TFs调控不同靶基因的表达，可能参与抵抗病原菌入侵反应，在抵御冬瓜枯萎病菌侵染中起着重要作用。
-
分别选取叶片4个DEGs与根部2个DEGs进行qRT-PCR以验证转录组结果的准确性。结果表明，除冬瓜叶片LOC120086723基因的相对表达量在某些处理间的变化趋势与RNA-seq结果有出入，其他的冬瓜叶片3个DEGs与南瓜根部2个DEGs的相对表达量表现与其对应的RNA-seq结果一致（图7），说明本研究转录组测序的结果可靠。
图 7 冬瓜叶片4个基因（A−D）南瓜根部2个基因（E−F）荧光定量实验结果
Figure 7. qRT-PCR verification of 4 genes from wax gourd in leaves（A-D）and in roots（E-F）
-
本研究将在嫁接冬瓜中上调，且在自接冬瓜中下调的差异表达基因定义为潜在抗病基因，分别在叶片与根部中筛选出了284个和2个潜在抗病基因。对叶片与根部的潜在抗病基因进行GO与KEGG分析，发现苯丙烷生物合成是嫁接冬瓜叶片与根部共有的代谢通路。苯丙烷代谢的限速酶是PAL，是连接植物新陈代谢中初生代谢和次生代谢的桥梁，在植物木质素、类黄酮、香豆酸酯类等次生代谢物的形成过程以及植物的抗逆调控中发挥重要作用[19]。因此，推测该代谢通路可能在防御尖孢镰刀菌侵染上扮演着重要角色。
转录因子在植物抗病起着重要的调控作用。逆境胁迫下转录因子可与结构基因上游启动子序列中的顺式作用元件特异性结合，激活或抑制特定抗病基因的表达，从而影响植物对病原体的抗性[20]。本研究发现，在病原菌侵染后，叶片与根部大量转录因子集中在MYB、b HLH、ERF与NAC等转录因子家族上，MYB、NAC、bHLH、ERF、bZIP等转录因子家族在植物抗病过程中发挥重要作用[21 − 23]。MYB转录因子几乎参与植株生理的各个活动过程，在植株抗逆，尤其是抗病中发挥重要的作用[24 − 25]，NAC转录因子也被发现在参与植物抗病[26 − 27]，且在次生壁的合成上发挥着重要的调控作用[28]。研究发现转录因子MYB与NAC通过结合PAL1、C4H、4CL1、C3'H1、HCT、CCoAOMT、CCR1、F5H1与CAD6等关键酶基因启动子中的AC元件，激活其表达并调控木质素合成[29]。本课题组前期研究发现，随着处理时间的延长，病原菌侵染下嫁接冬瓜的木质素含量呈现逐渐上升的趋势[11]，结合本研究推测，可能是由于大量的转录因子MYB与NAC结合了木质素关键酶基因的启动子中的AC原件，从而促进了木质素的持续合成。木质素的合成沉积有利于细胞壁形成更加坚硬与牢固的物理屏障，进一步减少病原菌的侵染扩散[11]。
bHLH基因家族是植物中第二大类转录因子家族，广泛参与植物生长发育、信号传导、代谢及逆境胁迫等过程[30 − 31]。bHLH转录因子在抵御各种细菌性与真菌性病害侵染上扮演重要的角色[32]。南瓜bHLH转录因子CmbHLH87在烟草植株中的异源表达诱导细胞发生了程序性死亡，促进了H2O2积累与病程相关蛋白（PRs）的表达，提高了对南瓜白粉病的抗性[33]。本研究发现，bHLH转录因子数量在叶片与根部都是排名第二，其可能在植株抵御尖孢镰刀菌侵染上发挥重要的作用。
ERF俗称乙烯应答因子（ethylene responsive factor），可以与乙烯诱导顺式作用元件GCC-box结合，在植物抵御植物逆境胁迫中发挥重要的作用[34 − 35]。研究发现，过表达EFR1基因的拟南芥植株对土传尖孢镰刀菌十字花科专化型菌株（Fusarium oxysporum f. sp. Conglutinans）、番茄尖孢镰刀菌番茄专化型菌株（Fusarium oxysporum f. sp. Lycopersici） 引起的枯萎病的抗性增强[36]。本研究发现，叶片与根部的ERF转录因子数量分别是146与315个，可能也在植株抵御病原菌侵染过程中发挥了重要的防御作用。
综上，MYB、b HLH、ERF与NAC等转录因子在嫁接冬瓜叶片与根部抵御病原菌侵染上发挥重要的作用，但这些转录因子是如何参与尖孢镰刀菌侵染下嫁接冬瓜叶片与根部的调节机制与功能尚未清楚，还需进一步深入研究。
植物抗病基因（R基因）是植物在与病原体的长期斗争中逐渐形成的，它们为植物提供了特定的受体和防御机制[37]。这些基因共有一些保守的结构域，根据其特定结构域被分成不同的类别，主要分为CNL、TNL, RLK、RLP、LYK、LYP和LECRK[38]。富含亮氨酸重复序列（LRRs）是大多数R基因蛋白的主要结构域，对识别特异性和抗性至关重要[39]。例如，多聚半乳糖醛酸酶抑制蛋白（polygalacturonase-inhibiting proteins，PGIPs）包含一个LRR结构域，它限制植物病原体的生长并诱导植物的防御反应[40]。有些受体蛋白也存在除LRR外的结构域，如LECRK，其含有一个凝集素样基序（lectin-like motifs，LECM）而不是LRR结构域，水稻Pid2蛋白就属于这一类别，其能诱导水稻对稻瘟病菌的抗性[41]。本研究通过DRAGO 3在叶片的差异表达基因中进一步预测出16个与抗病相关的基因，其中4个具有典型R基因结构域，分别是RLP类：LOC120072345与LOC120085476，LECRK类：LOC120090305与LOC120090267，因此推测嫁接冬瓜可能通过胞外模式识别受体识别病原菌以触发植物的免疫反应，但还需要对预测出的抗病基因进行进一步的分析与功能验证。
Transcriptome analysis of wax gourd grafted with pumpkin to improve resistance to Fusarium wilt
-
摘要: 为解析南瓜（Cucurbita spp.）嫁接提高冬瓜（Benincasa hispida（Thunb.）Cogn.）枯萎病抗性的分子机制，本研究以‘海砧1号’南瓜、‘铁柱2号’冬瓜作为材料，插接法嫁接或自接（自接是砧木与接穗均来自同一种材料）获得嫁接冬瓜（GW）、自接冬瓜（SW）和自接南瓜（SP）3种幼苗，将其种植在8.0×105 CFU·g−1尖孢镰刀菌接菌浓度土壤中，以0 CFU·g−1接菌浓度土壤为对照，培养5天后采集根部和叶片样品进行转录组测序，分析差异表达基因。结果表明，在叶片与根部中分别筛选出284个与2个潜在的抗病基因，对其进行GO与KEGG分析发现，苯丙烷生物合成是嫁接冬瓜叶片与根部共有的代谢通路；利用DRAGO 3对潜在抗病基因进行结构域预测，预测出16个与植物抗病相关的基因，主要与模式触发免疫(pattern-triggered immunity, PTI)途径有关；对叶片与根部中所有转录组序列进行转录因子(Transcription factor，TF)分析，分别鉴定出2 828个和5 426个差异表达TFs，其中共有的大量转录因子集中在MYB、b HLH、ERF与NAC等转录因子家族上。综上可推测，苯丙烷生物合成、MYB、b HLH、ERF与NAC等转录因子在防御尖孢镰刀菌侵染上发挥重要作用。Abstract: To investigate the molecular mechanisms by which pumpkin grafting enhances resistance to Fusarium wilt in wax gourd, pumpkin variety ‘Haizhen 1’ was used as the rootstock and wax gourd variety ‘Tiezhu 2’ as the scion in this experiment. Three types of seedling materials were produced through grafting: grafted wax gourd (GW) to pumpkin, self-grafted wax gourd (SW), and self-grafted pumpkin (SP). These graftings were cultivated in soil inoculated with Fusarium oxysporum at a concentration of 8.0 × 105 CFU·g−1, with a control group using soil at 0 CFU·g−1. After five days of cultivation, root and leaf samples were collected for transcriptome sequencing to analyze differentially expressed genes (DEGs). The results showed that 284 and 2 potential disease-resistant genes were identified in the leaves and roots, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that phenylpropane biosynthesis was a common metabolic pathway in both the leaves and roots of grafted wax gourd. Using DRAGO 3 for domain prediction of the potential disease-resistant genes, 16 disease resistance-related genes were predicted, mainly associated with the PTI (pattern-triggered immunity) pathway. Transcription factor (TF) analysis of all the transcriptome sequences in the leaves and roots showed that there were 2,828 and 5,426 differentially expressed TFs identified in the leaves and roots, respectively. A significant number of these TFs were concentrated in the MYB, bHLH, ERF, and NAC families. It is hence hypothesized that phenylpropane biosynthesis and transcription factors such as MYB, bHLH, ERF, and NAC play crucial roles in defending against Fusarium oxysporum infection.
-
Key words:
- wax gourd /
- pumpkin /
- grafting /
- Fusarium oxysporum /
- transcriptome sequencing
-
图 3 潜在抗病基因的GO功能分类（A：叶片，B：根部）
Fig. 3 GO functional classification of potential disease resistance genes in leaf（A）and root（B）
图 4 潜在抗病基因的KEGG富集分析（A：叶片，B：根部）
Fig. 4 KEGG enrichment analysis of potential disease resistance genes in leaf（A）and root（B）
图 5 苯丙烷生物合成通路图（A：叶片，B：根部）
Fig. 5 Map of phenylpropanoid biosynthesis pathway in leaf（A）and root（B）
图 6 尖孢镰刀菌侵染下响应的TFs（A：叶片，B：根部）
Fig. 6 TFs in response to Fusarium oxysporum infection in leaf（A）and root（B）
图 7 冬瓜叶片4个基因（A−D）南瓜根部2个基因（E−F）荧光定量实验结果
Fig. 7 qRT-PCR verification of 4 genes from wax gourd in leaves（A-D）and in roots（E-F）
表 1 叶片和根部基因qRT-PCR引物
Table 1 qRT-PCR primers for leaf and root genes
部位
Tissue基因
Gene上游引物序列
Forward primer sequence（5'→3'）下游引物序列
Reverse primer sequence（5'→3'）叶片
LeafLOC120076313 CGCTTGTAGCACAGCGAAAG TCGATTATCTTGGGCGGGAC LOC120067870 TGCAGACGGCATAGACAGTG TGCTATTCGGCCGTCTCTTC LOC120086723 TCACTGCTGATCCCCGAAAC GAAACTCAATGCTCGCCACC LOC120075868 GAAACGCATCTGGAAGTGCC ACGTAGGGTCCTGCTCATCA F-box protein CAAACTCCCCAATGCTCAGT CAAACAAGACCACCAGATGAAC 根部
RootLOC111431282 AACACCAAGCCAACTCGGAA AGAGGTTGCCGAAGGAGTTC LOC111434154 AGATGACCCGGAAGCCTTTG GTTGACTGGGTCAGCAAGGT Actin AGCCATCTCTCATCGGTAT CATGGTTGAACCACCACTG 
下载: 导出CSV
表 2 叶片测序数据统计表
Table 2 Statistical table of leaf sequencing data
处理
Treatment生物学重复
Biological repeat高质量测序读段
Clean reads高质量碱基数
Clean bases错误率
Error rate（%）Q20（%）
Q20（%）Q30（%）
Q30（%）GC含量
GC content（%）SWFL 1 43 179 818 6 458 447 311 0.0253 97.30 95.21 46.33 2 42 651 922 6 376 573 009 0.0245 97.70 95.82 45.16 3 42 875 394 6 406 161 037 0.0249 97.50 95.51 45.24 4 45 009 320 6 752 776 974 0.0263 96.83 94.48 45.98 GWFL 1 43 320 074 6 459 906 009 0.0241 97.92 96.16 45.05 2 43 053 652 6 427 466 383 0.0241 97.92 96.15 44.86 3 42 666 922 6 322 273 436 0.0241 97.91 96.17 45.65 4 42 812 064 6 373 731 835 0.0241 97.91 96.18 45.34 SWL 1 42 549 842 6 340 332 373 0.0241 97.91 96.16 45.45 2 43 429 956 6 458 514 075 0.0241 97.91 96.19 45.85 3 47 369 184 7 007 856 484 0.0239 97.98 96.35 45.87 4 43 195 524 6 390 480 048 0.0247 97.58 95.76 46.07 GWL 1 49 023 592 7 288 104 998 0.0243 97.83 96.02 46.81 2 40 544 988 6 041 735 660 0.0244 97.8 95.94 46.68 3 40 588 474 6 060 309 165 0.0243 97.84 96.02 47.10 4 41 032 682 6 117 046 724 0.0244 97.80 95.92 46.78
下载: 导出CSV
表 3 根部测序数据统计表
Table 3 Statistical table of root sequencing data
处理
Treatment生物学重复
Biological repeat高质量测序读段
Clean reads高质量碱基数
Clean bases错误率
Error rate（%）Q20（%）
Q20（%）Q30（%）
Q30（%）GC含量
GC content（%）SPFR 1 55 068 966 8 274 790 336 0.0253 97.95 93.88 46.00 2 41 816 552 6 288 523 427 0.0250 98.05 94.16 45.84 3 56 613 890 8 512 784 710 0.0252 98.00 93.97 45.69 4 53 390 314 8 004 266 521 0.0254 97.91 93.79 45.94 GWFR 1 49 987 502 7 494 885 674 0.0250 98.03 94.1 45.93 2 45 607 438 6 842 619 711 0.0251 97.99 94.01 46.14 3 42 957 918 6 456 667 868 0.0253 97.92 93.84 46.02 4 51 699 286 7 766 345 815 0.0256 97.83 93.57 45.95 SPR 1 47 556 512 7 125 387 830 0.0250 98.03 94.09 46.07 2 48 154 996 7 207 863 128 0.0260 97.64 93.17 45.96 3 53 055 568 7 969 322 838 0.0256 97.83 93.57 45.94 4 52 423 320 7 867 678 851 0.0255 97.86 93.70 45.87 GWR 1 50 296 086 7 558 923 882 0.0256 97.81 93.51 45.69 2 46 241 656 6 952 236 167 0.0257 97.76 93.42 45.92 3 45 286 002 6 803 395 798 0.0256 97.82 93.55 45.98 4 48 664 334 7 309 752 018 0.0254 97.90 93.71 45.93
下载: 导出CSV
表 4 叶片KEGG富集分析详情表
Table 4 Detailed table of leaf KEGG enrichment analysis
代谢通路描述
Pathway description通路编号
Pathway ID基因数量
Gene number基因
Gene基因功能描述
Gene description丙氨酸、天门冬氨酸和谷氨酸代谢
Alanine, aspartate and
glutamate metabolismmap00250 3 LOC120076313 Asparagine synthetase [glutamine-hydrolyzing]1-like LOC120089431 Asparagine synthetase [glutamine-hydrolyzing] 1-like, transcript variant X1 LOC120083554 Probable aldehyde dehydrogenase, transcript variant X1 苯丙烷生物合成
Phenylpropanoid biosynthesismap00940 5 LOC120067870 Peroxidase 2-like LOC120073918 Berberine bridge enzyme-like 15 LOC120081818 Peroxidase 11, transcript variant X1 LOC120085310 4-coumarate--CoA ligase 2-like LOC120091171 Caffeoylshikimate esterase-like, transcript variant X1 角质、亚蜡和蜡生物合成
Cutin, suberine and wax biosynthesismap00073 2 LOC120086723 Omega-hydroxypalmitate O-feruloyl transferase-like LOC120071639 Cytochrome P450 86B1-like 氮代谢
Nitrogen metabolismmap00910 2 LOC120074685 Uncharacterized LOC120074685 LOC120075868 Beta carbonic anhydrase 5, chloroplastic-like, transcript variant X1
下载: 导出CSV
表 5 根部KEGG富集分析详情表
Table 5 Detailed table of root KEGG enrichment analysis
代谢通路描述
Pathway description通路编号
Pathway ID基因数量
Gene number基因
Gene基因功能描述
Gene description芪类、二芳基庚烷类和姜酚生物合成 Stilbenoid, diarylheptan-oid and gingerol biosynthesis map00945 1 LOC111431282 Vinorine synthase-like 光合作用天线蛋白Photosynthesis-antenna proteins map00196 1 LOC111434154 Chlorophyll a-b binding protein of LHCII type 1-like 类黄酮生物合成 Flavonoid biosynthesis map00941 1 LOC111431282 Vinorine synthase-like 苯丙烷生物合成Phenylpropanoid biosynt-hesis map00940 1 LOC111431282 Vinorine synthase-like
下载: 导出CSV
表 6 叶片潜在抗病基因蛋白结构域预测
Table 6 Prediction of protein domain of potential disease resistance genes in leaves
基因 Gene 分类 Class 结构域 Domain 基因功能描述 Description LOC120087570 CK CC,TM,Kinase Contains a coiled-coil domain（CC）,a transmembrane
domain（TM）and a kinase domainLOC120090355 CL CC,LRR,TM Contains a coiled-coil domain,a leucine-rich repeat
domain（LRR）and a kinase domain（Kinase）LOC120076759 KIN TM,Kinase Contains a transmembrane domain and a kinase domain LOC120083806 KIN TM,Kinase LOC120082723 KIN TM,Kinase LOC120088593 KIN TM,Kinase LOC120089817 KIN TM,Kinase LOC120092573 KIN TM,Kinase LOC120072525 KIN TM,Kinase LOC120090305 LECRK LECM,TM,Kinase Lectin-receptor kinase, contains a lectin-like motifs（LECM）
instead of LRR,a transmembrane domain and a kinase domainLOC120090267 LECRK LECM,TM,Kinase LOC120090948 N NBS,TM Contains a nucleotide-binding site domain（NBS）domain and a
transmembrane domainLOC120085520 N NBS,TM LOC120072345 RLP LRR,TM Receptor like proteins, contains a leucine-rich repeat domain and
a kinase domainLOC120085476 RLP LRR,TM LOC120085949 TRAM TM Contains a transmembrane domain only
下载: 导出CSV
-
[1] YANG S G, CAI J S, WANG M, et al. The construction and analysis of wax gourd pangenome uncover fruit quality-related and resistance genes[J]. Scientia Horticulturae, 2023, 318: 112084. doi: 10.1016/j.scienta.2023.112084 [2] 谢大森, 江彪, 刘文睿, 等. 优质、抗病冬瓜多样化育种研究进展[J]. 广东农业科学, 2020, 47(11): 50 − 59. doi: 10.16768/j.issn.1004-874X.2020.11.006 [3] 何晓明, 彭庆务, 王敏, 等. 我国节瓜遗传育种研究进展[J]. 广东农业科学, 2021, 48(9): 1 − 11. doi: 10.16768/j.issn.1004-874X.2021.09.001 [4] 石博, 关峰, 张景云, 等. 冬瓜枯萎病研究进展[J]. 中国瓜菜, 2022, 35(11): 11 − 14. doi: 10.3969/j.issn.1673-2871.2022.11.002 [5] DITA M, BARQUERO M, HECK D, et al. Fusarium wilt of banana: current knowledge on epidemiology and research needs toward sustainable disease management[J]. Frontiers in Plant Science, 2018, 9: 1468. doi: 10.3389/fpls.2018.01468 [6] 袁飞, 刘子凡, 廖道龙, 等. 南瓜嫁接提高冬瓜枯萎病抗性的化感机制[J]. 中国瓜菜, 2021, 34(5): 26 − 29. doi: 10.3969/j.issn.1673-2871.2021.05.005 [7] 朱白婢, 胡艳平, 云天海, 等. 冬瓜砧木新品种“海砧1号”的选育[J]. 北方园艺, 2018(10): 208 − 210. doi: 10.11937/bfyy.20180051 [8] 廖道龙, 冯学杰, 朱白俾, 等. 冬瓜嫁接砧木评价及其产量与性状的相关性分析[J]. 南方农业学报, 2017, 48(8): 1434 − 1440. doi: 10.3969/j.issn.2095-1191.2017.08.16 [9] 黄项心兰, 廖道龙, 刘子凡, 等. 不同砧木材料嫁接冬瓜苗期氮营养特性综合评价[J]. 热带作物学报, 2022, 43(1): 137 − 144. doi: 10.3969/j.issn.1000-2561.2022.01.018 [10] 符厚隆, 廖道龙, 刘子凡, 等. 南瓜作砧木嫁接冬瓜根系分泌物对土壤微生态的影响[J]. 中国瓜菜, 2022, 35(6): 56 − 61. doi: 10.3969/j.issn.1673-2871.2022.06.009 [11] FU H L, FU J Y, ZHOU B, et al. Biochemical mechanisms preventing wilting under grafting: a case study on pumpkin rootstock grafting to wax gourd[J]. Frontiers in Plant Science, 2024, 15: 1331698. doi: 10.3389/fpls.2024.1331698 [12] SABRY S, ALI A Z, ABDEL-KADER D A, et al. Histopathological and biochemical aspects of grafted and non-grafted cucumber infected with stem rot caused by Fusarium spp.[J]. Saudi Journal of Biological Sciences, 2022, 29(3): 1770 − 1780. doi: 10.1016/j.sjbs.2021.10.053 [13] ATTIA E Z, FAROUK H M, ABDELMOHSEN U R, et al. Antimicrobial and extracellular oxidative enzyme activities of endophytic fungi isolated from alfalfa (Medicago sativa) assisted by metabolic profiling[J]. South African Journal of Botany, 2020, 134: 156 − 162. doi: 10.1016/j.sajb.2019.12.003 [14] WESTPHAL K R, HEIDELBACH S, ZEUNER E J, et al. The effects of different potato dextrose agar media on secondary metabolite production in Fusarium[J]. International Journal of Food Microbiology, 2021, 347: 109171. doi: 10.1016/j.ijfoodmicro.2021.109171 [15] WANG B B, YUAN J, ZHANG J, et al. Effects of novel bioorganic fertilizer produced by Bacillus amyloliquefaciens W19 on antagonism of Fusarium wilt of banana[J]. Biology and Fertility of Soils, 2013, 49(4): 435 − 446. doi: 10.1007/s00374-012-0739-5 [16] 张金平, 何晓明, 江彪, 等. 镰刀菌酸胁迫下节瓜实时荧光定量PCR分析内参基因的选择[J]. 广东农业科学, 2014, 41(10): 120 − 125. doi: 10.3969/j.issn.1004-874X.2014.10.029 [17] 张浩, 李海玉, 施松梅, 等. 2种黑籽南瓜响应枯萎病菌侵染的转录组学研究[J]. 中国农学通报, 2023, 39(18): 107 − 116. doi: 10.11924/j.issn.1000-6850.casb2023-0088 [18] 韩彦卿. 水稻抗稻曲病QTL定位及抗/感品种响应稻曲病菌侵染的表达谱分析[D]. 北京: 中国农业大学, 2015. [19] ZHANG C C, YAO X H, REN H D, et al. Genome-wide identification and characterization of the phenylalanine ammonia-lyase gene family in pecan (Carya illinoinensis)[J]. Scientia Horticulturae, 2022, 295: 110800. doi: 10.1016/j.scienta.2021.110800 [20] 谢迎秋, 孟蒙, 朱祯. 植物反式作用因子研究进展[J]. 高技术通讯, 2000, 10(2): 97 − 102. doi: 10.3321/j.issn:1002-0470.2000.02.024 [21] 马达. 拟南芥bHLH转录因子HFR1调控抗病性与发育的机制研究[D]. 长沙: 湖南农业大学, 2014. [22] 赵明奇, 刘晓洁, 梁玉青, 等. 新疆野苹果AP2/ERF转录因子家族鉴定与响应腐烂病的表达分析[J]. 西北植物学报, 2022, 42(6): 930 − 942. doi: 10.7606/j.issn.1000-4025.2022.06.0930 [23] 李铃仙. 小麦抗病基因的发现与应用[J]. 农业灾害研究, 2023, 13(10): 49 − 51. doi: 10.3969/j.issn.2095-3305.2023.10.017 [24] ZHANG Z J, CHEN J M, SU Y Y, et al. TaLHY, a 1R-MYB transcription factor, plays an important role in disease resistance against stripe rust fungus and ear heading in wheat[J]. PLoS One, 2015, 10(5): e0127723. doi: 10.1371/journal.pone.0127723 [25] NGUYEN N H, LEE H. MYB-related transcription factors function as regulators of the circadian clock and anthocyanin biosynthesis in Arabidopsis[J]. Plant Signaling & Behavior, 2016, 11(3): e1139278. doi: 10.1080/15592324.2016.1139278 [26] LIU Q, YAN S J, HUANG W J, et al. NAC transcription factor ONAC066 positively regulates disease resistance by suppressing the ABA signaling pathway in rice[J]. Plant Molecular Biology, 2018, 98(4): 289 − 302. doi: 10.1007/s11103-018-0768-z [27] PEROCHON A, KAHLA A, VRANIĆ M, et al. A wheat NAC interacts with an orphan protein and enhances resistance to Fusarium head blight disease[J]. Plant Biotechnology Journal, 2019, 17(10): 1892 − 1904. doi: 10.1111/pbi.13105 [28] ZHONG R Q, MCCARTHY R L, LEE C, et al. Dissection of the transcriptional program regulating secondary wall biosynthesis during wood formation in poplar[J]. Plant Physiology, 2011, 157(3): 1452 − 1468. doi: 10.1104/pp.111.181354 [29] VANHOLME R, DE MEESTER B, RALPH J, et al. Lignin biosynthesis and its integration into metabolism[J]. Current Opinion in Biotechnology, 2019, 56: 230 − 239. doi: 10.1016/j.copbio.2019.02.018 [30] 李陈莹, 王冉, 梁岩. 维管束木质化调控植物抗青枯病的研究进展[J]. 浙江大学学报(农业与生命科学版), 2023, 49(5): 633 − 643. doi: 10.3785/j.issn.1008-9209.2023.06.171 [31] WANG Y, JIANG C J, LI Y Y, et al. CsICE1 and CsCBF1: two transcription factors involved in cold responses in Camellia sinensis[J]. Plant Cell Reports, 2012, 31(1): 27 − 34. doi: 10.1007/s00299-011-1136-5 [32] XU X B, TIAN S P. Salicylic acid alleviated pathogen-induced oxidative stress in harvested sweet cherry fruit[J]. Postharvest Biology and Technology, 2008, 49(3): 379 − 385. doi: 10.1016/j.postharvbio.2008.02.003 [33] 刘茜, 王爱云, 荣玮. ERF转录因子在作物抗病基因工程中的研究进展[J]. 种子, 2014, 33(1): 48 − 53. doi: 10.3969/j.issn.1001-4705.2014.01.013 [34] VALENZUELA-RIFFO F, ZÚÑIGA P E, MORALES-QUINTANA L, et al. Priming of defense systems and upregulation of MYC2 and JAZ1 genes after Botrytis cinerea inoculation in methyl jasmonate-treated strawberry fruits[J]. Plants, 2020, 9(4): 447. doi: 10.3390/plants9040447 [35] 孙钧政, 李美玲, 唐金艳, 等. 水杨酸诱导果实采后抗病性机制研究进展[J]. 亚热带植物科学, 2021, 50(5): 413 − 420. doi: 10.3969/j.issn.1009-7791.2021.05.013 [36] 温晓丽, 罗维巍, 范翌婷, 等. 茉莉酸甲酯诱导的采后果蔬抗病性及其机制研究进展[J]. 食品安全质量检测学报, 2023, 14(12): 166 − 172. doi: 10.19812/j.cnki.jfsq11-5956/ts.2023.12.032 [37] ANDOLFO G, ERCOLANO M R. Plant innate immunity multicomponent model[J]. Frontiers in Plant Science, 2015, 6: 987. doi: 10.3389/fpls.2015.00987 [38] CALLE GARCÍA J, GUADAGNO A, PAYTUVI-GALLART A, et al. PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process[J]. Nucleic Acids Research, 2022, 50(D1): D1483 − D1490. doi: 10.1093/nar/gkab1087 [39] SHAN Y F, LI M Y, WANG R Z, et al. Evaluation of the early defoliation trait and identification of resistance genes through a comprehensive transcriptome analysis in pears[J]. Journal of Integrative Agriculture, 2023, 22(1): 120 − 138. doi: 10.1016/j.jia.2022.08.040 [40] DI C X, ZHANG M X, XU S J, et al. Role of poly-galacturonase inhibiting protein in plant defense[J]. Critical Reviews in Microbiology, 2006, 32(2): 91 − 100. doi: 10.1080/10408410600709834 [41] WANG J, QU B Y, DOU S J, et al. The E3 ligase OsPUB15 interacts with the receptor-like kinase PID2 and regulates plant cell death and innate immunity[J]. BMC Plant Biology, 2015, 15: 49. doi: 10.1186/s12870-015-0442-4 -
点击查看大图
计量
- 文章访问数: 13
- HTML全文浏览量: 8
- 被引次数: 0
