Root system architecture affects water and mineral uptake and is important for plant adaptation to fluctuating nutrient availability. Small signaling peptides and their receptors influence root traits associated with macronutrient uptake. In this study, genome-wide association analyses were performed using 2D images of agar plate-grown Medicago truncatula accessions to understand the impact of GOLVEN10 peptide (GLV10) treatment on three root traits: root tortuosity, lateral root (LR) branch angle, and the gravity setpoint angle (GSA). Upon GLV10 treatment, roots of wild-type M. truncatula Jemalong A17 and R108 accessions showed increased primary root coiling (or tortuosity), increased LR branch angle, and reduced GSA. We identified 88 significant single nucleotide polymorphisms (SNPs) associated with these traits in GLV10-treated plants, distinct from the 163 SNPs in untreated plants. Importantly, the ethylene regulatory pathway was implicated in root tortuosity and LR emergence relative to the primary root. Application of the ethylene precursor 1-aminocyclopropane-1-carboxylic acid reduced root sensitivity to GLV10, while the ethylene signaling mutant sickle was hypersensitive, indicating that GLV10 and ethylene pathways act antagonistically to control root tortuosity. These findings have implications for root gravitropic responses, and the ability of roots to penetrate deeper soil layers for nutrients and water.

Translating biological knowledge from Arabidopsis to crop species is important to advance agriculture and secure food production in the face of dwindling fertilizer resources and biotic and abiotic stresses. However, it is often not trivial to identify functional homologs (orthologs) of Arabidopsis genes in crops. Combining sequence and expression data can improve the correct prediction of orthologs. Here, we conducted a large-scale RNA sequencing based transcriptomics study of Arabidopsis, Medicago, Brachypodium, and Setaria grown side-by-side in phosphorus (P)-sufficient and P-limited conditions to generate comparable transcriptomics datasets. Comparison of top 200 P-limitation-induced genes in Arabidopsis revealed that ~80% of these genes have identifiable close homologs in the other three species but only ~50% retain their P-limitation response in the legume and grasses. Most of the hallmark genes of the P-starvation response were found conserved in all four species. This study reveals many known, novel, unannotated, conserved, and species-specific forms of regulation of the transcriptional P-starvation response. Identification and experimental verification of expressologs by independent RT-qPCR for P-limitation marker genes in Prunus showed the usefulness of comparative transcriptomics in pinpointing the functional orthologs in diverse crop species. This study provides an unprecedented resource for functional genomics and translational research to create P-efficient crops.

Background and aims Hakea prostrata (Proteaceae) is a highly phosphorus-use-efficient plant native to southwest Australia. It maintains a high photosynthetic rate at low leaf phosphorus (P) and exhibits delayed leaf greening, a convergent adaptation that increases nutrient-use efficiency. This study aimed to provide broad physiological and gene expression profiles across leaf development, uncovering pathways leading from young leaves as nutrient sinks to mature leaves as low-nutrient, energy-transducing sources. Methods To explore gene expression underlying delayed greening, we analysed a de novo transcriptome for H. prostrata across five stages of leaf development. Photosynthesis and respiration rates, and foliar pigment, P and nitrogen (N) concentrations were determined, including the division of P into five biochemical fractions. Key results Transcripts encoding functions associated with leaf structure generally decreased in abundance across leaf development, concomitant with decreases in foliar concentrations of 85% for anthocyanins, 90% for P and 70% for N. The expression of genes associated with photosynthetic function increased during or after leaf expansion, in parallel with increases in photosynthetic pigments and activity, much later in leaf development than in species that do not have delayed greening. As leaves developed, transcript abundance for cytosolic and mitochondrial ribosomal proteins generally declined, whilst transcripts for chloroplast ribosomal proteins increased. Conclusions There was a much longer temporal separation of leaf cell growth from chloroplast development in H. prostrata than is found in species that lack delayed greening. Transcriptome-guided analysis of leaf development in H. prostrata provided insight into delayed greening as a nutrient-saving strategy in severely phosphorus-impoverished landscapes.

