Transcriptome analysis in roots and leaves of wheat seedlings in response to low-phosphorus stress

Transcriptome analysis of Pennisetum americanum × Pennisetum purpureum and Pennisetum americanum leaves in response to high-phosphorus stress

Excessive phosphorus (P) levels can disrupt nutrient balance in plants, adversely affecting growth. The molecular responses of Pennisetum species to high phosphorus stress remain poorly understood. This study examined two Pennisetum species, Pennisetum americanum × Pennisetum purpureum and Pennisetum americanum, under varying P concentrations (200, 600 and 1000 µmol·L- 1 KH2PO4) to elucidate transcriptomic alterations under high-P conditions. Our findings revealed that P. americanum exhibited stronger adaption to high-P stress compared to P. americanum× P. purpureum. Both species showed an increase in plant height and leaf P content under elevated P levels, with P. americanum demonstrating greater height and higher P content than P. americanum× P. purpureum. Transcriptomic analysis identified significant up- and down-regulation of key genes (e.g. SAUR, GH3, AHP, PIF4, PYL, GST, GPX, GSR, CAT, SOD1, CHS, ANR, P5CS and PsbO) involved in plant hormone signal transduction, glutathione metabolism, peroxisomes, flavonoid biosynthesis, amino acid biosynthesis and photosynthesis pathways. Compared with P. americanum× P. purpureum, P. americanum has more key genes in the KEGG pathway, and some genes have higher expression levels. These results contribute valuable insights into the molecular mechanisms governing high-P stress in Pennisetum species and offer implications for broader plant stress research.

RNAseq of diverse spring wheat cultivars released during last 110 years

Here, we performed RNA-seq based expression analysis of root and leaf tissues of a set of 24 historical spring wheat cultivars representing 110 years of temporal genetic variations. This huge 130 tissues RNAseq dataset was initially used to study expression pattern of 97 genes regulating root growth and development in wheat. Root system architecture (RSA) is an important target for breeding stress-resilient and high-yielding wheat cultivars under climatic fluctuations. However, root transcriptome analysis is usually obscured due to challenges in root research due to their below ground presence. We also validated the dataset by performing correlation analysis between expression of RSA related genes in roots and leaves with 25 root traits analyzed under varying moisture conditions and 10 yield-related traits. The Pearson's correlation coefficients between root phenotypes and expression of root-specific genes varied from -0.72 to 0.78, and strong correlations with genes such as DRO1, TaMOR, ARF4, PIN1 was observed. The presented datasets have multiple uses such as a) studying the change in expression pattern of genes during time, b) differential expression of genes in two very important tissues of wheat i.e., leaf and roots, and c) studying customized expression of genes associated with important phenotypes in diverse wheat cultivars. The initial findings presented here provided key insights into understanding the transcriptomic basis of phenotypic variability of RSA in wheat cultivars.

Metabolomics reveals the impact of nitrogen combined with the zinc supply on zinc availability in calcareous soil via root exudates of winter wheat (Triticum aestivum)

A possible mechanism for the improved availability of zinc (Zn) in soil by combining nitrogen (N) with Zn supply was investigated based on the root exudates of winter wheat. N, Zn supply as well as their combination significantly regulated nine root exudates in winter wheat; in which, the secretion of cis-aconitic acid involving in the TCA cycle, C5-branched dibasic acid metabolism, glyoxylate and dicarboxylate metabolism and 2-oxocarboxylic acid metabolism was upregulated by N, Zn supply as well as their combination. N-Zn combination induced the activities of citrate synthase and cis-aconitase in roots and shoots of winter wheat thus to increase the concentrations of citric and aconitic acid; the decrease of isocitric acid concentrations in shoots indicated the inhibited conversion of aconitic acid to isocitric acid by N-Zn combination. It revealed a possible reason for the enhanced secretion of cis-aconitic acid by N-Zn combination. Exogenous addition of 10 μ plant-1 cis-aconitate significantly increased available Zn concentrations in soil and Zn concentrations in winter wheat under N-Zn combination. Thus, the N-Zn combination regulated the metabolism of cis-aconitic acid in winter wheat, thus enhancing the secretion of cis-aconitic acid to increase the bioavailability of Zn in soil.

Rice straw-derived smoke water promotes rice root growth under phosphorus deficiency by modulating oxidative stress and photosynthetic gene expression

Plant-derived smoke has been shown to promote plant growth and seed germination, but its roles and mechanisms in response to nutrient deficiency stress remain unclear. Plants respond to phosphorus (P) deficiency by undergoing morphological, physiological, and transcriptional changes in order to improve nutrient uptake efficiency. Here, we showed that rice straw-derived smoke water could promote root growth in rice (Oryza sativa cv. Nipponbare) grown under P-sufficient and P-deficient conditions. Transcriptome analysis of the root tissues identified 1309 genes up-regulated and 1311 genes down-regulated by smoke water under P-deficient conditions. The GO terms 'glutathione transferase activity' and 'photosynthesis-light reaction' were found to be significantly enriched among the genes that were up- and down-regulated by smoke water, respectively. Biochemical analysis showed that smoke water reduced P-deficient-induced accumulation of H2O2 and malondialdehyde (MDA), a lipid peroxidation marker, reduced sucrose contents, but increased Fe accumulation. Furthermore, smoke water suppressed the expression of strigolactone biosynthesis genes, which were strongly induced by P deficiency as an adaptive strategy to improve root P uptake. These results revealed a potential mechanism by which smoke water promotes root growth and interacts with P deficiency-induced transcriptional regulation to mitigate P deficiency stress in rice.

Comparative proteome analysis of phosphorus-responsive genotypes reveals the proteins differentially expressed under phosphorous starvation stress in rice

Phosphorus (P)-deficiency is one of the major nutrient constraints for global rice production. P-deficiency tolerance in rice involves complex regulatory mechanisms. To gain insights into the proteins involved in phosphorus acquisition and use efficiency in rice, proteome analysis of a high-yielding rice cultivar Pusa-44 and its near-isogenic line (NIL)-23 harboring a major phosphorous uptake (Pup1) QTL, grown under control and P-starvation stress, was performed. Comparative proteome profiling of shoot and root tissues from the plants grown hydroponically with P (16 ppm, +P) or without P (0 ppm, -P) resulted in the identification of 681 and 567 differentially expressed proteins (DEPs) in shoot of Pusa-44 and NIL-23, respectively. Similarly, 66 and 93 DEPs were identified in root of Pusa-44 and NIL-23, respectively. These P-starvation responsive DEPs were annotated to be involved in metabolic processes like photosynthesis, starch-, sucrose-, energy-metabolism, transcription factors (mainly ARF, ZFP, HD-ZIP, MYB), and phytohormone signaling. Comparative analysis of the expression patterns observed by proteome analysis with that reported at the transcriptome level indicated the Pup1 QTL-mediated post-transcriptional regulation plays an important role under -P stress. Thus, the present study describes the molecular aspect of the regulatory functions of Pup1 QTL under P-starvation stress in rice, which might help develop an efficient rice cultivar with enhanced P acquisition and assimilation for better performance in P-deficient soil.

