Research of Genetic Regulation of Bombax ceiba Prickle Based on Hormone，Transcriptome and Metabolome

doi:10.16420/j.issn.0513-353x.2024-0376

Acta Horticulturae Sinica ›› 2025, Vol. 52 ›› Issue (4): 933-946.doi: 10.16420/j.issn.0513-353x.2024-0376

• Genetic & Breeding·Germplasm Resources·Molecular Biology • Previous Articles Next Articles

Research of Genetic Regulation of Bombax ceiba Prickle Based on Hormone，Transcriptome and Metabolome

SHAO Yifan¹, ZHU Baoqing¹, WANG Tongxin¹, LIAO Jianhe², WU Fanhua³, YANG Siyi¹, FENG Jianhang¹, YU Xudong¹^,^*()

¹ Key Laboratory of Genetics and Germplasm Innovation of Tropical Forest Trees and Ornamental Plants（Ministry of Education）;Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province，School of Tropical Agriculture and Forestry, Hainan University，Danzhou，Hainan 571737，China
² School of Materials Science and Engineering，Hainan University，Haikou 570228，China
³ College of life sciences，Hainan University，Haikou 570228，China

Received:2025-01-08 Revised:2025-02-20 Online:2025-05-08 Published:2025-04-25
Contact: YU Xudong

Abstract

Abstract:

Key regulatory genes that create prickles in Bombax ceiba were explored via phytohormone assays，transcriptomics，and metabolomics by use of the stem bark from the upper less prickly zone and lower prickly zone. The phytohormone assay revealed that abscisic acid，jasmonic acid，and salicylic acid were higher in the less prickly zone than in the prickly zone，whereas gibberellin was higher in the more prickly zone than in the less prickly zone. Transcriptome sequencing revealed that in the growth hormone biosynthesis pathway，the expression of Aux/IAA protein genes and growth hormone response factor ARF genes was higher in the prickly zone than in the less prickly zone. In the cytokinin biosynthesis pathway，the expression of B-ARR genes was higher in the prickly zone than in the less prickly zone whereas，in the gibberellin biosynthesis pathway，the expression of the GID1，DELLA，and TF genes was higher in the less prickly zone than in the prickly area. Furthermore，the expression of the JAR1，JAZ，and MYC2 genes was higher in the less prickly zone than in the prickly zone in the jasmonic acid biosynthesis pathway. Overall，13 key genes（BRI1，BZR1/2，AUX1，AUX/IAA，ARF，DELLA，JAR1，TGA，PR-1，SIMKK，MYB62，MYB123，and WRKY28） that affected the development of cottonwood prickles were verified by qRT-PCR，which were similar to the results of high-throughput sequencing；thus，supporting the reliability of the RNA-Seq results. Metabolomics screened 2 187 metabolites and 92 differential metabolites and identified that the CCR and COMT genes in the phenylpropanoid biosynthesis pathway as well as FLS genes in the flavonoid biosynthesis pathway, were related to the formation of cottonwood prickles. In addition，the metabolites（2S）-2-isopropylmalic acid and （2R,3S）-3-isopropyl malate in the valine，leucine，and isoleucine metabolic pathways and quercetin in the flavonoid biosynthesis pathway are involved in cottonwood prickle formation.

Key words: Bombax ceiba, prickle, plant hormones, transcriptome, metabolome

SHAO Yifan, ZHU Baoqing, WANG Tongxin, LIAO Jianhe, WU Fanhua, YANG Siyi, FENG Jianhang, YU Xudong. Research of Genetic Regulation of Bombax ceiba Prickle Based on Hormone，Transcriptome and Metabolome[J]. Acta Horticulturae Sinica, 2025, 52(4): 933-946.

Figures/Tables 11

References 26

[1]	Boter M, Ruíz-Rivero O, Abdeen A, Prat S. 2004. Conserved MYC transcription factors play a key role in jasmonate signaling both in tomato and Arabidopsis. Genes & Development, 18 (13):1577-1591.
[2]	Fambrini M, Pugliesi C. 2019. The dynamic genetic-hormonal regulatory network controlling the trichome development in leaves. Plants-basel, 8 (8):253.
[3]	Feng Fayu, Wang Yi, Lu Bin, He Jiani, Yuan XiaoLong. 2022. Identification and expression analysis of Zanthoxylum armatum WRKY gene family based on the transcripome data. Molecular Plant Breeding, 20 (18):5984-5994. (in Chinese)
	冯发玉, 王毅, 陆斌, 何家妮, 原晓龙. 2022. 基于转录组数据竹叶花椒WRKY基因家族的鉴定及表达分析. 分子植物育种, 20 (18):5984-5994.