GOLVEN/ROOT MERISTEM GROWTH FACTOR family of signaling peptides have been shown to control root lateral organ number, density and positioning in plants, although the signaling pathways involved remain obscure. A diverse set of 171 Medicago truncatula HapMap accessions with variation in responses to the GOLVEN 10 peptide, GLV10, were used to identify 74 significant loci controlling seven traits related to nodule formation and root architecture. Importantly, a single nucleotide polymorphism (SNP) in the upstream region of the MtGLV10 peptide-inducible strigolactone receptor gene, MtDWARF14 was significantly associated with insensitivity of nodule density to GLV10, suggesting a link between strigolactone signaling and GLV10 responsiveness. Three independent d14 mutants of the D WARF14 gene were found to hypernodulate, while overexpression of the gene led to reduction in nodule number, phenocopying GLV10. A null mutant, mtd14-1 , remained sensitive to GLV10’s effect on nodule density. However, at the transcriptional level, the mutant failed to effectively induce the expression of the GOLVEN marker genes, MtPLETHORA3 and MtPINLIKES2 . Our study uncovers a hitherto unknown link between the strigolactone and GLV peptide signaling pathways using genotype x environment analysis of Medicago HapMap lines and provides a putative molecular mechanism for recovery from frost damage to fine roots. ### Competing Interest Statement The authors have declared no competing interest.

Translating biological knowledge from Arabidopsis to crop species is important to advance agriculture and secure food production in the face of dwindling fertilizer resources, biotic and abiotic stresses. However, it is often not trivial to identify functional homologs (orthologs) of Arabidopsis genes in crops. Combining sequence and expression data can improve the correct prediction of orthologs. Here, we conducted a large-scale RNA sequencing based transcriptomics study of Arabidopsis, Medicago, Brachypodium and Setaria grown side-by-side in Phosphorus (P)-sufficient and P-limited conditions to generate comparable transcriptomics datasets. Comparison of top 200 P-limitation induced genes in Arabidopsis revealed that ~80% of these genes have identifiable close homologs in the other three species but only ~50% retain their P-limitation response in the legume and grasses. Most of the hallmark genes of the P-starvation response were found conserved in all four species. This study reveals many known, novel, unannotated, conserved and species-specific regulations of the transcriptional P-starvation response. Identification and experimental verification of expressologs by independent RT-qPCR for P-limitation marker genes in Prunus showed the usefulness of comparative transcriptomics in pinpointing the functional orthologs in diverse crop species. This study provides an unprecedented resource for functional genomics and translational research to create P-efficient crops. ### Competing Interest Statement The authors have declared no competing interest.

The conservation of GOLVEN (GLV)/ROOT MERISTEM GROWTH FACTOR (RGF) peptide encoding genes across plant genomes capable of forming roots or root-like structures underscores their potential significance in the terrestrial adaptation of plants. This study investigates the function and role of GOLVEN peptide-coding genes in Medicago truncatula. Five out of fifteen GLV/RGF genes were notably upregulated during nodule organogenesis and were differentially responsive to nitrogen deficiency and auxin treatment. Specifically, the expression of MtGLV9 and MtGLV10 at nodule initiation sites was contingent upon the NODULE INCEPTION transcription factor. Overexpression of these five nodule-induced GLV genes in hairy roots of M. truncatula and application of their synthetic peptide analogues led to a decrease in nodule count by 25-50%. Uniquely, the GOLVEN10 peptide altered the positioning of the first formed lateral root and nodule on the primary root axis, an observation we term 'noduletaxis'; this decreased the length of the lateral organ formation zone on roots. Histological section of roots treated with synthetic GOLVEN10 peptide revealed an increased cell number within the root cortical cell layers without a corresponding increase in cell length, leading to an elongation of the root likely introducing a spatiotemporal delay in organ formation. At the transcription level, the GOLVEN10 peptide suppressed expression of microtubule-related genes and exerted its effects by changing expression of a large subset of Auxin responsive genes. These findings advance our understanding of the molecular mechanisms by which GOLVEN peptides modulate root morphology, nodule ontogeny, and interactions with key transcriptional pathways.

