A universal pipeline for mobile mRNA detection and insights into heterografting advantages under chilling stress

Abstract

Heterografting has long been used to enhance the chilling tolerance of temperature-sensitive crops, including watermelon, whose mechanism is known to involve bidirectional long-distance mRNA movements. Despite several studies reporting on mobile mRNA (mb-mRNA) profiles in plants, accurate identification of mb-mRNAs is challenging owing to an array of technical problems. Here, we developed a bioinformatical pipeline that took most of the known technical concerns into consideration and is considered to be a universal tool for mb-mRNA detection in heterografts. By applying this pipeline to a commercial watermelon–bottle gourd heterografting system, we detected 130 and 1144 mb-mRNAs upwardly and 167 and 1051 mb-mRNAs downwardly transmitted under normal and chilling-stress conditions, respectively. Quantitative real-time PCR indicated a high accuracy rate (88.2%) of mb-mRNA prediction with our pipeline. We further revealed that the mobility of mRNAs was not associated with their abundance. Functional annotation and classification implied that scions may convey the stress signal to the rootstock, subsequently triggering energy metabolism reprogramming and abscisic acid-mediated stress responses by upward movement of effective mRNAs, ultimately leading to enhanced chilling tolerance. This study provides a universal tool for mb-mRNA detection in plant heterografting systems and novel insights into heterografting advantages under chilling stress.

References

[1] Noh, J. et al. Effect of heat treatment around the fruit set region on growth and yield of watermelon [Citrullus lanatus (Thunb.) Matsum.and Nakai]. Physiol. Mol. Biol. Plants 19, 509-514 (2013).

[2] Xu, J., Zhang, M., Liu, G., Yang, X. & Hou, X. Comparative transcriptome profiling of chilling stress responsiveness in grafted watermelon seedlings. Plant Physiol. Biochem. 109, 561-570 (2016).

[3] Davis, A. R. et al. Cucurbit grafting. Crit. Rev. Plant Sci. 27, 50-74 (2008).

[4] Ding, M. et al. Physiological advantages of grafted watermelon (Citrullus lanatus) seedlings under low-temperature storage in darkness. Hortscience 46, 993-996 (2011).

[5] Petropoulos, S. et al. Fruit volatiles, quality, and yield of watermelon as affected by grafting. J. Agr. Sci. Technol. 16, 873-885 (2014).

[6] Kijima, J . Grafting experiment on watermelon. J. Okitsu Hort. Soc. 34, 57-71 (1938).

[7] Xu, P. et al. Population genomic analyses from low-coverage rad-seq data: a case study on the non-model cucurbit bottle gourd. Plant J. 77, 430-442 (2014).

[8] Liu, N. et al. Genome-wide identification and comparative analysis of grafting-responsive mRNA in watermelon grafted onto bottle gourd and squash rootstocks by high-throughput sequencing. Mol. Genet. Genomics 291, 621-633 (2016).

[9] Wang, L. et al. Generating homo- and heterografts between watermelon and bottle gourd for the study of cold-responsive MicroRNAs. J. Vis. Exp. 141, e58242 (2018).

[10] Thieme, C. J. et al. Endogenous Arabidopsis messenger RNAs transported to distant tissues. Nat. Plants 1, 15025-15032 (2015).

[11] Yang, Y. et al. Messenger RNA exchange between scions and rootstocks in grafted grapevines. BMC Plant Biol. 15, 251 (2015).

[12] Xia, C. et al. Elucidation of the mechanisms of long-distance mRNA movement in a Nicotiana benthamiana/tomato heterograft system. Plant Physiol. 177, 745-758 (2018).

[13] Zhang, Z. et al. Vascular-mediated signalling involved in early phosphate stress response in plants. Nat. Plants 2, 16033 (2016).

[14] Notaguchi, M., Higashiyama, T. & Suzuki, T. Identification of mRNAs that move over long distances using an RNA-Seq analysis of Arabidopsis/Nicotiana benthamiana heterografts. Plant Cell Physiol. 56, 311-321 (2015).

[15] Bremer, A. et al. Intrinsically disordered stress protein COR15A resides at the membrane surface during dehydration. Biophys. J. 113, 572-579 (2017).

[16] Bremer, A., Wolff, M., Thalhammer, A. & Hincha, D. K. Folding of intrinsically disordered plant LEA proteins is driven by glycerol‐induced crowding and the presence of membranes. FEBS J. 284, 919-936 (2017).

[17] Houde, M., Dallaire, S., N'Dong, D. & Sarhan, F. Overexpression of the acidic dehydrin WCOR410 improves freezing tolerance in transgenic strawberry leaves. Plant. Biotechnol. J. 2, 381-387 (2004).

[18] Liu, H. et al. Overexpression of ShDHN, a dehydrin gene from Solanum habrochaites enhances tolerance to multiple abiotic stresses in tomato. Plant Sci. 231, 198-211 (2015).

[19] Ruelland, E., Vaultier, M. N., Zachowski, A. & Hurry, V. Chapter 2 cold signalling and cold acclimation in plants. Adv. Bot. Res. 49, 35-150 (2009).

[20] Thomashow, M. F. Plant cold acclimation: freezing tolerance genes and regulatory mechanisms. Annu. Rev. Plant Physiol. Plant Mol. Biol. 50, 571-599 (1999).

[21] Martinière, A. et al. Homeostasis of plasma membrane viscosity in fluctuating temperatures. New Phytol 192, 328-337 (2011).

[22] Mustilli, A. C., Merlot, S., Vavasseur, A., Fenzi, F. & Giraudat, J. Arabidopsis ost1 protein kinase mediates the regulation of stomatal aperture by abscisic acid and acts upstream of reactive oxygen species production. Plant Cell. 14, 3089-3099 (2002).

[23] Ding, Y. et al. Ost1 kinase modulates freezing tolerance by enhancing ice1 stability in arabidopsis. Dev. Cell. 32, 278-289 (2015).

[24] Kim, G., LeBlanc, M. L., Wafula, E. K. & Westwood, J. H. Genomic-scale exchange of mRNA between a parasitic plant and its hosts. Science 345, 808-811 (2014).

[25] Mortazavi, A. et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5, 621 (2008).

[26] Calderwood, A., Kopriva, S. & Morris, R. J. Transcript abundance explains mRNA mobility data in Arabidopsis thaliana. Plant Cell. 28, 610-615 (2016).

[27] Paultre, D. S. G., Gustin, M. P., Molnar, A. & Oparka, K. J. Lost in transit: long-distance trafficking and phloem unloading of protein signals in arabidopsis homografts. Plant Cell. 28, 2016-2025 (2016).

[28] Shi, H. et al. Manipulation of arginase expression modulates abiotic stress tolerance in Arabidopsis: effect on arginine metabolism and ROS accumulation. J. Exp. Bot. 64, 1367-1379 (2013).

[29] Maguire, T. L., Collins, G. G. & Sedgley, M. A modified CTAB DNA extraction procedure for plants belonging to the family Proteaceae. Plant Mol. Biol. Rep. 12, 106-109 (1994).

[30] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔCT method . Methods 25, 402-408 (2001).

[31] Guo, S. et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions. Nat. Genet. 45, 51-58 (2013).

[32] Wang, Y. et al. GourdBase: a genome-centered multi-omics database for the bottle gourd (Lagenaria siceraria), an economically important cucurbit crop. Sci. Rep. 8, 3604 (2018).

[33] Eden, E. et al. Gorilla: a tool for discovery and visualization of enriched go terms in ranked gene lists. BMC Bioinformatics 10, 48-0 (2009).