[4]	Feng Fayu, Wang Yi, Wang Lijuan, Romania, Xu Lingwen, Yang Yanping, Lu Bin. 2021. Identification and expression analysis of MYB gene family in Zanthoxylum armatum. Economic forest research, 39 (4):148-157. (in Chinese)
	冯发玉, 王毅, 王丽娟, 罗玛妮娅, 徐令文, 杨彦萍, 陆斌. 2021. 竹叶花椒MYB基因家族的鉴定及其表达特性的分析. 经济林研究, 39 (4):148-157.
[5]	Gan Y, Kumimoto R, Liu C, Ratcliffe O, Yu H, Broun P. 2006. GLABROUS INFLORESCENCE STEMS modulates the regulation by gibberellins of epidermal differentiation and shoot maturation in Arabidopsis. Plant Cell, 18 (6):1383-1395.
[6]	Guan Shuwen. 2022. Study on the molecular mechanism of prickly development of Zanthoxylum armatum[M. D. Dissertation]. Kunming: Southwest Forestry University. (in Chinese)
	关淑文. 2022. 竹叶花椒的皮刺发育分子机理研究[硕士论文]. 昆明: 西南林业大学.
[7]	Hua B, Chang J, Xu Z, Han X, Xu M, Yang M, Yang C, Ye Z, Wu S. 2021. HOMEODOMAIN PROTEIN8 mediates jasmonate-triggered trichome elongation in tomato. New Phytologist, 230 (3):1063-1077.
[8]	Huang H, Liu B, Liu L, Song S. 2017. Jasmonate action in plant growth and development. Journal of Experimental Botany, 68 (6):1349-1359.
[9]	Li J, Wang X, Jiang R, Dong B, Fang S, Li Q, Lv Z, Chen W. 2021. Phytohormone-based regulation of trichome development. Frontiers in Plant Science,12:734776.
[10]	Li Xuhui, Chen Weiwei, Qi Yongwen, Lu Siqi. 2023. Application of maize ZmWRKY28 gene in regulating root hair length. CN116640773A,1-11. (in Chinese)
	李旭辉, 陈卫卫, 齐永文, 陆思奇. 2023. 玉米ZmWRKY28基因在调控根毛长度中的应用. CN116640773A,1-11.
[11]	Liu Jinqiu, Chen Kai, Zhang Zhenzhu, Chen Xiuling, Wang Aoxue. 2016. Effects of exogenous GA,MeJA,IAA,SA and KT on trichome formation in tomato. Acta Horticulturae Sinica, 43 (11):2151-2160. (in Chinese)
	刘金秋, 陈凯, 张珍珠, 陈秀玲, 王傲雪. 2016. 外施GA、MeJA、IAA、SA和KT对番茄表皮毛发生的作用. 园艺学报, 43 (11):2151-2160.
