PHOSPHATE1-mediated phosphate translocation from roots to shoots regulates floral transition in plants.

Phosphorus nutrition has been known for a long time to influence floral transition in plants, but the underlying mechanism is unclear. Arabidopsis phosphate transporter PHOSPHATE1 (PHO1) plays a critical role in phosphate translocation from roots to shoots, but whether and how it regulates floral transition is unknown. Here, we show that knockout mutation of PHO1 delays flowering under both long- and short-day conditions. The late flowering of pho1 mutants can be partially rescued by Pi supplementation in rosettes or shoot apices. Grafting assay indicates that the late flowering of pho1 mutants is a result of impaired phosphate translocation from roots to shoots. Knockout mutation of SPX1 and SPX2, two negative regulators of the phosphate starvation response, partially rescues the late flowering of pho1 mutants. PHO1 is epistatic to PHO2, a negative regulator of PHO1, in flowering time regulation. Loss of PHO1 represses the expression of some floral activators, including FT encoding florigen, and induces the expression of some floral repressors in shoots. Genetic analyses indicate that at least jasmonic acid signaling is partially responsible for the late flowering of pho1 mutants. In addition, we find that rice PHO1;2, the homolog of PHO1, plays a similar role in floral transition. These results suggest that PHO1 integrates phosphorus nutrition and flowering time, and could be used as a potential target in modulating phosphorus nutrition-mediated flowering time in plants.

Stomatal properties of Arabidopsis cauline and rice flag leaves and their contributions to seed production and grain yield.

Cauline leaves on the inflorescence stem of Arabidopsis thaliana may play important roles in supplying photosynthetic products to sinks, such as floral organs. Flag leaves in rice (Oryza sativa) have a higher photosynthetic capacity than other leaves, and are crucial for increasing grain yield. However, the detailed properties of stomata in cauline and flag leaves have not been investigated. In Arabidopsis, stomatal conductance and CO2 assimilation rate were higher in cauline leaves under white light than in rosette leaves, consistent with higher levels of plasma membrane (PM) H+-ATPase, a key enzyme for stomatal opening, in guard cells. Moreover, removal of cauline leaves significantly reduced the shoot biomass by approximately 20% and seed production by approximately 46%. In rice, higher stomatal density, stomatal conductance, and CO2 assimilation rate were observed in flag leaves than in fully expanded second leaves. Removal of the flag leaves significantly reduced grain yield by approximately 49%. Taken together, these results show that cauline and flag leaves have important roles in seed production and grain yield through enhanced stomatal conductance and CO2 assimilation rate.

Adsorptive decontamination of antibiotics from livestock wastewater by using alkaline-modified biochar.

Efficient abatement of antibiotics from livestock wastewater is in urgent demand, but still challenging. In this study, alkaline-modified biochar with larger surface area (130.520 m2 g-1) and pore volume (0.128 cm3 g-1) was fabricated and explored for the adsorption of different types of antibiotics from livestock wastewater. Batch adsorption experiments demonstrated that the adsorption process was mainly determined by chemisorption and was heterogeneous, which could be moderately affected by the variations of solution pH (3-10). Furthermore, the computational analysis based on density functional theory (DFT) indicated that the -OH groups on biochar surface could serve as the dominant active sites for antibiotics adsorption due to the strongest adsorption energies between antibiotics and -OH groups. In addition, the antibiotics removal was also evaluated in multi-pollutants system, where biochar performed synergistic adsorption towards Zn2+/Cu2+ and antibiotics. Overall, these findings not only deepen our understandings on the adsorption mechanism between biochar and antibiotics, but also promote the application of biochar in the remediation of livestock wastewater.

Analyses of Botrytis cinerea-responsive LrWRKY genes from Lilium regale reveal distinct roles of two LrWRKY transcription factors in mediating responses to B. cinerea.

