High transcriptome plasticity drives phosphate starvation responses in tomato

Abstract

Tomato is an important vegetable crop and fluctuating available soil phosphate (Pi) level elicits several morpho-physiological responses driven by underlying molecular responses. Therefore, understanding these molecular responses at the gene and isoform levels has become critical in the quest for developing crops with improved Pi use efficiency. A quantitative time-series RNA-seq analysis was performed to decipher the global transcriptomic changes that accompany Pi starvation in tomato. Apart from changes in the expression levels of genes, there were also alterations in the expression of alternatively-spliced transcripts. Physiological responses such as anthocyanin accumulation, reactive oxygen species generation and cell death are obvious 7 days after Pi deprivation accompanied with the maximum amount of transcriptional change in the genome making it an important stage for in-depth study while studying Pi stress responses (PSR). Our study demonstrates that transcriptomic changes under Pi deficiency are dynamic and complex in tomato. Overall, our study dwells on the dynamism of the transcriptome in eliciting a response to adapt to low Pi stress and lays it bare. Findings from this study will prove to be an invaluable resource for researchers using tomato as a model for understanding nutrient deficiency.

Keywords

Tomato / Pi starvation / Transcriptome / Gene expression / Alternative splicing / Co-expression network

References

[1] AmesBN. NeufeldE, GinsburgV. Assay of inorganic phosphate, total phosphate and phosphatases. Methods in Enzymology, 1966 New York Academic Press 115-118

[2] ApelK, HirtH. Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Annu Rev Plant Biol, 2004, 55: 373-399

[3] BaldwinJC, KarthikeyanAS, CaoA, RaghothamaKG. Biochemical and molecular analysis of LePS2;1: a phosphate starvation induced protein phosphatase gene from tomato. Planta, 2008, 228(2):273-280

[4] BaldwinJC, KarthikeyanAS, RaghothamaKG. LEPS2, a phosphorus starvation-induced novel acid phosphatase from tomato. Plant Physiol, 2001, 125(2):728-737

[5] BustosR, CastrilloG, LinharesF, PugaMI, RubioV, Pérez-PérezJ, SolanoR, LeyvaA, Paz-AresJ. A central regulatory system largely controls transcriptional activation and repression responses to phosphate starvation in Arabidopsis. PLoS Genet, 2010, 6(9

[6] CalixtoCPG, GuoW, JamesAB, TzioutziouNA, EntizneJC, PanterPE, KnightH, NimmoHG, ZhangR, BrownJWS. Rapid and dynamic alternative splicing impacts the Arabidopsis cold response transcriptome. Plant Cell, 2018, 30 7):1424-1444 http://www.plantcell.org/lookup/doi/10.1105/tpc.18

[7] CastrilloG, Lima TeixeiraPJP, ParedesSH, LawTF, de LorenzoL, FeltcherME, FinkelOM, BreakfieldNW, MieczkowskiP, JonesCD, Paz-AresJ, DanglJL. Root microbiota drive direct integration of phosphate stress and immunity. Nature, 2017, 543: 513-518

[8] ChardinC, KrappA, SchenkST, HirtH, ColcombetJ. Review: mitogen-activated protein kinases in nutritional signaling in Arabidopsis. Plant Sci, 2017, 260: 101-108

[9] ChiouT-J, LinS-I. Signaling network in sensing phosphate availability in plants. Annu Rev Plant Biol, 2011, 62(1):185-206

[10] CuiP, LinQ, DingF, et al.. A comparison between ribo-minus RNA-sequencing and polyA-selected RNA-sequencing. Genomics, 2010, 96: 259-265

[11] DesnosT. Root branching responses to phosphate and nitrate. Curr Opin Plant Biol, 2008, 11: 82-87

[12] DobinA, DavisCA, SchlesingerF, DrenkowJ, ZaleskiC, JhaS, BatutP, ChaissonM, GingerasTR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics, 2013, 29(1):15-21

[13] DongC, HeF, BerkowitzO, LiuJ, CaoP, TangM, ShiH, WangW, LiQ, ShenZ, WhelanJ, ZhengL. Alternative splicing plays a critical role in maintaining mineral nutrient homeostasis in Rice (Oryza sativa). Plant Cell, 2018, 30(10):2267-2285