[12]	Liu X, He X, Liu Z, Wu P, Tang N, Chen Z, Zhang W, Rao S, Cheng S, Luo C, Xu F. 2022. Transcriptome mining of genes in Zanthoxylum armatum revealed ZaMYB 86 as a negative regulator of prickly development. Genomics, 114 (3):110374.
[13]	Pandey S, Goel R, Bhardwaj A, Asif M, H Sawant, S V, Misra P. 2018. Transcriptome analysis provides insight into prickle development and its link to defense and secondary metabolism in Solanum viarum dunal. Scientific Reports, 11 (1):16098.
[14]	Qi T, Huang H, Wu D, Yan J, Qi Y, Song S, Xie D. 2014. Arabidopsis DELLA and JAZ proteins bind the WD-repeat/bHLH/MYB complex to modulate gibberellin and jasmonate signaling synergy. Plant Cell, 26 (3):1118-1133.
[15]	Ren Z, Li Z, Miao Q, Yang Y, Deng W, Hao Y. 2011. The auxin receptor homologue in Solanum lycopersicum stimulates tomato fruit set and leaf morphogenesis. Journal of Experimental Botany, 62 (8):2815-2826.
[16]	Tang N, Cao Z, Wu P, Liu Y, Lou J, Hu Y, Sun X, Si S, Chen Z. 2023. Comparative transcriptome analysis reveals hormone,transcriptional and epigenetic regulation involved in prickle formation in Zanthoxylum armatum. Gene,871:147434.
[17]	Traw M B, Bergelson J. 2003. Interactive effects of jasmonic acid,salicylic acid,and gibberellin on induction of trichomes in Arabidopsis. Plant Physiology, 133 (3):1367-1375.
[18]	Wang T, Jia Q, Wang W, Hussain S, Ahmed S, Adnan Zhou, D X, Ni Z, Wang S. 2019. GCN5 modulates trichome initiation in Arabidopsis by manipulating histone acetylation of core trichome initiation regulator genes. Plant Cell Reports, 38 (6):755-765.
[19]	Werner S, Bartrina I, Schmülling T. 2021. Cytokinin regulates vegetative phase change in Arabidopsis thaliana through the miR172/TOE1-TOE 2 module. Nature Communications, 12 (1):5816.
[20]	Wu Jinshan, Yu Xudong, Xu Zekang, Wang Yachen, Gao Haojie, Zhao Ying. 2017. Photosynthetic characteristics and carbon sink capacity of Hainan Bombax ceiba in different months. Tropical Agriculture in China,(3):67-71. (in Chinese)
	吴金山, 于旭东, 许泽康, 王亚沉, 高豪杰, 赵莹. 2017. 海南木棉不同月份光合特性及碳汇能力研究. 中国热带农业,(3):67-71.
[21]	Wu Yansheng, Yan Wenpei, Zhang Xiaoli, Niu Weitao, Zhang Pengfei, Wu Lina, Li Min, Xing Cuijuan, Wang Senghu. 2021. Discrimination on the types and functions of plant thorns. Bulletin of Biology, 56 (4):13-15. (in Chinese)
	武延生, 严文培, 张晓丽, 牛伟涛, 张鹏飞, 武丽娜, 李敏, 邢翠娟, 王僧虎. 2021. 植物刺的类型与功能辨析. 生物学通报, 56 (4):13-15.
[22]	Xia X, Hu Q, Li W, Chen Y, Han L, Miao T, Wu W, Li X, Huang G. 2018. Cotton(Gossypium hirsutum) JAZ3 and SLR1 function in jasmonate and gibberellin mediated epidermal cell differentiation and elongation. Plant Cell Tissue and Organ Culture, 133 (2):1-14.
[23]	Xia Xiaocong. 2017. Study on the function of Gossypium hirsutum GhJAZ3 and GhGAI 1 in JA and GA mediated trichome differentiation and elongation of transgenic Arabidopsis thaliana[M. D. Dissertation]. Wuhan: Central China Normal University. (in Chinese)
	夏小丛. 2017. 棉花(Gossypium hirsutum)GhJAZ3和GhGAI1在JA和GA介导的转基因拟南芥表皮毛分化和伸长中的功能研究[硕士论文]. 武汉: 华中师范大学.
[24]	Yao Kexin. 2023. Functional study of the BsMYB62 transcription factor in Begonia semperflorens[M. D. Dissertation]. Zhengzhou: Henan Agricultural University. (in Chinese)
	姚珂心. 2023. 四季秋海棠BsMYB62转录因子的功能研究[硕士论文]. 郑州: 河南农业大学.
[25]	Yu N, Cai W J, Wang S C, Shan C M, Wang L J, Chen X Y. 2010. Temporal control of trichome distribution by microRNA156-targeted SPL genes in Arabidopsis thaliana. Plant Cell, 22 (7):2322-2335.
[26]	Zhuang Yuan. 2015. Cloning and functional analysis of chalcone synthase gene and Cinnamoyl-CoA reductase gene in Pinus koraiensis[M. D. Dissertation]. Changchun: Jilin University. (in Chinese)
	庄园. 2015. 红松查尔酮合酶基因和肉桂酰辅酶A还原酶基因的克隆及功能分析[硕士论文]. 长春: 吉林大学.