Three PHO2-like genes encoding putative ubiquitin-conjugating E2 enzymes of Medicago truncatula were characterized for potential roles in phosphorous (P) homeostasis and symbiotic nitrogen fixation (SNF). All three genes, MtPHO2A, B and C, contain miR399-binding sites characteristic of PHO2 genes in other plant species. Distinct spatiotemporal expression patterns and responsiveness of gene expression to P- and N-deprivation in roots and shoots indicated potential roles, especially for MtPHO2B, in P and N homeostasis. Phenotypic analysis of pho2 mutants revealed that MtPHO2B is integral to Pi homeostasis, affecting Pi allocation during plant growth under nutrient-replete conditions, while MtPHO2C had a limited role in controlling Pi homeostasis. Genetic analysis also revealed a connection between Pi allocation, plant growth and SNF performance. Under N-limited, SNF conditions, Pi allocation to different organs was dependent on MtPHO2B and, to a lesser extent, MtPHO2C and MtPHO2A. MtPHO2A also affected Pi homeostasis associated with nodule formation. Thus, MtPHO2 genes play roles in systemic and localized, i.e., nodule, P homeostasis affecting SNF.

Abstract To cope with environmental stresses and ensure maximal reproductive success, plants have developed strategies to adjust the timing of their transition to reproductive growth. This has a substantial impact on the stress resilience of crops and ultimately on agricultural productivity. Here, we report a previously uncharacterized, plant-specific gene family designated as Regulator of Flowering and Stress (RFS). Overexpression of the BdRFS gene in Brachypodium distachyon delayed flowering, increased biomass accumulation, and promoted drought tolerance, whereas clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated knockout mutants exhibited opposite phenotypes. A double T-DNA insertional mutant in the two Arabidopsis (Arabidopsis thaliana) homologs replicated the effects on flowering and water deprivation seen in the B. distachyon CRISPR knockout lines, highlighting the functional conservation of the family between monocots and dicots. Lipid analysis of B. distachyon and Arabidopsis revealed that digalactosyldiacylglycerol (DGDG) and phosphatidylcholine (PC) contents were significantly, and reciprocally, altered in overexpressor and knockout mutants. Importantly, alteration of C16:0-containing PC, a Flowering Locus T-interacting lipid, associated with flowering phenotype, with elevated levels corresponding to earlier flowering. Co-immunoprecipitation analysis suggested that BdRFS interacts with phospholipase Dα1 as well as several other abscisic acid-related proteins. Furthermore, reduction of C18:3 fatty acids in DGDG corresponded with reduced jasmonic acid metabolites in CRISPR mutants. Collectively, we suggest that stress-inducible RFS proteins represent a regulatory component of lipid metabolism that impacts several agronomic traits of biotechnological importance.

Stomata are crucial for gas exchange and water evaporation, and environmental stimuli influence their density (SD) and size (SS). Although genes and mechanisms underlying stomatal development have been elucidated, stress-responsive regulators of SD and SS are less well-known. Previous studies have shown that the stress-inducible Brachypodium RFS (REGULATOR OF FLOWERING AND STRESS, BdRFS) gene affects heading time and enhances drought tolerance by reducing leaf water loss. Here, we report that overexpression lines (OXs) of BdRFS have reduced SD and increased SS, regardless of soil water status. Furthermore, biomass and plant water content of OXs were significantly increased compared to wild type. CRISPR/Cas9-mediated BdRFS knockout mutant (KO) exhibited the opposite stomatal characteristics and biomass changes. Reverse transcription-quantitative polymerase chain reaction analysis revealed that expression of BdICE1 was reversely altered in OXs and KO, pointing to a potential cause for the observed changes in stomatal phenotypes. Stomatal and transcriptional changes were not observed in the Arabidopsis rfs double mutant. Taken together, RFS is a novel regulator of SD and SS and is a promising candidate for genetic engineering of climate-resilient crops.