Transcriptomics Insights into Phosphorus Stress Response of Myriophyllum aquaticum

Through excellent absorption and transformation, the macrophyte Myriophyllum (M.) aquaticum can considerably remove phosphorus from wastewater. The results of changes in growth rate, chlorophyll content, and roots number and length showed that M. aquaticum could cope better with high phosphorus stress compared with low phosphorus stress. Transcriptome and differentially expressed genes (DEGs) analyses revealed that, when exposed to phosphorus stresses at various concentrations, the roots were more active than the leaves, with more DEGs regulated. M. aquaticum also showed different gene expression and pathway regulatory patterns when exposed to low phosphorus and high phosphorus stresses. M. aquaticum's capacity to cope with phosphorus stress was maybe due to its improved ability to regulate metabolic pathways such as photosynthesis, oxidative stress reduction, phosphorus metabolism, signal transduction, secondary metabolites biosynthesis, and energy metabolism. In general, M. aquaticum has a complex and interconnected regulatory network that deals efficiently with phosphorus stress to varying degrees. This is the first time that the mechanisms of M. aquaticum in sustaining phosphorus stress have been fully examined at the transcriptome level using high-throughput sequencing analysis, which may indicate the direction of follow-up research and have some guiding value for its future applications.

Transcriptome analysis provides insights into the response of Lotus corniculatus roots to low-phosphorus stress

INTRODUCTION

A lack of soil phosphorus (P) is a principal factor restricting the normal growth of Lotus corniculatus in the karst area of Guizhou Province, China, but the response mechanism of L. corniculatus under low-phosphorus stress remains unclear.

METHODS

Therefore, we treated two selected L. corniculatus lines (low-P-intolerant line 08518 and low-P-tolerant line 01549) from 13 L. corniculatus lines with normal phosphorus (0.5 mmol/L KH₂PO₄, NP) and low phosphorus (0.005 mmol/L KH₂PO₄, LP) concentrations to study changes in morphological, physiological and transcriptome data under low-phosphorus stress.

RESULTS

The low-P-tolerant line 01549 exhibited better performance under low-phosphorus stress. Compared with the NP treatment, all root morphological indicators of the low-P-tolerant line 01549 increased, and those of the low-P-intolerant line 08518 decreased under low-P stress. Compared with the NP treatment, acid phosphatase (ACP), catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activities, and the malondialdehyde (MDA), soluble sugar (SS), soluble protein (SP) and proline (Pro) contents of the two L. corniculatus lines increased under low-P stress. A transcriptome analysis of L. corniculatus showed that a total of 656 and 2243 differentially expressed genes (DEGs) were identified in line 01549 and line 08518, respectively. Meanwhile, the main pathways, such as carbohydrate metabolism, acid phosphatases, phosphate transporters and biosynthesis of secondary metabolites, as well as related genes were also screened by performing a KEGG enrichment analysis.

DISCUSSION

The findings provide an essential point of reference for studying the physiological and molecular mechanism of the response to low-P stress in L. corniculatus.

Integrative physiological, transcriptome and metabolome analysis reveals the involvement of carbon and flavonoid biosynthesis in low phosphorus tolerance in cotton

Phosphorus (P) is an essential nutrient controlling plant growth and development through the regulation of basic metabolic processes; however, the molecular details of these pathways remain largely unknown. In this study, physiological, transcriptome, and metabolome analysis were compared for two cotton genotypes with different low P tolerance under P starvation and resupply. The results showed that the glucose, fructose, sucrose, and starch contents increased by 18.2%, 20.4%, 20.2%, and 14.3% in the roots and 18.3%, 23.3%, 11.0%, and 13.6% in the shoot of Jimian169 than DES926, respectively. Moreover, the activities of enzymes related to carbon and phosphorus metabolism were higher in the roots and shoots of Jimian169 than DES926. In addition, transcriptome analysis revealed that the number of differentially expressed genes (DEGs) was higher in both roots (830) and shoots (730) under P starvation and the DEGs drastically reduced upon P resupply. The KEGG analysis indicated that DEGs were mainly enriched in phenylpropanoid biosynthesis, carbon metabolism, and photosynthesis. The metabolome analysis showed the enrichment of phenylpropanoid, organic acids and derivatives, and lipids in all the pairs at a given time point. The combined transcriptome and metabolome analysis revealed that carbon metabolism and flavonoid biosynthesis are involved in the P starvation response in cotton. Moreover, co-expression network analysis identified 3 hub genes in the roots and shoots that regulate the pathways involved in the P starvation response. This study provides the foundation for understanding the mechanisms of low P tolerance and the hub genes as a potential target for the development of low P tolerant genotypes.

Integrated transcriptomic and metabolomic analyses reveal key metabolic pathways in response to potassium deficiency in coconut (Cocos nucifera L.) seedlings

