留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

转录组差异解析文心兰跳芽现象及其调控

兰寒俏 刘乐 黎维诗 郝代成 陈泰臻 陆玲 唐敏强 凌鹏

兰寒俏, 刘乐, 黎维诗, 等. 转录组差异解析文心兰跳芽现象及其调控[J]. 热带生物学报(中英文), 2026, 17(3): 532−540 doi:  10.15886/j.cnki.rdswxb.20240197
引用本文: 兰寒俏, 刘乐, 黎维诗, 等. 转录组差异解析文心兰跳芽现象及其调控[J]. 热带生物学报(中英文), 2026, 17(3): 532−540 doi:  10.15886/j.cnki.rdswxb.20240197
Lan Hanqiao, Liu Le, Li Weishi, Hao Daicheng, Chen Taizhen, Lu Ling, Tang Minqiang, Ling Peng. Transcriptome differences resolve the phenomenon of bud jumping in Oncidium and its regulation[J]. Journal of Tropical Biology, 2026, 17(3): 532-540. doi: 10.15886/j.cnki.rdswxb.20240197
Citation: Lan Hanqiao, Liu Le, Li Weishi, Hao Daicheng, Chen Taizhen, Lu Ling, Tang Minqiang, Ling Peng. Transcriptome differences resolve the phenomenon of bud jumping in Oncidium and its regulation[J]. Journal of Tropical Biology, 2026, 17(3): 532-540. doi: 10.15886/j.cnki.rdswxb.20240197

转录组差异解析文心兰跳芽现象及其调控

DOI: 10.15886/j.cnki.rdswxb.20240197 CSTR: 32425.14.j.cnki.rdswxb.20240197
基金项目: 国家自然科学基金项目(32201624);海南大学生态文明协同创新中心项目(XTCX2022STC10);海南大学启动科研基金项目(KYQD(ZR)-21039);海南大学横向自然科学项目(RH2200000309)
详细信息
    第一作者:

    兰寒俏(1998—),女,海南大学热带农林学院风景园林2022级硕士研究生。E-mail:hanqiao3414@foxmail.com

    通信作者:

    唐敏强(1988—),男,博士。研究方向:热带植物生物信息与遗传进化。E-mail:tangminqiang@hainanu.edu.cn

    凌鹏(1965—),男,教授,博士生导师。研究方向:园林观赏植物、果树资源遗传育种及应用。E-mail:993037@hainanu.edu.cn

  • 中图分类号: S682.31

