Pi-starvation induced transcriptional changes in barley revealed by a comprehensive RNA-Seq and degradome analyses

Abstract

Background: Small RNAs (sRNAs) are 18–24 nt regulatory elements which are responsible for plant development regulation and participate in many plant stress responses. Insufficient inorganic phosphate (Pi) concentration triggers plant responses to balance the internal Pi level. Results: In this study, we describe Pi-starvation-responsive small RNAs and transcriptome changes in barley (Hordeum vulgare L.) using Next-Generation Sequencing (NGS) data derived from three different types of NGS libraries: (i) small RNAs, (ii) degraded RNAs, and (iii) functional mRNAs. We find that differentially and significantly expressed miRNAs (DEMs, p-value < 0.05) are represented by 162 (44.88 % of total differentially expressed small RNAs) molecules in shoot and 138 (7.14 %) in root; mainly various miR399 and miR827 isomiRs. The remaining small RNAs (i.e., those without perfect match to reference sequences deposited in miRBase) are considered as differentially expressed other sRNAs (DESs, Bonferroni correction). In roots, a more abundant and diverse set of other sRNAs (1796 unique sequences, 0.13 % from total unique reads obtained under low-Pi) contributes more to the compensation of low-Pi stress than that in shoots (199 unique sequences, 0.01 %). More than 80 % of differentially expressed other sRNAs are upregulated in both organs. Additionally, in barley shoots, upregulation of small RNAs is accompanied by strong induction of two nucleases (S1/P1 endonuclease and 3’-5’ exonuclease). This suggests that most small RNAs may be generated upon endonucleolytic cleavage to increase the internal Pi pool. Transcriptomic profiling of Pi-starved barley shoots identify 98 differentially expressed genes (DEGs). A majority of the DEGs possess characteristic Pi-responsive cis-regulatory elements (P1BS and/or PHO element), located mostly in the proximal promoter regions. GO analysis shows that the discovered DEGs primarily alter plant defense, plant stress response, nutrient mobilization, or pathways involved in the gathering and recycling of phosphorus from organic pools.Conclusions: Our results provide comprehensive data to demonstrate complex responses at the RNA level in barley to maintain Pi homeostasis and indicate that barley adapts to Pi scarcity through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.