[14] DongJ, MaG, SuiL, WeiM, SatheeshV, ZhangR, GeS, LiJ, ZhangT-E, WittwerC, JessenHJ, ZhangH, AnGY, ChaoDY, LiuD, LeiM. Inositol pyrophosphate InsP₈ acts as an intracellular phosphate signal in Arabidopsis. Mol Plant, 2019, 2019(12):1463-1473

[15] DuQ, WangK, XuC, ZouC, XieC, XuY, et al.. Strand-specific RNA-Seq transcriptome analysis of genotypes with and without low-phosphorus tolerance provides novel insights into phosphorus-use efficiency in maize. BMC Plant Biol, 2016, 16(1):1-12

[16] FAO (2015) Food and Agriculture Organization of the United Nations. In: FAO Statistical Yearbook 2015. World Food and Agriculture, United Nations. Available at http://faostat.fao.org

[17] FinnRD, ClementsJ, EddySR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res, 2011, 39: W29-37

[18] GeC, GeorgievA, ÖhmanA, WieslanderÅ, KellyAA. Tryptophan residues promote membrane association for a plant lipid glycosyltransferase involved in phosphate stress. J Biol Chem, 2011, 286(8):6669-6684

[19] GuZ, EilsR, SchlesnerM. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics., 2016, 32(18):2847-2849

[20] GuoW, CalixtoCPG, BrownJWS, ZhangR. TSIS: an R package to infer alternative splicing isoform switches for time-series data. Bioinformatics, 2017, 33(20):3308-3310

[21] HinsingerP. Bioavailability of soil inorganic P in the rhizosphere as affected by root induced chemical changes: a review. Plant and Soil, 2001, 237(2):173-195

[22] HoehenwarterW, MönchgesangS, NeumannS, MajovskyP, AbelS, MüllerJ. Comparative expression profiling reveals a role of the root apoplast in local phosphate response. BMC Plant Biol, 2016, 16(1):106

[23] HorstWJ, KamhM, JibrinJM, ChudeVO. Agronomic measures for increasing P availability to crops. Plant and Soil, 2001, 237(2):211-223

[24] Huerta-CepasJ, SerraF, BorkP. ETE 3: reconstruction, analysis, and visualization of phylogenomic data. Mol Biol Evol, 2016, 33(6):1635-1638

[25] JammesF, SongC, ShinD, et al.. MAP kinases MPK9 and MPK12 are preferentially expressed in guard cells and positively regulate ROS-mediated ABA signaling. Proc Natl Acad Sci U S A, 2009, 106: 20520-20525

[26] JinJ, TianF, YangD-C, et al.. PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res, 2017, 45: D1040-D1045

[27] LangfelderP, HorvathS. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics, 2008, 9(1):559

[28] LiP, LiYJ, ZhangFJ, ZhangGZ, JiangXY, YuHM, HouBK. The Arabidopsis UDP-glycosyltransferases UGT79B2 and UGT79B3, contribute to cold, salt and drought stress tolerance via modulating anthocyanin accumulation. Plant J, 2017, 89(1):85-103

[29] LiQ, XiaoG, ZhuYX. Single-nucleotide resolution mapping of the Gossypium raimondii transcriptome reveals a new mechanism for alternative splicing of introns. Mol Plant, 2014, 7: 829-840

[30] LiW, LinW-D, RayP, LanP. Schmidt W (2013) genome-wide detection of condition-sensitive alternative splicing in Arabidopsis roots. Plant Physiol, 2013, 162(3):1750-1763

[31] LiuC, MuchhalUS, RaghothamaKG. Differential expression of TPSI1, a phosphate starvation-induced gene in tomato. Plant Mol Biol, 1997, 33(5):867-874

[32] LiuZ, QinJ, TianX, XuS, WangY, LiH, WangX, PengH, YaoY, HuZ, NiZ, XinM, SunQ. Global profiling of alternative splicing landscape responsive to drought, heat and their combination in wheat (Triticum aestivum L.). Plant Biotechnol J, 2018, 16(3):714-726

[33] López-ArredondoDL, Leyva-GonzálezMA, González-MoralesSI, López-BucioJ, Herrera-EstrellaL. Phosphate nutrition: improving low-phosphate tolerance in crops. Annu Rev Plant Biol, 2014, 65(1):95-123