Transcriptome differences resolve the phenomenon of bud jumping in Oncidium and its regulation

  • 摘要: 文心兰(Oncidium hybridum)切花生产中常出现不定芽跳跃性发育打断正常花芽分化生长周期现象,由此造成切花产量下降,为解决此问题,以文心兰‘博大一号’的花芽和营养芽为材料,利用高通量技术进行转录测序。结果获得127 452 717个高质量序列,共37.36 Gb,鉴定出7 671个差异表达基因(DEGs)。COG功能分类主要在信号转导途径和碳水化合物运输和代谢途径,DEGs显著富集的KEGG通路有淀粉和蔗糖代谢途径、植物激素信号转导等途径,其中植物激素信号转导途径获得78个DEGs,差异最显著的涉及生长素、细胞分裂素、水杨酸、赤霉素。筛选到影响开花的相关转录因子和基因有13个,包括MADS1AP2FLK等。结果在一定程度上探析了生长素、细胞分裂素等激素及这些开花基因对文心兰花芽和营养芽分化的影响,为进一步深入研究文心兰生产过程中不定芽跳跃性发育现象形成的机制提供理论基础,便于后续提高文心兰生产效率和质量、优化文心兰栽培管理技术,这对推动中国兰花育种、种苗生产和产业发展具有十分重要的现实意义。
  • 图  2  文心兰花芽及营养芽样本示意图

    注:1.第一片子叶;2.第二片子叶;3.第3片子叶;4.第4片子叶。

    Fig.  2  Diagram of flower bud and vegetative bud samples of Oncidium

    Note: 1. First cotyledon; 2. Second cotyledon; 3. Third cotyledon; 4. Fourth cotyledon.

    图  3  营养芽 VS 花芽的DEGs COG功能分类图

    Fig.  3  COG functional classification diagram of DEGs between vegetative bud vs flower bud

    图  4  营养芽 VS 花芽的DEGs KEGG通路富集条形图

    Fig.  4  Barplot graph of KEGG pathway enrichment for DEGs between vegetative bud vs flower bud

    图  5  营养芽 VS 花芽显著差异表达基因激素信号转导途径及相关热图

    Fig.  5  Heatmap of significantly differentially expressed genes involved in hormone signaling pathways between vegetative bud vs flower bud

    图  6  qRT-PCR数据统计

    注:*表示差异,P≤0.05;**表示差异显著,P≤0.005;****表示差异极显著,P≤0.000 1。

    Fig.  6  Statistics of qRT-PCR data

    Note: * indicates difference at P ≤ 0.05; ** indicates a significant difference at P ≤ 0.005; **** indicates a highly significant difference at P ≤ 0.0001.

    表  1  关键基因引物序列

    Table  1  The primer sequences of key genes.

    基因
    Genes
    正向引物(5′—3′)
    Forward primers (5′—3′)
    反向引物(5′—3′)
    Reverse primers (5′—3′)
    Actin GCAACATTGTTCTTAGCGGAGGCT TCTTCATGCTGCTTGGTGCAAGTG
    AUX/IAA22D AACAAAGCTCGACCTCCAAACG GTGTTGCGCTACAATGGTGTTG
    SAUR-6B ATGGGTTCATCGTCTCCGTCTC AACAGCACGATCCGCAGAACG
    MADS1 TGTGAGTGTGCTTGGTTGCTATC TTGATTGAAGACCCGAAAGATCCC
    PIF4 AGGAGACTATGAATCTGCCGATGG AAACTTGAACTTTGCTTACCCTTTCTG
    FLK GCAGAATCTACAGATGGAGCACAG GAGGCGGCACGTACTGAGG
    下载: 导出CSV

    表  2  文心兰营养芽 VS 花芽的DEGs注释统计

    Table  2  Annotation and statistics of DEGs between vegetative buds and flower buds in Oncidium.

    数据库
    Database
    COG GO KEGG KOG NR Pfam Swiss-Prot eggNOG
    数量
    Number
    2 662 3 237 2 894 3 893 7 004 5 891 6 194 7 227
      注:共有7420个DEGs,至少被一个数据库注释(各数据库注释数量有重叠,不计为总和)。
      Note: A total of 7420 DEGs were annotated by at least one database (the number of annotations across databases overlaps and is not included in the total count).
    下载: 导出CSV

    表  3  开花相关的差异表达转录因子及基因统计

    Table  3  Statistics of differentially expressed transcription factors and genes related to flowering

    转录因子/基因
    transcription factors/genes
    表达趋势
    Expression trend
    转录因子/基因
    transcription factors/genes
    表达趋势
    Expression trend
    转录因子/基因
    transcription factors/genes
    表达趋势
    Expression trend
    OhGene045464(Forever
    young flower protein)
    上up OhGene033882(AP2) 下down OhGene046073(MADS16) 下down
    OhGene015200(FLK) 下down OhGene015557(AP2) 下down MSTRG.8385(MADS2) 上up
    OhGene004083(CAULIFLOWER/MADS-box) 上up OhGene033883(AP2) 下down OhGene004083(MADS1) 上up
    OhGene004080(CAULIFLOWER/MADS-box) 上up MSTRG.52270(MADS2) 上up OhGene004080(MADS1) 上up
    OhGene021037(AP2) 下down
    下载: 导出CSV
  • [1] 张正梁, 尹华. 文心兰分株繁殖和控花栽培[J]. 中国花卉园艺, 2003(12): 26−27.
    [2] Chin D C, Hsieh C C, Lin H Y, et al. A low glutathione redox state couples with a decreased ascorbate redox ratio to accelerate flowering in Oncidium orchid [J]. Plant & Cell Physiology, 2016, 57(2): 423−436.
    [3] 张海良. 兰花芽分化的分子机理和相关功能基因克隆 [D]. 广州: 华南农业大学, 2018.
    [4] 董晓宇, 张晶, 符真珠, 等. 蝴蝶兰腋芽增殖过程中的转录组特性[J]. 基因组学与应用生物学, 2018, 37(3): 1265−1270.
    [5] Li Z, Xiao W, Chen H, et al. Transcriptome analysis reveals endogenous hormone changes during spike development in Phalaenopsis [J]. International Journal of Molecular Sciences, 2022, 23(18): 10461.
    [6] Fan Z, Li J, Li X, et al. Genome-wide transcriptome profiling provides insights into floral bud development of summer-flowering Camellia azalea [J]. Scientific Reports, 2015, 5: 9729.
    [7] 石玉波. 百子莲花芽分化过程中比较转录组分析及开花相关基因的克隆 [D]. 哈尔滨: 东北林业大学, 2014.
    [8] Pochamreddy M, Haim D, Halon E, et al. Alternate bearing in ‘hass’ avocado: fruit load-induced changes in bud auxin homeostasis are associated with flowering repression [J]. Journal of Experimental Botany, 2024, 75(18): 5717−5733.
    [9] Chandler J W. The hormonal regulation of flower development [J]. Journal of Plant Growth Regulation, 2011, 30(2): 242−254.
    [10] Li Y, Zhang B, Yu H. Molecular genetic insights into orchid reproductive development [J]. Journal of Experimental Botany, 2022, 73(7): 1841−1852.
    [11] 史梅容, 舒文波, 邱明萱, 等. 兰花花器官及成花基因调控研究进展[J]. 中国农业大学学报, 2023, 28(7): 57−67. https://doi.org/10.11841/j.issn.1007-4333.2023.07.05 doi:  10.11841/j.issn.1007-4333.2023.07.05
    [12] 罗远华, 王振波, 黄敏玲, 等. 文心兰不同生育期茎叶生理指标的动态变化[J]. 福建农业学报, 2017, 32(7): 719−723.
    [13] Mouhu K, Kurokura T, Koskela E A, et al. The Fragaria vesca homolog of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 represses flowering and promotes vegetative growth [J]. The Plant Cell, 2013, 25(9): 3296−3310. https://doi.org/10.1105/tpc.113.115055 doi:  10.1105/tpc.113.115055
    [14] Feng J, Wang Y, Ge W, et al. Regulatory mechanism of the miR172e-LbrAP2 module during the vegetative growth phase transition in Lilium [J]. Planta, 2023, 259(1): 26.
    [15] 史绍林. 红松营养生长与生殖生长转换中植物激素动态研究 [D]. 哈尔滨: 东北林业大学, 2020.
    [16] Chen S, Zhou Y, Chen Y, et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor [J]. Bioinformatics, 2018, 34(17): 884−890. https://doi.org/10.1093/bioinformatics/bty560 doi:  10.1093/bioinformatics/bty560
    [17] 金洲, 卢山, 江俊浩, 等. 园艺植物花芽分化影响因素及机理研究进展[J]. 园艺学报, 2023, 50(5): 1151−1164.
    [18] 李川, 王瑞娴, 崔鸣. 花魔芋花芽和叶芽的转录组分析[J]. 分子植物育种, 2020, 18(16): 5315−5322.
    [19] Pei L, Gao Y, Feng L, et al. Phenolic acids and flavonoids play important roles in flower bud differentiation in Mikania micrantha: transcriptomics and metabolomics [J]. International Journal of Molecular Sciences, 2023, 24(23): 16550. https://doi.org/10.3390/ijms242316550 doi:  10.3390/ijms242316550
    [20] Liu W, Feng Y, Yu S, et al. The flavonoid biosynthesis network in plants [J]. International Journal of Molecular Sciences, 2021, 22(23): 12824. https://doi.org/10.3390/ijms222312824 doi:  10.3390/ijms222312824
    [21] Li W, Xu P, Qian C, et al. The combined analysis of the transcriptome and metabolome revealed the possible mechanism of flower bud formation in Amorphophallus bulbifer [J]. Agronomy, 2024, 14(3): 519. https://doi.org/10.3390/agronomy14030519 doi:  10.3390/agronomy14030519
    [22] 石长双. 马尾松短枝腋芽萌发关键基因的挖掘 [D]. 贵阳: 贵州大学, 2020.
    [23] 邹礼平, 潘铖, 王梦馨, 等. 激素调控植物成花机理研究进展[J]. 遗传, 2020, 42(8): 739−751.
    [24] Chao W S, Doğramaci M, Horvath D P, et al. Phytohormone balance and stress-related cellular responses are involved in the transition from bud to shoot growth in leafy spurge [J]. BMC Plant Biology, 2016, 16: 47. https://doi.org/10.1186/s12870-016-0735-2 doi:  10.1186/s12870-016-0735-2
    [25] Tong N, Shu Q, Wang B, et al. Histology, physiology, and transcriptomic and metabolomic profiling reveal the developmental dynamics of annual shoots in tree peonies (Paeonia suffruticosa Andr.)[J]. Horticulture Research, 2023, 10(9): uhad152. https://doi.org/10.1093/hr/uhad152 doi:  10.1093/hr/uhad152
    [26] Dharmasiri N, Dharmasiri S, Estelle M. The F-box protein TIR1 is an auxin receptor [J]. Nature, 2005, 435(7041): 441−445.
    [27] Chapman E J, Estelle M. Mechanism of auxin-regulated gene expression in plants [J]. Annual Review of Genetics, 2009, 43: 265−285. https://doi.org/10.1146/annurev-genet-102108-134148 doi:  10.1146/annurev-genet-102108-134148
    [28] Swarup R, Parry G, Graham N, et al. Auxin cross-talk: integration of signalling pathways to control plant development [J]. Plant Molecular Biology, 2002, 49(3/4): 411−426.
    [29] Werner T, Motyka V, Laucou V, et al. Cytokinin-deficient transgenic Arabidopsis plants show multiple developmental alterations indicating opposite functions of cytokinins in the regulation of shoot and root meristem activity [J]. The Plant Cell, 2003, 15(11): 2532−2550. https://doi.org/10.1105/tpc.014928 doi:  10.1105/tpc.014928
    [30] Zhang Y J, Nie C R, Guo W J, et al. Exploring flowering genes in Phalaenopsis through transcriptome analysis and critical gene validation of hormone signal transduction pathway [J]. Russian Journal of Plant Physiology, 2023, 70(3): 25. https://doi.org/10.1134/S1021443722601938 doi:  10.1134/S1021443722601938
    [31] 徐东东, 东琳, 邵丽, 等. TGA转录因子在调控植物逆境应答和生长发育中的作用研究进展[J]. 植物生理学报, 2024, 60(7): 1079−1086.
    [32] Amasino R. Seasonal and developmental timing of flowering [J]. The Plant Journal, 2010, 61(6): 1001−1013. https://doi.org/10.1111/j.1365-313X.2010.04148.x doi:  10.1111/j.1365-313X.2010.04148.x
    [33] Lyons R, Rusu A, Stiller J, et al. Investigating the association between flowering time and defense in the Arabidopsis thaliana-Fusarium oxysporum interaction [J]. PLoS One, 2015, 10(6): e0127699. https://doi.org/10.1371/journal.pone.0127699 doi:  10.1371/journal.pone.0127699
    [34] Lim M H, Kim J, Kim Y S, et al. A new Arabidopsis gene, FLK encodes an RNA binding protein with K homology motifs and regulates flowering time via FLOWERING LOCUS C [J]. The Plant Cell, 2004, 16(3): 731−740. https://doi.org/10.1105/tpc.019331 doi:  10.1105/tpc.019331
    [35] Xie W, Ding C, Hu H, et al. Molecular events of rice AP2/ERF transcription factors [J]. International Journal of Molecular Sciences, 2022, 23(19): 12013. https://doi.org/10.3390/ijms231912013 doi:  10.3390/ijms231912013
    [36] Mehrnia M, Balazadeh S, Zanor M I, et al. EBE, an AP2/ERF transcription factor highly expressed in proliferating cells, affects shoot architecture in Arabidopsis [J]. Plant Physiology, 2013, 162(2): 842−857. https://doi.org/10.1104/pp.113.214049 doi:  10.1104/pp.113.214049
    [37] 李元元, 王鲁, 苏振刚, 等. MADS-box基因控制植物成花的分子机理[J]. 基因组学与应用生物学, 2010, 29(6): 1122−1132. https://doi.org/10.3969/gab.029.001122 doi:  10.3969/gab.029.001122
    [38] Thiruvengadam M, Chung I M, Yang C H. Overexpression of Oncidium MADS box (OMADS1) gene promotes early flowering in transgenic orchid (Oncidium Gower Ramsey)[J]. Acta Physiologiae Plantarum, 2012, 34(4): 1295−1302. https://doi.org/10.1007/s11738-012-0926-x doi:  10.1007/s11738-012-0926-x
    [39] Hsu H F, Hsieh W P, Chen M K, et al. C/D class MADS box genes from two monocots, orchid (Oncidium Gower Ramsey) and lily (Lilium longiflorum), exhibit different effects on floral transition and formation in Arabidopsis thaliana [J]. Plant and Cell Physiology, 2010, 51(6): 1029−1045. https://doi.org/10.1093/pcp/pcq052 doi:  10.1093/pcp/pcq052
  • [1] 孟丽华, 何子婷, 张育华, 李轲祎, 吴金烽, 吴奕彤, 罗敏, 胡平.  基于mt COⅠ序列的广西松材线虫遗传多样性和遗传分化的初探 . 热带生物学报(中英文), 2026, 17(1): 101-107. doi: 10.15886/j.cnki.rdswxb.20240135
    [2] 范弋豪, 张叶, 叶蓓蕾, 黎维诗, 郝代成, 凌鹏.  高温胁迫下文心兰的生理响应机制 . 热带生物学报(中英文), 2026, 17(): 1-11. doi: 10.15886/j.cnki.rdswxb.20260060
    [3] 杨敏, 黄婧, 曹毓杰, 李东栋, 郑育声.  基于转录组的蛇皮果雌雄株差异基因的筛选和性别鉴定 . 热带生物学报(中英文), 2025, 16(1): 13-20. doi: 10.15886/j.cnki.rdswxb.20240110
    [4] 李佳雪, 丁一, 王猛, 李涛, 郭攀阳, 刘成立, 韦双双, 黄家权, 李洪立, 胡文斌, 汤华.  火龙果HubHLH基因家族的全基因组分析及其对冬季补光诱导开花的表达响应 . 热带生物学报(中英文), 2024, 15(2): 198-209. doi: 10.15886/j.cnki.rdswxb.20220108
    [5] 全飞, 兰国玉, 魏亚情, 李明美, 孙树晴, 杜昊楠.  基于宏基因组测序的橡胶林病原微生物的分析 . 热带生物学报(中英文), 2022, 13(1): 27-35. doi: 10.15886/j.cnki.rdswxb.2022.01.005
    [6] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才.  转录组荔枝Dof基因家族的鉴定及其表达 . 热带生物学报(中英文), 2021, 12(1): 7-14. doi: 10.15886/j.cnki.rdswxb.2021.01.002
    [7] 汪玉玲, 宋希强, 郭向阳, 王洪星, 郭梨锦.  文心兰炭疽病生防潜力菌的鉴定及其抑菌性的测定 . 热带生物学报(中英文), 2021, 12(2): 236-243. doi: 10.15886/j.cnki.rdswxb.2021.02.013
    [8] 宋娜, 张刘宁颖, 曹敏, 郭文雅, 洪雨慧, 吴金山, 陈银华, 于晓惠.  木薯ERF转录因子调控的靶基因筛选与表达分析 . 热带生物学报(中英文), 2020, 11(2): 177-189. doi: 10.15886/j.cnki.rdswxb.2020.02.008
    [9] 闫冰玉, 巩笑笑, 谭玉荣, 王丹, 高璇, 张恒, 李双江, 王鹏, 潘英文, 刘进平.  文心兰OnRR10基因的克隆及表达分析 . 热带生物学报(中英文), 2020, 11(3): 288-295. doi: 10.15886/j.cnki.rdswxb.2020.03.005
    [10] 闫冰玉, 巩笑笑, 谭玉荣, 王鹏, 李双江, 潘英文, 刘进平.  文心兰OnRR9基因克隆及表达 . 热带生物学报(中英文), 2018, 9(4): 393-400. doi: 10.15886/j.cnki.rdswxb.2018.04.005
    [11] 程裕轩, 杨翠萍, 田晓岩, 石乐松, 胡进, 闫冰玉, 巩笑笑, 庄玉粉, 潘英文, 刘进平.  在文心兰不同开花阶段施用保鲜剂对切花衰老的影响 . 热带生物学报(中英文), 2017, 8(1): 48-51,57. doi: 10.15886/j.cnki.rdswxb.2017.01.008
    [12] 石乐松, 孙慧, 田晓岩, 吴繁花, 蔡秀清, 庄玉粉, 刘进平.  文心兰南茜种花发育进程与分期 . 热带生物学报(中英文), 2016, 7(4): 440-443,449. doi: 10.15886/j.cnki.rdswxb.2016.04.006
    [13] 黄浩宇, 肖云骁, 石乐松, 田晓岩, 刘进平.  文心兰切花花被不同开放阶段酶活性的变化 . 热带生物学报(中英文), 2016, 7(2): 199-203. doi: 10.15886/j.cnki.rdswxb.2016.02.010
    [14] 勾昕, 胡薇薇, 范亚飞, 熊兵, 石乐松, 田晓岩, 刘进平.  文心兰切花不同开放阶段花被的生理生化变化 . 热带生物学报(中英文), 2016, 7(1): 70-75. doi: 10.15886/j.cnki.rdswxb.2016.01.012
    [15] 范玉龙, 姜达, 李明容, 黄惜, 夏志辉.  适合于高通量测序的槟榔总RNA的提取方法 . 热带生物学报(中英文), 2015, 6(1): 87-90. doi: 10.15886/j.cnki.rdswxb.2015.01.013
    [16] 张贺, 漆艳香, 刘晓妹, 张欣, 蒲金基, 张辉强, 谢艺贤.  粉蕉枯萎病病原菌1号小种致病力的分化 . 热带生物学报(中英文), 2014, 5(2): 136-141,146. doi: 10.15886/j.cnki.rdswxb.2014.02.010
    [17] 赵羿鸾, 黄琴, 贾贤, 王颖, 吴坤鑫, 陈雄庭.  文心兰染色体的制片技术 . 热带生物学报(中英文), 2013, 4(2): 198-202. doi: 10.15886/j.cnki.rdswxb.2013.02.016
    [18] 陈兴汉, 刘晓春, 蒙子宁, 张勇, 林浩然, 叶卫.  尼罗罗非鱼仔鱼个体发育及性腺分化的发育 . 热带生物学报(中英文), 2012, 3(3): 197-203. doi: 10.15886/j.cnki.rdswxb.2012.03.017
    [19] 蔡波, 陈施明, 敖苏, 周卫川, 林明光.  文心兰设施栽培中蜗牛和蛞蝓为害规律及综合防治 . 热带生物学报(中英文), 2012, 3(4): 329-332. doi: 10.15886/j.cnki.rdswxb.2012.04.003
    [20] 柯海丽, 龙建勇, 黎维诗, 凌绪柏.  不同的遮光率对文心兰切花生长及品质的影响 . 热带生物学报(中英文), 2011, 2(3): 264-266. doi: 10.15886/j.cnki.rdswxb.2011.03.017
  • 加载中
  • 图(6) / 表 (3)
    计量
    • 文章访问数:  806
    • HTML全文浏览量:  479
    • PDF下载量:  28
    • 被引次数: 0
    出版历程
    • 收稿日期:  2024-12-24
    • 录用日期:  2025-03-21
    • 修回日期:  2025-03-03
    • 网络出版日期:  2026-05-20
    • 刊出日期:  2026-05-25

    转录组差异解析文心兰跳芽现象及其调控

    DOI: 10.15886/j.cnki.rdswxb.20240197
      基金项目:  国家自然科学基金项目(32201624);海南大学生态文明协同创新中心项目(XTCX2022STC10);海南大学启动科研基金项目(KYQD(ZR)-21039);海南大学横向自然科学项目(RH2200000309)
      作者简介:

      兰寒俏(1998—),女,海南大学热带农林学院风景园林2022级硕士研究生。E-mail:hanqiao3414@foxmail.com

      通讯作者: 唐敏强(1988—),男,博士。研究方向:热带植物生物信息与遗传进化。E-mail:tangminqiang@hainanu.edu.cn凌鹏(1965—),男,教授,博士生导师。研究方向:园林观赏植物、果树资源遗传育种及应用。E-mail:993037@hainanu.edu.cn
    • 中图分类号: S682.31

    摘要: 文心兰(Oncidium hybridum)切花生产中常出现不定芽跳跃性发育打断正常花芽分化生长周期现象,由此造成切花产量下降,为解决此问题,以文心兰‘博大一号’的花芽和营养芽为材料,利用高通量技术进行转录测序。结果获得127 452 717个高质量序列,共37.36 Gb,鉴定出7 671个差异表达基因(DEGs)。COG功能分类主要在信号转导途径和碳水化合物运输和代谢途径,DEGs显著富集的KEGG通路有淀粉和蔗糖代谢途径、植物激素信号转导等途径,其中植物激素信号转导途径获得78个DEGs,差异最显著的涉及生长素、细胞分裂素、水杨酸、赤霉素。筛选到影响开花的相关转录因子和基因有13个,包括MADS1AP2FLK等。结果在一定程度上探析了生长素、细胞分裂素等激素及这些开花基因对文心兰花芽和营养芽分化的影响,为进一步深入研究文心兰生产过程中不定芽跳跃性发育现象形成的机制提供理论基础,便于后续提高文心兰生产效率和质量、优化文心兰栽培管理技术,这对推动中国兰花育种、种苗生产和产业发展具有十分重要的现实意义。

    English Abstract

    兰寒俏, 刘乐, 黎维诗, 等. 转录组差异解析文心兰跳芽现象及其调控[J]. 热带生物学报(中英文), 2026, 17(3): 532−540 doi:  10.15886/j.cnki.rdswxb.20240197
    引用本文: 兰寒俏, 刘乐, 黎维诗, 等. 转录组差异解析文心兰跳芽现象及其调控[J]. 热带生物学报(中英文), 2026, 17(3): 532−540 doi:  10.15886/j.cnki.rdswxb.20240197
    Lan Hanqiao, Liu Le, Li Weishi, Hao Daicheng, Chen Taizhen, Lu Ling, Tang Minqiang, Ling Peng. Transcriptome differences resolve the phenomenon of bud jumping in Oncidium and its regulation[J]. Journal of Tropical Biology, 2026, 17(3): 532-540. doi: 10.15886/j.cnki.rdswxb.20240197
    Citation: Lan Hanqiao, Liu Le, Li Weishi, Hao Daicheng, Chen Taizhen, Lu Ling, Tang Minqiang, Ling Peng. Transcriptome differences resolve the phenomenon of bud jumping in Oncidium and its regulation[J]. Journal of Tropical Biology, 2026, 17(3): 532-540. doi: 10.15886/j.cnki.rdswxb.20240197
    • 文心兰(Oncidium hybridum)为多年生草本,其生长史从幼苗期发育到假鳞茎发育再到假鳞茎成熟萌发不定芽阶段,紧接抽苔阶段,由不定芽发育成花芽并分化成花序和各花器官直至开花阶段,待花衰败后,由基部的不定芽萌发阶段过渡到发育成营养芽阶段(图1[1-2]。前3个阶段(i~iii)为营养生长阶段,最后阶段(vi)的营养芽会继续发育成另一个新植株。若在假鳞茎成熟阶段(iii)时不定芽分化出的不是花芽,而直接跳跃发育到营养芽阶段(vi),即为跳芽现象。

      图  1  文心兰生长史

      Figure 1.  Growth history of Oncidium hybridum

      利用现代分子手段及生物信息分析方法研究关于兰花营养芽、花芽分化发育的基因表达调控途径越来越常见。在兰花快繁中,用转录测序初步确定与兰花芽分化相关的基因[3],蝴蝶兰(Phalaenopsis spp.)腋芽增殖能力相关基因[4]。Li等[5]研究内源激素代谢和信号通路对蝴蝶兰腋芽发育至成花发育的分子调控机制。Fan等[6]认为杜鹃叶山茶(Camellia azalea)的花芽发育与植物激素途径和植物代谢等多种途径调控有关。在东方百子莲(Agapanthus praecox ssp. orientalis)的成花过程中生长素信号发挥着重要调控作用[7],也有研究发现鳄梨(Persea americana) 高果实负荷(heavy fruit load,HFL) 在芽中产生生长素信号,诱导产生花抑制因子 PaTFL1,从而抑制花序发育[8]。除了生长素、细胞分裂素对花芽分化过程的调控研究外[9],在兰花中,开花启动子相关基因能促进开花转变,而抑制开花启动的相关基因则以拮抗的形式抑制开花[10-11]。而在文心兰芽分化研究中,假鳞茎内糖类物质的变化对子代假鳞茎侧芽的萌发影响很大,假鳞茎成熟阶段淀粉会大量转化为供应芽萌发的能量[12]

      文心兰作为一种有着较高观赏价值的兰花切花和盆花,具有重要的市场价值和潜力。海南省气候资源优越,加上政策支持、技术加持等条件,为兰花产业的发展增添了强劲助力。在海南博大兰花科技有限公司东山兰花基地的文心兰切花生产过程中发现了文心兰的几个品种(‘博大一号’、‘柠檬绿’)中均出现跳芽现象,该现象发生过多会导致一整个产花周期无花可产,总体切花产量下降,而在相同条件下其他品种的文心兰却少有跳芽现象。在触发新生营养芽和花芽分化过程中,因某些开花相关基因受到抑制或延迟[13],可能导致花芽分化失败,直接进入营养芽阶段;或营养生长相关基因被过度激活,而直接促进营养芽发育[14];又或者复杂的激素信号通路、基因网络调控也会打破营养生长与生殖生长的平衡[15],导致该现象的发生。目前,虽有前人指出文心兰跳芽的现象,但并未将跳芽现象进行概念化归纳,同时从生产实际中探索解决跳芽现象的方法仍很少,对于文心兰花芽/营养芽分化的分子机制少见报道。

      本研究探究文心兰花芽及营养芽分化的分子机制,以期初步了解调控文心兰花芽和营养芽生长发育的关键因素,为减少文心兰切花生产过程中不定芽跳跃性发育现象提供理论基础,以提高文心兰生产效率和质量、优化文心兰栽培管理技术,对推动中国兰花育种、种苗生产和产业发展具有十分重要的现实意义。

      • 检测样品从基地同一批次生长健壮的文心兰‘博大一号’(Oncidium hybridum 'Boda NO1')植株中随机选择处于抽苔阶段的当年生假鳞茎基部第1或2子叶内长出约10 mm花芽的植株共3株;随机选择处于营养芽阶段的当年生假鳞茎基部第3、4子叶内长出约10 mm营养芽的植株共3株,取其花芽、营养芽(图2),剥去叶(鞘)片,称取约2~3 g,放入冻存管中,共计6个样,用于转录组测序实验。

        图  2  文心兰花芽及营养芽样本示意图

        Figure 2.  Diagram of flower bud and vegetative bud samples of Oncidium

      • 首先,使用天根植物多糖多酚试剂盒提供的标准提取方法提取 RNA,然后对提取的 RNA 进行浓度、纯度、完整性检测,以确保其符合测序质量标准。其次,构建 cDNA 文库,并在 Illumina NovaSeq 6000 平台上对 cDNA 文库进行测序。最后,对转录组数据进行质量控制,以消除原始转录组数据中的冗余序列,从而获得高质量的纯净读数。

      • 使用FASTP v0.19.3[16]对所有样本数据进行质量控制和优化评估(合格数据定义为 Q30>85%),而后将高质量序列与参考基因组进行序列比对,以营养芽为对照组,将差异倍数(Fold Change)≥2且PDR<0.05作为筛选标准。将筛选出的DEGs与GO、COG、KEGG等数据库进行序列比对,获得这些基因的注释信息。

      • 提取转录组样本花芽、营养芽的RNA,检测合格的RNA用于后续实验。使用HiScript Ⅱ 1st Strand cDNA Synthesis Kit(+gDNA wiper)反转录试剂盒对RNA进行反转录,反转后的cDNA用于后续实验。使用Primer v6.0设计qRT-PCR特异性引物,引物序列见表1。使用并按照ChamQ SYBR qPCR Master Mix试剂盒说明书进行荧光定量PCR,选择Actin基因作为内参基因。数据结果分析使用2−ΔΔCT法计算相对表达量。

        表 1  关键基因引物序列

        Table 1.  The primer sequences of key genes.

        基因
        Genes
        正向引物(5′—3′)
        Forward primers (5′—3′)
        反向引物(5′—3′)
        Reverse primers (5′—3′)
        Actin GCAACATTGTTCTTAGCGGAGGCT TCTTCATGCTGCTTGGTGCAAGTG
        AUX/IAA22D AACAAAGCTCGACCTCCAAACG GTGTTGCGCTACAATGGTGTTG
        SAUR-6B ATGGGTTCATCGTCTCCGTCTC AACAGCACGATCCGCAGAACG
        MADS1 TGTGAGTGTGCTTGGTTGCTATC TTGATTGAAGACCCGAAAGATCCC
        PIF4 AGGAGACTATGAATCTGCCGATGG AAACTTGAACTTTGCTTACCCTTTCTG
        FLK GCAGAATCTACAGATGGAGCACAG GAGGCGGCACGTACTGAGG
      • 纯化后的RNA经Agilent 2100仪器检测,结果显示RIN值在8.4~9.5范围内,样品质量满足建库测序要求。通过清洗数据后,得到高质量序列127 452 717个,共37.36 Gb。各样品Q30均在94.62%以上,表明测序数据可靠,所构建文库符合要求,可用于后续分析。

        通过筛选得到7 671个DEGs,其中上调表达为3 695个,下调表达为3 976个。为更清晰地了解这些DEGs的情况,对其进行COG、GO、KEGG等数据库的功能注释。从表2中可知,共有7 420个DEGs被赋予注释信息,注释最多的数据库为eggNOG,为7 227个。注释较少的KEGG数据库注释到2 894个DEGs,COG数据库注释到2 662个。而一些未注释到的基因可能是未发现的新基因,也可能与参考基因组有关。

        表 2  文心兰营养芽 VS 花芽的DEGs注释统计

        Table 2.  Annotation and statistics of DEGs between vegetative buds and flower buds in Oncidium.

        数据库
        Database
        COG GO KEGG KOG NR Pfam Swiss-Prot eggNOG
        数量
        Number
        2 662 3 237 2 894 3 893 7 004 5 891 6 194 7 227
          注:共有7420个DEGs,至少被一个数据库注释(各数据库注释数量有重叠,不计为总和)。
          Note: A total of 7420 DEGs were annotated by at least one database (the number of annotations across databases overlaps and is not included in the total count).
      • 为进一步了解DEGs的直系同源信息,在COG数据库中(图3),DEGs分类到执行信号传导机制的数目最多,而后是碳水化合物运输和代谢,其次是一般功能预测和次生代谢产物的生物合成、运输和分解代谢。这些结果可以看出信号传导机制和碳水化合物是不定芽向花芽分化或营养芽分化的关键因素。

        图  3  营养芽 VS 花芽的DEGs COG功能分类图

        Figure 3.  COG functional classification diagram of DEGs between vegetative bud vs flower bud

      • 展示显著性q值最小的前20个通路,从图4可看出,DEGs显著富集在氨基酸、黄酮类化合物、苯丙烷等生物合成途径,淀粉和蔗糖、色氨酸等代谢途径,以及植物激素信号转导途径。其中,植物激素信号转导途径、氨基酸生物合成途径、苯丙烷生物合成途径富集的DEG数量最多。淀粉和蔗糖代谢、苯丙烷生物合成等途径参与植物次生代谢产物合成,色氨酸和植物激素信号转导途径等参与调控激素合成等植物生长发育代谢生理过程。

        图  4  营养芽 VS 花芽的DEGs KEGG通路富集条形图

        Figure 4.  Barplot graph of KEGG pathway enrichment for DEGs between vegetative bud vs flower bud

      • 植物激素信号转导途径是环境信息处理功能中的重要途径,在植物激素信号转导途径中,共筛选出78个DEGs,上调表达有33个,涉及生长素(Aux)、赤霉素(GA)、油菜素甾醇(BR)、脱落酸(ABA)、水杨酸(SA)、茉莉酸(JA);下调表达有45个,涉及Aux、GA、BR、ABA、SA、JA、细胞分裂素(CTK)。经分析,筛选出差异表达较为显著的4个信号通路(图5),在Aux信号转导途径中,SAURTIR、AUX/IAA22D基因及其同源基因表达量较高,后两者均上调表达,GH3基因也上调表达。在SA 信号途径中,TGA1TGA10 上调表达,TGA2下调表达,NPR5TGA5 则有上调或下调表达。在CTK信号转导途径中,AHP显著下调表达,赤霉素信号转导途径中仅有PIF4基因呈上调表达。

        图  5  营养芽 VS 花芽显著差异表达基因激素信号转导途径及相关热图

        Figure 5.  Heatmap of significantly differentially expressed genes involved in hormone signaling pathways between vegetative bud vs flower bud

      • 植物开花是一个多方因素参与调控的复杂过程,研究发现多种基因或转录因子与开花密切相关。本研究筛选出13个相关的基因及转录因子(表3),其中,上调表达7个,下调表达6个。影响花发育的相关关键基因有APETALA2AP2)、Flowering locus K homology domain isoform X1FLK)同源基因均下调表达,Forever young flower proteinCAULIFLOWER/MADS-box、MADS1、MADS2等6个MADS-box基因家族则上调表达。这些转录因子和基因的表达对文心兰芽分化起重要调控作用。

        表 3  开花相关的差异表达转录因子及基因统计

        Table 3.  Statistics of differentially expressed transcription factors and genes related to flowering

        转录因子/基因
        transcription factors/genes
        表达趋势
        Expression trend
        转录因子/基因
        transcription factors/genes
        表达趋势
        Expression trend
        转录因子/基因
        transcription factors/genes
        表达趋势
        Expression trend
        OhGene045464(Forever
        young flower protein)
        上up OhGene033882(AP2) 下down OhGene046073(MADS16) 下down
        OhGene015200(FLK) 下down OhGene015557(AP2) 下down MSTRG.8385(MADS2) 上up
        OhGene004083(CAULIFLOWER/MADS-box) 上up OhGene033883(AP2) 下down OhGene004083(MADS1) 上up
        OhGene004080(CAULIFLOWER/MADS-box) 上up MSTRG.52270(MADS2) 上up OhGene004080(MADS1) 上up
        OhGene021037(AP2) 下down
      • 为进一步验证转录组数据的可用性,用qRT-PCR实验对5个关键基因OhGene009120AUX/IAA22D)、OhGene011730SAUR-6B)、OhGene004080MADS1)、OhGene001497PIF4)、OhGene015200FLK)进行定量检测(图6),其变化趋势与转录组数据趋势一致。

        图  6  qRT-PCR数据统计

        Figure 6.  Statistics of qRT-PCR data

      • 芽分化是在环境条件和自身遗传等多种因素相互调控、共同作用下完成的。有研究表明,碳水化合物代谢和植物内源激素调控影响着花芽分化过程[17],在花魔芋(Amorphophallus konjac)的这一阶段过程中,糖代谢及内源激素扮演着至关重要的角色[18]。百子莲芽的KEGG通路富集在植物激素信号转导、淀粉和蔗糖代谢、苯丙烷生物合成途径上[7]。苯丙烷、黄酮类化合物是促进植物开花的次生活性代谢物质,这些物质生物合成途径在植物生长中发挥重要作用[19-20],本研究中均显著富集到这些途径,与魔芋(Amorphophallus bulbifer)的花芽和叶芽KEGG富集研究结果一致[21]。花芽和营养芽分化在植物激素信号转导及一些碳水化合物等代谢上有所差异,这与石长双[22]认为腋芽萌发与植物激素信号转导途径、糖类生物合成代谢途径有关的研究结果相近。

      • 植物感受外部环境变化信号,通过植物激素信号转导途径调控多种激素的协同或拮抗作用以调整自身的生长和分化进程[23]。从营养生长到开花的过程,植物内源激素尤其重要。根据激素调节假说,Aux、CTK、GA等是营养生长转向生殖生长的重要调控因子[24]AUX/IAA22D、SAUR、GH3同源基因均上调表达,而ARFs则多数下调表达,两者的相互作用可能影响生长素信号对芽分化的调控,推测花芽分化受其影响较大。TIR1与生长素分子结合,并通过促进AUX/IAA蛋白降解来调控ARF活性,进而调控植物生长。生长素信号输出受ARF-AUX/IAA相互作用的调节,在下游基因GH3作用下,IAA与氨基酸结合,同时生长素依赖极性运输向芽点运输,可能更有利于茎尖分生组织内花芽分化的调节[25-27]。此外,本研究中的所有CTK信号均下调表达,说明CTK在对营养芽分化有正向调节作用,而在花芽分化中受限。CTK与Aux也能相互作用于植物生长,特别是芽的发育和分化,CTK信号促进生长素转运,从而引发自身的降解或糖基化[28],CTK调节分生组织活性,并在生殖阶段限制发育器官中的细胞形成[29]。赤霉素在芽发育过程中起着重要作用,对于开花诱导的作用已经被前人在拟南芥上深入研究过。在本研究中,PIF4在花芽中高表达,PIF4可能是花芽分化的关键诱导基因,与PIF4可能是直接响应赤霉素促进蝴蝶兰成花的关键基因的结果一致[30]TGA10为上调表达,在花芽分化中起重要作用。TGA1TGA2TGA5广泛参与植株的抗病反应,并与SA信号途径中关键因子NPR1相互作用,而TGA10对植物花器官发育起关键作用[31]。此外,ABA、JA等一些激素会通过生物钟和光周期相关的开花途径参与开花调节,以响应非生物胁迫[32-33]

      • 本研究中FLK、AP2基因均下调表达,推测可能是这些基因在文心兰花期抑制着花芽分化,而导致不定芽发育为营养芽。在拟南芥中,FLK影响FLC等开花基因的表达,延迟花芽分化[34]AP2的过表达有时候会抑制开花[35],有时会直接参与侧芽的调控[36]MADS1MADS2在花芽中均高表达,对花芽发育分化起关键作用。MADS-box基因家族在花发育过程中起重要作用[37]MADS1能使转基因文心兰更早开花[38],来自兰花 (Oncidium Gower Ramsey)的MADS2能使转基因拟南芥提前或中度提前开花[39]。此外,Forever young flower protein、CAULIFLOWER/MADS-box也在花芽中高表达,这两种基因或许也促进花芽的发育,但仍需进一步实验验证。

      • 本研究对文心兰‘博大一号’的花芽和营养芽转录测序分析,结果在一定程度上探析了生长素、细胞分裂素、水杨酸等激素及一些影响开花的相关基因MADS1、AP2、FLK等对文心兰花芽和营养芽分化的影响,为进一步深入研究文心兰生产过程中不定芽跳跃性发育现象形成的机制提供理论基础,便于提高文心兰生产效率和质量、优化文心兰栽培管理技术,对推动中国兰花产业发展具有十分重要的现实意义。

    参考文献 (39)

    目录

      /

      返回文章
      返回