Transcriptome Analysis of Zygophyllum xanthoxylum Adaptation Strategies to Phosphate Stress

Soil phosphate (Pi) deficiency is a global issue and a major constraint on plant growth. Plants typically acclimatize to low Pi by enhancing their P utilization and/or P acquisition efficiencies; however, different species have variable preferred strategies. RNA sequencing analysis was performed on the shoots and roots of Zygophyllum xanthoxylum, under 1 day and 10 days of Pi stress, to investigate their adaptation strategies to P deprivation. A total of 364,614 unigenes and 9,270 differentially expressed genes (DEGs) were obtained via transcriptome sequencing. An analysis of the DEGs revealed that under the 10D treatment, anthocyanin synthesis genes were upregulated under Pi stress, whereas gibberellin, ethylene, and cytokinins synthesis genes were upregulated, and abscisic acid synthesis genes were downregulated. Genes related to organic acid synthesis, encoding for purple acid phosphatases (APase) and nucleases (RNase) were upregulated under the 1D and 10D treatments, respectively. Furthermore, genes associated with Pi transport were induced by Pi stress. Zygophyllum xanthoxylum has special P adaptation strategies, the variation trends of genes involved in external P mobilization and acquisition, which were different from that of most other species; however, the expression levels of organophosphorus mobilization related genes, such as APases and RNases, were significantly increased. Meanwhile, PHT2s and TPTs, which distributed Pi to effective sites (e.g., chloroplast), played critical roles in the maintenance of photosynthesis. We speculated that these were economic and energy saving strategies, and there are critical adaptive mechanisms that Z. xanthoxylum employs to cope with deficits in Pi.

The ACID PHOSPHATASE 1 regulates Pi stress adaptation by maintaining intracellular Pi homeostasis

The concentration and homeostasis of intracellular phosphate (Pi) are crucial for sustaining cell metabolism and growth. During short-term Pi starvation, intracellular Pi is maintained relatively constant at the expense of vacuolar Pi. After the vacuolar stored Pi is exhausted, the plant cells induce the synthesis of intracellular acid phosphatase (APase) to recycle Pi from expendable organic phosphate (Po). In this study, the expression, enzymatic activity and subcellular localization of ACID PHOSPHATASE 1 (OsACP1) were determined. OsACP1 expression is specifically induced in almost all cell types of leaves and roots under Pi stress conditions. OsACP1 encodes an acid phosphatase with broad Po substrates and localizes in the endoplasmic reticulum (ER) and Golgi apparatus (GA). Phylogenic analysis demonstrates that OsACP1 has a similar structure with human acid phosphatase PHOSPHO1. Overexpression or mutation of OsACP1 affected Po degradation and utilization, which further influenced plant growth and productivity under both Pi-sufficient and Pi-deficient conditions. Moreover, overexpression of OsACP1 significantly affected intracellular Pi homeostasis and Pi starvation signalling. We concluded that OsACP1 is an active acid phosphatase that regulates rice growth under Pi stress conditions by recycling Pi from Po in the ER and GA.

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

Background Small RNAs (sRNAs) are 20–30 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) RNA-Seq 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, Bonferroni adjusted p-value < 0.05) are represented by 15 molecules in shoot and 13 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, p-value Bonferroni correction < 0.05). In roots, a more abundant and diverse set of other sRNAs (DESs, 1796 unique sequences, 0.13% from the average of the unique small RNA expressed under low-Pi) contributes more to the compensation of low-Pi stress than that in shoots (DESs, 199 unique sequences, 0.01%). More than 80% of differentially expressed other sRNAs are up-regulated in both organs. Additionally, in barley shoots, up-regulation 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 nucleolytic cleavage to increase the internal Pi pool. Transcriptomic profiling of Pi-starved barley shoots identifies 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-starvation through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.

Genome-wide analysis of haloacid dehalogenase genes reveals their function in phosphate starvation responses in rice