KEY MESSAGE: Botrytis cinerea induced expression of 15 LrWRKY genes; overexpression of LrWRKY39 and LrWRKY41a increased resistance and susceptibility, respectively, to B. cinerea in a manner related to SA and JA signaling. WRKY transcription factors (TFs), a large family, play important roles in coping with biotic stresses. Lilium regale Wilson is a lily species with strong resistance to fungi and viruses; however, functional characterization of LrWRKY TFs remains very limited. Here, a total of 25 LrWRKY members were identified from the L. regale transcriptome, and 15 LrWRKY genes were significantly induced by Botrytis cinerea. Based on their structural features, B. cinerea-responsive LrWRKY genes could be classified into six subgroups (Groups I, IIa-d, and III), and sequence alignment showed that 12 LrWRKY proteins have a well-conserved WRKYGQK domain, while 3 LrWRKYs have a variant sequence (WRKYGKK or WRMYEQK). Quantitative RT-PCR analysis revealed tissue-specific expression of B. cinerea-responsive LrWRKY genes and their expression profiles in response to defense-related hormones salicylic acid (SA), methyl jasmonate (MeJA) and hydrogen peroxide. LrWRKY39 and LrWRKY41a, which encode two LrWRKY TFs with different three-dimensional (3D) models of the WRKY domain, were cloned, and both proteins were targeted to the nucleus. Overexpression of LrWRKY39 and LrWRKY41a in Arabidopsis thaliana increased the resistance and susceptibility to B. cinerea, respectively, compared to the wild type. Similar results were also observed in tobacco and lily (L. longiflorum 'Snow Queen') by transient transformation analyses. Their distinct roles may be related to changes in the transcript levels of SA-/JA-responsive genes. Our results provide new insights into B. cinerea-responsive LrWRKY members and the biological functions of two different 3D models of LrWRKY TFs in defense responses to B. cinerea infection.

Overexpression of a Plasma Membrane H+-ATPase Gene OSA1 Stimulates the Uptake of Primary Macronutrients in Rice Roots.

Plasma membrane (PM) H+-ATPase is a master enzyme involved in various plant physiological processes, such as stomatal movements in leaves and nutrient uptake and transport in roots. Overexpression of Oryza sativa PM H+-ATPase 1 (OSA1) has been known to increase NH4+ uptake in rice roots. Although electrophysiological and pharmacological experiments have shown that the transport of many substances is dependent on the proton motive force provided by PM H+-ATPase, the exact role of PM H+-ATPase on the uptake of nutrients in plant roots, especially for the primary macronutrients N, P, and K, is still largely unknown. Here, we used OSA1 overexpression lines (OSA1-oxs) and gene-knockout osa1 mutants to investigate the effect of modulation of PM H+-ATPase on the absorption of N, P, and K nutrients through the use of a nutrient-exhaustive method and noninvasive microtest technology (NMT) in rice roots. Our results showed that under different concentrations of P and K, the uptake rates of P and K were enhanced in OSA1-oxs; by contrast, the uptake rates of P and K were significantly reduced in roots of osa1 mutants when compared with wild-type. In addition, the net influx rates of NH4+ and K+, as well as the efflux rate of H+, were enhanced in OSA1-oxs and suppressed in osa1 mutants under low concentration conditions. In summary, this study indicated that overexpression of OSA1 stimulated the uptake rate of N, P, and K and promoted flux rates of cations (i.e., H+, NH4+, and K+) in rice roots. These results may provide a novel insight into improving the coordinated utilization of macronutrients in crop plants.

Proteomic Analysis Demonstrates a Molecular Dialog Between Trichoderma guizhouense NJAU 4742 and Cucumber (Cucumis sativus L.) Roots: Role in Promoting Plant Growth.

Trichoderma is a genus of filamentous fungi that play notable roles in stimulating plant growth after colonizing the root surface. However, the key proteins and molecular mechanisms governing this stimulation have not been completely elucidated. In this study, Trichoderma guizhouense NJAU 4742 was investigated in a hydroponic culture system after interacting with cucumber roots. The total proteins of the fungus were characterized, and the key metabolic pathways along with related genes were analyzed through proteomic and transcriptomic analyses. The roles played by the regulated proteins during the interaction between plants and NJAU 4742 were further examined. The intracellular or extracellular proteins from NJAU 4742 and extracellular proteins from cucumber were quantified, and the high-abundance proteins were determined which were primarily involved in the shikimate pathway (tryptophan, tyrosine, and phenylalanine metabolism, auxin biosynthesis, and secondary metabolite synthesis). Moreover, 15N-KNO3 labeling analysis indicated that NJAU 4742 had a strong ability to convert nitrogenous amino acids, nitrate, nitrile, and amines into ammonia. The auxin synthesis and ammonification metabolism pathways of NJAU 4742 significantly contributed to plant growth. The results of this study demonstrated the crucial metabolic pathways involved in the interactions between Trichoderma spp. and plants.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A mycorrhiza-specific H+ -ATPase is essential for arbuscule development and symbiotic phosphate and nitrogen uptake.