Potassium ions (K⁺) are important for plant growth and crop yield. However, the effects of K⁺ deficiency on the biomass of coconut seedlings and the mechanism by which K⁺ deficiency regulates plant growth remain largely unknown. Therefore, in this study, we compared the physiological, transcriptome, and metabolite profiles of coconut seedling leaves under K⁺-deficient and K⁺-sufficient conditions using pot hydroponic experiments, RNA-sequencing, and metabolomics technologies. K⁺ deficiency stress significantly reduced the plant height, biomass, and soil and plant analyzer development value, as well as K content, soluble protein, crude fat, and soluble sugar contents of coconut seedlings. Under K⁺ deficiency, the leaf malondialdehyde content of coconut seedlings were significantly increased, whereas the proline (Pro) content was significantly reduced. Superoxide dismutase, peroxidase, and catalase activities were significantly reduced. The contents of endogenous hormones such as auxin, gibberellin, and zeatin were significantly decreased, whereas abscisic acid content was significantly increased. RNA-sequencing revealed that compared to the control, there were 1003 differentially expressed genes (DEGs) in the leaves of coconut seedlings under K⁺ deficiency. Gene Ontology analysis revealed that these DEGs were mainly related to "integral component of membrane," "plasma membrane," "nucleus", "transcription factor activity," "sequence-specific DNA binding," and "protein kinase activity." Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that the DEGs were mainly involved in "MAPK signaling pathway-plant," "plant hormone signal transduction," "starch and sucrose metabolism," "plant-pathogen interaction," "ABC transporters," and "glycerophospholipid metabolism." Metabolomic analysis showed that metabolites related to fatty acids, lipidol, amines, organic acids, amino acids, and flavonoids were generally down-regulated in coconut seedlings under K⁺ deficiency, whereas metabolites related to phenolic acids, nucleic acids, sugars, and alkaloids were mostly up-regulated. Therefore, coconut seedlings respond to K⁺ deficiency stress by regulating signal transduction pathways, primary and secondary metabolism, and plant-pathogen interaction. These results confirm the importance of K⁺ for coconut production, and provide a more in-depth understanding of the response of coconut seedlings to K⁺ deficiency and a basis for improving K⁺ utilization efficiency in coconut trees.

Comprehensive analysis of aberrantly expressed circRNAs, mRNAs and lncRNAs in patients with nasopharyngeal carcinoma

BACKGROUND

The location of nasopharyngeal cancer is hidden, so it is difficult to diagnose at an early stage. In this study, we aimed to investigate the expression profiles of circRNAs, mRNAs and IncRNAs and to provide some basis for further studies.

METHODS

Expression profiles of circRNAs, mRNAs, and lncRNAs were analyzed using microarray techniques. The differentially expressed ncRNA was calculated by bioinformatics.

RESULTS

A total of 3048 circRNAs, 2179 lncRNAs, and 2015 mRNAs were detected to be significantly differentially expressed in NPC. The most upregulated circRNAs, lncRNAs, and mRNAs were hsa-circ-0067562, NONHSAT232922.1, and HOXB13, respectively. And, the most downregulated circRNAs, lncRNAs, and mRNAs were hsa_circ_0078837, lnc-TTC8-4:3, and LTF, respectively. The number of upregulated DE lncRNAs was more than twice than those downregulated. Our data showed that 80.44% of pairs of lncRNAs and cis-mRNAs demonstrated positive correlations. For lncRNAs and trans-mRNAs pairs, 53.7% of pairs showed positive correlation. LncRNA-mediated cis regulation is a prevalent regulatory mode in the development of nasopharyngeal carcinoma. CR1, LRMP and SORBS2 are predicted to be mediated not only by cis-acting lncRNA modes of action, but also by trans-acting lncRNA mechanisms. Additionally, we constructed a diagnostic prediction model with a high sensitivity and specificity.

CONCLUSION

Our study characterized the landscape of circRNAs, mRNAs and lncRNAs in NPC tissue and provided novel insights into the molecular mechanisms of NPC.

Protoplast Isolation for Plant Single-Cell RNA-seq

The growth and development of plants depends on diversified gene expression in different cell types. Compared to traditional bulk RNA sequencing, droplet-based single-cell RNA sequencing (scRNA-seq) allows for transcriptome profiling of individual cells within heterogeneous tissues. scRNA-seq provides a high-resolution atlas of cellular characterization and vastly improves our understandings of the interactions between individual cells and the microenvironment. However, the difficulty in protoplast isolation has limited the application of single-cell sequencing technology in plant research. Here we describe a high-efficiency protoplast isolation protocol for scRNA-seq.

CRISPR-Cas9-mediated editing of BADH2 gene triggered fragrance revolution in rice

Fragrance is one of the most important quality traits for breeding in rice. The natural aroma substance 2-acetyl-1-pyrroline (2-AP) is a key fragrance compound among over 200 volatiles identified in fragrant rice. In addition to rice, there are other plant species that contain a germplasm that naturally produces a fragrant aroma. These other plant species all have lower activity levels of the enzyme BETAINE ALDEHYDE DEHYDROGENASE 2 (BADH2). Therefore, improving fragrance efficiency has been a focus of intensive research. Recent studies have engineered BADH2 gene, which is responsible for fragrance trait in non-fragrant cultivars of rice, using CRISPR-Cas9. Although engineering rice BADH2 can be useful for upregulating 2-AP, there are still a lot of restrictions on how it can be applied in practice. In this review article, we discuss the recent developments in BADH2 editing and propose potential future strategies to effectively target BADH2 for transcriptional regulation, with the goal of producing a better fragrance.

Physiological and Transcriptomic Analyses Uncover the Reason for the Inhibition of Photosynthesis by Phosphate Deficiency in Cucumis melo L

Phosphate (Pi) deficiency is a common phenomenon in agricultural production and limits plant growth. Recent work showed that long-term Pi deficiency caused the inhibition of photosynthesis and inefficient electron transport. However, the underlying mechanisms are still unknown. In this study, we used the physiological, histochemical, and transcriptomic methods to investigate the effect of low-Pi stress on photosynthetic gas exchange parameters, cell membrane lipid, chloroplast ultrastructure, and transcriptional regulation of key genes in melon seedlings. The results showed that Pi deficiency significantly downregulated the expression of aquaporin genes, induced an increase in ABA levels, and reduced the water content and free water content of melon leaves, which caused physiological drought in melon leaves. Therefore, gas exchange was disturbed. Pi deficiency also reduced the phospholipid contents in leaf cell membranes, caused the peroxidation of membrane lipids, and destroyed the ultrastructure of chloroplasts. The transcriptomic analysis showed that 822 differentially expressed genes (DEGs) were upregulated and 1254 downregulated by Pi deficiency in leaves. GO and KEGG enrichment analysis showed that DEGs significantly enriched in chloroplast thylakoid membrane composition (GO:0009535), photosynthesis-antenna proteins (map00196), and photosynthesis pathways (map00195) were downregulated by Pi deficiency. It indicated that Pi deficiency regulated photosynthesis-related genes at the transcriptional level, thereby affecting the histochemical properties and physiological functions, and consequently causing the reduced light assimilation ability and photosynthesis efficiency. It enriches the mechanism of photosynthesis inhibition by Pi deficiency.

Differentially reset transcriptomes and genome bias response orchestrate wheat response to phosphate deficiency