引物 Primer	正引物序列（5′-3′） Forward primer	反引物序列（5′-3′） Reverse primer	基因ID Gene ID
BRI1	GGAGGTCATAAACAGCCC	TCCTTCAAGTCCCCAATC	evm.model.Scaffold196.5
BZR1/2	ACCTCTATCTTCCCCTAC	CATCAGACTCATCGCATT	evm.model.Scaffold21.213
AUX1	ACCCGAGAAGAGGGAAAA	AGAAGGAATAGGGCAGCG	evm.model.Scaffold500.291
AUX/IAA	GCGAACTCACCAAAAACA	AATAGAGCGAGCAGCAAC	evm.model.Scaffold262.58
ARF	TGTGGAATCAGACGAAGT	TCCCGATAGAGTAAATGG	evm.model.Scaffold181.29
DELLA	ATGGAGGTAAAAGGAACG	CATCAGCAAAATAGTGGG	evm.model.Scaffold132.216
JAR1	GCTCTGGCACTACTCAA	CAGCAAACAAACCACCTT	evm.model.Scaffold221.79
TGA	AGATTGGTATGTGGGGAG	TCGTCAGCATCACTTCCT	evm.model.Scaffold127.460
PR-1	CCCAATCTATTCCACCTC	ATCGGGCACATTTCTCAA	evm.model.Scaffold2351.33
SIMKK	TCCCAGAAACGGAGAGAC	AGTAGAAGGAGTAGGCGG	evm.model.Scaffold66.339
MYB62	GGGAGAACGGAAAATGAG	TGGATGAACAAGAGGATG	evm.model.Scaffold144.195
MYB123	CGTTACAGCCACCACCAA	CATCAAGCACTGCGAATC	evm.model.Scaffold498.25
WRKY28	TCAGTTTCTTCCTCCTCT	TCCTCCACCTGTATCCAT	evm.model.Scaffold46.3

引物 Primer	正引物序列（5′-3′） Forward primer	反引物序列（5′-3′） Reverse primer	基因ID Gene ID
BRI1	GGAGGTCATAAACAGCCC	TCCTTCAAGTCCCCAATC	evm.model.Scaffold196.5
BZR1/2	ACCTCTATCTTCCCCTAC	CATCAGACTCATCGCATT	evm.model.Scaffold21.213
AUX1	ACCCGAGAAGAGGGAAAA	AGAAGGAATAGGGCAGCG	evm.model.Scaffold500.291
AUX/IAA	GCGAACTCACCAAAAACA	AATAGAGCGAGCAGCAAC	evm.model.Scaffold262.58
ARF	TGTGGAATCAGACGAAGT	TCCCGATAGAGTAAATGG	evm.model.Scaffold181.29
DELLA	ATGGAGGTAAAAGGAACG	CATCAGCAAAATAGTGGG	evm.model.Scaffold132.216
JAR1	GCTCTGGCACTACTCAA	CAGCAAACAAACCACCTT	evm.model.Scaffold221.79
TGA	AGATTGGTATGTGGGGAG	TCGTCAGCATCACTTCCT	evm.model.Scaffold127.460
PR-1	CCCAATCTATTCCACCTC	ATCGGGCACATTTCTCAA	evm.model.Scaffold2351.33
SIMKK	TCCCAGAAACGGAGAGAC	AGTAGAAGGAGTAGGCGG	evm.model.Scaffold66.339
MYB62	GGGAGAACGGAAAATGAG	TGGATGAACAAGAGGATG	evm.model.Scaffold144.195
MYB123	CGTTACAGCCACCACCAA	CATCAAGCACTGCGAATC	evm.model.Scaffold498.25
WRKY28	TCAGTTTCTTCCTCCTCT	TCCTCCACCTGTATCCAT	evm.model.Scaffold46.3