Most land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi to enhance uptake of mineral nutrients, particularly phosphate (Pi) and nitrogen (N), from the soil. It is established that transport of Pi from interfacial apoplast into plant cells depends on the H+ gradient generated by the H+ -ATPase located on the periarbuscular membrane (PAM); however, little evidence regarding the potential link between mycorrhizal N transport and H+ -ATPase activity is available to date. Here, we report that a PAM-localized tomato H+ -ATPase, SlHA8, is indispensable for arbuscule development and mycorrhizal P and N uptake. Knockout of SlHA8 resulted in truncated arbuscule morphology, reduced shoot P and N accumulation, and decreased H+ -ATPase activity and acidification of apoplastic spaces in arbusculated cells. Overexpression of SlHA8 in tomato promoted both P and N uptake, and increased total colonization level, but did not affect arbuscule morphology. Heterogeneous expression of SlHA8 in the rice osha1 mutant could fully complement its defects in arbuscule development and mycorrhizal P and N uptake. Our results propose a pivotal role of the SlHA8 in energizing both the symbiotic P and N transport, and highlight the evolutionary conservation of the AM-specific H+ -ATPase orthologs in maintaining AM symbiosis across different mycorrhizal plant species.

Integrated analyses of miRNAome and transcriptome reveal zinc deficiency responses in rice seedlings.

BACKGROUND: Zinc (Zn) deficiency is one of the most widespread soil constraints affecting rice productivity, but the molecular mechanisms underlying the regulation of Zn deficiency response is still limited. Here, we aim to understand the molecular mechanisms of Zn deficiency response by integrating the analyses of the global miRNA and mRNA expression profiles under Zn deficiency and resupply in rice seedlings by integrating Illumina's high-throughput small RNA sequencing and transcriptome sequencing. RESULTS: The transcriptome sequencing identified 360 genes that were differentially expressed in the shoots and roots of Zn-deficient rice seedlings, and 97 of them were recovered after Zn resupply. A total of 68 miRNAs were identified to be differentially expressed under Zn deficiency and/or Zn resupply. The integrated analyses of miRNAome and transcriptome data showed that 12 differentially expressed genes are the potential target genes of 10 Zn-responsive miRNAs such as miR171g-5p, miR397b-5p, miR398a-5p and miR528-5p. Some miRNA genes and differentially expressed genes were selected for validation by quantitative RT-PCR, and their expressions were similar to that of the sequencing results. CONCLUSION: These results provide insights into miRNA-mediated regulatory pathways in Zn deficiency response, and provide candidate genes for genetic improvement of Zn deficiency tolerance in rice.

Boron Alleviates Aluminum Toxicity by Promoting Root Alkalization in Transition Zone via Polar Auxin Transport.

Boron (B) alleviates aluminum (Al) toxicity in higher plants; however, the underlying mechanisms behind this phenomenon remain unknown. Here, we used bromocresol green pH indicator, noninvasive microtest, and microelectrode ion flux estimation techniques to demonstrate that B promotes root surface pH gradients in pea (Pisum sativum) roots, leading to alkalization in the root transition zone and acidification in the elongation zone, while Al inhibits these pH gradients. B significantly decreased Al accumulation in the transition zone (∼1.0-2.5 mm from the apex) of lateral roots, thereby alleviating Al-induced inhibition of root elongation. Net indole acetic acid (IAA) efflux detected by an IAA-sensitive platinum microelectrode showed that polar auxin transport, which peaked in the root transition zone, was inhibited by Al toxicity, while it was partially recovered by B. Electrophysiological experiments using the Arabidopsis (Arabidopsis thaliana) auxin transporter mutants (auxin resistant1-7; pin-formed2 [pin2]) and the specific polar auxin transporter inhibitor1-naphthylphthalamic acid showed that PIN2-based polar auxin transport is involved in root surface alkalization in the transition zone. Our results suggest that B promotes polar auxin transport driven by the auxin efflux transporter PIN2 and leads to the downstream regulation of the plasma membrane-H+-ATPase, resulting in elevated root surface pH, which is essential to decrease Al accumulation in this Al-targeted apical root zone. These findings provide a mechanistic explanation for the role of exogenous B in alleviation of Al accumulation and toxicity in plants.