[34] López-BucioJ, Cruz-RamírezA, Herrera-EstrellaL. The role of nutrient availability in regulating root architecture. Curr Opin Plant Biol, 2003, 6(3):280-287

[35] Lopez-BucioJ, Hernández-AbreuE, Sánchez-CalderónL, Nieto-JacoboMF, SimpsonJ, Herrera-EstrellaL. Phosphate availability alters architecture and causes changes in hormone sensitivity in the Arabidopsis root system. Plant Physiol, 2002, 129(1):244-256

[36] LoveMI, HuberW, AndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014, 15(12):550

[37] LynchJP, BrownKM. LynchJP, DeikmanJ. Regulation of root architecture by phosphorus availability. Phosphorus in plant biology: regulatory roles in molecular, cellular, organismic, and ecosystem processes, 1998 Maryland Amer. Soc. Plant Physiol 148-156

[38] MarquezY, BrownJW, SimpsonC, BartaA, KalynaM. Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res, 2012, 22(6):1184-1195

[39] MarrsKA, AlfenitoMR, LloydAM, WalbotV. A glutathione S-transferase involved in vacuolar transfer encoded by the maize gene Bronze-2. Nature, 1995, 375(6530):397-400

[40] MasonMJ, FanG, PlathK, ZhouQ, HorvathS. Signed weighted gene co-expression network analysis of transcriptional regulation in murine embryonic stem cells. Signed weighted gene co-expression network analysis of transcriptional regulation in murine embryonicstem cells BMC Genomics, 2009, 10(1):327

[41] MissonJ, RaghothamaKG, JainA, JouhetJ, BlockMA, BlignyR, OrtetP, CreffA, SomervilleS, RollandN, DoumasP, NacryP, Herrerra-EstrellaL, NussaumeL, ThibaudM-C. A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Natl Acad Sci U S A, 2005, 102: 11934

[42] MiuraK, RusA, SharkhuuA, et al.. The Arabidopsis SUMO E3 ligase SIZ1 controls phosphate deficiency responses. Proc Natl Acad Sci U S A, 2005, 102: 7760-7765

[43] MuellerLA, GoodmanCD, SiladyRA, WalbotV. AN9, a petunia glutathione S-transferase required for anthocyanin sequestration, is a flavonoid-binding protein. Plant Physiol, 2000, 123(4):1561-1570

[44] MüllerJ, ToevT, HeistersM, TellerJ, MooreKL, HauseG, DineshDC, BürstenbinderK, AbelS. Iron-dependent callose deposition adjusts root meristem maintenance to phosphate availability. Dev Cell, 2015, 33: 216-230

[45] MüllerR, MorantM, JarmerH, NilssonL, NielsenTH. Genome-wide analysis of the arabidopsis leaf transcriptome reveals interaction of phosphate and sugar metabolism. Plant Physiol, 2007, 143(1):156-171

[46] MuneerS, JeongBR. Proteomic analysis provides new insights in phosphorus homeostasis subjected to pi ( inorganic phosphate ) starvation in tomato plants (Solanum lycopersicum L.). PLoS one, 2015, 10(7):e0134103

[47] NishidaS, KakeiY, ShimadaY, FujiwaraT. Genome-wide analysis of specific alterations in transcript structure and accumulation caused by nutrient deficiencies in Arabidopsis thaliana. Plant J, 2017, 91(4):741-753

[48] PerteaM, PerteaGM, AntonescuCM, et al.. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol, 2015, 33: 290-295

[49] Pfaff J, Denton AK, Usadel B, Pfaff C (2020) Phosphate starvation causes different stress responses in the lipid metabolism of tomato leaves and roots. Biochim Biophys Acta Mol Cell Biol Lipids 1865:158763. https://doi.org/10.1016/j.bbalip.2020.158763

[50] PitzschkeA, HirtH. Mitogen-activated protein kinases and reactive oxygen species signaling in plants. Plant Physiol, 2006, 141(2):351-356

[51] PugaI, Rojas-TrianaM, de LorenzoL, LeyvaA, RubioV, Paz-aresJ. Novel signals in the regulation of pi starvation responses in plants : facts and promises. Curr Opin Plant Bio, 2017, 39: 40-49

