Time-course analysis and transcriptomic identification of key response strategies of Nelumbo nucifera to complete submergence

Abstract

Water submergence is an environmental stress with detrimental effects on plant growth and survival. As a wetland plant species, lotus (Nelumbo nucifera) is widely cultivated in flood-prone lowlands throughout Asian countries, but little is known about its endurance and acclimation mechanisms to complete submergence. Here, we performed a time-course submergence experiment and an RNA- sequencing transcriptome analysis of the two lotus varieties “Qiuxing” and “China Antique”. Both varieties showed low submergence tolerance, with a median lethal time of approximately 10 days. Differentially expressed gene (DEG) analysis and weighted gene co- expression network analysis (WGCNA) identified a number of key genes putatively involved in lotus submergence responses. Lotus plants under complete submergence developed thinned leaves and elongated petioles containing a high density of aerenchyma. All four lotus submergence-responsive ERF-VII genes and gene sets corresponding to the low oxygen “escape” strategy (LOES) were elevated. In addition, a number of lotus innate immunity genes were rapidly induced by submergence, probably to confer resistance to possible pathogen infections. Our data also reveal the probable involvement of jasmonic acid in the modulation of lotus submergence responses, although to a lesser extent than the gaseous hormone ethylene. These results suggest that lotus plants primarily use the LOES strategy to cope with complex submergence-induced stresses, and they will be valuable for understanding the molecular basis underlying plant submergence acclimation.

References

[1] Food and Agriculture Organization of the United Nations. The impact of disasters and crisis on agriculture and food security . Rome: FAO; 2017: 168.

[2] Hattori Y, Nagai K, Ashikari M . Rice growth adapting to Deepwa- ter. Curr Opin Plant Biol. 2011; 14: 100-5.

[3] Fukao T, Barrera-Figueroa BE, Juntawong P et al. Submer- gence and waterlogging stress in plants: a review highlighting research opportunities and understudied aspects. Front Plant Sci. 2019; 10.

[4] Gibbs J, Greenway H . Review: mechanisms of anoxia tolerance in plants. I. Growth, survival and anaerobic catabolism. Funct Plant Biol. 2003; 30: 1-47.

[5] Fukao T, Bailey-Serres J . Plant responses to hypoxia-is survival a balancing act? Trends Plant Sci. 2004; 9: 449-56.

[6] Vervuren PJA, Bolm CWPM et al. Extreme flooding events on the Rhine and the survival and distribution of riparian plant species. J Ecol. 2003; 91: 135-46.

[7] Bailey-Serres J, Voesenek LACJ . Flooding stress: acclimations and genetic diversity. Annu Rev Plant Biol. 2008; 59: 313-39.

[8] Greenway H, Setter TL. Is there anaerobic metabolism in sub- merged rice plants? a view point. In: Singh VP et al. , eds. Phys- iology of Stress Tolerance in Rice: Proceedings of the International Conference on Stress Physiology of Rice.Faizabad, U.P., India: IRRI, 1996, 11-30.

[9] Ram PC, Singh AK, Singh BB et al. Environmental characteriza- tion of floodwater in eastern India: relevance to submergence tolerance of lowland rice. Exp Agric. 1999; 35: 141-52.

[10] Crawford RMM . Oxygen availability as an ecological limit to plant distribution. Adv Ecol Res. 1992; 23: 93-185.

[11] Fukao T, Yeung E, Bailey-Serres J . The submergence tolerance regulator SUB1A mediates crosstalk between submergence and drought tolerance in rice. Plant Cell. 2011; 23: 412-27.

[12] Xu K, Xu X, Fukao T et al. Sub1A is an ethylene-response-factor- like gene that confers submergence tolerance to rice. Nature. 2006; 442: 705-8.

[13] Singh N, Dang TTM, Vergara GV et al. Molecular marker survey and expression analyses of the rice submergence-tolerance gene SUB1A . Theor Appl Genet. 2010; 121: 1441-53.

[14] Perata P, Voesenek LACJ . Submergence tolerance in rice requires Sub1A, an ethylene-response-factor-like gene . Trends Plant Sci. 2007; 12: 43-6.

[15] Schmitz AJ, Folsom JJ, Jikamaru Y et al. SUB1A-mediated submer- gence tolerance response in rice involves differential regulation of the brassinosteroid pathway. New Phytol. 2013; 198: 1060-70.

[16] Hattori Y, Nagai K, Furukawa S et al. The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature. 2009; 460: 1026-30.

[17] Ayano M, Kani T, Kojima M et al. Gibberellin biosynthesis and signal transduction is essential for internode elongation in Deepwater rice. Plant Cell Environ. 2014; 37: 2313-24.

[18] Voesenek LACJ, Rijnders JHGM, Peeters AJM et al. Plant hormones regulate fast shoot elongation under water: from genes to com- munities. Ecology. 2004; 85: 16-27.