Phosphorus (P) limitation in the majority of world soils is a major constraint for plant growth and crop productivity. RNA sequencing was used to discover novel P-responsive gene transcripts (PRGT) in leaves and roots of Arabidopsis. Hisat StringTie and Cufflinks TopHat transcript assembler were used to analyze reads and identify 1,074 PRGTs with a >5-fold altered abundance during P-limitation. Interestingly, 60% of these transcripts were not previously reported. Among the novel PRGT, 106 were from unannotated genes, and some were among the most P-responsive, including At2g36727 which encodes a novel microRNA. Annotated novel PRGTs encode for transcription factors, microRNAs, small signaling peptides, long non-coding RNAs, defense-related proteins, and transporters, along with proteins involved in many biological processes. We identified several genes that undergo alternative splicing during P-limitation, including a novel miR399 resistant splice variant of PHOSPHATE2 ( PHO2 . 2 ). Several novel P-responsive genes were regulated by PHOSPHATE STARVATION RESPONSE1 ( PHR1 ), PHR1-LIKE 1 ( PHL1 ) and PHO2 . We discovered that P-limited plants show increased resistance to pathogens and drought stress mediated by PHR1-PHL1. Identification of novel P-responsive transcripts and the discovery of the influence of P-limitation on biotic and abiotic stress adds a significant component to our understanding of plant P-signaling. Highlight Phosphorus limitation elicits the expression of several novel genes including many previously unannotated genes, noncoding RNAs, small peptides and alternatively spliced RNAs, and leads to enhanced disease and drought tolerance.

The root system of a plant provides vital functions including resource uptake, storage, and anchorage in soil. The uptake of macro-nutrients like nitrogen (N), phosphorus (P), potassium (K), and sulphur (S) from the soil is critical for plant growth and development. Small signaling peptide (SSP) hormones are best known as potent regulators of plant growth and development with a few also known to have specialized roles in macronutrient utilization. Here we describe a high throughput phenotyping platform for testing SSP effects on root uptake of multiple nutrients. The SSP, CEP1 (C-TERMINALLY ENCODED PEPTIDE) enhanced nitrate uptake rate per unit root length in Medicago truncatula plants deprived of N in the high-affinity transport range. Single structural variants of M. truncatula and Arabidopsis thaliana specific CEP1 peptides, MtCEP1D1:hyp4,11 and AtCEP1:hyp4,11, enhanced uptake not only of nitrate, but also phosphate and sulfate in both model plant species. Transcriptome analysis of Medicago roots treated with different MtCEP1 encoded peptide domains revealed that hundreds of genes respond to these peptides, including several nitrate transporters and a sulfate transporter that may mediate the uptake of these macronutrients downstream of CEP1 signaling. Likewise, several putative signaling pathway genes including LEUCINE-RICH REPEAT RECPTOR-LIKE KINASES and Myb domain containing transcription factors, were induced in roots by CEP1 treatment. Thus, a scalable method has been developed for screening synthetic peptides of potential use in agriculture, with CEP1 shown to be one such peptide.