[52] PugaMI, MateosI, CharukesiR, WangZ, Franco-ZorrillaJM, de LorenzoL, IrigoyenML, MasieroS, BustosR, RodriguezJ, LeyvaA, RubioV, SommerH, Paz-AresJ. SPX1 is a phosphate-dependent inhibitor of PHOSPHATE STARVATION RESPONSE 1 in Arabidopsis. Proc Natl Acad Sci, 2014, 111(41):14947-14952

[53] Qi W, Manfield IW, Muench SP, Baker A (2017) AtSPX1 affects the AtPHR1–DNA-binding equilibrium by binding monomeric AtPHR1 in solution. Biochem J 474:3675–3687. 474:3675-3687. https://doi.org/10.1042/bcj20170522

[54] RaghothamaKG. Phosphate acquisition. Annual Rev Plant Physiol Plant Mol Biol, 1999, 50(1):665-693

[55] RawatV, AbdelsamadA, PietzenukB, SeymourDK, KoenigD, WeigelD, PecinkaA, SchneebergerK. Improving the annotation of Arabidopsis lyrata using RNA-Seq data. PLoS One, 2015, 10(9

[56] RiedMK, WildR, ZhuJ, PipercevicJ, SturmK, BrogerL, HarmelRK, AbriataLA, HothornLA, FiedlerD, HillerS, HothornM. Inositol pyrophosphates promote the interaction of SPX domains with the coiled-coil motif of PHR transcription factors to regulate plant phosphate homeostasis. Nat Commun, 2021, 12 1):384

[57] RubioV, LinharesF, SolanoR, MartinAC, IglesiasJ, LeyvaA, et al.. A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Genes Dev, 2001, 15(16):2122-2133

[58] SeccoD, JabnouneM, WalkerH, ShouH, WuP, PoirierY, WhelanJ. Spatio-temporal transcript profiling of Rice roots and shoots in response to phosphate starvation and recovery. Plant Cell, 2013, 25(11):4285-4304

[59] ShenY, ZhouZ, WangZ, LiW, FangC, WuM, MaY, LiuT, KongL-A, PengD-L, TianZ. Global dissection of alternative splicing in paleopolyploid soybean. Plant Cell, 2014, 26(3):996-1008

[60] SiebersM, DörmannP, HölzlG. Membrane remodelling in phosphorus-deficient plants. Annu Plant Rev, 2015 Hoboken Wiley 237-263

[61] SrivastavaR, AkashParida AP, ChauhanPK, KumarR. Identification, structure analysis, and transcript profiling of purple acid phosphatases under Pi deficiency in tomato (Solanum lycopersicum L.) and its wild relatives. Int J Biol Macromol, 2020, 165: 2253-2266

[62] StaigerD, BrownJW. Alternative splicing at the intersection of biological timing, development, and stress responses. Plant Cell, 2013, 25(10):3640-3656

[63] Suen PK, Zhang S, Sun SS-M (2015) Molecular characterization of a tomato purple acid phosphatase during seed germination and seedling growth under phosphate stress. Plant Cell Rep 34:981–992. https://doi.org/10.1007/s00299-015-1759-z

[64] TeigeM, ScheiklE, EulgemT, et al.. The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis. Mol Cell, 2004, 15: 141-152

[65] ThatcherSR, DanilevskayaON, MengX, BeattyM, Zastrow-HayesG, HarrisC, van AllenB, HabbenJ, LiB. Genome-wide analysis of alternative splicing during development and drought stress in maize. Plant Physiol, 2016, 170(1):586-599

[66] TianP, LinZ, LinD, DongS, HuangJ, HuangT. The pattern of DNA methylation alteration, and its association with the changes of gene expression and alternative splicing during phosphate starvation in tomato. Plant J, 2021, 108(3):841-858

[67] TyburskiJ, DunajskaK, TretynA. Reactive oxygen species localization in roots of Arabidopsis thaliana seedlings grown under phosphate deficiency. Plant Growth Regul, 2009, 59(1):27-36

[68] Vitting-SeerupK, SandelinA. The landscape of isoform switches in human cancers. Mol Cancer Res American Association for Cancer Res, 2017, 15(9):1206-1220

[69] WangY, GarvinDF, KochianLV. Rapid Induction of Regulatory and Transporter Genes in Response to Phosphorus , Potassium , and Iron Deficiencies in Tomato Roots . Evidence for Cross Talk and Root / Rhizosphere-Mediated Signals. Plant Physiol, 2002, 130(3):1361-1370