[19] Mergemann H, Sauter M . Ethylene induces epidermal cell death at the site of adventitious root emergence in rice. Plant Physiol. 2000; 124: 609-14.

[20] Steffens B, Sauter M . Epidermal cell death in rice (Oryza sativa L.) is regulated by ethylene, gibberellin and abscisic acid . Plant Physiol. 2005; 139: 713-21.

[21] Steffens B, Wang J, Sauter M . Interactions between ethylene, gibberellin and abscisic acid regulate emergence and growth rate of adventitious roots in deep water rice. Planta. 2006; 223: 604-12.

[22] Steffens B, Geske T, Sauter M . Aerenchyma formation in the rice stem and its promotion by H₂O₂ . New Phytol. 2011; 190: 369-78.

[23] Nakamura M, Noguchi K . Tolerant mechenism to O₂ deficiency under submergence conditions in plants . J Plant Res. 2020; 133: 343-71.

[24] Mignolli F, Todaro JS, Vidoz ML . Internal aeration and respiration of submerged tomato hypocotyls are enhanced by ethylene- mediated aerenchyma formation and hypertrophy. Physiol Plant. 2020; 169: 49-63.

[25] Bartoli G, Forino LM, Durante M et al. A lysigenic programmed cell death-dependent process shapes schizogenously formed aerenchyma in the stems of the waterweed Egeria densa . Ann Bot. 2015; 116: 91-9.

[26] Drew MC . Oxygen deficiency and root metabolism: injury and acclimation under hypoxia and anoxia. Annu Rev Plant Physiol Plant Mol Biol. 1997; 48: 223-50.

[27] Geigenberger P. Response of plant metabolism to too little oxy- gen. Curr Opin Plant Biol. 2003; 6: 247-56.

[28] Bailey-Serres J, Chang R . Sensing and signalling in response to oxygen deprivation in plants and other organisms. Ann Bot. 2005; 96: 507-18.

[29] Paudel KR, Panth N . Phytochemical profile and biological activ- ity of Nelumbo nucifera . Evid Based Complement Alternat Med. 2015; 2015: 1-16.

[30] Deng X, Zhu L, Fang T et al. Analysis of isoquinoline alkaloid composition and wound-induced variation in Nelumbo using HPLC-MS/MS . J Agric Food Chem. 2016; 64: 1130-6.

[31] Lin Z, Zhang C, Cao D et al. The latest studies on lotus (Nelumbo nucifera) - an emerging horticultural model plant . Int J Mol Sci. 2019; 20: 3680.

[32] Matthews PGD, Seymour RS . Stomata actively regulate internal aeration of the sacred lotus Nelumbo nucifera . Plant Cell Environ. 2014; 37: 402-13.

[33] Ymamoto M, Nishikawa N, Mayama H et al. Theoretical explana- tion of the lotus effect: superhydrophobic property changes by removal of nanostructures from the surface of a lotus leaf. Lang Des. 2015; 31: 7355-63.

[34] Colmer TD, Pedersen O . Underwater photosynthesis and respi- ration in leaves of submerged wetland plants: gas films improve CO₂ and O₂ exchange . New Phytol. 2008; 177: 918-26.

[35] Pedersen O, Rich SM, Colmer TD . Surviving floods: leaf gas films improve O₂ and CO₂ exchange root aeration, and growth of completely submerged rice . Plant J. 2009; 58: 147-56.

[36] Wang B, Jin Q, Zhang S et al. Genome-wide transcriptional analysis of submerged lotus reveals cooperative regulation and gene responses. Scientific reports. 2018; 8: 9187.

[37] Jin Q, Xu Y, Mattson N et al. Identification of submergence- responsive microRNAs and their targets reveals complex miRNA-mediated regulatory networks in lotus (Nelumbo nucifera Gaertn) . Front Plant Sci. 2019; 8: 6.

[38] Verbruggen N, Hermans C . Proline accumulation in plants: a review. Amino Acids. 2008; 35: 753-9.

[39] Steele MR, Gitelson AA, Rundquist DC . A comparison of two techniques for nondestructive measurement of chlorophyll con- tent in grapevine leaves. Agron J. 2008; 100: 779-82.

[40] Deng X, Zhao L, Fang T et al. Investigation of benzylisoquinoline alkaloid biosynthetic pathway and its transcriptional regulation in lotus. Hortic Res. 2018; 5: 29.

[41] Ming R, Van Buren R, Liu Y et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.) . Genome Biol. 2013; 14: R41.

[42] Muhammad R, He G, Yang G et al. AP2/ERF transctiption factor in rice: genome-wide canvas and synthenic relationships between monocots and eudicots. Evol Bioinforma. 2012; 8: 321-55.

[43] Li J, Chory J . A putative leucine-rich repeat receptor kinase involved in brassinosteroid signal transduction. Cell. 1997; 90: 929-38.