SUMMARYGLV/RGF peptide encoding genes can be identified in genomes of all plants that can form roots or root-like structures suggesting they were essential for transition of plants to land.In Medicago truncatula, five of fifteen GOLVEN(GLV)/ROOT MERISTEM GROWTH FACTOR (RGF) peptide coding genes were induced during nodule organogenesis and to a varying extent under nitrogen deficiency and auxin treatment. Expression of MtGLV9 and MtGLV10 at nodule initiation sites was dependent on the transcription factor NODULE INCEPTION.Overexpression of all five nodule-induced GLV genes in M. truncatula hairy roots as well as application of the corresponding synthetic peptides resulted in a 25-50% reduction in nodule number indicating GOLVENs are negative regulators of nodule organogenesis.The peptide GOLVEN10 shifted the position of the first formed lateral root (rhizotaxis) as well as the first formed nodule along the longitudinal primary root axis, a phenomenon we term ‘nodulotaxis’, thereby reducing the absolute length of the zone of lateral organ formation on roots.Application of synthetic GOLVEN10 peptide caused an increase in cell number but not cell length in each root cortical cell layer causing an increase in root length and a consequent spatiotemporal delay in formation of the first lateral organ.Plain Language SummaryNodule positioning is an understudied trait, yet it determines the length of the root that can support nodule formation and consequently the total number of functional nodules formed. We identify for the first time, genetic factors called GOLVEN peptides that alter nodule and lateral root positioning on the primary root along with several other traits including nodule organ initiation and root architecture.

Root hairs (RHs) function in nutrient and water acquisition, root metabolite exudation, soil anchorage and plant-microbe interactions. Longer or more abundant RHs are potential breeding traits for developing crops that are more resource-use efficient and can improve soil health. While many genes are known to promote RH elongation, relatively little is known about genes and mechanisms that constrain RH growth. Here we demonstrate that a DOMAIN OF UNKNOWN FUNCTION 506 (DUF506) protein, AT3G25240, negatively regulates Arabidopsis thaliana RH growth. The AT3G25240 gene is strongly and specifically induced during phosphorus (P)-limitation. Mutants of this gene, which we call REPRESSOR OF EXCESSIVE ROOT HAIR ELONGATION 1 (RXR1), have much longer RHs, higher phosphate content and seedling biomass, while overexpression of the gene exhibits opposite phenotypes. Co-immunoprecipitation, pull-down and bimolecular fluorescence complementation (BiFC) analyses reveal that RXR1 physically interacts with a RabD2c GTPase in nucleus, and a rabd2c mutant phenocopies the rxr1 mutant. Furthermore, N-terminal variable region of RXR1 is crucial for inhibiting RH growth. Overexpression of a Brachypodium distachyon RXR1 homolog results in repression of RH elongation in Brachypodium. Taken together, our results reveal a novel DUF506-GTPase module with a prominent role in repression of plant RH elongation especially under P stress.

The perennial native switchgrass adapts better than other plant species do to marginal soils with low plant-available nutrients, including those with low phosphorus (P) content. Switchgrass roots and their associated microorganisms can alter the pools of available P throughout the whole soil profile making predictions of P availability in situ challenging. Plant P homeostasis makes monitoring of P limitation via measurements of plant P content alone difficult to interpret. To address these challenges, we developed a machine-learning model trained with high accuracy using the leaf tissue chemical profile, rather than P content. By applying this learned model in field trials across two sites with contrasting extractable soil P, we observed that actual plant available P in soil was more similar than expected, suggesting that adaptations occurred to alleviate the apparent P constraint. These adaptations come at a metabolic cost to the plant that have consequences for feedstock chemical components and quality. We observed that other biochemical signatures of P limitation, such as decreased cellulose-to-lignin ratios, were apparent, indicating re-allocation of carbon resources may have contributed to increased P acquisition. Plant P allocation strategies also differed across sites, and these differences were correlated with the subsequent year's biomass yields.