Phosphorus (P) is an essential macronutrient for all organisms. Phosphate (Pi) deficiency reduces grain yield and quality in wheat. Understanding how wheat responds to Pi deficiency at the global transcriptional level remains limited. We revisited the available RNA-seq transcriptome from Pi-starved wheat roots and shoots subjected to Pi starvation. Genome-wide transcriptome resetting was observed under Pi starvation, with a total of 917 and 2338 genes being differentially expressed in roots and shoots, respectively. Chromosomal distribution analysis of the gene triplets and differentially expressed genes (DEGs) revealed that the D genome displayed genome induction bias and, specifically, the chromosome 2D might be a key contributor to Pi-limiting triggered gene expression response. Alterations in multiple metabolic pathways pertaining to secondary metabolites, transcription factors and Pi uptake-related genes were evidenced. This study provides genomic insight and the dynamic landscape of the transcriptional changes contributing to the hexaploid wheat during Pi starvation. The outcomes of this study and the follow-up experiments have the potential to assist the development of Pi-efficient wheat cultivars.

Local and systemic responses conferring acclimation of Brassica napus roots to low phosphorus conditions

Due to the non-uniform distribution of inorganic phosphate (Pi) in the soil, plants modify their root architecture to improve acquisition of this nutrient. In this study, a split-root system was employed to assess the nature of local and systemic signals that modulate root architecture of Brassica napus grown with non-uniform Pi availability. Lateral root (LR) growth was regulated systemically by non-uniform Pi distribution, by increasing the second-order LR (2°LR) density in compartments with high Pi supply but decreasing it in compartments with low Pi availability. Transcriptomic profiling identified groups of genes regulated, both locally and systemically, by Pi starvation. The number of systemically induced genes was greater than the number of genes locally induced, and included genes related to abscisic acid (ABA) and jasmonic acid (JA) signalling pathways, reactive oxygen species (ROS) metabolism, sucrose, and starch metabolism. Physiological studies confirmed the involvement of ABA, JA, sugars, and ROS in the systemic Pi starvation response. Our results reveal the mechanistic basis of local and systemic responses of B. napus to Pi starvation and provide new insights into the molecular and physiological basis of root plasticity.

Transcriptome unveiled the gene expression patterns of root architecture in drought-tolerant and sensitive wheat genotypes

Drought is a big challenge for agricultural production. Root attributes are the important target traits for breeding high-yielding sustainable wheat varieties against ever changing climatic conditions. However, the transcriptomic of wheat concerning root architecture remained obscure. Here, we explored RNA-Seq based transcriptome to dissect putative genes involved in root system variations in naturally occurring six genotypes (drought-tolerant and sensitive) of wheat. Global RNA-Seq based root transcriptome analysis revealed single nucleotide polymorphisms (SNPs) variations and differentially expressed genes. Putative 56 SNPs were identified related to 15 genes involved in root architecture. Enrichment of these genes using GO terms demonstrated that differentially expressed genes (DEGs) are divided into sub-categories implicated in molecular functions, cellular components and biological processes. The KEGG analysis of DEGs in each comparison of genotype include metabolic, biosynthesis of secondary metabolites, microbial metabolism in diverse environments and biosynthesis of antibiotics. A deeper insight into DEGs unveiled various pathways involved in drought response and positive gravitropism. These genes belong to various transcription factor families such as DOF, C3H, MYB, and NAC involved in root developmental and stress-related pathways. Local White and UZ-11-CWA-8, which are drought-tolerant genotypes, harbor over-representation of most of DEGs or transcription factors. Notably, a microtubule-associated protein MAPRE1 belonging to RP/EB family recruited in positive gravitropism was enriched. Real-time PCR analysis revealed expression of MAPRE1 and PAL genes is consistent with RNA-seq data. The presented data and genetic resources seem valuable for providing genes involved in the root system architecture of drought-tolerant and susceptible genotypes.

Transcriptome analysis of a near-isogenic line and its recurrent parent reveals the role of Pup1 QTL in phosphorus deficiency tolerance of rice at tillering stage

Phosphorus (P) is essential for cellular processes like respiration, photosynthesis, biosynthesis of membrane phospholipids, etc. To cope with P deficiency stress, plants adopt reprograming of the expression of genes involved in different metabolic/signaling pathways for survival, growth, and development. Plants use transcriptional, post-transcriptional, and/or post-translational machinery to achieve P homeostasis. Several transcription factors (TFs), miRNAs, and P transporters play important roles in P deficiency tolerance; however, the underlying mechanisms responsible for P deficiency tolerance remain poorly understood. Studies on P starvation/deficiency responses in plants at early (seedling) stage of growth have been reported but only a few of them focused on molecular responses of the plant at advanced (tillering or reproductive) stage of growth. To decipher the strategies adopted by rice at tillering stage under P deficiency stress, a pair of contrasting genotypes [Pusa-44 (a high-yielding, P deficiency sensitive cultivar) and its near-isogenic line (NIL-23, P deficiency tolerant) for Pup1 QTL] was used for morphophysiological, biochemical, and molecular analyses. Comparative analyses of shoot and root tissues from 45-day-old plants grown hydroponically under P sufficient (16 ppm) or P deficient (4 ppm) medium confirmed some of the known morphophysiological responses. Moreover, RNA-seq analysis revealed the important roles of phosphate transporters, TFs, auxin-responsive proteins, modulation in the cell wall, fatty acid metabolism, and chromatin architecture/epigenetic modifications in providing P deficiency tolerance to NIL-23, which were brought in due to the introgression of the Pup1 QTL in Pusa-44. This study provides insights into the molecular functions of Pup1 for P deficiency tolerance, which might be utilized to improve P-use efficiency of rice for better productivity in P deficient soils. KEY MESSAGE: Introgression of Pup1 QTL in high-yielding rice cultivar modulates mainly phosphate transporters, TFs, auxin-responsive proteins, cell wall structure, fatty acid metabolism, and chromatin architecture/epigenetic modifications at tillering stage of growth under phosphorus deficiency stress.

Role of Wheat Phosphorus Starvation Tolerance 1 Genes in Phosphorus Acquisition and Root Architecture