[44] Tran LSP, Nakashima K, Sakuma Y et al. Co-expression of the stress-inducible zinc finger homeodomain ZFHD1 and NAC transcription factors enhances expression of the ERD1 gene in Arabidopsis . Plant J. 2007; 49: 46-63.

[45] Yang SF, Hoffman NE . Ethylene biosynthesis and its regulation in higher plants. Annu Rev Plant Physiol. 1984; 35: 155-89.

[46] Kende H. Ethylene biosynthesis. Annu Rev Plant Physiol Plant Mol Biol. 1993; 44: 283-307.

[47] Lorenzo O, Piqueras R, Sanchez-Serrano JJ et al. ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jas- monate pathways in plant defense. Plant Cell. 2003; 15: 165-78.

[48] Mao JL, Miao ZQ, Wang Z et al. Arabidopsis ERF1 mediates cross-talk beween ethylene and auxin biosynthesis during pri- mary root elongation by regulating ASA1 expression. PLoS Genet. 2016; e1005760.

[49] Cho HT, Kende H . Expression of expansin genes is correlated with growth in deep water rice . Plant Cell. 1997; 9: 1661-71.

[50] Vreeburg RAM, Benschop JJ, Peeters AJM et al. Ethylene regulates fast apoplastic acidification and expansin a transcription dur- ing submergence-induced petiole elongation in Rumex palustris . Plant J. 2005; 43: 597-610.

[51] Gibbs DJ, Md Isa N, Movahedi M et al. Nitric oxide sensing in plants is mediated by proteolytic control of group VII ERF transcription factors. Mol Cell. 2014; 53: 369-79.

[52] van Loon LC, van Strien EA . The families of pathogenesis-related proteins, their activities, and comparative analysis of PR-1 type proteins. Physiol Mol Plant Pathol. 1999; 55: 85-97.

[53] Kende H, van der Knaap E, Cho HT . Deepwater rice: a model plant to study stem elongation. Plant Physiol. 1998; 118: 1105-10.

[54] Fukao T, Bailey-Serres J . Submergence tolerance conferred by Sub1A is mediated by SLR1 and SLRL1 restriction of gibberellin responses in rice. Proc Natl Acad Sci U S A. 2008; 105: 16814-9.

[55] Zhu GH, Ye NH, Zhang JH . Glucose-induced delay of seed germi- nation in rice is mediated by the suppression of ABA catabolism rather than an enhancement of ABA biosynthesis. Plant Cell Physiol. 2009; 50: 644-51.

[56] Li J, Xiong Y, Li Y et al. Comprehensive analysis and functional studies of WRKY transcriptional factors in Nelumbo nucifera . Int J Mol Sci. 2019; 20: 5006.

[57] Mommer L, Pons TL, Visser EJW . Photosynthetic consequences of phenotypic plasticity in response to submergence: Rumex palustris as a case study. J Exp Bot. 2006; 57: 283-90.

[58] Vashisht D, Hesselink A, Pierik R et al. Natural variation of submergence tolerance among Arabidopsis thaliana accessions . New Phytol. 2011; 190: 299-310.

[59] Lee Y, Kende H . Expressions of a-expansins is correlated with intermodal elongation in Deepwater rice. Plant Physiol. 2001; 127: 645-54.

[60] Muhanmmad R, He G, Yang G et al. AP2/ERF transcription factor in rice: genome-wide canvas and syntenic relationships between monocots and eudicots. Evol Bioinforma. 2012; 8: 321-55.

[61] Bailey-Serres J, Fukao T, Gibbs DJ et al. Making sense of low oxygen sensing. Trends Plant Sci. 2012; 17: 129-38.

[62] Mommer L, Pons TL, Wolters-Arts M et al. Submergence-induced morphological, anatomical, and biochemical responses in a ter- restrial species affects gas diffusion resisance and photosyn- thetic performance. Plant Physiol. 2005; 139: 497-508.

[63] Porra RJ . The chequered history of the development and use of simultaneous equations for the accurate determination of chlorophylls a and b . Photosynth Res. 2002; 73: 149-56.

[64] Ramírez V, Coego A, Lopez A et al. Drought tolerance in Arabidop- sis is controlled by the OCP3 disease resistance regulator . Plant J. 2009; 58: 578-91.

[65] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[66] Pertea M, Pertea GM, Antonescu CM . String tie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5.

[67] Florea L, Song L, Salzberg SL . Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues. F1000Res. 2013; 2: 188.

[68] Wang LK, Feng ZK, Wang X et al. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010; 26: 136-8.

[69] Sturn A, Quackenbush J, Trajanoski Z . Genesis: cluster analysis of microarray data. Bioinformatics. 2002; 18: 207-8.

[70] Xie C, Mao X, Huang J et al. KOBAS 2.0: a web server for anno- tation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011; 39: W316-22.

[71] Shannon P, Markiel A, Owen O et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13: 2498-504.

[72] Livak KJ, Schmittgen TD . Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001; 25: 402-8.