Domain of Unknown Function 506 proteins are ubiquitous in plants. The phosphorus (P) stress-inducible REPRESSOR OF EXCESSIVE ROOT HAIR GROWTH1 (AtRXR1) gene encodes the first characterized DUF506. AtRXR1 inhibits root hair elongation by interacting with RabD2c GTPase. However, functions of other P-responsive DUF506 genes are still missing. Here, we selected two additional P-inducible DUF506 genes for further investigation. The expression of both genes was induced by auxin. Under P-stress, At3g07350 gene expressed ubiquitously in seedlings, whereas At1g62420 (AtRXR3) expression was strongest in roots. AtRXR3 overexpressors and knockouts had shorter and longer root hairs, respectively. A functional AtRXR3-green fluorescent protein fusion localized to root epidermal cells. Chromatin immunoprecipitation and quantitative reverse-transcriptase-polymerase chain reaction revealed that AtRXR3 was transcriptionally activated by RSL4. Bimolecular fluorescence complementation and calmodulin (CaM)-binding assays showed that AtRXR3 interacted with CaM in the presence of Ca2+ . Moreover, cytosolic Ca2+ ([Ca2+ ]cyt ) oscillations in root hairs of rxr3 mutants exhibited elevated frequencies and dampened amplitudes compared to those of wild type. Thus, AtRXR3 is another DUF506 protein that attenuates P-limitation-induced root hair growth through mechanisms that involve RSL4 and interaction with CaM to modulate tip-focused [Ca2+ ]cyt oscillations.

The HapMap (haplotype map) projects have produced valuable genetic resources in life science research communities, allowing researchers to investigate sequence variations and conduct genome-wide association study (GWAS) analyses. A typical HapMap project may require sequencing hundreds, even thousands, of individual lines or accessions within a species. Due to limitations in current sequencing technology, the genotype values for some accessions cannot be clearly called. Additionally, allelic heterozygosity can be very high in some lines, causing genetic and sometimes phenotypic segregation in their descendants. Genetic and phenotypic segregation degrades the original accession’s specificity and makes it difficult to distinguish one accession from another. Therefore, it is vitally important to determine and validate HapMap accessions before one conducts a GWAS analysis. However, to the best of our knowledge, there are no prior methodologies or tools that can readily distinguish or validate multiple accessions in a HapMap population. We devised a bioinformatics approach to distinguish multiple HapMap accessions using only a minimum number of genetic markers. First, we assign each candidate marker with a distinguishing score (DS), which measures its capability in distinguishing accessions. The DS score prioritizes those markers with higher percentages of homozygous genotypes (allele combinations), as they can be stably passed on to offspring. Next, we apply the “set-partitioning” concept to select optimal markers by recursively partitioning accession sets. Subsequently, we build a hierarchical decision tree in which a specific path represents the selected markers and the homogenous genotypes that can be used to distinguish one accession from others in the HapMap population. Based on these algorithms, we developed a web tool named MAD-HiDTree (Multiple Accession Distinguishment-Hierarchical Decision Tree), designed to analyze a user-input genotype matrix and construct a hierarchical decision tree. Using genetic marker data extracted from the Medicago truncatula HapMap population, we successfully constructed hierarchical decision trees by which the original 262 M. truncatula accessions could be efficiently distinguished. PCR experiments verified our proposed method, confirming that MAD-HiDTree can be used for the identification of a specific accession. MAD-HiDTree was developed in C/C++ in Linux. Both the source code and test data are publicly available at https://bioinfo.noble.org/MAD-HiDTree/.

Virus-induced gene silencing (VIGS) is a rapid and powerful method to evaluate gene function, especially for species like hexaploid wheat that have large, redundant genomes and are difficult and time-consuming to transform. The Brome mosaic virus (BMV)-based VIGS vector is widely used in monocotyledonous species but not wheat. Here we report the establishment of a simple and effective VIGS procedure in bread wheat using BMVCP5, the most recently improved BMV silencing vector, and wheat genes PHYTOENE DESATURASE ( TaPDS ) and PHOSPHATE2 ( TaPHO2 ) as targets. Time-course experiments revealed that smaller inserts (~100 nucleotides, nt) were more stable in BMVCP5 and conferred higher silencing efficiency and longer silencing duration, compared with larger inserts. When using a 100-nt insert and a novel coleoptile inoculation method, BMVCP5 induced extensive silencing of TaPDS transcript and a visible bleaching phenotype in the 2nd to 5th systemically-infected leaves from nine to at least 28 days post inoculation (dpi). For TaPHO2 , the ability of BMVCP5 to simultaneously silence all three homoeologs was demonstrated. To investigate the feasibility of BMV VIGS in wheat roots, ectopically expressed enhanced GREEN FLUORESCENT PROTEIN ( eGFP ) in a transgenic wheat line was targeted for silencing. Silencing of eGFP fluorescence was observed in both the maturation and elongation zones of roots. BMVCP5 mediated significant silencing of eGFP and TaPHO2 mRNA expression in roots at 14 and 21 dpi, and TaPHO2 silencing led to the doubling of inorganic phosphate concentration in the 2nd through 4th systemic leaves. All 54 wheat cultivars screened were susceptible to BMV infection. BMVCP5-mediated TaPDS silencing resulted in the expected bleaching phenotype in all eight cultivars examined, and decreased TaPDS transcript was detected in all three cultivars examined. This BMVCP5 VIGS technology may serve as a rapid and effective functional genomics tool for high-throughput gene function studies in aerial and root tissues and in many wheat cultivars.