样品 Sample	总读取数Total Reads	比对数 Reads mapped	唯一比对数 Unique mapped	多重比对数 Multi mapped	Read1比对数 Read1 mapped	Read2比对数 Read2 mapped	正链比对数 ‘+’mapped	负链比对数 ‘-’mapped
LP-1	43 929 340	29 341 819（66.79%）	28 260 900（64.33%）	10 809 19（2.46%）	14 203 787（32.33%）	14 057 113（32.00%）	14 115 917（32.13%）	14 144 983（32.20%）
LP-2	47 777 116	32 339 797（67.69%）	31 122 901（65.14%）	12 168 96（2.55%）	15 637 980（32.73%）	15 484 921（32.41%）	15 545 737（32.54%）	15 577 164（32.60%）
LP-3	46 733 462	30 526 628（65.32%）	29 400 369（62.91%）	11 262 59（2.41%）	14 773 962（31.61%）	14 626 407（31.30%）	14 684 755（31.42%）	14 715 614（31.49%）
P-1	45 162 140	28 053 596（62.12%）	26 901 920（59.57%）	11 516 76（2.55%）	13 517 914（29.93%）	13 384 006（29.64%）	13 440 797（29.76%）	13 461 123（29.81%）
P-2	44 355 554	28 766 857（64.86%）	27 560 164（62.13%）	12 066 93（2.72%）	13 849 789（31.22%）	13 710 375（30.91%）	13 770 449（31.05%）	13 789 715（31.09%）
P-3	43 934 068	27 565 859（62.74%）	26 383 331（60.05%）	11 825 28（2.69%）	13 135 656（29.90%）	13 247 675（30.15%）	13 184 883（30.01%）	13 198 448（30.04%）

样品 Sample	总读取数Total Reads	比对数 Reads mapped	唯一比对数 Unique mapped	多重比对数 Multi mapped	Read1比对数 Read1 mapped	Read2比对数 Read2 mapped	正链比对数 ‘+’mapped	负链比对数 ‘-’mapped
LP-1	43 929 340	29 341 819（66.79%）	28 260 900（64.33%）	10 809 19（2.46%）	14 203 787（32.33%）	14 057 113（32.00%）	14 115 917（32.13%）	14 144 983（32.20%）
LP-2	47 777 116	32 339 797（67.69%）	31 122 901（65.14%）	12 168 96（2.55%）	15 637 980（32.73%）	15 484 921（32.41%）	15 545 737（32.54%）	15 577 164（32.60%）
LP-3	46 733 462	30 526 628（65.32%）	29 400 369（62.91%）	11 262 59（2.41%）	14 773 962（31.61%）	14 626 407（31.30%）	14 684 755（31.42%）	14 715 614（31.49%）
P-1	45 162 140	28 053 596（62.12%）	26 901 920（59.57%）	11 516 76（2.55%）	13 517 914（29.93%）	13 384 006（29.64%）	13 440 797（29.76%）	13 461 123（29.81%）
P-2	44 355 554	28 766 857（64.86%）	27 560 164（62.13%）	12 066 93（2.72%）	13 849 789（31.22%）	13 710 375（30.91%）	13 770 449（31.05%）	13 789 715（31.09%）
P-3	43 934 068	27 565 859（62.74%）	26 383 331（60.05%）	11 825 28（2.69%）	13 135 656（29.90%）	13 247 675（30.15%）	13 184 883（30.01%）	13 198 448（30.04%）

Research of Genetic Regulation of Bombax ceiba Prickle Based on Hormone，Transcriptome and Metabolome

RichHTML

PDF

Knowledge

Abstract

Cite this article

share this article

Figures/Tables 11

References 26

Related Articles 15

Recommended Articles

Metrics

Comments