The wheat plant requires elevated phosphorus levels for its normal growth and yield, but continuously depleting non-renewable phosphorus reserves in the soil is one of the biggest challenges in agricultural production worldwide. The Phosphorus Starvation Tolerance 1 (PSTOL1) gene has been reported to play a key role in efficient P uptake, deeper rooting, and high yield in rice. However, the function of the PSTOL1 gene in wheat is still unclear. In this study, a total of 22 PSTOL1 orthologs were identified in the wheat genome, and found that wheat PSTOL1 orthologs are unevenly distributed on chromosomes, and these genes were under strong purifying selection. Under different phosphorus regimes, wheat PSTOL1 genes showed differential expression patterns in different tissues. These results strengthen the classification of Pakistan-13 as a P-efficient cultivar and Shafaq-06 as a P-inefficient cultivar. Phenotypic characterization demonstrated that Pakistan-13 wheat cultivar has significantly increased P uptake, root length, root volume, and root surface area compared to Shafaq-06. Some wheat PSTOL1 orthologs are co-localized with phosphorus starvation’s related quantitative trait loci (QTLs), suggesting their potential role in phosphorus use efficiency. Altogether, these results highlight the role of the wheat PSTOL1 genes in wheat P uptake, root architecture, and efficient plant growth. This comprehensive study will be helpful for devising sustainable strategies for wheat crop production and adaptation to phosphorus insufficiency.

Impacts, Tolerance, Adaptation, and Mitigation of Heat Stress on Wheat under Changing Climates

Heat stress (HS) is one of the major abiotic stresses affecting the production and quality of wheat. Rising temperatures are particularly threatening to wheat production. A detailed overview of morpho-physio-biochemical responses of wheat to HS is critical to identify various tolerance mechanisms and their use in identifying strategies to safeguard wheat production under changing climates. The development of thermotolerant wheat cultivars using conventional or molecular breeding and transgenic approaches is promising. Over the last decade, different omics approaches have revolutionized the way plant breeders and biotechnologists investigate underlying stress tolerance mechanisms and cellular homeostasis. Therefore, developing genomics, transcriptomics, proteomics, and metabolomics data sets and a deeper understanding of HS tolerance mechanisms of different wheat cultivars are needed. The most reliable method to improve plant resilience to HS must include agronomic management strategies, such as the adoption of climate-smart cultivation practices and use of osmoprotectants and cultured soil microbes. However, looking at the complex nature of HS, the adoption of a holistic approach integrating outcomes of breeding, physiological, agronomical, and biotechnological options is required. Our review aims to provide insights concerning morpho-physiological and molecular impacts, tolerance mechanisms, and adaptation strategies of HS in wheat. This review will help scientific communities in the identification, development, and promotion of thermotolerant wheat cultivars and management strategies to minimize negative impacts of HS.

Improving phosphate use efficiency in the aquatic crop watercress (Nasturtium officinale)

Watercress is a nutrient-dense leafy green crop, traditionally grown in aquatic outdoor systems and increasingly seen as well-suited for indoor hydroponic systems. However, there is concern that this crop has a detrimental impact on the environment through direct phosphate additions causing environmental pollution. Phosphate-based fertilisers are supplied to enhanced crop yield, but their use may contribute to eutrophication of waterways downstream of traditional watercress farms. One option is to develop a more phosphate use efficient (PUE) crop. This review identifies the key traits for this aquatic crop (the ideotype), for future selection, marker development and breeding. Traits identified as important for PUE are (i) increased root surface area through prolific root branching and adventitious root formation, (ii) aerenchyma formation and root hair growth. Functional genomic traits for improved PUE are (iii) efficacious phosphate remobilisation and scavenging strategies and (iv) the use of alternative metabolic pathways. Key genomic targets for this aquatic crop are identified as: PHT phosphate transporter genes, global transcriptional regulators such as those of the SPX family and genes involved in galactolipid and sulfolipid biosynthesis such as MGD2/3, PECP1, PSR2, PLDζ1/2 and SQD2. Breeding for enhanced PUE in watercress will be accelerated by improved molecular genetic resources such as a full reference genome sequence that is currently in development.

Transfer Learning for High-Dimensional Linear Regression: Prediction, Estimation and Minimax Optimality

This paper considers estimation and prediction of a high-dimensional linear regression in the setting of transfer learning where, in addition to observations from the target model, auxiliary samples from different but possibly related regression models are available. When the set of informative auxiliary studies is known, an estimator and a predictor are proposed and their optimality is established. The optimal rates of convergence for prediction and estimation are faster than the corresponding rates without using the auxiliary samples. This implies that knowledge from the informative auxiliary samples can be transferred to improve the learning performance of the target problem. When the set of informative auxiliary samples is unknown, we propose a data-driven procedure for transfer learning, called Trans-Lasso, and show its robustness to non-informative auxiliary samples and its efficiency in knowledge transfer. The proposed procedures are demonstrated in numerical studies and are applied to a dataset concerning the associations among gene expressions. It is shown that Trans-Lasso leads to improved performance in gene expression prediction in a target tissue by incorporating data from multiple different tissues as auxiliary samples.

Integrative analysis of the metabolome and transcriptome reveal the phosphate deficiency response pathways of alfalfa

Understanding the mechanisms underlying the responses to inorganic phosphate (Pi) deficiency in alfalfa will help enhance Pi acquisition efficiency and the sustainable use of phosphorous resources. Integrated global metabolomic and transcriptomic analyses of mid-vegetative alfalfa seedlings under 12-day Pi deficiency were conducted. Limited seedling growth were found, including 13.24%, 16.85% and 33.36% decreases in height, root length and photosynthesis, and a 24.10% increase in root-to-shoot ratio on day 12. A total of 322 and 448 differentially abundant metabolites and 1199 and 1061 differentially expressed genes were identified in roots and shoots. Increased (>3.68-fold) inorganic phosphate transporter 1;4 and SPX proteins levels in the roots (>2.15-fold) and shoots (>2.50-fold) were related to Pi absorption and translocation. The levels of phospholipids and Pi-binding carbohydrates and nucleosides were decreased, while those of phosphatases and pyrophosphatases in whole seedlings were induced under reduced Pi. In addition, nitrogen assimilation was affected by inhibiting high-affinity nitrate transporters (NRT2.1 and NRT3.1), and nitrate reductase. Increased delphinidin-3-glucoside might contribute to the gray-green leaves induced by Pi limitation. Stress-induced MYB, WRKY and ERF transcription factors were identified. The responses of alfalfa to Pi deficiency were summarized as local systemic signaling pathways, including root growth, stress-related responses consisting of enzymatic and nonenzymatic systems, and hormone signaling and systemic signaling pathways including Pi recycling and Pi sensing in the whole plant, as well as Pi recovery, and nitrate and metal absorption in the roots. This study provides important information on the molecular mechanism of the response to Pi deficiency in alfalfa.