Overexpression of rice aquaporin OsPIP1;2 improves yield by enhancing mesophyll CO2 conductance and phloem sucrose transport.

Aquaporins are involved in CO2 transport from the leaf intercellular air space to the chloroplast, which contributes to CO2 assimilation. However, the mechanism of CO2 transport by rice (Oryza sativa L.) aquaporins is unknown. Here, we investigated the function of the aquaporin OsPIP1;2 in CO2 diffusion-associated photosynthesis and phloem sucrose transport. Moreover, the grain yield of rice lines overexpressing OsPIP1;2 was determined. OsPIP1;2 was localized to the plasma membrane and the relative expression of OsPIP1;2 was approximately 5-fold higher in leaves in the presence of an elevated CO2 concentration. Overexpression of OsPIP1;2 increased mesophyll conductance by approximately 150% compared with wild-type (WT) rice. The OsPIP1;2-overexpressing lines had higher biomass than the WT, possibly due to increased phloem sucrose transport. In addition, the grain yield of OsPIP1;2-overexpressing lines was approximately 25% higher than that of the WT in three-season field experiments, due to the increased numbers of effective tillers and spikelets per panicle. Our results suggest that OsPIP1;2 modulates rice growth and grain yield by facilitating leaf CO2 diffusion, which increases both the net CO2 assimilation rate and sucrose transport.

Transcriptome profiles of soybean leaves and roots in response to zinc deficiency.

Zinc (Zn) deficiency is a widespread agricultural problem in arable soils of the whole world. However, the molecular mechanisms underlying Zn-deficiency response are largely unknown. Here, we analyzed the transcriptomic profilings of soybean leaves and roots in response to Zn deficiency through Illumina's high-throughput RNA sequencing in order to understand the molecular basis of Zn-deficiency response in the plants. A total of 614 and 1011 gene loci were found to be differentially expressed in leaves and roots, respectively, and 88 loci were commonly found in both leaves and roots. Twelve differentially expressed genes (DEGs) were randomly selected for validation by quantitative reverse transcription polymerase chain reaction, and their fold changes were similar to those of RNA-seq. Gene ontology enrichment analysis showed that ion transport, nicotianamine (NA) biosynthetic process and queuosine biosynthetic process were enriched in the upregulated genes, while oxidation-reduction process and defense response were enriched in the downregulated genes. Among the DEGs, 20 DEGs are potentially involved in Zn homeostasis, including seven ZRT, IRT-related protein (ZIP) transporter genes, three NA synthase genes, and seven metallothionein genes; 40 DEGs are possibly involved in diverse hormonal signals such as auxin, cytokinin, ethylene and gibberellin; nine DEGs are putatively involved in calcium signaling; 85 DEGs are putative transcription factor genes. Nine DEGs were found to contain zinc-deficiency-response element in their promoter regions. These results could provide comprehensive insights into the soybean response to Zn deficiency and will be helpful for further elucidation of the molecular mechanisms of Zn-deficiency response and Zn-deficiency tolerance in plants.

Early Transcriptomic Response to Phosphate Deprivation in Soybean Leaves as Revealed by RNA-Sequencing.

Low phosphate (Pi) availability is an important limiting factor affecting soybean production. However, the underlying molecular mechanisms responsible for low Pi stress response and tolerance remain largely unknown, especially for the early signaling events under low Pi stress. Here, a genome-wide transcriptomic analysis in soybean leaves treated with a short-term Pi-deprivation (24 h) was performed through high-throughput RNA sequencing (RNA-seq) technology. A total of 533 loci were found to be differentially expressed in response to Pi deprivation, including 36 mis-annotated loci and 32 novel loci. Among the differentially expressed genes (DEGs), 303 were induced and 230 were repressed by Pi deprivation. To validate the reliability of the RNA-seq data, 18 DEGs were randomly selected and analyzed by quantitative RT-PCR (reverse transcription polymerase chain reaction), which exhibited similar fold changes with RNA-seq. Enrichment analyses showed that 29 GO (Gene Ontology) terms and 8 KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were significantly enriched in the up-regulated DEGs and 25 GO terms and 16 KEGG pathways were significantly enriched in the down-regulated DEGs. Some DEGs potentially involved in Pi sensing and signaling were up-regulated by short-term Pi deprivation, including five SPX-containing genes. Some DEGs possibly associated with water and nutrient uptake, hormonal and calcium signaling, protein phosphorylation and dephosphorylation and cell wall modification were affected at the early stage of Pi deprivation. The cis-elements of PHO (phosphatase) element, PHO-like element and P responsive element were present more frequently in promoter regions of up-regulated DEGs compared to that of randomly-selected genes in the soybean genome. Our transcriptomic data showed an intricate network containing transporters, transcription factors, kinases and phosphatases, hormone and calcium signaling components is involved in plant responses to early Pi deprivation.