[70] WangZ, RuanW, ShiJ, ZhangL, XiangD, YangC, LiC, WuZ, LiuY, YuY, ShouH, MoX, MaoC, WuP. Rice SPX1 and SPX2 inhibit phosphate starvation responses through interacting with PHR2 in a phosphate-dependent manner. Proc Natl Acad Sci, 2014, 111(41):14953-14958

[71] WildR, GerasimaiteR, JungJY, TruffaultV, PavlovicI, SchmidtA, SaiardiA, JessenHJ, PoirierY, HothornM, MayerA. Control of eukaryotic phosphate homeostasis by inositol polyphosphate sensor domains. Science, 2016, 352(6288):986-990

[72] WuP, ShouH, XuG, LianX. Improvement of phosphorus efficiency in rice on the basis of understanding phosphate signaling and homeostasis. Curr Opin Plant Biol, 2013, 16(2):205-212

[73] XieY-J, XuS, HanB, WuM-Z, YuanX-X, HanY, GuQ, XuDK, YangQ, ShenWB. Evidence of Arabidopsis salt acclimation induced by up-regulation of HY1 and the regulatory role of RbohD-derived reactive oxygen species synthesis. Plant J, 2011, 66(2):280-292

[74] XingY, JiaW, ZhangJ. AtMKK1 mediates ABA-induced CAT1 expression and H2O2 production via AtMPK6-coupled signaling in Arabidopsis. Plant J, 2008, 54: 440-451

[75] YuB, XuC, BenningC. Arabidopsis disrupted in SQD2 encoding sulfolipid synthase is impaired in phosphate-limited growth. Proc Natl Acad Sci U S A National Academy of Sciences, 2002, 99(8):5732-5737

[76] Yuan H, Liu D (2008) Signaling components involved in plant responses to phosphate starvation. J Integr Plant Biol 50:849-859. https://onlinelibrary.wiley.com/doi/10.1111/j.1744-7909.2008.00709.x.

[77] ZhaiL, SunC, FengY, LiD, ChaiX, WangL, SunQ, ZhangG, LiY, WuT, ZhangX, XuX, WangY, HanZ. AtROP6 is involved in reactive oxygen species signaling in response to iron-deficiency stress in Arabidopsis thaliana. FEBS Lett, 2018, 592(20):3446-3459

[78] ZhangR, CalixtoCPG, MarquezY, VenhuizenP, TzioutziouNA, GuoW, SpensleyM, EntizneJC, LewandowskaD, ten HaveS, Frei dit FreyN, HirtH, JamesAB, NimmoHG, BartaA, KalynaM, BrownJWS. A high quality Arabidopsis transcriptome for accurate transcript-level analysis of alternative splicing. Nucleic Acids Res, 2017, 45(9):5061-5073

[79] ZhangY, WangY, WangE, WuX, ZhengQ, HanY, LinW, LiuZ, LinW. SlPHL1, a MYB-CC transcription factor identified from tomato, positively regulates the phosphate starvation response. Physiol Plant, 2021, 173(3):1063-1077

[80] ZhangY, ZhouY, ChenS, et al.. Gibberellins play dual roles in response to phosphate starvation of tomato seedlings, negatively in shoots but positively in roots. J Plant Physiol, 2019, 234–235: 145-153

[81] Zhao P, You Q, Lei M (2019) A CRISPR/Cas9 deletion into the phosphate transporter SlPHO1;1 reveals its role in phosphate nutrition of tomato seedlings. Physiol Plantarum 167: 556–563. https://onlinelibrary.wiley.com/doi/full/10.1111/ppl.12897

[82] ZhouJ, HuQ, XiaoX, et al.. Mechanism of phosphatesensing and signaling revealed by rice SPX1-PHR2 complex structure. Nat Commun, 2021, 12: 7040

[83] ZhuJ, LauK, PuschmannR, HarmelRK, ZhangY, PriesV, GauglerP, BrogerL, DuttaAK, JessenHJ, SchaafG, FernieAR, HothornLA, FiedlerD, HothornM. Two bifunctional inositol pyrophosphate kinases/phosphatases control plant phosphate homeostasis. ELife, 2019, 8