[1]	YUAN Juanwei, JIA Li, WANG Han, YAN Congsheng, ZHANG Qi’an, YU Feifei, GAN Defang, JIANG Haikun. Screening of Resistant Germplasm and Identifying Candidate Gene for Phytophthora capsici Resistance in Pepper Germplasm Resources [J]. Acta Horticulturae Sinica, 2025, 52(4): 857-871.
[2]	TANG Lingli, HE Yuhua, WANG Qingtao, AN Lulu, XU Yongyang, ZHANG Jian, KONG Weihu, HU Keyun, ZHAO Guangwei. Hormonal and Transcriptome Analysis of Ripe Fruits of Non-Climacteric and Climacteric Melon [J]. Acta Horticulturae Sinica, 2025, 52(4): 883-896.
[3]	LIN Ruxue, WANG Fei, ZHANG Yanjie, MA Li, LIU Xiaofeng, LI Shuran, LIU Yalong, WANG Chen, JIANG Shuling, OU Chunqing. Transcriptome Screening and Functional Analysis of Dwarf Related Genes in Pear [J]. Acta Horticulturae Sinica, 2025, 52(3): 545-560.
[4]	KUANG Meimei, LI Li, MA Jianwei, LIU Yuan, JIANG Hongfei, LEI Rui, MAN Yuping, WANG Yifan, HUANG Bo, WANG Yanchang, LIU Shibiao. Screening of Kiwifruit Canker Resistance-Related Genes Based on the Transcriptome Sequencing Analysis of Hybrids from Actinidia chinensis [J]. Acta Horticulturae Sinica, 2024, 51(8): 1743-1757.
[5]	YAO Fengping, WANG Yanbin, QIN Yuchuan, WANG Liling. Analysis of the Difference of Flavonoids in Dendropanax dentiger Leaves with Different Drying Processes Based on Untargeted Metabolomics [J]. Acta Horticulturae Sinica, 2024, 51(5): 1069-1082.
[6]	ZHANG Wenhao, ZHANG Hui, LIU Yuting, WANG Yan, ZHANG Yingying, WANG Xinman, WANG Quanhua, ZHU Weimin, YANG Xuedong. Preliminary Transcriptome Analysis in Two Tomato Fruits Materials with Different Sugar Content [J]. Acta Horticulturae Sinica, 2024, 51(2): 281-294.
[7]	WANG Hong, MIAO Chen, DING Xiaotao, SHEN Haibin, ZHU Weimin. Growth，Development and Transcriptome Analysis of Tomato Cultivars with Different Low Light Tolerances Under Low Light Condition [J]. Acta Horticulturae Sinica, 2024, 51(2): 335-345.
[8]	BU Wenxuan, YAO Yiping, HUANG Yu, YANG Xingyu, LUO Xiaoning, ZHANG Minhuan, LEI Weiqun, WANG Zheng, TIAN Jianing, CHEN Lujie, QIN Liping. Transcriptome Analysis of the Response of Tree Peony‘Hu Hong’Under High Temperature Stress [J]. Acta Horticulturae Sinica, 2024, 51(12): 2800-2816.
[9]	LU Yanqing, LIN Yanjin, LU Xinkun. Cell Wall Material Metabolism and the Response to High Temperature and Water Deficiency Stresses in Pericarp are Associated with Fruit Cracking of‘Duwei’Pummelo [J]. Acta Horticulturae Sinica, 2023, 50(8): 1747-1768.
[10]	RUAN Ruoxin, LUO Huifeng, ZHANG Chen, HUANG Kangkang, XI Dujun, PEI Jiabo, XING Mengyun, LIU Hui. Transcriptome Analysis of Flower Buds of Sweet Cherry Cultivars with Different Chilling Requirements During Dormancy Stages in Hangzhou [J]. Acta Horticulturae Sinica, 2023, 50(6): 1187-1202.
[11]	XU Yue, LI Zixiong, CHEN Jie, SUN Liang. Transcriptional Bases of the Regulation of 2,4-D on Tomato Fruit Shape [J]. Acta Horticulturae Sinica, 2023, 50(4): 802-814.
[12]	WANG Zehan, YU Wentao, WANG Pengjie, LIU Caiguo, FAN Xiaojing, GU Mengya, CAI Chunping, WANG Pan, YE Naixing. WGCNA Analysis of Differentially Expressed Genes in Floral Organs of Tea Germplasms with Ovary-glabrous and Ovary-trichome [J]. Acta Horticulturae Sinica, 2023, 50(3): 620-634.
[13]	JIANG Yu, TU Xunliang, HE Junrong. Analysis of Differential Expression Genes in Leaves of Leaf Color Mutant of Chinese Orchid [J]. Acta Horticulturae Sinica, 2023, 50(2): 371-381.
[14]	XIE Qian, ZHANG Shiyan, JIANG Lai, DING Mingyue, LIU Lingling, WU Rujian, CHEN Qingxi. Analysis of Canarium album Transcriptome SSR Information and Molecular Marker Development and Application [J]. Acta Horticulturae Sinica, 2023, 50(11): 2350-2364.
[15]	ZHAO Xueyan, WANG Qi, WANG Li, WANG Fangyuan, WANG Qing, LI Yan. Comparative Transcriptome Analysis of Differential Expression in Different Tissues of Corydalis yanhusuo [J]. Acta Horticulturae Sinica, 2023, 50(1): 177-187.