The Nitrification Inhibitor Methyl 3-(4-Hydroxyphenyl)Propionate Modulates Root Development by Interfering with Auxin Signaling via the NO/ROS Pathway.

Methyl 3-(4-hydroxyphenyl)propionate (MHPP) is a root exudate that functions as a nitrification inhibitor and as a modulator of the root system architecture (RSA) by inhibiting primary root (PR) elongation and promoting lateral root formation. However, the mechanism underlying MHPP-mediated modulation of the RSA remains unclear. Here, we report that MHPP inhibits PR elongation in Arabidopsis (Arabidopsis thaliana) by elevating the levels of auxin expression and signaling. MHPP induces an increase in auxin levels by up-regulating auxin biosynthesis, altering the expression of auxin carriers, and promoting the degradation of the auxin/indole-3-acetic acid family of transcriptional repressors. We found that MHPP-induced nitric oxide (NO) production promoted reactive oxygen species (ROS) accumulation in root tips. Suppressing the accumulation of NO or ROS alleviated the inhibitory effect of MHPP on PR elongation by weakening auxin responses and perception and by affecting meristematic cell division potential. Genetic analysis supported the phenotype described above. Taken together, our results indicate that MHPP modulates RSA remodeling via the NO/ROS-mediated auxin response pathway in Arabidopsis. Our study also revealed that MHPP significantly induced the accumulation of glucosinolates in roots, suggesting the diverse functions of MHPP in modulating plant growth, development, and stress tolerance in plants.

Arabidopsis plasma membrane H+-ATPase genes AHA2 and AHA7 have distinct and overlapping roles in the modulation of root tip H+ efflux in response to low-phosphorus stress.

Phosphorus deficiency in soil is one of the major limiting factors for plant growth. Plasma membrane H+-ATPase (PM H+-ATPase) plays an important role in the plant response to low-phosphorus stress (LP). However, few details are known regarding the action of PM H+-ATPase in mediating root proton (H+) flux and root growth under LP. In this study, we investigated the involvement and function of different Arabidopsis PM H+-ATPase genes in root H+ flux in response to LP. First, we examined the expressions of all Arabidopsis PM H+-ATPase gene family members (AHA1-AHA11) under LP. Expression of AHA2 and AHA7 in roots was enhanced under this condition. When the two genes were deficient in their respective Arabidopsis mutant plants, root growth and responses of the mutants to LP were highly inhibited compared with the wild-type plant. AHA2-deficient plants exhibited reduced primary root elongation and lower H+ efflux in the root elongation zone. AHA7-deficient plants exhibited reduced root hair density and lower H+ efflux in the root hair zone. The modulation of H+ efflux by AHA2 or AHA7 was affected by the action of 14-3-3 proteins and/or auxin regulatory pathways in the context of root growth and response to LP. Our results suggest that under LP conditions, AHA2 acts mainly to modulate primary root elongation by mediating H+ efflux in the root elongation zone, whereas AHA7 plays an important role in root hair formation by mediating H+ efflux in the root hair zone.

Plasma membrane H+-ATPases in mineral nutrition and crop improvement.

Plasma membrane H+-ATPases (PMAs) pump H+ out of the cytoplasm by consuming ATP to generate a membrane potential and proton motive force for the transmembrane transport of nutrients into and out of plant cells. PMAs are involved in nutrient acquisition by regulating root growth, nutrient uptake, and translocation, as well as the establishment of symbiosis with arbuscular mycorrhizas. Under nutrient stresses, PMAs are activated to pump more H+ and promote organic anion excretion, thus improving nutrient availability in the rhizosphere. Herein we review recent progress in the physiological functions and the underlying molecular mechanisms of PMAs in the efficient acquisition and utilization of various nutrients in plants. We also discuss perspectives for the application of PMAs in improving crop production and quality.

Spatial-temporal analysis of zinc homeostasis reveals the response mechanisms to acute zinc deficiency in Sorghum bicolor.