Physiological and molecular response of colored wheat seedlings against phosphate deficiency is linked to accumulation of distinct anthocyanins

Anthocyanin rich colored wheat with additional health benefits has created interest among breeders, consumers and policy makers to address the prevailing malnutrition in the vulnerable population. Researchers are exploring how colored wheat could perform under different nutrient conditions for the maintenance of growth and development. The present study was aimed to investigate the differential response of phosphorous (Pi) deficiency at the seedling stage using hydroponics. Our results showed that Pi-deficiency triggered typical response in the wheat along with the changes in the plant root morphology, total biomass, micronutrient concentration and distinct anthocyanin accumulation. Our physiological and biochemical data revealed that these parameters were positively altered under stress in the colored wheat and the adaptation followed the trend of white < blue <purple < black. Our results also confirmed that stress induced accumulation of distinct anthocyanins including derivatives of cyanidin, delphinidin and peonidin in a genotype dependent manner. Differential expression pattern visualized for the transcripts encoding phosphate transporters, anthocyanin biosynthesis, putative transporters and regulators may be one of the underlying factors. Altogether, our data showed that the black wheat genotype with highest anthocyanin content could able to adapt better with the P stress. This study will help in identifying suitable colored wheat adapting the stress condition and have potential for influence on the future agricultural cultivation practices.

Allyl Isothiocyanate in the Volatiles of Brassica juncea Inhibits the Growth of Root Rot Pathogens of Panax notoginseng by Inducing the Accumulation of ROS

The cultivation of Panax notoginseng is often seriously hindered by root rot disease caused by the accumulation of soil-borne pathogens. Here, the inhibitory activity of Brassica juncea volatiles on P. notoginseng root rot pathogens was assessed and compounds in volatiles were identified. Furthermore, the antimicrobial activity and mechanism of allyl isothiocyanate (AITC) were deciphered by integrated transcriptome and metabolome analyses. The volatiles of B. juncea showed dose-dependent antimicrobial activity against root rot pathogens. AITC, identified as the main volatile compound, not only significantly inhibited pathogen growth in vitro but also suppressed root rot disease in the field. Integrated transcriptomic and metabolomics analysis revealed that AITC inhibited Fusarium solani by interfering with energy production and induced the accumulation of ROS by decreasing the content of glutathione (GSH). In summary, B. juncea releases AITC to inhibit soil-borne pathogens and could be used as a rotation crop or soil fumigant to alleviate root rot disease.

A new biomimetic nanozyme of hemin/graphdiyne oxide with superior peroxidase-like activity for colorimetric bioassays

Graphdiyne oxide (GDYO) is a novel type of two-dimensional carbon allotrope nanomaterial consisting of a large conjugated system and excellent chemical stability. To date, application of GDYO as a nanozyme in biosensing has been rarely reported. In this study, a novel ultrasensitive colorimetric bioassay was constructed using a hemin/GDYO nanocomposite (H/GDYO) as a new nanozyme with superior peroxidase-like activity for the detection of H₂O₂ and glucose. It was discovered that H/GDYO exhibited 6-fold higher peroxidase-like activity than pure hemin. Catalytic kinetic analysis showed that H/GDYO had a much higher affinity for H₂O₂ and glucose than that of hemin. The designed colorimetric bioassay displayed excellent sensitivity for H₂O₂ and glucose detection with a wide linear range of 0.015-0.5 mM and 0.1-10 mM, respectively, while the limit of detection (LOD) was as low as 4.39 μM and 38 μM, respectively. Moreover, it was successfully applied for the analysis of H₂O₂ in milk and glucose in real human serum samples with acceptable recoveries. Importantly, the developed colorimetric bioassay shows good agreement with the results obtained from a commercial blood glucose meter. We believe that the proposed method could provide a promising prospect for medical diagnosis and biotechnology.

Temperature Dynamics in Rat Brains Exposed to Near-Field Waveguide Outputs at 2.8 GHz

Biological effects in the microwave band of the radiofrequency (RF) spectrum are thermally mediated. For acute high-power microwave exposures, these effects will depend on transient time-temperature histories within the tissue. In this article, we summarize the transient temperature response of rats exposed to RF energy emanating from an open-ended rectangular waveguide. These exposures produced specific absorption rates of approximately 36 and 203 W/kg in the whole body and brain, respectively. We then use the experimentally measured thermal data to infer the baseline perfusion rate in the brain and modify a custom thermal modeling tool based upon these findings. Finally, we compare multi-physics simulations of rat brain temperature against empirical measurements in both live and euthanized subjects and find close agreement between model and experimentation. This research revealed that baseline brain perfusion rates in rat subjects could be larger than previously assumed in the RF thermal modeling literature, and plays a significant role in the transient thermal response to high-power microwave exposures. © 2021 Bioelectromagnetics Society.

Convergent evolution of gene regulatory networks underlying plant adaptations to dry environments

Plants transitioned from an aquatic to a terrestrial lifestyle during their evolution. On land, fluctuations on water availability in the environment became one of the major problems they encountered. The appearance of morpho-physiological adaptations to cope with and tolerate water loss from the cells was undeniably useful to survive on dry land. Some of these adaptations, such as carbon concentrating mechanisms (CCMs), desiccation tolerance (DT) and root impermeabilization, appeared in multiple plant lineages. Despite being crucial for evolution on land, it has been unclear how these adaptations convergently evolved in the various plant lineages. Recent advances on whole genome and transcriptome sequencing are revealing that co-option of genes and gene regulatory networks (GRNs) is a common feature underlying the convergent evolution of these adaptations. In this review, we address how the study of CCMs and DT has provided insight into convergent evolution of GRNs underlying plant adaptation to dry environments, and how these insights could be applied to currently emerging understanding of evolution of root impermeabilization through different barrier cell types. We discuss examples of co-option, conservation and innovation of genes and GRNs at the cell, tissue and organ levels revealed by recent phylogenomic (comparative genomic) and comparative transcriptomic studies.

Transcriptome and metabolome analyses revealed the response mechanism of apple to different phosphorus stresses

