Abstract
Nutrient stress is the most important environmental stress that limits plant growth and development. Although recent evidence highlights the vital functions of long non-coding RNAs (lncRNA) in response to single nutrient stress in some model plants, a comprehensive investigation of the effect of lncRNAs in response to nutrient stress has not been performed in Arabidopsis thaliana. Here, we presented the identification and characterization of lncRNAs under seven nutrient stress conditions. The expression pattern analysis revealed that aberrant expression of lncRNAs is a stress-specific manner under nutrient stress conditions and that lncRNAs are more sensitive to nutrient stress than protein-coding genes (PCGs). Moreover, competing endogenous RNA (ceRNA) network and lncRNA-mRNA co-expression network (CEN) were constructed to explore the potential function of these lncRNAs under nutrient stress conditions. We further combined different expressed lncRNAs with ceRNA network and CEN to select key lncRNAs in response to nutrient stress. Together, our study provides important information for further insights into the role of lncRNAs in response to stress in plants.
Similar content being viewed by others
References
Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11:Artn R106. https://doi.org/10.1186/Gb-2010-11-10-R106
Anders S, Pyl PT, Huber W (2015) HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31:166–169. https://doi.org/10.1093/bioinformatics/btu638
Argueso CT, Ferreira FJ, Kieber JJ (2009) Environmental perception avenues: the interaction of cytokinin and environmental response pathways. Plant Cell Environ 32:1147–1160
Assenov Y, Ramirez F, Schelhorn SE, Lengauer T, Albrecht M (2008) Computing topological parameters of biological networks. Bioinformatics 24:282–284. https://doi.org/10.1093/bioinformatics/btm554
Ausin I et al (2012) INVOLVED IN DE NOVO 2-containing complex involved in RNA-directed DNA methylation in Arabidopsis. Proc Natl Acad Sci USA 109:8374–8381. https://doi.org/10.1073/pnas.1206638109
Bader GD, Hogue CW (2003) An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics 4:2
Bernal M, Casero D, Singh V, Wilson GT, Grande A, Yang H, Dodani SC, Pellegrini M, Huijser P, Connolly EL, Merchant SS, Kramer U (2012) Transcriptome sequencing identifies SPL7-regulated copper acquisition genes FRO4/FRO5 and the copper dependence of iron homeostasis in Arabidopsis. Plant cell 24:738-761. DOI 10.1105/tpc.111.090431
Bonnet E, He Y, Billiau K, Van de Peer Y (2010) TAPIR, a web server for the prediction of plant microRNA targets, including target mimics. Bioinformatics 26:1566–1568. https://doi.org/10.1093/bioinformatics/btq233
Chekanova JA (2015) Long non-coding RNAs and their functions in plants. Curr Opin Plant Biol 27:207–216. https://doi.org/10.1016/j.pbi.2015.08.003
Chen M, Wang C, Bao H, Chen H, Wang Y (2016) Genome-wide identification and characterization of novel lncRNAs in Populus under nitrogen deficiency. Mol Genet Genom : MGG 291:1663–1680. https://doi.org/10.1007/s00438-016-1210-3
Cheng CY, Krishnakumar V, Chan AP, Thibaud-Nissen F, Schobel S, Town CD (2017) Araport11: a complete reannotation of the Arabidopsis thaliana reference genome. Plant J 89:789–804. https://doi.org/10.1111/tpj.13415
Chung PJ, Jung H, Jeong D-H, Ha S-H, Do Choi Y, Kim J-K (2016) Transcriptome profiling of drought responsive noncoding RNAs and their target genes in rice. BMC Genomics 17:563
Cruz de Carvalho MH, Sun HX, Bowler C, Chua NH (2016) Noncoding and coding transcriptome responses of a marine diatom to phosphate fluctuations. New Phytol 210:497–510. https://doi.org/10.1111/nph.13787
Dai XB, Zhao PX (2011) psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res 39:W155–W159. https://doi.org/10.1093/nar/gkr319
Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, Lagarde J, Veeravalli L, Ruan X, Ruan Y, Lassmann T, Carninci P, Brown JB, Lipovich L, Gonzalez JM, Thomas M, Davis CA, Shiekhattar R, Gingeras TR, Hubbard TJ, Notredame C, Harrow J, Guigo R (2012) The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res 22:1775–1789
Di C et al (2014) Characterization of stress-responsive lncRNAs in Arabidopsis thaliana by integrating expression, epigenetic and structural features. Plant J 80:848–861. https://doi.org/10.1111/tpj.12679
Fang S et al (2018) NONCODEV5: a comprehensive annotation database for long non-coding. RNAs Nucleic Acids Res 46:D308–D314. https://doi.org/10.1093/nar/gkx1107
Franco-Zorrilla JM et al (2007) Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet 39:1033–1037. https://doi.org/10.1038/ng2079
Fu X-D (2014) Non-coding RNA: a new frontier in regulatory biology. Natl Sci Rev 1:190–204
Haak DC, Fukao T, Grene R, Hua Z, Ivanov R, Perrella G, Li S (2017) Multilevel Regulation of Abiotic Stress Responses in Plants. Front Plant Sci 8:1564. https://doi.org/10.3389/fpls.2017.01564
Hsieh TF, Ibarra CA, Silva P, Zemach A, Eshed-Williams L, Fischer RL, Zilberman D (2009) Genome-Wide Demethylation of Arabidopsis Endosperm. Science 324:1451–1454. https://doi.org/10.1126/science.1172417
Hu L, Xu ZY, Hu BQ, Lu ZJ (2017) COME: a robust coding potential calculation tool for lncRNA identification and characterization based on multiple features. Nucleic Acids Res 45:Artn E2. https://doi.org/10.1093/Nar/Gkw798
Jiang C-J, Shimono M, Sugano S, Kojima M, Liu X, Inoue H, Sakakibara H, Takatsuji H (2013) Cytokinins act synergistically with salicylic acid to activate defense gene expression in rice. Mol Plant-Microbe Interact 26:287–296
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL (2013) TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36. https://doi.org/10.1186/gb-2013-14-4-r36
Kozomara A, Birgaoanu M, Griffiths-Jones S (2019) miRBase: from microRNA sequences to function. Nucleic Acids Res 47:D155–D162. https://doi.org/10.1093/nar/gky1141
Kumar S, Verma S, Trivedi PK (2017) Involvement of small RNAs in phosphorus and sulfur sensing, signaling and stress: current update. Front Plant Sci 8:Artn 285. https://doi.org/10.3389/Fpls.2017.00285
Lan P, Li W, Schmidt W (2012) Complementary proteome and transcriptome profiling in phosphate-deficient Arabidopsis roots reveals multiple levels of gene regulation. Mol Cell Proteomics 11:1156-1166. DOI 10.1074/mcp.M112.020461
Li S et al (2017) Genome-wide identification and functional prediction of cold and/or drought-responsive lncRNAs in cassava. Sci Rep 7:45981
Li W, Lin WD, Ray P, Lan P, Schmidt W (2013) Genome-wide detection of condition-sensitive alternative splicing in Arabidopsis roots. Plant Physiol. 162:1750-1763. DOI 10.1104/pp.113.217778
Liang G, Ai Q, Yu DQ (2015) Uncovering miRNAs involved in crosstalk between nutrient deficiencies in Arabidopsis. Sci Rep-Uk 5:Artn 11813. https://doi.org/10.1038/Srep11813
Liu J et al (2012) Genome-wide analysis uncovers regulation of long Intergenic noncoding RNAs in Arabidopsis. Plant Cell 24:4333–4345. https://doi.org/10.1105/tpc.112.102855
Liu X, Hao L, Li D, Zhu L, Hu S (2015) Long non-coding RNAs and their biological roles in plants. Genomics Proteomics Bioinformatics 13:137–147. https://doi.org/10.1016/j.gpb.2015.02.003
Lu Q, Ren S, Lu M, Zhang Y, Zhu D, Zhang X, Li T (2013) Computational prediction of associations between long non-coding RNAs and proteins. BMC Genomics 14:651
Lv Y et al (2016) Genome-wide identification and functional prediction of nitrogen-responsive intergenic and intronic long non-coding RNAs in maize (Zea mays L.). BMC Genomics 17:350. https://doi.org/10.1186/s12864-016-2650-1
Maere S, Heymans K, Kuiper M (2005) BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21:3448–3449. https://doi.org/10.1093/bioinformatics/bti551
Meng J, Liu D, Luan YS (2015) Inferring plant microRNA functional similarity using a weighted protein-protein interaction network. BMC bioinformatics 16:Artn 360. https://doi.org/10.1186/S12859-015-0789-4
Ning S et al (2016) Construction of a lncRNA-mediated feed-forward loop network reveals global topological features and prognostic motifs in human cancers. Oncotarget. https://doi.org/10.18632/oncotarget.10004
Niu Y, Ahammed GJ, Tang C, Guo L, Yu J (2016) Physiological and Transcriptome Responses to Combinations of Elevated CO2 and Magnesium in Arabidopsis thaliana. PLoS One 11:e0149301. DOI 10.1371/journal.pone.0149301
Paytuvi Gallart A, Hermoso Pulido A, Anzar Martinez de Lagran I, Sanseverino W, Aiese Cigliano R (2016) GREENC: a Wiki-based database of plant lncRNAs. Nucleic Acids Res 44:D1161–D1166. https://doi.org/10.1093/nar/gkv1215
Qian W et al (2012a) A histone acetyltransferase regulates active DNA demethylation in Arabidopsis. Science 336:1445–1448
Qian WQ, Miki D, Zhang H, Liu Y, Zhang X, Tang K, Kan Y, la H, Li X, Li S, Zhu X, Shi X, Zhang K, Pontes O, Chen X, Liu R, Gong Z, Zhu JK (2012b) A Histone Acetyltransferase Regulates Active DNA Demethylation in Arabidopsis. Science 336:1445–1448. https://doi.org/10.1126/science.1219416
Roudier F et al (2011) Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. EMBO J 30:1928–1938. https://doi.org/10.1038/emboj.2011.103
Saeed AI et al (2003) TM4: a free, open-source system for microarray data management and analysis. Biotechniques 34:374−+
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP (2011) A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell 146:353–358. https://doi.org/10.1016/j.cell.2011.07.014
Sani E, Herzyk P, Perrella G, Colot V, Amtmann A (2013) Hyperosmotic priming of Arabidopsis seedlings establishes a long-term somatic memory accompanied by specific changes of the epigenome. Genome Biol. 14:R59. DOI 10.1186/gb-2013-14-6-r59
Shafiq S, Li J, Sun Q (2016) Functions of plants long non-coding RNAs. Biochim Biophys Acta 1859:155–162. https://doi.org/10.1016/j.bbagrm.2015.06.009
Shannon P et al (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13:2498–2504. https://doi.org/10.1101/gr.1239303
Shumayla SS, Taneja M, Tyagi S, Singh K, Upadhyay SK (2017) Survey of high throughput RNA-Seq data reveals potential roles for lncRNAs during development and stress response in bread wheat. Front Plant Sci 8:1019. https://doi.org/10.3389/fpls.2017.01019
Srivastava S, Srivastava AK, Suprasanna P, D'Souza SF (2013) Identification and profiling of arsenic stress-induced microRNAs in Brassica juncea. J Exp Bot 64:303–315. https://doi.org/10.1093/jxb/ers333
Stauffer E, Maizel A (2014) Post-transcriptional regulation in root development. Wiley Interdiscip Rev: RNA 5:679–696
Sun L et al (2013) Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res 41:ARTN e166. https://doi.org/10.1093/nar/gkt646
Supek F, Bosnjak M, Skunca N, Smuc T (2011) REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One 6:e21800. https://doi.org/10.1371/journal.pone.0021800
Szczesniak MW, Bryzghalov O, Ciomborowska-Basheer J, Makalowska I (2019) CANTATAdb 2.0: expanding the collection of plant long noncoding RNAs. Methods Mol Biol 1933:415–429. https://doi.org/10.1007/978-1-4939-9045-0_26
Tran L-SP, Urao T, Qin F, Maruyama K, Kakimoto T, Shinozaki K, Yamaguchi-Shinozaki K (2007) Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in Arabidopsis. Proc Natl Acad Sci 104:20623–20628
Trapnell C et al (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28:511–515. https://doi.org/10.1038/nbt.1621
Verma V, Ravindran P, Kumar PP (2016) Plant hormone-mediated regulation of stress responses. BMC Plant Biol 16:86. https://doi.org/10.1186/s12870-016-0771-y
Vidal EA, Moyano TC, Krouk G, Katari MS, Tanurdzic M, McCombie WR, Coruzzi GM, Gutierrez RA (2013) Integrated RNA-seq and sRNA-seq analysis identifies novel nitrate-responsive genes in Arabidopsis thaliana roots. BMC Genomics 14:701. DOI 10.1186/1471-2164-14-701
Villegas VE, Zaphiropoulos PG (2015) Neighboring gene regulation by antisense long non-coding RNAs. Int J Mol Sci 16:3251–3266
Wang L, Park HJ, Dasari S, Wang SQ, Kocher JP, Li W (2013) CPAT: coding-potential assessment tool using an alignment-free logistic regression model. Nucleic Acids Res 41:ARTN e74. https://doi.org/10.1093/nar/gkt006
Wang J, Meng X, Dobrovolskaya OB, Orlov YL, Chen M (2017) Non-coding RNAs and their roles in stress response in plants genomics. Proteomics Bioinformatics 15:301–312. https://doi.org/10.1016/j.gpb.2017.01.007
Winter D, Vinegar B, Nahal H, Ammar R, Wilson GV, Provart NJ (2007) An "electronic fluorescent pictograph" browser for exploring and analyzing large-scale biological data sets. PLoS One 2:e718. https://doi.org/10.1371/journal.pone.0000718
Xu J et al (2015) The mRNA related ceRNA-ceRNA landscape and significance across 20 major cancer types. Nucleic Acids Res 43:8169–8182. https://doi.org/10.1093/nar/gkv853
Xu XW, Zhou XH, Wang RR, Peng WL, An Y, Chen LL (2016) Functional analysis of long intergenic non-coding RNAs in phosphate-starved rice using competing endogenous. RNA network Sci rep-Uk 6:Artn 20715. https://doi.org/10.1038/Srep20715
Yamada M (2017) Functions of long intergenic non-coding (linc) RNAs in plants. J Plant Res 130:67–73. https://doi.org/10.1007/s10265-016-0894-0
Yanai I et al (2005) Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification. Bioinformatics 21:650–659. https://doi.org/10.1093/bioinformatics/bti042
Yuan C, Wang J, Harrison AP, Meng X, Chen D, Chen M (2015) Genome-wide view of natural antisense transcripts in Arabidopsis thaliana. DNA Res 22:233–243. https://doi.org/10.1093/dnares/dsv008
Yuan C, Meng X, Li X, Illing N, Ingle RA, Wang J, Chen M (2016a) PceRBase: a database of plant competing endogenous. RNA Nucleic Acids Res 45:D1009–D1014
Yuan JP, Zhang Y, Dong JS, Sun YZ, Lim BL, Liu D, Lu ZJ (2016b) Systematic characterization of novel lncRNAs responding to phosphate starvation in Arabidopsis thaliana. BMC Genomics 17:Artn 655. https://doi.org/10.1186/S12864-016-2929-2
Zeng HQ, Wang GP, Hu XY, Wang HZ, Du LQ ZYY (2014) Role of microRNAs in plant responses to nutrient stress. Plant Soil 374:1005–1021. https://doi.org/10.1007/s11104-013-1907-6
Zhang Z et al (2010) PMRD: plant microRNA database. Nucleic Acids Res 38:D806–D813. https://doi.org/10.1093/nar/gkp818
Zhang Y-C et al (2014) Genome-wide screening and functional analysis identify a large number of long noncoding RNAs involved in the sexual reproduction of rice. Genome Biol 15:512
Zhou M et al (2016) Characterization of long non-coding RNA-associated ceRNA network to reveal potential prognostic lncRNA biomarkers in human ovarian cancer. Oncotarget 7:12598–12611. https://doi.org/10.18632/oncotarget.7181
Acknowledgements
We are grateful to the other members of Ming Chen’s laboratory for helpful discussions and helpful comments. The authors are also grateful for the critical comments of anonymous reviewers.
Funding
This work was supported by the National Natural Sciences Foundation of China (No. 31771477).
Author information
Authors and Affiliations
Contributions
J.W. and Q.C. designed the study. W.W., Y.C., and Y. Z. collected the data and performed the analysis. J.W. drafted the manuscript. G.G. and M. C. completed the final version of the manuscript. All co-authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supporting Information
Supplementary Fig. S1
Heat map of differentially expressed lncRNAs fold change value pattern in different nutrient stress conditions. (PNG 2098 kb)
Supplementary Fig. S2
The ceRNA network under nutrient stress conditions. Blue nodes denote mRNAs, red nodes denote lncRNAs and orange nodes denote miRNAs that are represented as circles. (PNG 4901 kb)
Supplementary Fig. S3
Visualization of genome region for four isoforms of XLOC_014118. Blue bars refer to protein coding genes on the forward strand (Gene (+)) and reverse strand (Gene (-)). Red bars refer to lncRNA isoforms. Green bars refer to A. thaliana lncRNAs stored in NONCODE database. (PNG 2842 kb)
Supplementary Fig. S4
The DNA methylation pattern of XLOC_014118 and XLOC_001771. Red box presents the DNA methylation in 5kb upstream of TSS regions. (PNG 5069 kb)
ESM 1
Supplementary Table S1. The gff file of identified lncRNAs. Supplementary Table S2. Nutrient-responsive miRNAs in A. thaliana. Supplementary Table S3. List of 360 key lncRNAs in this study. (XLSX 282 kb)
Rights and permissions
About this article
Cite this article
Wang, J., Chen, Q., Wu, W. et al. Genome-wide analysis of long non-coding RNAs responsive to multiple nutrient stresses in Arabidopsis thaliana. Funct Integr Genomics 21, 17–30 (2021). https://doi.org/10.1007/s10142-020-00758-5
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10142-020-00758-5