Search
2024 Volume 4
Article Contents
ARTICLE   Open Access    

Alternative splicing of the FLOWERING LOCUS C-like gene MaMADS33 is associated with endodormancy in mulberry

More Information
  • Received: 28 September 2023
    Revised: 02 August 2024
    Accepted: 14 August 2024
    Published online: 04 September 2024
    Forestry Research  4 Article number: e029 (2024)  |  Cite this article
  • Alternative splicing (AS) is an important post-transcriptional process that generates multiple mRNA isoforms. FLOWERING LOCUS C (FLC) is a pivotal gene in both the vernalization and autonomous pathways of flowering plants, and MaMADS33 is one of the FLC homologs in white mulberry (Morus alba). Recent studies have revealed that MaMADS33 is involved in endodormancy, but the underlying molecular mechanism remains to be characterized. Here, a comparison of MaMADS33 expression among three mulberry cultivars with different degrees of dormancy revealed a positive association between MaMADS33 expression and dormancy. Further 3' and 5' rapid amplification of cDNA ends (RACE) analyses led to identifying four MaMADS33 isoforms derived from AS and designated MaMADS33-AS14. Analysis of their coding potential revealed that MaMADS33-AS1 was a long non-coding RNA. Expression profiling and splicing-efficiency analyses showed that cold stress during endodormancy induced AS of MaMADS33, resulting in a predominance of truncated isoforms, especially MaMADS33-AS1. MaMADS33-AS2 expression was upregulated during both endodormancy and ecodormancy, whereas MaMADS33-AS3 and MaMADS33-AS4 were endodormancy-associated isoforms that were upregulated during endodormancy and then downregulated during ecodormancy. MaMADS33-AS4 was used as bait for a yeast two-hybrid screen because its gene expression was higher than that of MaMADS33-AS3, and mulberry winter-accumulating 18 kDa protein (MaWAP18) was identified as an MaMADS33-AS4 interaction partner. The interaction between MaWAP18 and MaMADS33-AS4 was confirmed by a bimolecular fluorescence complementation assay. These findings offer insight into the role of FLC homologs in the endodormancy of woody plants.
  • 加载中
  • Supplemental Table S1 Primer sequences used in this study.
    Supplemental Table S2 Clone statistics of MaMADS33-AS3 and MaMADS33-AS4.
    Supplemental Table S3 MaMADS33-AS4 yeast two-hybrid screening results.
  • [1]

    Nilsson O. 2022. Winter dormancy in trees. Current Biology 32:R630−R634

    doi: 10.1016/j.cub.2022.04.011

    CrossRef   Google Scholar

    [2]

    Lang GA, Early JD, Martin GC, Darnell RL. 1987. Endodormancy, paradormancy, and ecodormancy: physiological terminology and classification for dormancy research. HortScience 22:371−77

    doi: 10.21273/HORTSCI.22.3.371

    CrossRef   Google Scholar

    [3]

    Luna V, Reinoso H, Lorenzo E, Bottini R, Abdala G. 1991. Dormancy in peach ( Prunus persica L.) flower buds. Trees 5:244−46

    doi: 10.1007/BF00227532

    CrossRef   Google Scholar

    [4]

    Atkinson CJ, Brennan RM, Jones HG. 2013. Declining chilling and its impact on temperate perennial crops. Environmental and Experimental Botany 91:48−62

    doi: 10.1016/j.envexpbot.2013.02.004

    CrossRef   Google Scholar

    [5]

    Campoy JA, Ruiz D, Egea J . 2010. Effects of shading and thidiazuron + oil treatment on dormancy breaking, blooming and fruit set in apricot in a warm-winter climate. Scientia Horticulturae 125:203−10

    doi: 10.1016/j.scienta.2010.03.029

    CrossRef   Google Scholar

    [6]

    Guillamón JG, Dicenta F, Sánchez-Pérez R. 2021. Advancing endodormancy release in temperate fruit trees using agrochemical treatments. Frontiers in Plant Science 12:812621

    doi: 10.3389/fpls.2021.812621

    CrossRef   Google Scholar

    [7]

    Ito A, Sakaue T, Fujimaru O, Iwatani A, Ikeda T, et al. 2018. Comparative phenology of dormant Japanese pear (Pyrus pyrifolia) flower buds: a possible cause of 'flowering disorder'. Tree Physiology 38:825−39

    doi: 10.1093/treephys/tpx169

    CrossRef   Google Scholar

    [8]

    Tominaga A, Ito A, Sugiura T, Yamane H. 2021. How is global warming affecting fruit tree blooming? "Flowering (dormancy) disorder" in Japanese pear (Pyrus pyrifolia) as a case study. Frontiers in Plant Science 12:787638

    doi: 10.3389/fpls.2021.787638

    CrossRef   Google Scholar

    [9]

    Pagter M, Kjær KH. 2022. Winter warming stimulates vegetative growth and alters fruit quality of blackcurrant (Ribes nigrum). International Journal of Biometeorology 66:1391−401

    doi: 10.1007/s00484-022-02284-4

    CrossRef   Google Scholar

    [10]

    Olsen JE. 2010. Light and temperature sensing and signaling in induction of bud dormancy in woody plants. Plant Molecular Biology 73:37−47

    doi: 10.1007/s11103-010-9620-9

    CrossRef   Google Scholar

    [11]

    Bielenberg D, Wang Y, Fan S, Reighard GL, Scorza R, Abbott A. 2004. A deletion affecting several gene candidates is present in the Evergrowing peach mutant. Journal of Heredity 95:436−44

    doi: 10.1093/jhered/esh057

    CrossRef   Google Scholar

    [12]

    Bielenberg DG, Wang Y, Li Z, Zhebentyayeva T, Fan S, et al. 2008. Sequencing and annotation of the evergrowing locus in peach [Prunus persica (L.) Batsch] reveals a cluster of six MADS-box transcription factors as candidate genes for regulation of terminal bud formation. Tree Genetics & Genomes 4:495−507

    doi: 10.1007/s11295-007-0126-9

    CrossRef   Google Scholar

    [13]

    Jiménez S, Lawton-Rauh AL, Reighard GL, Abbott AG, Bielenberg DG. 2009. Phylogenetic analysis and molecular evolution of the dormancy associated MADS-box genes from peach. BMC Plant Biology 9:81

    doi: 10.1186/1471-2229-9-81

    CrossRef   Google Scholar

    [14]

    Li Z, Reighard GL, Abbott AG, Bielenberg DG. 2009. Dormancy-associated MADS genes from the EVG locus of peach [Prunus persica (L.) Batsch] have distinct seasonal and photoperiodic expression patterns. Journal of Experimental Botany 60:3521−30

    doi: 10.1093/jxb/erp195

    CrossRef   Google Scholar

    [15]

    Horvath D, Sung S, Kim D, Chao W, Anderson J. 2010. Characterization, expression and function of DORMANCY ASSOCIATED MADS-BOX genes from leafy spurge. Plant Molecular Biology 73:169−79

    doi: 10.1007/s11103-009-9596-5

    CrossRef   Google Scholar

    [16]

    Sasaki R, Yamane H, Ooka T, Jotatsu H, Kitamura Y, et al. 2011. Functional and expressional analyses of PmDAM genes associated with endodormancy in Japanese apricot. Plant Physiology 157:485−97

    doi: 10.1104/pp.111.181982

    CrossRef   Google Scholar

    [17]

    Wu RM, Walton EF, Richardson AC, Wood M, Hellens RP, et al. 2012. Conservation and divergence of four kiwifruit SVP-like MADS-box genes suggest distinct roles in kiwifruit bud dormancy and flowering. Journal of Experimental Botany 63:797−807

    doi: 10.1093/jxb/err304

    CrossRef   Google Scholar

    [18]

    Yamane H, Kashiwa Y, Ooka T, Tao R, Yonemori K. 2008. Suppression subtractive hybridization and differential screening reveals endodormancy-associated expression of an SVP/AGL24-type MADS-box gene in lateral vegetative buds of Japanese apricot. Journal of the American Society for Horticultural Science 133:708−16

    doi: 10.21273/JASHS.133.5.708

    CrossRef   Google Scholar

    [19]

    Saito T, Bai S, Ito A, Sakamoto D, Saito T, et al. 2013. Expression and genomic structure of the dormancy-associated MADS box genes MADS13 in Japanese pears (Pyrus pyrifolia Nakai) that differ in their chilling requirement for endodormancy release. Tree Physiology 33:654−67

    doi: 10.1093/treephys/tpt037

    CrossRef   Google Scholar

    [20]

    Wu R, Tomes S, Karunairetnam S, Tustin SD, Hellens RP, et al. 2017. SVP-like MADS box genes control dormancy and budbreak in apple. Frontiers in Plant Science 8:477

    doi: 10.3389/fpls.2017.00477

    CrossRef   Google Scholar

    [21]

    Rinne PL, Welling A, Vahala J, Ripel L, Ruonala R, et al. 2011. Chilling of dormant buds hyperinduces FLOWERING LOCUS T and recruits GA-inducible 1,3-β-glucanases to reopen signal conduits and release dormancy in Populus. The Plant Cell 23:130−46

    doi: 10.1105/tpc.110.081307

    CrossRef   Google Scholar

    [22]

    Mohamed R, Wang CT, Ma C, Shevchenko O, Dye SJ, et al. 2010. Populus CEN/TFL1 regulates first onset of flowering, axillary meristem identity and dormancy release in Populus. The Plant Journal 62:674−88

    doi: 10.1111/j.1365-313X.2010.04185.x

    CrossRef   Google Scholar

    [23]

    Voogd C, Wang T, Varkonyi-Gasic E. 2015. Functional and expression analyses of kiwifruit SOC1-like genes suggest that they may not have a role in the transition to flowering but may affect the duration of dormancy. Journal of Experimental Botany 66:4699−710

    doi: 10.1093/jxb/erv234

    CrossRef   Google Scholar

    [24]

    Li J, Xu Y, Niu Q, He L, Teng Y, et al. 2018. Abscisic acid (ABA) promotes the induction and maintenance of pear (Pyrus pyrifolia white pear group) flower bud endodormancy. International Journal of Molecular Sciences 19:310

    doi: 10.3390/ijms19010310

    CrossRef   Google Scholar

    [25]

    Li S, Wang Q, Wen B, Zhang R, Jing X, et al. 2021. Endodormancy release can be modulated by the GA4-GID1c-DELLA2 module in peach leaf buds. Frontiers in Plant Science 12:713514

    doi: 10.3389/fpls.2021.713514

    CrossRef   Google Scholar

    [26]

    Andrés F, Coupland G. 2012. The genetic basis of flowering responses to seasonal cues. Nature Reviews Genetics 13:627−39

    doi: 10.1038/nrg3291

    CrossRef   Google Scholar

    [27]

    Zhang J, Li Z, Mei L, Yao J, Hu C. 2009. PtFLC homolog from trifoliate orange (Poncirus trifoliata) is regulated by alternative splicing and experiences seasonal fluctuation in expression level. Planta 229:847−59

    doi: 10.1007/s00425-008-0885-z

    CrossRef   Google Scholar

    [28]

    Kumar G, Arya P, Gupta K, Randhawa V, Acharya V, et al. 2016. Comparative phylogenetic analysis and transcriptional profiling of MADS-box gene family identified DAM and FLC-like genes in apple (Malus × domestica). Scientific Reports 6:20695

    doi: 10.1038/srep20695

    CrossRef   Google Scholar

    [29]

    Porto DD, Bruneau M, Perini P, Anzanello R, Renou JP, et al. 2015. Transcription profiling of the chilling requirement for bud break in apples: a putative role for FLC-like genes. Journal of Experimental Botany 66:2659−72

    doi: 10.1093/jxb/erv061

    CrossRef   Google Scholar

    [30]

    Miotto Y, Tessele C, Czermainski A, Porto D, Falavigna V, et al. 2019. Spring is coming: genetic analyses of the bud break bate locus reveal candidate genes from the cold perception pathway to dormancy release in apple (Malus × domestica Borkh.). Frontiers in Plant Science 10:33

    doi: 10.3389/fpls.2019.00033

    CrossRef   Google Scholar

    [31]

    Nishiyama S, Matsushita MC, Yamane H, Honda C, Okada K, et al. 2021. Functional and expressional analyses of apple FLC-like in relation to dormancy progress and flower bud development. Tree Physiology 41:562−70

    doi: 10.1093/treephys/tpz111

    CrossRef   Google Scholar

    [32]

    Voogd C, Brian LA, Wu R, Wang T, Allan AC, et al. 2022. A MADS-box gene with similarity to FLC is induced by cold and correlated with epigenetic changes to control budbreak in kiwifruit. New Phytologist 233:2111−26

    doi: 10.1111/nph.17916

    CrossRef   Google Scholar

    [33]

    James AB, Syed NH, Bordage S, Marshall J, Nimmo GA, et al. 2012. Alternative splicing mediates responses of the Arabidopsis circadian clock to temperature changes. The Plant Cell 24:961−81

    doi: 10.1105/tpc.111.093948

    CrossRef   Google Scholar

    [34]

    Wright CJ, Smith CWJ, Jiggins CD. 2022. Alternative splicing as a source of phenotypic diversity. Nature Reviews Genetics 23:697−710

    doi: 10.1038/s41576-022-00514-4

    CrossRef   Google Scholar

    [35]

    Gallegos J. 2018. Alternative splicing plays a major role in plant response to cold temperatures. The Plant Cell 30:1378−79

    doi: 10.1105/tpc.18.00430

    CrossRef   Google Scholar

    [36]

    Seo PJ, Park MJ, Park CM. 2013. Alternative splicing of transcription factors in plant responses to low temperature stress: mechanisms and functions. Planta 237:1415−24

    doi: 10.1007/s00425-013-1882-4

    CrossRef   Google Scholar

    [37]

    John S, Olas JJ, Mueller-Roeber B. 2021. Regulation of alternative splicing in response to temperature variation in plants. Journal of Experimental Botany 72:6150−63

    doi: 10.1093/jxb/erab232

    CrossRef   Google Scholar

    [38]

    Li J, Yan X, Ahmad M, Yu W, Song Z, et al. 2021. Alternative splicing of the dormancy-associated MADS-box transcription factor gene PpDAM1 is associated with flower bud dormancy in 'Dangshansu' pear (Pyrus pyrifolia white pear group). Plant Physiology and Biochemistry 166:1096−108

    doi: 10.1016/j.plaphy.2021.07.017

    CrossRef   Google Scholar

    [39]

    Hao J, Gao Y, Xue J, Yang Y, Yin J, et al. 2022. Phytochemicals, pharmacological effects and molecular mechanisms of mulberry. Foods 11:1170

    doi: 10.3390/foods11081170

    CrossRef   Google Scholar

    [40]

    Hasegawa K, Tsuboi A. 1960. The effect of low temperature on the breaking of rest for winter bud in mulberry tree. The Journal of Sericultural Science of Japan 29:63−68

    doi: 10.11416/kontyushigen1930.29.63

    CrossRef   Google Scholar

    [41]

    Zhu Z, Yu C, Dong Z, Mo R, Zhang C, et al. 2024. Phylogeny and fungal community structures of Helotiales associated with sclerotial disease of mulberry fruits in China. Plant Disease 108:502−12

    doi: 10.1094/PDIS-02-23-0223-RE

    CrossRef   Google Scholar

    [42]

    Lü R, Zhao A, Yu J, Wang C, Liu C, et al. 2017. Biological and epidemiological characteristics of the pathogen of hypertrophy sorosis scleroteniosis, Ciboria shiraiana. Wei Sheng Wu Xue Bao 57:388−98

    doi: 10.13343/j.cnki.wsxb.20160258

    CrossRef   Google Scholar

    [43]

    Luo Y, Li H, Xiang Z, He N. 2018. Identification of Morus notabilis MADS-box genes and elucidation of the roles of MnMADS33 during endodormancy. Scientific Reports 8:5860

    doi: 10.1038/s41598-018-23985-0

    CrossRef   Google Scholar

    [44]

    Qi X, Shuai Q, Chen H, Fan L, Zeng Q, et al. 2014. Cloning and expression analyses of the anthocyanin biosynthetic genes in mulberry plants. Molecular Genetics and Genomics 289:783−93

    doi: 10.1007/s00438-014-0851-3

    CrossRef   Google Scholar

    [45]

    He N, Zhang C, Qi X, Zhao S, Tao Y, et al. 2013. Draft genome sequence of the mulberry tree Morus notabilis. Nature Communications 4:2445

    doi: 10.1038/ncomms3445

    CrossRef   Google Scholar

    [46]

    Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, et al. 2003. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Research 31:3497−500

    doi: 10.1093/nar/gkg500

    CrossRef   Google Scholar

    [47]

    Walton W. 1959. Results of test drilling and aquifer tests in Jerome, Lincoln, and Minidoka counties, Idaho. Report. pp. 58−108. doi:10.3133/ofr58108

    [48]

    El-Sese A, Mohamed A. 2003. Chilling, heat requirements and hormonal control in relation to bud dormancy in red roomy and thompson seedless grape cultivars, (Vitis vinifera L.). Assiut Journal of Agricultural Sciences 34:221−36

    Google Scholar

    [49]

    Marquardt S, Raitskin O, Wu Z, Liu F, Sun Q, et al. 2014. Functional consequences of splicing of the antisense transcript COOLAIR on FLC transcription. Molecular Cell 54:156−65

    doi: 10.1016/j.molcel.2014.03.026

    CrossRef   Google Scholar

    [50]

    Xiong F, Ren JJ, Yu Q, Wang YY, Lu CC, et al. 2019. AtU2AF65b functions in abscisic acid mediated flowering via regulating the precursor messenger RNA splicing of ABI5 and FLC in Arabidopsis. New Phytologist 223:277−92

    doi: 10.1111/nph.15756

    CrossRef   Google Scholar

    [51]

    Kang Y, Yang D, Kong L, Hou M, Meng Y, et al. 2017. CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Research 45:W12−W16

    doi: 10.1093/nar/gkx428

    CrossRef   Google Scholar

    [52]

    Wang W, Liu D, Chen D, Cheng Y, Zhang X, et al. 2019. MicroRNA414c affects salt tolerance of cotton by regulating reactive oxygen species metabolism under salinity stress. RNA Biology 16:362−75

    doi: 10.1080/15476286.2019.1574163

    CrossRef   Google Scholar

    [53]

    Zhu Z, An F, Feng Y, Li P, Xue L, et al. 2011. Derepression of ethylene-stabilized transcription factors (EIN3/EIL1) mediates jasmonate and ethylene signaling synergy in Arabidopsis. Proceedings of the National Academy of Sciences of the United States of America 108:12539−44

    doi: 10.1073/pnas.1103959108

    CrossRef   Google Scholar

    [54]

    Li T, Qi X, Zeng Q, Xiang Z, He N. 2014. MorusDB: a resource for mulberry genomics and genome biology. Database 2014:bau054

    doi: 10.1093/database/bau054

    CrossRef   Google Scholar

    [55]

    Horvath DP, Anderson JV, Chao WS, Foley ME. 2003. Knowing when to grow: signals regulating bud dormancy. Trends in Plant Science 8:534−40

    doi: 10.1016/j.tplants.2003.09.013

    CrossRef   Google Scholar

    [56]

    Amasino R. 2004. Vernalization, competence, and the epigenetic memory of winter. The Plant Cell 16:2553−59

    doi: 10.1105/tpc.104.161070

    CrossRef   Google Scholar

    [57]

    Purvis ON, Gregory FG. 1952. Studies in vernalisation XII. The reversibility by high temperature of the vernalised condition in petkus winter rye. Annals of Botany 16:1−21

    doi: 10.1093/oxfordjournals.aob.a083297

    CrossRef   Google Scholar

    [58]

    Erez A, Couvillon GA, Hendershott CH. 1979. The effect of cycle length on chilling negation by high temperatures in dormant peach leaf buds. Journal of the American Society for Horticultural Science 104:573−76

    doi: 10.21273/JASHS.104.4.573

    CrossRef   Google Scholar

    [59]

    Chouard P. 1960. Vernalization and its relations to dormancy. Annual Review of Plant Physiology 11:191−238

    doi: 10.1146/annurev.pp.11.060160.001203

    CrossRef   Google Scholar

    [60]

    Horvath D. 2009. Common mechanisms regulate flowering and dormancy. Plant Science 177:523−31

    doi: 10.1016/j.plantsci.2009.09.002

    CrossRef   Google Scholar

    [61]

    Brunner AM, Evans LM, Hsu CY, Sheng X. 2014. Vernalization and the chilling requirement to exit bud dormancy: shared or separate regulation? Frontiers in Plant Science 5:732

    doi: 10.3389/fpls.2014.00732

    CrossRef   Google Scholar

    [62]

    Song J, Irwin J, Dean C. 2013. Remembering the prolonged cold of winter. Current Biology 23:R807−R811

    doi: 10.1016/j.cub.2013.07.027

    CrossRef   Google Scholar

    [63]

    Chiang GCK, Barua D, Kramer EM, Amasino RM, Donohue K. 2009. Major flowering time gene, FLOWERING LOCUS C, regulates seed germination in Arabidopsis thaliana. Proceedings of the National Academy of Sciences of the United States of America 106:11661−66

    doi: 10.1073/pnas.090136710

    CrossRef   Google Scholar

    [64]

    Blair L, Auge G, Donohue K. 2017. Effect of FLOWERING LOCUS C on seed germination depends on dormancy. Functional Plant Biology 44:493−506

    doi: 10.1071/FP16368

    CrossRef   Google Scholar

    [65]

    Chen M, Penfield S. 2018. Feedback regulation of COOLAIR expression controls seed dormancy and flowering time. Science 360:1014−17

    doi: 10.1126/science.aar7361

    CrossRef   Google Scholar

    [66]

    Liu Y, Dreni L, Zhang H, Zhang X, Li N, et al. 2022. A tea plant (Camellia sinensis) FLOWERING LOCUS C-like gene, CsFLC1, is correlated to bud dormancy and triggers early flowering in Arabidopsis. International Journal of Molecular Sciences 23:15711

    doi: 10.3390/ijms232415711

    CrossRef   Google Scholar

    [67]

    Chen W, Tamada Y, Yamane H, Matsushita M, Osako Y, et al. 2022. H3K4me3 plays a key role in establishing permissive chromatin states during bud dormancy and bud break in apple. The Plant Journal 111:1015−31

    doi: 10.1111/tpj.15868

    CrossRef   Google Scholar

    [68]

    Calixto CPG, Tzioutziou NA, James AB, Hornyik C, Guo W, et al. 2019. Cold-dependent expression and alternative splicing of Arabidopsis long non-coding RNAs. Frontiers in Plant Science 10:235

    doi: 10.3389/fpls.2019.00235

    CrossRef   Google Scholar

    [69]

    Leviatan N, Alkan N, Leshkowitz D, Fluhr R. 2013. Genome-wide survey of cold stress regulated alternative splicing in Arabidopsis thaliana with tiling microarray. PLoS One 8:e66511

    doi: 10.1371/journal.pone.0066511

    CrossRef   Google Scholar

    [70]

    Liu JJ, Ekramoddoullah AKM. 2006. The family 10 of plant pathogenesis-related proteins: their structure, regulation, and function in response to biotic and abiotic stresses. Physiological and Molecular Plant Pathology 68:3−13

    doi: 10.1016/j.pmpp.2006.06.004

    CrossRef   Google Scholar

    [71]

    Lopes NDS, Santos AS, Silva de Novais DP, Pirovani CP, Micheli F. 2023. Pathogenesis-related protein 10 in resistance to biotic stress: progress in elucidating functions, regulation and modes of action. Frontiers in Plant Science 14:1193873

    doi: 10.3389/fpls.2023.1193873

    CrossRef   Google Scholar

    [72]

    Desouky AF, Ahmed AHH, Stützel H, Jacobsen HJ, Pao YC, et al. 2021. Enhanced abiotic stress tolerance of Vicia faba L. plants heterologously expressing the PR10a gene from potato. Plants 10:173

    doi: 10.3390/plants10010173

    CrossRef   Google Scholar

    [73]

    Ukaji N, Kuwabara C, Takezawa D, Arakawa K, Fujikawa S. 2004. Accumulation of pathogenesis-related (PR) 10/Bet v 1 protein homologues in mulberry (Morus bombycis Koidz.) tree during winter. Plant, Cell & Environment 27:1112−21

    doi: 10.1111/j.1365-3040.2004.01216.x

    CrossRef   Google Scholar

    [74]

    Pnueli L, Hallak-Herr E, Rozenberg M, Cohen M, Goloubinoff P, et al. 2002. Molecular and biochemical mechanisms associated with dormancy and drought tolerance in the desert legume Retama raetam. The Plant Journal 31:319−30

    doi: 10.1046/j.1365-313X.2002.01364.x

    CrossRef   Google Scholar

    [75]

    Ekramoddoullah AKM, Taylor D, Hawkins BJ. 1995. Characterization of a fall protein of sugar pine and detection of its homologue associated with frost hardiness of western white pine needles. Canadian Journal of Forest Research 25:1137−47

    doi: 10.1139/x95-126

    CrossRef   Google Scholar

    [76]

    Agurla S, Gahir S, Munemasa S, Murata Y, Raghavendra AS. 2018. Mechanism of stomatal closure in plants exposed to drought and cold stress. In Survival Strategies in Extreme Cold and Desiccation, vol 1081, eds Iwaya-Inoue M, Sakurai M, Uemura M. pp. 215−32. doi: 10.1007/978-981-13-1244-1_12

    [77]

    Yu J, Cang J, Lu Q, Fan B, Xu Q, et al. 2020. ABA enhanced cold tolerance of wheat 'dn1' via increasing ROS scavenging system. Plant Signaling & Behavior 15:1780403

    doi: 10.1080/15592324.2020.1780403

    CrossRef   Google Scholar

  • Cite this article

    Luo Y, Liu H, Han Y, Li W, Wei W, et al. 2024. Alternative splicing of the FLOWERING LOCUS C-like gene MaMADS33 is associated with endodormancy in mulberry. Forestry Research 4: e029 doi: 10.48130/forres-0024-0027
    Luo Y, Liu H, Han Y, Li W, Wei W, et al. 2024. Alternative splicing of the FLOWERING LOCUS C-like gene MaMADS33 is associated with endodormancy in mulberry. Forestry Research 4: e029 doi: 10.48130/forres-0024-0027