Phosphorus (P) is an important element in numerous metabolic reactions and signalling pathways, but the molecular details of these pathways remain largely unknown. In this study, physiological, transcriptome and metabolite analyses of apple leaves and roots were compared under different P conditions. The results showed that different P stresses influenced phenotypic characteristics, soil plant analytical development (SPAD) values and the contents of flavonoids and anthocyanins in apple seedlings. The contents of hydrogen peroxide (H₂O₂) and malondialdehyde (MDA) and the activities of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), acid phosphatase (ACP) and purple acid phosphatase (PAP) were also affected by different P stresses. In addition, RNA sequencing (RNA-seq) was used to characterize the influence of different P stresses on apple seedlings. Compared with control apple plants, there were 1246 and 1183 differentially expressed genes (DEGs) in leaves and roots under the low-P treatment and 60 and 1030 DEGs in leaves and roots under the high-P treatment, respectively. Gene Ontology (GO) analysis indicated that apple trees might change their responses to metabolic processes, cell proliferation, regulation of biological processes, reactive oxygen species metabolic processes and flavonoid metabolic processes under P stress. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis further indicated that DEGs act on the mitogen-activated protein kinase (MAPK) signalling pathway, flavonoid biosynthesis, phenylpropanoid biosynthesis, and ATP-binding cassette (ABC) transporters. The metabolome analysis revealed that the levels of most amino acids and their derivatives, organic acids and flavonoids in roots treated with low-P stress were higher than those in roots of apple seedlings under control growth conditions. Apple seedlings regulate the flavonoid pathway to respond to different phosphorus environments. The results provide a framework for understanding the metabolic processes underlying different P responses and provide a foundation for improving the utilization efficiency of P in apple trees.

ANOVA-particle swarm optimization-based feature selection and gradient boosting machine classifier for improved protein-protein interaction prediction

Feature fusion and selection strategies have been applied to improve accuracy in the prediction of protein-protein interaction (PPI). In this paper, an embedded feature selection framework is developed by integrating a cost function based on analysis of variance (ANOVA) with the particle swarm optimization (PSO), termed AVPSO. Initially, the features of the protein sequences extracted using pseudo-amino acid composition (PseAAC), conjoint triad composition, and local descriptor are fused. Then, AVPSO is employed to select the optimal set of features. The light gradient boosting machine (LGBM) classifier is used to predict the PPIs using the optimal feature subset. On the five-fold cross-validation analysis, the proposed model (AVPSO-LGBM) achieved an average accuracy of 97.12% and 95.09%, respectively, on the intraspecies PPI datasets Saccharomyces cerevisiae and Helicobacter pylori. On the interspecies, PPI datasets of the Human-Bacillus and Human-Yersinia, an average accuracy of 95.20% and 93.44%, are achieved. Results obtained on independent test datasets, and network datasets show that the prediction accuracy of the AVPSO-LGBM is better than the existing methods, demonstrating its generalization ability. The improved prediction performance obtained by the proposed model makes it a reliable and effective PPI prediction model.

Transcriptomic Analysis of Salt-Stress-Responsive Genes in Barley Roots and Leaves

Barley is characterized by a rich genetic diversity, making it an important model for studies of salinity response with great potential for crop improvement. Moreover, salt stress severely affects barley growth and development, leading to substantial yield loss. Leaf and root transcriptomes of a salt-tolerant Tunisian landrace (Boulifa) exposed to 2, 8, and 24 h salt stress were compared with pre-exposure plants to identify candidate genes and pathways underlying barley's response. Expression of 3585 genes was upregulated and 5586 downregulated in leaves, while expression of 13,200 genes was upregulated and 10,575 downregulated in roots. Regulation of gene expression was severely impacted in roots, highlighting the complexity of salt stress response mechanisms in this tissue. Functional analyses in both tissues indicated that response to salt stress is mainly achieved through sensing and signaling pathways, strong transcriptional reprograming, hormone osmolyte and ion homeostasis stabilization, increased reactive oxygen scavenging, and activation of transport and photosynthesis systems. A number of candidate genes involved in hormone and kinase signaling pathways, as well as several transcription factor families and transporters, were identified. This study provides valuable information on early salt-stress-responsive genes in roots and leaves of barley and identifies several important players in salt tolerance.

Endothelial cells from umbilical cord of women affected by gestational diabetes: A suitable in vitro model to study mechanisms of early vascular senescence in diabetes

Human umbilical cord endothelial cells (HUVECs) obtained from women affected by gestational diabetes (GD-HUVECs) display durable pro-atherogenic modifications and might be considered a valid in vitro model for studying chronic hyperglycemia effects on early endothelial senescence. Here, we demonstrated that GD- compared to C-HUVECs (controls) exhibited oxidative stress, altered both mitochondrial membrane potential and antioxidant response, significant increase of senescent cells characterized by a reduced NAD-dependent deacetylase sirtuin-1 (SIRT1) activity together with an increase in cyclin-dependent kinase inhibitor-2A (P16), cyclin-dependent kinase inhibitor-1 (P21), and tumor protein p53 (P53) acetylation. This was associated with the p300 activation, and its silencing significantly reduced the GD-HUVECs increased protein levels of P300 and Ac-P53 thus indicating a persistent endothelial senescence via SIRT1/P300/P53/P21 pathway. Overall, our data suggest that GD-HUVECs can represent an "endothelial hyperglycemic memory" model to investigate in vitro the early endothelium senescence in cells chronically exposed to hyperglycemia in vivo.

Characterization of contrasting rice (Oryza sativa L.) genotypes reveals the Pi-efficient schema for phosphate starvation tolerance

BACKGROUND

Phosphorus (P), being one of the essential components of nucleic acids, cell membranes and enzymes, indispensable for diverse cellular processes like photosynthesis/carbohydrate metabolism, energy production, redox homeostasis and signaling. Crop yield is severely affected due to Phosphate (Pi) deficiency; and to cope with Pi-deficiency, plants have evolved several strategies. Some rice genotypes are compatible with low Pi availability, whereas others are sensitive to Pi deficiency. However, the underlying molecular mechanism for low Pi tolerance remains largely unexplored.

RESULT

Several studies were carried out to understand Pi-deficiency responses in rice at seedling stage, but few of them targeted molecular aspects/responses of Pi-starvation at the advanced stage of growth. To delineate the molecular mechanisms for low Pi tolerance, a pair of contrasting rice (Oryza sativa L.) genotypes [viz. Pusa-44 (Pi-deficiency sensitive) and its near isogenic line (NIL-23, Pi-deficiency tolerant) harboring Phosphorus uptake 1 (Pup1) QTL from an aus landrace Kasalath] were used. Comparative morphological, physiological, and biochemical analyses confirmed some of the well-known findings. Transcriptome analysis of shoot and root tissues from 45-day-old rice plants grown hydroponically under P-sufficient (16 ppm Pi) or P-starved (0 ppm Pi) medium revealed that Pi-starvation stress causes global transcriptional reprogramming affecting several transcription factors, signaling pathways and other regulatory genes. We could identify several significantly up-regulated genes in roots of NIL-23 under Pi-starvation which might be responsible for the Pi starvation tolerance. Pathway enrichment analysis indicated significant role of certain phosphatases, transporters, transcription factors, carbohydrate metabolism, hormone-signaling, and epigenetic processes in improving P-starvation stress tolerance in NIL-23.