The HAD superfamily is named after the halogenated acid dehalogenase found in bacteria, which hydrolyses a diverse range of organic phosphate substrates. Although certain studies have shown the involvement of HAD genes in Pi starvation responses, systematic classification and bioinformatics analysis of the HAD superfamily in plants is lacking. In this study, 41 and 40 HAD genes were identified by genomic searching in rice and Arabidopsis, respectively. According to sequence similarity, these proteins are divided into three major groups and seven subgroups. Conserved motif analysis indicates that the majority of the identified HAD proteins contain phosphatase domains. A further structural analysis showed that HAD proteins have four conserved motifs and specified cap domains. Fewer HAD genes show collinearity relationships in both rice and Arabidopsis, which is consistent with the large variations in the HAD genes. Among the 41 HAD genes of rice, the promoters of 28 genes contain Pi-responsive cis-elements. Mining of transcriptome data and qRT-PCR results showed that at least the expression of 17 HAD genes was induced by Pi starvation in shoots or roots. These HAD proteins are predicted to be involved in intracellular or extracellular Po recycling under Pi stress conditions in plants.

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

Background: Small RNAs (sRNAs) are 20–30 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) RNA-Seq 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, p-value Bonferroni correction < 0.05). In roots, a more abundant and diverse set of other sRNAs (DESs, 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 (DESs, 199 unique sequences, 0.01 %). More than 80 % of differentially expressed other sRNAs are up-regulated in both organs. Additionally, in barley shoots, up-regulation 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 identifies 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-starvation through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.

Differential Gene Expression Responding to Low Phosphate Stress in Leaves and Roots of Maize by cDNA-SRAP

Phosphate (Pi) deficiency in soil can have severe impacts on the growth, development, and production of maize worldwide. In this study, a cDNA-sequence-related amplified polymorphism (cDNA-SRAP) transcript profiling technique was used to evaluate the gene expression in leaves and roots of maize under Pi stress for seven days. A total of 2494 differentially expressed fragments (DEFs) were identified in response to Pi starvation with 1202 and 1292 DEFs in leaves and roots, respectively, using a total of 60 primer pairs in the cDNA-SRAP analysis. These DEFs were categorized into 13 differential gene expression patterns. Results of sequencing and functional analysis showed that 63 DEFs (33 in leaves and 30 in roots) were annotated to a total of 54 genes involved in diverse groups of biological pathways, including metabolism, photosynthesis, signal transduction, transcription, transport, cellular processes, genetic information, and organismal system. This study demonstrated that (1) the cDNA-SRAP transcriptomic profiling technique is a powerful method to analyze differential gene expression in maize showing advantageous features among several transcriptomic methods; (2) maize undergoes a complex adaptive process in response to low Pi stress; and (3) a total of seven differentially expressed genes were identified in response to low Pi stress in leaves or roots of maize and could be used in the genetic modification of maize.

Integrated mRNA and miRNA expression profiling analyses of Pinus massoniana roots and shoots in response to phosphate deficiency

Background Masson pine ( Pinus massoniana ) is primarily present in subtropical and tropical areas of China, which are severely deficient in inorganic phosphate (Pi). Although some studies identified transcriptomic and proteomic responses to Pi deficiency in Masson pine seedlings, different tissues, especially the roots, that exhibit primary responses to low-Pi stress, have not been well studied. To shed further light on the complex responses of Masson pine to Pi deficiency, a spatiotemporal experiment was performed to identify differentially expressed mRNAs and miRNAs under Pi-deficient conditions. Results Spatiotemporal analyses of 72 RNA sequencing libraries provided a comprehensive overview of the dynamic responses of Masson pine to low-Pi stress. Differentially expressed gene analysis revealed several high-affinity phosphate transporters and nitrate transporters, reflecting the crosstalk between nitrate and Pi homeostasis in plants. The MYB family was the most abundant transcription factor family identified. miRNA differential expression analysis identified several families that were associated with Pi deficiency, such as miR399. In addition, some other families, including novel miRNA families in Masson pine, were dramatically changed in response to Pi starvation. GO and KEGG analyses of these mRNAs and targets of miRNAs indicated that metabolic processes were most enriched under Pi deficiency. Conclusions This study provided abundant spatiotemporal transcriptomic information to functionally dissect the response of Masson pine seedlings to Pi deficiency, which will aid in further elucidation of the biological regulatory mechanisms of pines in response to low-Pi stress.

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

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.