DUF506 proteins are omnipresent in higher plants. Phosphorus (P) stress-inducible AtRXR1/REPRESSOR OF EXCESSIVE ROOT HAIR ELONGATION 1 gene, as the first functionally characterized DUF506 gene, is proved to inhibit root hair elongation through interaction of RXR2/RabD2c GTPase. However, the knowledge of other P-responsive DUF506 genes is still limited. Here, we identify four additional P-inducible DUF506 genes and select two of candidates for functional investigation. Expression analysis results reveal that both of candidates are induced by auxin. At3g07350, the duplicated pair of RXR1, expresses ubiquitously in seedlings under P-stress, whereas At1g62420 is mainly induced in roots. Overexpression and knockout mutants of At1g62420, which is called RXR3, exhibit shorter or longer root hair length, respectively. Cellular imaging results demonstrate RXR3 localizes in root epidermal cells. ChIP, synthetic peptide treatment and qRT-PCR assay results indicate RXR3 is transcriptionally activated by RSL4 or RALF1. BiFC and CaM-binding assay suggest that RXR3 interacts with various CaMs in presence of Ca2+. Moreover, the frequencies of [Ca2+]cyt in rxr3 mutants are approximately 20% higher compared to that of wild type. Taken together, our results illustrate a divergent mechanism by which RSL4-directed RXR3 interacts with calmodulin and functions in parallel of RXR1, to prevent root hair excessive growth.

A growing number of small secreted peptides (SSPs) in plants are recognized as important regulatory molecules with roles in processes such as growth, development, reproduction, stress tolerance, and pathogen defense. Recent discoveries further implicate SSPs in regulating root nodule development, which is of particular significance for legumes. SSP-coding genes are frequently overlooked, because genome annotation pipelines generally ignore small open reading frames, which are those most likely to encode SSPs. Also, SSP-coding small open reading frames are often expressed at low levels or only under specific conditions, and thus are underrepresented in non-tissue-targeted or non-condition-optimized RNA-sequencing projects. We previously identified 4,439 SSP-encoding genes in the model legume Medicago truncatula To support systematic characterization and annotation of these putative SSP-encoding genes, we developed the M. truncatula Small Secreted Peptide Database (MtSSPdb; https://mtsspdb.noble.org/). MtSSPdb currently hosts (1) a compendium of M. truncatula SSP candidates with putative function and family annotations; (2) a large-scale M. truncatula RNA-sequencing-based gene expression atlas integrated with various analytical tools, including differential expression, coexpression, and pathway enrichment analyses; (3) an online plant SSP prediction tool capable of analyzing protein sequences at the genome scale using the same protocol as for the identification of SSP genes; and (4) information about a library of synthetic peptides and root and nodule phenotyping data from synthetic peptide screens in planta. These datasets and analytical tools make MtSSPdb a unique and valuable resource for the plant research community. MtSSPdb also has the potential to become the most complete database of SSPs in plants.

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.