Author summary Plants thrive in highly heterogenous soils. How they compute a multitude of contrasting stimuli and mount an adaptive response without a centralized information processing unit is an intriguing question. For instance, below ground, roots can sense and respond to the single or multiple nutrient stresses, and adjust its growth rate accordingly. Nevertheless, the genetic architecture of root growth responses under single and combined stress remains poorly understood. To fill this gap in our understanding about such crucial phenomenon for plant survival, we explored the natural variation of root growth rate (RGR) in Arabidopsis grown under single and combined nutritional stress, including deficiencies of iron (-Fe), zinc (-Zn), phosphate and iron (-P-Fe) and phosphate and zinc (-P-Zn). Our GWAS revealed distinct genetic architectures underlying root growth responses to single or combined nutrient stresses. By integrating GWAS and coexpression networks, we identified and validated genes controlling the variation of root growth to combined nutrient-deficiency, namely VARIANT IN METHYLATION 1, FORMIN-LIKE-PROTEIN-6 and VOLTAGE-DEPENDENT ANION-SELECTIVE CHANNEL PROTEIN 3. Our findings provide a framework to accelerate future research aiming at better understanding how plants sense and respond to multiple environmental inputs, and promise to help designing new agronomical and biotechnological strategies to improve root growth.

Senna tora is a widely used medicinal plant. Its health benefits have been attributed to the large quantity of anthraquinones, but how they are made in plants remains a mystery. To identify the genes responsible for plant anthraquinone biosynthesis, we reveal the genome sequence of S. tora at the chromosome level with 526 Mb (96%) assembled into 13 chromosomes. Comparison among related plant species shows that a chalcone synthase-like (CHS-L) gene family has lineage-specifically and rapidly expanded in S. tora. Combining genomics, transcriptomics, metabolomics, and biochemistry, we identify a CHS-L gene contributing to the biosynthesis of anthraquinones. The S. tora reference genome will accelerate the discovery of biologically active anthraquinone biosynthesis pathways in medicinal plants. Anthraquinones are aromatic polyketides and have been used for treating various diseases, but the biosynthetic pathway is unclear. Here, the authors assemble the genome of an anthraquinone-producing medicinal plant Senna tora and show the evidences that CHS-like genes may be involved in anthraquinone biosynthesis.

Linkage mapping has been widely used to identify quantitative trait loci (QTL) in many plants and usually requires a time-consuming and labor-intensive fine mapping process to find the causal gene underlying the QTL. Previously, we described QTG-Finder, a machine-learning algorithm to rationally prioritize candidate causal genes in QTLs. While it showed good performance, QTG-Finder could only be used in Arabidopsis and rice because of the limited number of known causal genes in other species. Here we tested the feasibility of enabling QTG-Finder to work on species that have few or no known causal genes by using orthologs of known causal genes as the training set. The model trained with orthologs could recall about 64% of Arabidopsis and 83% of rice causal genes when the top 20% ranked genes were considered, which is similar to the performance of models trained with known causal genes. The average precision was 0.027 for Arabidopsis and 0.029 for rice. We further extended the algorithm to include polymorphisms in conserved non-coding sequences and gene presence/absence variation as additional features. Using this algorithm, QTG-Finder2, we trained and cross-validatedSorghum bicolorandSetaria viridismodels. TheS. bicolormodel was validated by causal genes curated from the literature and could recall 70% of causal genes when the top 20% ranked genes were considered. In addition, we applied theS. viridismodel and public transcriptome data to prioritize a plant height QTL and identified 13 candidate genes. QTL-Finder2 can accelerate the discovery of causal genes in any plant species and facilitate agricultural trait improvement.

In this Perspective, Diane Dickel and colleagues review recent progress and opportunities in applying single-cell sequencing and microfluidics methods to plants. The authors highlight the need for new tools developed with plants in mind, and advocate for the creation of a centralized, open-access database to house plant single-cell data.

We highlight the achievements of four pioneering Black plant scientists to raise awareness of the importance of diversity, equity, and inclusion in science. Their stories come alive at Historically Black Colleges and Universities through exhibits of science and art and classroom activities (https://www.plantcellatlas.org/pca-art-exhibit.html).

Temporal dynamics of gene expression underpin responses to internal and environmental stimuli. In eukaryotes, regulation of gene induction includes changing chromatin states at target genes and recruiting the transcriptional machinery that includes transcription factors. As one of the most potent defense compounds in Arabidopsis thaliana, camalexin can be rapidly induced by bacterial and fungal infections. Though several transcription factors controlling camalexin biosynthesis genes have been characterized, how the rapid activation of genes in this pathway upon a pathogen signal is enabled remains unknown. By combining publicly available epigenomic data with in vivo chromatin modification mapping, we found that camalexin biosynthesis genes are marked with two epigenetic modifications with opposite effects on gene expression, trimethylation of lysine 27 of histone 3 (H3K27me3) (repression) and acetylation of lysine 18 of histone 3 (H3K18ac) (activation), to form a previously uncharacterized type of bivalent chromatin. Mutants with reduced H3K27me3 or H3K18ac suggested that both modifications were required to determine the timing of gene expression and metabolite accumulation at an early stage of the stress response. Our study indicates that the H3K27me3-H3K18ac bivalent chromatin, which we name as kairostat, plays an important role in controlling the timely induction of gene expression upon stress stimuli in plants.

Despite the ecological and agronomical importance of seed germination, how seeds integrate environmental signals to trigger germination remains enigmatic. Recently we reported that a protein called FLOE1 is involved in sensing and responding to water availability during germination. Here, we present a live-imaging protocol to assess the subcellular localization of a protein of interest during imbibition of desiccated Arabidopsis thaliana seeds with the goal of understanding protein dynamics during the early stages of water uptake. For complete details on the use and execution of this profile, please refer to Dorone etal. (2021).

Iron deficiency hampers photosynthesis and is associated with chlorosis. We recently showed that iron deficiency-induced chlorosis depends on phosphorus availability. How plants integrate these cues to control chlorophyll accumulation is unknown. Here, we show that iron limitation downregulates photosynthesis genes in a phosphorus-dependent manner. Using transcriptomics and genome-wide association analysis, we identify two genes, PHT4;4 encoding a chloroplastic ascorbate transporter and bZIP58, encoding a nuclear transcription factor, which prevent the downregulation of photosynthesis genes leading to the stay-green phenotype under iron-phosphorus deficiency. Joint limitation of these nutrients induces ascorbate accumulation by activating expression of an ascorbate biosynthesis gene, VTC4, which requires bZIP58. Furthermore, we demonstrate that chloroplastic ascorbate transport prevents the downregulation of photosynthesis genes under iron-phosphorus combined deficiency through modulation of ROS homeostasis. Our study uncovers a ROS-mediated chloroplastic retrograde signaling pathway to adapt photosynthesis to nutrient availability.