Figures(7)  /  Tables(1)

Article Metrics

Article views(954) PDF downloads(124)

ARTICLE   Open Access    

Alternative splicing of the FLOWERING LOCUS C-like gene MaMADS33 is associated with endodormancy in mulberry

Forestry Research  4 Article number: e029  (2024)  |  Cite this article

Abstract: Alternative splicing (AS) is an important post-transcriptional process that generates multiple mRNA isoforms. FLOWERING LOCUS C (FLC) is a pivotal gene in both the vernalization and autonomous pathways of flowering plants, and MaMADS33 is one of the FLC homologs in white mulberry (Morus alba). Recent studies have revealed that MaMADS33 is involved in endodormancy, but the underlying molecular mechanism remains to be characterized. Here, a comparison of MaMADS33 expression among three mulberry cultivars with different degrees of dormancy revealed a positive association between MaMADS33 expression and dormancy. Further 3' and 5' rapid amplification of cDNA ends (RACE) analyses led to identifying four MaMADS33 isoforms derived from AS and designated MaMADS33-AS14. Analysis of their coding potential revealed that MaMADS33-AS1 was a long non-coding RNA. Expression profiling and splicing-efficiency analyses showed that cold stress during endodormancy induced AS of MaMADS33, resulting in a predominance of truncated isoforms, especially MaMADS33-AS1. MaMADS33-AS2 expression was upregulated during both endodormancy and ecodormancy, whereas MaMADS33-AS3 and MaMADS33-AS4 were endodormancy-associated isoforms that were upregulated during endodormancy and then downregulated during ecodormancy. MaMADS33-AS4 was used as bait for a yeast two-hybrid screen because its gene expression was higher than that of MaMADS33-AS3, and mulberry winter-accumulating 18 kDa protein (MaWAP18) was identified as an MaMADS33-AS4 interaction partner. The interaction between MaWAP18 and MaMADS33-AS4 was confirmed by a bimolecular fluorescence complementation assay. These findings offer insight into the role of FLC homologs in the endodormancy of woody plants.

    • Bud dormancy of perennials is a survival strategy for adaptation to the strikingly different environmental conditions of summer and winter[1]. Based on the specific conditions that trigger dormancy, it can be divided into endodormancy, aradormancy, and ecodormancy[2]. Endodormancy is the deep dormancy of plants that is established through low temperatures or short photoperiods in autumn and winter, and the release of endodormancy requires a specific period of chilling. Aradormancy refers to the dormancy of lateral buds that results from apical dominance, whereas ecodormancy refers to a temporary halt in growth caused by adverse environmental conditions[2]. Endodormancy is intimately linked to fruit set and quality, as a lack of sufficient chilling exposure can lead to abortion of reproductive whorls, low bud burst, and non-uniform bloom[35]. Given the ongoing trend of global warming, the mechanism of endodormancy has garnered increasing attention[69].

      Low temperatures and short photoperiods are two crucial factors required for the induction of endodormancy in plants[10], but the underlying mechanisms are not fully understood. For a long time, research on endodormancy in woody plants has concentrated primarily on external phenotypic manifestations and biochemical components. The first dormancy-associated MADS-box (DAM) genes were identified through localization and cloning of the evergrowing locus in peach[11,12]. This discovery significantly expedited the analysis of the molecular mechanisms underlying endodormancy. The DAM genes are homologs of the flowering genes SHORT VEGETATIVE PHASE (SVP) and AGAMOUS LIKE24 (AGL24) in Arabidopsis[13,14]. In recent years, the involvement of the DAM genes in the endodormancy process of perennial plants has been confirmed at the nucleic acid level[1520]. Furthermore, other flowering genes, including FLOWERING LOCUS T/TERMINAL FLOWER 1 (FT/TFL1)[21,22] and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1)[23], as well as regulatory genes associated with plant hormones such as ABA and GA[24,25], also play pivotal roles in regulating the dormancy of woody plants.

      The MADS-box family gene FLOWERING LOCUS C (FLC) serves as the central integrator of the vernalization pathway in Arabidopsis[26]. Numerous studies have demonstrated that FLC also plays an important role in the regulation of endodormancy in perennial plants. In Poncirus trifoliata, expression of PtFLC is upregulated during autumn and winter but downregulated in spring, coinciding with the release from dormancy[27]. Similarly, two FLC homologs in apple (Malus domestica), MdMADS135, and MdMADS136, exhibit expression patterns comparable to those of PtFLC[28]. Notably, the expression of MdoFLC is strongly correlated with endodormancy, and its chromosome position is linked to the timing of bud break[29,30]. Ectopic expression in Arabidopsis revealed that MdoFLC functions as a growth inhibitory factor during endodormancy[31]. Transcriptome studies in buds of kiwifruit (Actinidia chinensis) during endodormancy showed that AcFLCL expression peaks during the accumulation of chilling hours, and overexpression of AcFLCL promotes budbreak[32]. However, the precise role and mechanism by which FLC homologs regulate endodormancy in perennial plants remain to be clarified.

      Alternative splicing (AS) is a post-transcriptional process that generates multiple mRNA isoforms through seven major types of splicing events[33,34]. AS plays major roles in responses to environmental cues, including cold temperatures[3537]. Research in pear (Pyrus pyrifolia) showed that AS of DAM genes was involved in dormancy[38]. Similarly, two sense AS isoforms of AcFLCL have been identified in kiwifruit, with the longer isoform induced by chilling and the shorter isoform detected before, and after chilling treatment[32].

      White mulberry (Morus alba) is a perennial woody plant whose fruit is rich in anthocyanins and has numerous medicinal and health benefits[39]. A chilling period is essential for breaking endodormancy in mulberry[40]. Mulberry sclerotial disease, a widespread fungal infection, poses a significant threat to the production of mulberry[41]; it typically erupts from January to April, peaking in mid-March[42]. Staggering the flowering times of mulberry trees to avoid this infection period could be an effective strategy for disease control.

      Our previous study demonstrated that the mulberry FLC homolog MnMADS33 is upregulated by cold treatment under field and artificially controlled conditions, and its downregulation is associated with the release from endodormancy[43]. In this study, the role of MaMADS33 in mulberry endodormancy was further investigated. Comparative expression analysis revealed a positive association between MaMADS33 expression and mulberry endodormancy. Using 3' RACE and 5' RACE techniques, we identified four splicing isoforms of MaMADS33 during endodormancy. Cloning and expression analyses showed that cold stress during endodormancy induced AS of the MaMADS33 gene, resulting in a predominance of truncated isoforms, especially the long noncoding RNA MaMADS33-AS1. Expression profiling indicated that MaMADS33-AS4 expression was positively associated with mulberry endodormancy. Furthermore, winter-accumulating 18 kDa protein (MaWAP18) was identified as a MaMADS33-AS4 interaction partner through yeast two-hybrid screening and bimolecular fluorescence complementation (BiFC) assays. Collectively, the present findings suggest that MaMADS33 plays a crucial role as an endodormancy regulator in mulberry.

    • Flower buds were collected from the diploid M. alba varieties 'Jinqiang63' (JQ63), 'Lunjiao109' (LJ109), and 'Zhenzhubai' (ZZB) growing at the National Mulberry Breed Improvement Center of the Southwest University of China (29°49'18"N, 106°25'45"E). These distinct mulberry varieties hail from various geographic regions within China. LJ109 was sourced from Foshan City in Guangdong Province, JQ63 is a local cultivar from Dazu District in Chongqing, and ZZB originated in Linqing City in Shandong Province. The flower buds of LJ109, JQ63, and ZZB break in January, early February, and late February, respectively, demonstrating the progressive increase in dormancy duration of these three varieties.

      For comparative expression analysis, flower buds of LJ109, JQ63, and ZZB were collected at six time points from autumn through spring: 16 October 2017, 6 November 2017, 27 November 2017, 21 December 2017, 12 January 2018, and 25 February 2018. For other analyses, flower buds of JQ63 were collected at five time points: 23 October 2020, 26 November 2020, 25 December 2020, 23 January 2021, and 24 February 2021. At each time point, buds were harvested from three replicate groups, each consisting of three to four independent plants of the same genotype. Before molecular analyses, the buds were processed by removing their scales, immediately freezing them in liquid nitrogen, and storing them at −80°C. Daily maximum temperature data were retrieved from a reliable online source (http://lishi.tianqi.com/).

    • Flower buds frozen at −80 °C were ground in liquid nitrogen in an RNA-free environment. Total RNA was isolated using the TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. First-strand complementary DNA (cDNA) was synthesized using a PrimeScript RT reagent Kit with gDNA Eraser (Perfect Real Time) (Takara, Dalian, China), followed by second-strand cDNA synthesis. qRT–PCR was performed in a volume of 20 μL using the ABI7500 Fast Real-time PCR system (Applied Biosystems, Foster City, CA, USA) and SYBR Green I Master Mix (Takara, Dalian, China). The mulberry MaRPL15 gene (Morus024083) served as the internal control[44]. Relative gene expression was calculated using the formula 2−[Ct(target gene) − Ct(control gene)]. RT–PCR was performed using primers specific to the open reading frame (ORF) of each AS isoform. PCR products were run on 1% agarose gels and verified by sequencing. The primers used for qRT–PCR and RT–PCR are listed in Supplemental Table S1.

    • MaMADS33 transcripts were amplified by 3' RACE and 5' RACE using the 3'-Full RACE Core Set with PrimeScript RTase (Takara, Dalian, China) and the FirstChoice RLM-RACE Kit (Invitrogen Life Technologies, Carlsbad, CA, USA), respectively, according to the manufacturer's instructions. The template used for these amplifications was a cDNA pool derived from flower buds of JQ63 collected on 23 October 2020, 26 November 2020, 25 December 2020, 23 January 2021, and 24 February 2021. The primers are listed in Supplemental Table S1. The gene structure of MaMADS33 was visualized using the Gene Structure Display Server 2.0 (GSDS) with the Morus notabilis genome[45]. A multiple alignment of the MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4 protein sequences were obtained using ClustalW[46].

    • Mulberry endodormancy was defined as described by Walton[47,48]. In brief, buds were considered to be endodormant if fewer than 50% of them reached the green tip stage after exposure to suitable environmental conditions. Shoot cuttings with two buds were immersed in distilled water at room temperature (22–25 °C), and the water was changed every 2 d. The bud-break percentage was recorded after four weeks.

    • Splicing efficiency was measured as described previously[49]. Unspliced primers were designed to span intron–exon junctions, and spliced primers were designed to cross exon–exon junctions. All PCR products were verified by sequencing. Only primers that produced specific and single bands were used for further analysis. qRT–PCR was then performed to quantify the levels of spliced and unspliced transcripts of MaMADS33. MaRPL15 was used as an internal control. The splicing efficiency of individual introns was calculated as the spliced/unspliced ratio[50]. Data are presented as the means of three biological replicates, and specific primers used for this analysis are listed in Supplemental Table S1.

    • Coding Potential Calculator 2 (CPC2) was used to analyze the coding potential of various MaMADS33 isoforms[51]. The ORFs of MaMADS33-AS1 and MaWAP18 without termination codons were cloned into the pZYGC-GFP vector using the pEASY-Basic Seamless Cloning and Assembly Kit (TransGen Biotech, Beijing, China) to produce MaMADS33-AS1-GFP and MaWAP18-GFP fusion constructs driven by the 35S promoter. Primers used are presented in Supplemental Table S1. Plasmids containing the target fragments were then transformed into Escherichia coli competent Trans-T1 cells (TransGen Biotech, Beijing, China). The accuracy of the cloned sequences was verified by sequencing. pZYGC-GFP was used as a positive control plasmid. The correctly sequenced plasmids were transformed into Agrobacterium tumefaciens strain GV3101 by the freeze-thaw method. GV3101 cells harboring various expression plasmids were then individually transformed into fully expanded leaves of six-week-old Nicotiana benthamiana plants as described previously[52]. The leaf epidermal cell layers were observed at 2–3 d post-inoculation with an FV1200 laser scanning confocal microscope (Olympus, Tokyo, Japan). At least three independent biological replicates were performed.

    • Y2H screening was performed with assistance from Oebiotech (Shanghai, China). A normalized cDNA library for Y2H screening was generated by pooling RNA extracted from dormant flower buds, mature mulberry fruits, and flower buds at the floral formation stage collected in December 2016, May 2017, and June 2017, respectively. The Y2H library screen was performed using pGBKT7-MaMADS33-AS4 as the bait in the yeast strain AH109. To further validate the MaMADS33-AS4–MaWAP18 interaction identified in the Y2H screen, a Y2H assay was performed as described previously[53]. MaMADS33-AS4-BD served as the bait construct, and MaWAP18-AD served as the prey. pGBKT7-Lam and pGADT7-T were used as negative controls and pGBKT7-53 and pGADT7-T as positive controls. Bait and prey constructs were co-transformed into yeast strain AH109. Transformants were selected on SD-Leu-Trp (SD-LW) medium, gradient-diluted, and plated on SD-Leu-Trp-His-Ade (QX) medium. Pictures were taken after 4 d of incubation at 30 °C and β-galactosidase activity was quantified using the Yeast β-Galactosidase Assay Kit (Thermo Fisher Scientific, MA, USA) according to the manufacturer's instructions. All Y2H experiments were performed with three independent biological replicates. Primers used for this analysis are presented in Supplemental Table S1. Gene annotations were obtained from the Mulberry Genome Database, MorusDB (https://morus.biodb.org/index)[54].

    • BiFC assays were performed to confirm the protein–protein interaction identified above. The ORFs of MaMADS33-AS4 and MaWAP18 without stop codons were obtained by PCR amplification with gene-specific primers (Supplemental Table S1). MaMADS33-AS4 and MaWAP18 were then independently fused with the N-terminal and C-terminal halves of the yellow fluorescent protein (YFP). A BiFC assay was performed in the leaves of six-week-old N. benthamiana plants using the A. tumefaciens transformation and injection protocols described above. Leaf epidermal cell layers were observed at 2–3 d post-inoculation with an FV1200 laser scanning confocal microscope (Olympus, Tokyo, Japan). Three independent biological replicates were performed.

    • Differences between samples were analyzed by one-way analysis of variance (ANOVA) in SPSS version 19 (IBM Corp., Armonk, NY, USA). Means were compared using Duncan's multiple range test (p < 0.05).

    • The relationship between MaMADS33 and dormancy was investigated by comparing MaMADS33 expression in flower buds from three mulberry varieties with different dormancy durations (LJ109, JQ63, and ZZB) using primers targeting the conserved amino acid sequence of the MADS-box domain (primer I, Supplemental Table S1). MaMADS33 expression peaked at different times in the three varieties: early November in LJ109, late November in JQ63, and January in ZZB (Fig. 1). Notably, the relative timing of these peaks was consistent with the relative timing of bud break in the three varieties (earliest in LJ109 and latest in ZZB). Expression levels of MaMADS33 also differed among the mulberry varieties. In October, November, and December, MaMADS33 expression was highest in ZZB, moderate in JQ63, and lowest in LJ109, again consistent with their relative durations of dormancy. In January, MaMADS33 expression remained high in ZZB but was low in LJ109 and JQ63. By February, no unopened buds remained on LJ109, and MaMADS33 expression was low in both JQ63 and ZZB.

      Figure 1. 

      MaMADS33 expression is positively associated with dormancy in mulberry. MaMADS33 expression in the mulberry varieties 'Lunjiao109' (LJ109), 'Jinqiang63' (JQ63), and 'Zhenzhubai' (ZZB) was measured throughout dormancy by qRT–PCR using primers targeting the conserved sequence of the MADS-box domain. MaRPL15 was the internal control gene for normalization of the expression data (n = 3; mean ± measurement range). Significant differences are indicated by different lowercase letters (ANOVA and Duncan's multiple range test; p < 0.05).

    • To understand the possible roles of AS during dormancy, MaMADS33 transcripts were next examined in flower buds of JQ63, which exhibited an intermediate dormancy duration. Using 5' RACE and 3' RACE with specific primers targeting the conserved MADS-box domain sequence (Fig. 2a & b), four distinct complete ORFs of MaMADS33: MaMADS33-AS1, MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4, with lengths of 327, 555, 633, and 675 bp, respectively were identified. MaMADS33-AS1 and MaMADS33-AS2 were truncated isoforms that arose from alternative last exons; MaMADS33-AS1 terminated in the second intron and MaMADS33-AS2 in the sixth intron (Fig. 2c). MaMADS33-AS3 was derived from the skipping of the seventh exon. Finally, MaMADS33-AS4 represented the complete sequence, including all eight exons.

      Figure 2. 

      Alternative splicing (AS) of mulberry MaMADS33 from endodormancy through ecodormancy. (a) 5' RACE, and (b) 3' RACE methods were used to amplify MaMADS33 transcripts in cDNA pools from JQ63 flower buds collected on five dates from October 2020 through February 2021. Primers were designed to target the conserved sequence of the MADS-box domain. Bands corresponding to MaMADS33-AS1, MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4 are indicated by arrows. (c) Schematic of MaMADS33 AS isoforms. Exons are represented by green boxes and introns by lines.

    • To investigate the functions of the four MaMADS33 AS isoforms, their coding potential was examined and it was found that all contained a complete ORF. However, the coding potential of MaMADS33-AS1 was 0.23, whereas those of MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4 were all greater than 0.95. Based on their coding potentials, MaMADS33-AS1 was classified as a noncoding sequence, and the other three isoforms were classified as coding sequences (Table 1). The low coding potential of MaMADS33-AS1 was further validated by transient expression of a MaMADS33-AS1-GFP construct in N. benthamiana leaves, which resulted in no GFP signal (Fig. 3a). These findings confirmed that MaMADS33-AS1 was a long noncoding RNA.

      Table 1.  Coding potential of the four MaMADS33 AS isoforms were predicted using Coding Potential Calculator 2 (CPC2).

      ID Label Coding probability Peptide length (aa) ORF
      integrity
      MaMADS33-AS1 noncoding 0.232883 109 complete
      MaMADS33-AS2 coding 0.977475 185 complete
      MaMADS33-AS3 coding 0.969709 211 complete
      MaMADS33-AS4 coding 0.985029 225 complete

      Figure 3. 

      Assessment of MaMADS33-AS1 coding potential and alignment of amino acid sequences. (a) The coding potential of MaMADS33-AS1 was examined by transient expression in Nicotiana benthamiana leaves. Scale bar = 20 μm. (b) Amino acid sequences of MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4. The MADS-box and K-domains are indicated by red and black boxes, respectively. The I-domain and C-terminal domain are marked above the sequences.

      MaMADS33-AS2, MaMADS3-AS3, and MaMADS03-AS4 were composed of 185, 211, and 225 amino acids, respectively (Table 1). Alignment of their amino acid sequences revealed that MaMADS33-AS2 lacked the C-terminal domain, MaMADS33-AS3 lacked 14 amino acids, mainly in the C-terminal domain, and MaMADS33-AS4 contained all the characteristic domains of MADS family proteins, including the MADS-box, I-domain, K-domain, and C-terminal domain (Fig. 3b).

    • To investigate the potential functions of MaMADS33 in endodormancy, the endodormancy stages of JQ63 flower buds were first characterized(Fig. 4a). From October to December, temperatures gradually decreased (Fig. 4b), and the bud-break percentage remained below 50%, indicating that the flower buds were in an endodormant state. Next, the expression of the four MaMADS33 AS isoforms in bud tissues were measured by RT–PCR with specific primers during endodormancy and ecodormancy (Fig. 4c). Transcript levels of MaMADS33-AS1 remained low in October, increased significantly from November to December, and then decreased markedly in February. By contrast, expression of MaMADS33-AS2 increased gradually after October and remained high throughout ecodormancy. Expression of MaMASD33-AS3 and MaMADS33-AS4 was relatively high throughout the endodormancy period, peaking around November, then decreased significantly during ecodormancy. To further analyze the expression of these isoforms during endodormancy, nucleic acid bands of MaMASD33-AS3 and MaMADS33-AS4 were counted by cloning and sequencing. A total of 44 clones were sequenced, 16 (36.36%) of which were MaMASD33-AS3 and 28 (63.64%) of which were MaMADS33-AS4 (Supplemental Table S2). MaMADS33-AS4 was therefore the predominant isoform in the mixed nucleic acid bands during endodormancy. Expression of these two long protein-coding mRNAs, particularly MaMADS33-AS4, thus exhibited a positive association with endodormancy in mulberry flower buds.

      Figure 4. 

      Expression profiles of four MaMADS33 AS isoforms during dormancy. Dormancy stages of mulberry flower buds. Flower buds of JQ63 were collected from October 2020 through February 2021, and the (a) bud break percentage (n = 3, mean ± measurement range), and (b) temperature were recorded. (c) Expression profiles of four MaMADS33 AS isoforms in flower buds measured by RT–PCR with MaRPL15 as the reference gene. PCR products were separated on 1% agarose gels.

    • To gain a deeper understanding of the relationship between dormancy and MaMADS33 expression patterns, the splicing efficiency of intron 2 was analyzed using qRT–PCR data. A schematic representation of the MaMADS33 gene and the primers used for amplification is provided in Fig. 5a. The amplification products generated by primer I, which targeted the conserved MADS-box domain, corresponded to all MaMADS33 transcripts (Fig. 5b). Overall transcription of MaMADS33 was high during endodormancy, then decreased significantly during ecodormancy. The amplification products of primer II, which specifically targeted transcripts from which intron 2 had been spliced, corresponded to MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4. Notably, the abundance of these spliced transcripts decreased after October (Fig. 5c). The amplification products generated by primer III represented unspliced pre-mRNA containing intron 2 (Fig. 5d). To quantify the splicing efficiency of intron 2 in the isoforms MaMADS33-AS2, MaMADS33-AS3, and MaMADS33-AS4, the spliced/unspliced ratio was calculated based on qRT–PCR data obtained with primer II vs primer III. The results revealed a significant decrease in splicing efficiency during endodormancy, followed by a gradual increase during ecodormancy (Fig. 5e).

      Figure 5. 

      Splicing efficiency of MaMADS33 intron 2. (a) Schematic of the MaMADS33 gene. The locations of primers I, II, and III are indicated with arrows. (b) Primer I, located in the conserved MADS-box domain, targeted all spliced and unspliced isoforms of MaMADS33. (c) Primer II targeted mRNA from which intron 2 had been spliced, corresponding to transcripts of MaMADS-AS2, MaMADS-AS3, and MaMADS-AS4. (d) Primer III targeted unspliced pre-mRNA for intron 2. (e) Splicing efficiency was calculated as the spliced/unspliced ratio. Relative expression was measured by qRT–PCR using MaRPL15 as a reference gene (n = 3; mean ± measurement range). Flower buds of JQ63 were collected from October 2020 through February 2021. Significant differences are indicated by different lowercase letters (ANOVA and Duncan's multiple range test; p < 0.05).

    • To further investigate the potential mechanism by which MaMADS33 functions in endodormancy, a Y2H screen of pooled RNA from dormant flower buds, mature fruits, and flower buds at the floral formation stage was performed. All interacting proteins identified are listed in Supplemental Table S3; among them was the winter-accumulating protein MaWAP18. We next performed a Y2H and β-galactosidase activity assay to further examine the interaction between MaMADS33-AS4 and MaWAP18 (Fig. 6a). The interaction between MaMADS33-AS4 and MaWAP18 was strong when MaMADS33-AS4 served as the bait and MaWAP18 as the prey. To confirm this interaction in vivo, a BiFC assay in N. benthamiana was performed. Transient co-expression of MaMADS33-AS4-YFPN and MaWAP18-YFPC in tobacco leaves resulted in a distinct yellow fluorescence signal in the guard cells, conclusively demonstrating the interaction between MaMADS33-AS4 and MaWAP18 in planta (Fig. 6b). To determine whether the expression of MaWAP18 was also associated with endodormancy, its expression profile was analyzed through time. Expression of MaWAP18 was high during endodormancy and low during ecodormancy (Fig. 7a). Subcellular localization analysis suggested that MaWAP18 was localized primarily in the cell membrane and nucleus, although it was also detected in the guard cells of N. benthamiana leaves (Fig. 7b).

      Figure 6. 

      Interaction between MaMADS33-AS4 and MaWAP18. (a) Interaction of MaMADS33-AS4 and MaWAP18 in yeast. pGBKT7-p53 was mated with pGADT7-T as a positive control, and pGBKT7-Lam was mated with pGADT7-T as a negative control. Yeast was diluted 1, 10, and 100 fold before plating onto quadruple dropout (QDO) medium. Corresponding measurements of β-galactosidase activity are shown. Three independent experiments were performed with similar results. (b) BiFC assay in 4-week-old Agrobacterium-infiltrated N. benthamiana leaves. MaMADS33-AS4 and MaWAP18 were independently fused to the N-terminal and C-terminal halves of yellow fluorescent protein (YFP), respectively. Images of YFP fluorescence were obtained using a confocal microscope. Scale bar = 10 μm. Three independent experiments were performed with similar results.

      Figure 7. 

      Expression and subcellular localization analyses of MaWAP18. (a) Relative expression was measured by qRT–PCR using MaRPL15 as the reference gene (n = 3; mean ± measurement range). Flower buds of JQ63 were collected from October 2020 through February 2021. (b) Subcellular localization of MaWAP18. Images of 4-week-old Agrobacterium-infiltrated N. benthamiana leaves expressing the MaWAP18-GFP fusion protein driven by the CaMV 35S promoter were obtained under green fluorescence, merged light, and visible light. 35S:GFP was used as a positive control. Scale bar = 10 μm. Significant differences are indicated by different lowercase letters (ANOVA and Duncan's multiple range test; p < 0.05).

    • Endodormancy of woody plants and vernalization of annual plants involve similar biological processes induced by cold or light, including growth cessation, slowed development or stagnation, and release, enabling plants to survive adverse winter conditions[55,56]. The duration of effective cold exposure during winter has a direct effect on the release from endodormancy and the vernalization process[57,58]. Some researchers have reported that endodormancy and vernalization share common mechanisms[5961]. During the vernalization process, the regulation of the FLC gene can be divided into several stages. Before vernalization, FLC is activated to establish the conditions necessary for vernalization. As vernalization proceeds, expression of the long noncoding RNAs COOLAIR and COLDAIR suppresses FLC expression. Changes in chromosome conformation also contribute to the decrease in FLC expression. Finally, after vernalization and before the floral transition, epigenetic modifications and high expression of VERNALIZATION INSENSITIVE 3 (VIN3) maintain FLC in a silenced state[26,62].

      FLC is also involved in Arabidopsis seed dormancy through mechanisms similar to those of vernalization[6365], and the MdoFLC gene in apple co-localizes with a bud-break QTL[29]. Previous studies in perennial plants have reported upregulation of FLC genes during endodormancy and downregulation during ecodormancy; examples include PtFLC in trifoliate orange[27], MdoFLC in apple[29], and CsFLC1 in Camellia sinensis[66]. Recently, similar expression patterns were documented for MnMADS33 (i.e., MaMADS33-AS4) in mulberry using primers that targeted conserved domain sequences[43]. However, previous expression profiles were based solely on quantitative detection using single primers and may therefore not accurately reflect gene expression, especially for genes with AS. Therefore, in this study, MaMADS33 expression in mulberry was systematically analyzed using both qRT–PCR and RT–PCR with multiple primers. Combined with measurements of mulberry endodormancy, these expression data indicated that MaMADS33 expression was positively associated with endodormancy of flower buds, and this was particularly true for the longer MaMADS33-AS3 and MaMADS33-AS4 isoforms. The long noncoding RNA MaMADS33-AS1 was also identified, which accumulated during the late stage of endodormancy and in the early stage of ecodormancy. These findings further support the proposed role of FLC homologs as important regulators of endodormancy.

    • A brief period of chilling in autumn triggers the onset of bud dormancy in apple, and a longer duration of chilling ultimately results in dormancy release[67]. Low temperatures are known to trigger extensive and rapid alterations in the RNA isoforms produced in response to environmental fluctuations in Arabidopsis[68]. For the mulberry variety JQ63, temperatures gradually decreased during endodormancy and rose during ecodormancy (Fig. 4a), and the splicing efficiency of MaMADS33 intron 2 increased with increasing temperature (Figs 4b & 5e). Previous studies of Arabidopsis have shown that the majority of cold-regulated AS events introduce premature termination codons into transcripts[69]. Consistent with this observation, MaMADS33-AS1 and MaMADS33-AS2 were truncated versions of MaMADS33 that arose from alternative last exons. We therefore suggest that chilling temperatures induce AS of MaMADS33, resulting in the production of shorter isoforms, particularly MaMADS33-AS1. However, further study will be required to determine the role and mechanism by which MaMADS33-AS1 participates in dormancy.

    • Pathogenesis-related protein 10 (PR10) plays a crucial role in plant responses to biotic and abiotic stresses[7072]. MaWAP18, a PR10 homolog accumulates during winter in mulberry (Morus bombycis Koidz.)[73], and purified WAP18 from M. bombycis exhibits cryoprotective activity towards lactate dehydrogenase in vitro[73], suggesting that it has a role in freezing tolerance. A similar pattern of PR10-like protein accumulation during dormancy was reported in Retama raetam[74], and cold acclimation triggered the accumulation of a PR10 homolog in white pine[75]. It was also found that the expression of MaWAP18 was high during endodormancy (Fig. 7a). The subcellular localization assay demonstrated that MaWAP18 was localized to the cell membrane, nucleus, and guard cells in N. benthamiana leaves, displaying a localization pattern similar to that observed for MaMADS33-AS4[43]. Notably, an interaction between MaMADS33-AS4 and MaWAP18 was observed in stomata (Fig. 6b). Low temperatures are known to promote ABA accumulation, leading to stomatal closure[76], and studies have shown that application of exogenous ABA can improve plant cold tolerance[77]. In our previous report, MaMADS33-AS4 expression was upregulated by both exogenous ABA and cold treatment[43]. Taken together, these findings suggest that MaMADS33-AS4, induced by ABA, may interact with MaWAP18 at the stomata of flower buds, thereby contributing to cold tolerance, although this possibility remains to be tested.

    • Four isoforms of MaMADS33 resulting from AS were identified. Expression profiling revealed that the overall abundance of MaMADS33 transcripts increased during endodormancy. Expression levels of MaMADS33-AS3 and MaMADS33-AS4 showed a positive association with endodormancy in flower buds. The splicing efficiency of MaMADS33 intron 2 decreased during endodormancy, leading to the accumulation of the truncated long noncoding RNA MaMADS33-AS1. Furthermore, MaWAP18, which has been suggested to play a role in the acquisition of freezing tolerance, was identified as an interaction partner of MaMADS33-AS4. These findings shed light on the molecular mechanisms of endodormancy in woody plants and contribute to a better understanding of the roles of FLC homologs in endodormancy.

    • The authors confirm contribution to the paper as follows: study conception and design: He N, Luo Y; data collection: Luo Y, Liu H, Han Y, Li W, Wei W; analysis and interpretation of results: He N, Luo Y; draft manuscript preparation: He N, Luo Y, Han Y. All authors reviewed the results and approved the final version of the manuscript.

    • The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

      • This work was funded by the National Natural Science Foundation of China (Grant No. 32001328) and the National Key Research and Development Program of China (Grant No. 2022YFD1201602). We thank Xianzhi Meng, Xin Ding, Maoshuang Ran, and Yan Zhu for assistance with laser scanning confocal microscopy measurements.

      • The authors declare that they have no conflict of interest.

      • Copyright: © 2024 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (7)  Table (1) References (77)
  • About this article
    Cite this article
    Luo Y, Liu H, Han Y, Li W, Wei W, et al. 2024. Alternative splicing of the FLOWERING LOCUS C-like gene MaMADS33 is associated with endodormancy in mulberry. Forestry Research 4: e029 doi: 10.48130/forres-0024-0027
    Luo Y, Liu H, Han Y, Li W, Wei W, et al. 2024. Alternative splicing of the FLOWERING LOCUS C-like gene MaMADS33 is associated with endodormancy in mulberry. Forestry Research 4: e029 doi: 10.48130/forres-0024-0027

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return