CONCLUSION

We report the important candidate mechanisms for Pi acquisition/solubilization, recycling, remobilization/transport, sensing/signalling, genetic/epigenetic regulation, and cell wall structural changes to be responsible for P-starvation tolerance in NIL-23. The study provides some of the novel information useful for improving phosphorus-use efficiency in rice cultivars.

Assisted estimation of gene expression graphical models

In the study of gene expression data, network analysis has played a uniquely important role. To accommodate the high dimensionality and low sample size and generate interpretable results, regularized estimation is usually conducted in the construction of gene expression Gaussian Graphical Models (GGM). Here we use GeO-GGM to represent gene-expression-only GGM. Gene expressions are regulated by regulators. gene-expression-regulator GGMs (GeR-GGMs), which accommodate gene expressions as well as their regulators, have been constructed accordingly. In practical data analysis, with a "lack of information" caused by the large number of model parameters, limited sample size, and weak signals, the construction of both GeO-GGMs and GeR-GGMs is often unsatisfactory. In this article, we recognize that with the regulation between gene expressions and regulators, the sparsity structures of a GeO-GGM and its GeR-GGM counterpart can satisfy a hierarchy. Accordingly, we propose a joint estimation which reinforces the hierarchical structure and use the construction of a GeO-GGM to assist that of its GeR-GGM counterpart and vice versa. Consistency properties are rigorously established, and an effective computational algorithm is developed. In simulation, the assisted construction outperforms the separation construction of GeO-GGM and GeR-GGM. Two The Cancer Genome Atlas data sets are analyzed, leading to findings different from the direct competitors.

A comparison of shared patterns of differential gene expression and gene ontologies in response to water-stress in roots and leaves of four diverse genotypes of Lolium and Festuca spp. temperate pasture grasses

Ryegrasses (Lolium spp.) and fescues (Festuca spp.) are closely related and widely cultivated perennial forage grasses. As such, resilience in the face of abiotic stresses is an important component of their traits. We have compared patterns of differentially expressed genes (DEGs) in roots and leaves of two perennial ryegrass genotypes and a single genotype of each of a festulolium (predominantly Italian ryegrass) and meadow fescue with the onset of water stress, focussing on overall patterns of DEGs and gene ontology terms (GOs) shared by all four genotypes. Plants were established in a growing medium of vermiculite watered with nutrient solution. Leaf and root material were sampled at 35% (saturation) and, as the medium dried, at 15%, 5% and 1% estimated water contents (EWCs) and RNA extracted. Differential gene expression was evaluated comparing the EWC sampling points from RNAseq data using a combination of analysis methods. For all genotypes, the greatest numbers of DEGs were identified in the 35/1 and 5/1 comparisons in both leaves and roots. In total, 566 leaf and 643 root DEGs were common to all 4 genotypes, though a third of these leaf DEGs were not regulated in the same up/down direction in all 4 genotypes. For roots, the equivalent figure was 1% of the DEGs. GO terms shared by all four genotypes were often enriched by both up- and down-regulated DEGs in the leaf, whereas generally, only by either up- or down-regulated DEGs in the root. Overall, up-regulated leaf DEGs tended to be more genotype-specific than down-regulated leaf DEGs or root DEGs and were also associated with fewer GOs. On average, only 5–15% of the DEGs enriching common GO terms were shared by all 4 genotypes, suggesting considerable variation in DEGs between related genotypes in enacting similar biological processes.

Population structure and genome-wide association studies in bread wheat for phosphorus efficiency traits using 35 K Wheat Breeder's Affymetrix array

Soil bioavailability of phosphorus (P) is a major concern for crop productivity worldwide. As phosphatic fertilizers are a non-renewable resource associated with economic and environmental issues so, the sustainable option is to develop P use efficient crop varieties. We phenotyped 82 diverse wheat (Triticum aestivum L.) accessions in soil and hydroponics at low and sufficient P. To identify the genic regions for P efficiency traits, the accessions were genotyped using the 35 K-SNP array and genome-wide association study (GWAS) was performed. The high-quality SNPs across the genomes were evenly distributed with polymorphic information content values varying between 0.090 and 0.375. Structure analysis revealed three subpopulations (C1, C2, C3) and the phenotypic responses of these subpopulations were assessed for P efficiency traits. The C2 subpopulation showed the highest genetic variance and heritability values for numerous agronomically important traits as well as strong correlation under both P levels in soil and hydroponics. GWAS revealed 78 marker-trait associations (MTAs) but only 35 MTAs passed Bonferroni Correction. A total of 297 candidate genes were identified for these MTAs and their annotation suggested their involvement in several biological process. Out of 35, nine (9) MTAs were controlling polygenic trait (two controlling four traits, one controlling three traits and six controlling two traits). These multi-trait MTAs (each controlling two or more than two correlated traits) could be utilized for improving bread wheat to tolerate low P stress through marker-assisted selection (MAS).

Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation

Knowing how switchgrass (Panicum virgatum L.) responds and adapts to phosphorus (P)-limitation will aid efforts to optimize P acquisition and use in this species for sustainable biomass production. This integrative study investigated the impacts of mild, moderate, and severe P-stress on genome transcription and whole-plant metabolism, physiology and development in switchgrass. P-limitation reduced overall plant growth, increased root/shoot ratio, increased root branching at moderate P-stress, and decreased root diameter with increased density and length of root hairs at severe P-stress. RNA-seq analysis revealed thousands of genes that were differentially expressed under moderate and severe P-stress in roots and/or shoots compared to P-replete plants, with many stress-induced genes involved in transcriptional and other forms of regulation, primary and secondary metabolism, transport, and other processes involved in P-acquisition and homeostasis. Amongst the latter were multiple miRNA399 genes and putative targets of these. Metabolite profiling showed that levels of most sugars and sugar alcohols decreased with increasing P stress, while organic and amino acids increased under mild and moderate P-stress in shoots and roots, although this trend reversed under severe P-stress, especially in shoots.