Zinc (Zn) is an essential micronutrient in plants. The activity of copper/zinc superoxide dismutase (CSD) and carbonic anhydrase (CA) correlate with differences in Zn efficiency in plants; therefore, it is reasonable to hypothesize the existence of a Zn economy model that saves Zn for these essential Zn proteins during Zn deficiency. However, up to this point, direct evidence for the idea that CSD and/or CA might be priorities for Zn delivery has been lacking. Here, we investigated the spatial-temporal effects of acute Zn depletion and resupply by integrating physiological studies and molecular analyses using hydroponically grown Sorghum. The elevated expression of miR398 repressed CSD expression in roots, whereas the reduced expression of miR528 resulted in a relatively stable level of CSD expression in Sorghum leaves under Zn depletion. Spatial-temporal analysis after Zn resupply to previously depleted plants revealed that the expression and activity of CA were the first to recover after Zn addition, whereas the recovery of the activities of CSD and alcohol dehydrogenase (ADH) was delayed, suggesting that CA receives priority in Zn delivery over CSD and ADH. Our results also indicate that microRNAs (miRNAs) are important regulators of the response of Zn deficiency in plants.

Measurement of Total Phosphorus and Polyphosphate in Chlamydomonas reinhardtii.

Phosphorus is an essential nutrient for plants. Green algae usually store excess P as polyphosphate (polyP) in the vacuoles. PolyP, a linear chain of three to hundreds of phosphate residues linked by phosphoanhydride bonds, is important for cell growth. Based on the previous method of polyP purification with silica gel columns (Werner et al., 2005; Canadell et al., 2016) in yeast cells, we developed a protocol to purify and determine the total P and polyP in Chlamydomonas reinhardtii by a quick, simplified, and quantitative method. We use hydrochloric acid or nitric acid to digest polyP or total P in dried cells and analyze P content using the malachite green colorimetric method. This method may be applied to other microalgae.

Aquaporins and CO2 diffusion across biological membrane.

Despite the physiological significance of effective CO2 diffusion across biological membranes, the underlying mechanism behind this process is not yet resolved. Particularly debatable is the existence of CO2-permeable aquaporins. The lipophilic characteristic of CO2 should, according to Overton's rule, result in a rapid flux across lipid bilayers. However, experimental evidence of limited membrane permeability poses a challenge to this idea of free diffusion. In this review, we summarized recent progress with regard to CO2 diffusion, and discussed the physiological effects of altered aquaporin expression, the molecular mechanisms of CO2 transport via aquaporins, and the function of sterols and other membrane proteins in CO2 permeability. In addition, we highlight the existing limits in measuring CO2 permeability and end up with perspectives on resolving such argument either by determining the atomic resolution structure of CO2 permeable aquaporins or by developing new methods for measuring permeability.

Comparative transcriptome analyses under individual and combined nutrient starvations provide insights into N/P/K interactions in rice.

Crops often suffer from simultaneous limitations of multiple nutrients in soils, including nitrogen (N), phosphorus (P) and potassium (K), which are three major macronutrients essential for ensuring growth and yield. Although plant responses to individual N, P, and K deficiency have been well documented, our understanding of the responses to combined nutrient deficiencies and the crosstalk between nutrient starvation responses is still limited. Here, we compared the physiological responses in rice under seven kinds of single and multiple low nutrient stress of N, P and K, and used RNA sequencing approaches to compare their transcriptome changes. A total of 13,000 genes were found to be differentially expressed under all these single and multiple low N/P/K stresses, and 66 and 174 of them were shared by all these stresses in roots and shoots, respectively. Functional enrichment analyses of the DEGs showed that a group of biological and metabolic processes were shared by these low N/P/K stresses. Comparative analyses indicated that DEGs under multiple low nutrient stress was not the simple summation of single nutrient stress. N was found to be the predominant factor affecting the transcriptome under combined nutrient stress. N, P, or K availability exhibited massive influences on the transcriptomic responses to starvation of other nutrients. Many genes involved in nutrient transport, hormone signaling, and transcriptional regulation were commonly responsive to low N/P/K stresses. Some transcription factors were predicted to regulate the expression of genes that are commonly responsive to N, P, and K starvations. These results revealed the interactions between N, P, and K starvation responses, and will be helpful for further elucidation of the molecular mechanisms underlying nutrient interactions.