Identification of bZIP Transcription Factors That Regulate the Development of Leaf Epidermal Cells in Arabidopsis thaliana by Single-Cell RNA Sequencing.

Epidermal cells are the main avenue for signal and material exchange between plants and the environment. Leaf epidermal cells primarily include pavement cells, guard cells, and trichome cells. The development and distribution of different epidermal cells are tightly regulated by a complex transcriptional regulatory network mediated by phytohormones, including jasmonic acid, and transcription factors. How the fate of leaf epidermal cells is determined, however, is still largely unknown due to the diversity of cell types and the complexity of their regulation. Here, we characterized the transcriptional profiles of epidermal cells in 3-day-old true leaves of Arabidopsis thaliana using single-cell RNA sequencing. We identified two genes encoding BASIC LEUCINE-ZIPPER (bZIP) transcription factors, namely bZIP25 and bZIP53, which are highly expressed in pavement cells and early-stage meristemoid cells. Densities of pavement cells and trichome cells were found to increase and decrease, respectively, in bzip25 and bzip53 mutants, compared with wild-type plants. This trend was more pronounced in the presence of jasmonic acid, suggesting that these transcription factors regulate the development of trichome cells and pavement cells in response to jasmonic acid.

Breaking new ground: Decoding the root's molecular circuits to penetrate compacted soil.

Penetration capacity is a central root phenotype for enhancing rooting into compacted soils and improving abiotic stress tolerance in crop plants, but the molecular mechanisms are unclear. In this issue of Developmental Cell, Xu et al. demonstrate the vital role of the FER-PIF3-PIEZO circuitry in controlling primary root penetrability into compacted substrates.

The Shoot Apical Meristem: An Evolutionary Molding of Higher Plants.

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and other meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperms are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom and determines the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

A non-invasive method to predict drought survival in Arabidopsis using quantum yield under light conditions.

BACKGROUND: Survival rate (SR) is frequently used to compare drought tolerance among plant genotypes. While a variety of techniques for evaluating the stress status of plants under drought stress conditions have been developed, determining the critical point for the recovery irrigation to evaluate plant SR often relies directly on a qualitative inspection by the researcher or on the employment of complex and invasive techniques that invalidate the subsequent use of the tested individuals. RESULTS: Here, we present a simple, instantaneous, and non-invasive method to estimate the survival probability of Arabidopsis thaliana plants after severe drought treatments. The quantum yield (QY), or efficiency of photosystem II, was monitored in darkness (Fv/Fm) and light (Fv'/Fm') conditions in the last phase of the drought treatment before recovery irrigation. We found a high correlation between a plant's Fv'/Fm' value before recovery irrigation and its survival phenotype seven days after, allowing us to establish a threshold between alive and dead plants in a calibration stage. This correlation was maintained in the Arabidopsis accessions Col-0, Ler-0, C24, and Kondara under the same conditions. Fv'/Fm' was then applied as a survival predictor to compare the drought tolerance of transgenic lines overexpressing the transcription factors ATAF1 and PLATZ1 with the Col-0 control. CONCLUSIONS: The results obtained in this work demonstrate that the chlorophyll a fluorescence parameter Fv'/Fm' can be used as a survival predictor that gives a numerical estimate of the Arabidopsis drought SR before recovery irrigation. The procedure employed to get the Fv'/Fm' measurements is fast, non-destructive, and requires inexpensive and easy-to-handle equipment. Fv'/Fm' as a survival predictor can be used to offer an overview of the photosynthetic state of the tested plants and determine more accurately the best timing for rewatering to assess the SR, especially when the symptoms of severe dehydration between genotypes are not contrasting enough to identify a difference visually.

Pseudomonas aeruginosa LasI-dependent plant growth promotion requires the host nitrate transceptor AtNRT1.1/CHL1 and the nitrate reductases NIA1 and NIA2.

MAIN CONCLUSION: In P. aeruginosa, mutation of the gene encoding N-acyl-L-homoserine lactone synthase LasI drives defense and plant growth promotion, and this latter trait requires adequate nitrate nutrition. Cross-kingdom communication with bacteria is crucial for plant growth and productivity. Here, we show a strong induction of genes for nitrate uptake and assimilation in Arabidopsis seedlings co-cultivated with P. aeruginosa WT (PAO1) or ΔlasI mutants defective on the synthesis of the quorum-sensing signaling molecule N-(3-oxododecanoyl)-L-homoserine lactone. Along with differential induction of defense-related genes, the change from plant growth repression to growth promotion upon bacterial QS disruption, correlated with upregulation of the dual-affinity nitrate transceptor CHL1/AtNRT1/NPF6.3 and the nitrate reductases NIA1 and NIA2. CHL1-GUS was induced in Arabidopsis primary root tips after transfer onto P. aeruginosa ΔlasI streaks at low and high N availability, whereas this bacterium required high concentrations of nitrogen to potentiate root and shoot biomass production and to improve root branching. Arabidopsis chl1-5 and chl1-12 mutants and double mutants in NIA1 and NIA2 nitrate reductases showed compromised growth under low nitrogen availability and failed to mount an effective growth promotion and root branching response even at high NH4NO3. WT P. aeruginosa PAO1 and P. aeruginosa ΔlasI mutant promoted the accumulation of nitric oxide (NO) in roots of both the WT and nia1nia2 double mutants, whereas NO donors SNP or SNAP did not improve growth or root branching in nia1nia2 double mutants with or without bacterial cocultivation. Thus, inoculation of Arabidopsis roots with P. aeruginosa drives gene expression for improved nitrogen acquisition and this macronutrient is critical for the plant growth-promoting effects upon disruption of the LasI quorum-sensing system.

Multi-omic analyses reveal the unique properties of chia (Salvia hispanica) seed metabolism.

Chia (Salvia hispanica) is an emerging crop considered a functional food containing important substances with multiple potential applications. However, the molecular basis of some relevant chia traits, such as seed mucilage and polyphenol content, remains to be discovered. This study generates an improved chromosome-level reference of the chia genome, resolving some highly repetitive regions, describing methylation patterns, and refining genome annotation. Transcriptomic analysis shows that seeds exhibit a unique expression pattern compared to other organs and tissues. Thus, a metabolic and proteomic approach is implemented to study seed composition and seed-produced mucilage. The chia genome exhibits a significant expansion in mucilage synthesis genes (compared to Arabidopsis), and gene network analysis reveals potential regulators controlling seed mucilage production. Rosmarinic acid, a compound with enormous therapeutic potential, was classified as the most abundant polyphenol in seeds, and candidate genes for its complex pathway are described. Overall, this study provides important insights into the molecular basis for the unique characteristics of chia seeds.

Development of High-Quality Nuclei Isolation to Study Plant Root-Microbe Interaction for Single-Nuclei Transcriptomic Sequencing in Soybean.

Single-nucleus RNA sequencing (sNucRNA-seq) is an emerging technology that has been rapidly adopted and demonstrated to be a powerful tool for detailed characterization of each cell- and sub cell-types in complex tissues of higher eukaryotes. sNucRNA-seq has also been used to dissect cell-type-specific transcriptional responses to environmental or developmental signals. In plants, this technology is being utilized to identify cell-type-specific trajectories for the study of several tissue types and important traits, including the single-cell dissection of the genetic determinants regulating plant-microbe interactions. The isolation of high-quality nuclei is one of the prerequisite steps to obtain high-quality sNucRNA-seq results. Although nuclei isolation from several plant tissues is well established, this process is highly troublesome when plant tissues are associated with beneficial or pathogenic microbes. For example, root tissues colonized with rhizobium bacteria (nodules), leaf tissue infected with bacterial or fungal pathogens, or roots infected with nematodes pose critical challenges to the isolation of high-quality nuclei and use for downstream application. Therefore, isolation of microbe-free, high-quality nuclei from plant tissues are necessary to avoid clogging or interference with the microfluidic channel (e.g., 10× Genomics) or particle-templated emulsion that are used in sNucRNA-seq platforms. Here, we developed a simple, effective, and efficient method to isolate high-quality nuclei from soybean roots and root nodules, followed by washing out bacterial contamination. This protocol has been designed to be easily implemented into any lab environment, and it can also be scaled up for use with multiple samples and applicable to a variety of samples with the presence of microbes. We validated this protocol by successfully generating a barcoded library using the 10× Genomics microfluidic platform from tissue subjected to this procedure. This workflow was developed to provide an accessible alternative to instrument-based approaches (e.g., fluorescent cell sorting) and will expand the ability of researchers to perform experiments such as sNucRNA-seq and sNucATAC-seq on inherently heterogeneous plant tissue samples.

Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis.

Lonicera macranthoides (LM) and L. japonica (LJ) are medicinal plants widely used in treating viral diseases, such as COVID-19. Although the two species are morphologically similar, their secondary metabolite profiles are significantly different. Here, metabolomics analysis showed that LM contained ~86.01 mg/g hederagenin-based saponins, 2000-fold higher than LJ. To gain molecular insights into its secondary metabolite production, a chromosome-level genome of LM was constructed, comprising 9 pseudo-chromosomes with 40 097 protein-encoding genes. Genome evolution analysis showed that LM and LJ were diverged 1.30-2.27 million years ago (MYA). The two plant species experienced a common whole-genome duplication event that occurred ∼53.9-55.2 MYA before speciation. Genes involved in hederagenin-based saponin biosynthesis were arranged in clusters on the chromosomes of LM and they were more highly expressed in LM than in LJ. Among them, oleanolic acid synthase (OAS) and UDP-glycosyltransferase 73 (UGT73) families were much more highly expressed in LM than in LJ. Specifically, LmOAS1 was identified to effectively catalyse the C-28 oxidation of ß-Amyrin to form oleanolic acid, the precursor of hederagenin-based saponin. LmUGT73P1 was identified to catalyse cauloside A to produce α-hederin. We further identified the key amino acid residues of LmOAS1 and LmUGT73P1 for their enzymatic activities. Additionally, comparing with collinear genes in LJ, LmOAS1 and LmUGT73P1 had an interesting phenomenon of 'neighbourhood replication' in LM genome. Collectively, the genomic resource and candidate genes reported here set the foundation to fully reveal the genome evolution of the Lonicera genus and hederagenin-based saponin biosynthetic pathway.

Bacterial community and root endodermis: a complementary relationship.

There are feedforward and feedback loops along the microbiota-root-shoot axis to maintain plant growth or defense under environmental stresses. Here, we highlight a reciprocal interaction between the endodermis and the plant-bacterial community, which stabilizes the diffusion barriers to maintain nutrient homeostasis under nutritional stress.

Back to primary endosymbiosis: from plastids to artificial photosynthetic life-forms.

Throughout evolution, only two known primary photosynthetic endosymbiosis occurred, which originated the Archaeplastida and the Paulinella spp. Fundamental questions regarding primary endosymbiosis remain unsolved, but may now be addressed with the recent development of chimeric photosynthetic life-form. Cournoyer et al. could establish artificial photosynthetic endosymbiosis between yeast and cyanobacteria.

Engineering microalgae for water phosphorus recovery to close the phosphorus cycle.

As a finite and non-renewable resource, phosphorus (P) is essential to all life and crucial for crop growth and food production. The boosted agricultural use and associated loss of P to the aquatic environment are increasing environmental pollution, harming ecosystems, and threatening future global food security. Thus, recovering and reusing P from water bodies is urgently needed to close the P cycle. As a natural, eco-friendly, and sustainable reclamation strategy, microalgae-based biological P recovery is considered a promising solution. However, the low P-accumulation capacity and P-removal efficiency of algal bioreactors restrict its application. Herein, it is demonstrated that manipulating genes involved in cellular P accumulation and signalling could triple the Chlamydomonas P-storage capacity to ~7% of dry biomass, which is the highest P concentration in plants to date. Furthermore, the engineered algae could recover P from wastewater almost three times faster than the unengineered one, which could be directly used as a P fertilizer. Thus, engineering genes involved in cellular P accumulation and signalling in microalgae could be a promising strategy to enhance P uptake and accumulation, which have the potential to accelerate the application of algae for P recovery from the water body and closing the P cycle.

Assembly and annotation of the Gossypium barbadense L. 'Pima-S6' genome raise questions about the chromosome structure and gene content of Gossypium barbadense genomes.

BACKGROUND: Gossypium barbadense L. Pima cotton is known for its resistance to Fusarium wilt and for producing fibers of superior quality highly prized in the textile market. We report a high-quality genome assembly and annotation of Pima-S6 cotton and its comparison at the chromosome and protein level to other ten Gossypium published genome assemblies. RESULTS: Synteny and orthogroup analyses revealed important differences on chromosome structure and annotated proteins content between our Pima-S6 and other publicly available G. barbadense assemblies, and across Gossypium assemblies in general. Detailed synteny analyses revealed chromosomal rearrangements between Pima-S6 and other Pima genomes on several chromosomes, with three major inversions in chromosomes A09, A13 and D05, raising questions about the true chromosome structure of Gossypium barbadense genomes. CONCLUSION: Analyses of the re-assembled and re-annotated genome of the close relative G. barbadense Pima 3-79 using our Pima-S6 assembly suggest that contig placement of some recent G. barbadense assemblies might have been unduly influenced by the use of the G. hirsutum TM-1 genome as the anchoring reference. The Pima-S6 reference genome provides a valuable genomic resource and offers new insights on genomic structure, and can serve as G. barbadense genome reference for future assemblies and further support FOV4-related studies and breeding efforts.

Strigolactones positively regulate Verticillium wilt resistance in cotton via crosstalk with other hormones.

Verticillium wilt caused by Verticillium dahliae is a serious vascular disease in cotton (Gossypium spp.). V. dahliae induces the expression of the CAROTENOID CLEAVAGE DIOXYGENASE 7 (GauCCD7) gene involved in strigolactone (SL) biosynthesis in Gossypium australe, suggesting a role for SLs in Verticillium wilt resistance. We found that the SL analog rac-GR24 enhanced while the SL biosynthesis inhibitor TIS108 decreased cotton resistance to Verticillium wilt. Knock-down of GbCCD7 and GbCCD8b genes in island cotton (Gossypium barbadense) decreased resistance, whereas overexpression of GbCCD8b in upland cotton (Gossypium hirsutum) increased resistance to Verticillium wilt. Additionally, Arabidopsis (Arabidopsis thaliana) SL mutants defective in CCD7 and CCD8 putative orthologs were susceptible, whereas both Arabidopsis GbCCD7- and GbCCD8b-overexpressing plants were more resistant to Verticillium wilt than wild-type (WT) plants. Transcriptome analyses showed that several genes related to the jasmonic acid (JA)- and abscisic acid (ABA)-signaling pathways, such as MYELOCYTOMATOSIS 2 (GbMYC2) and ABA-INSENSITIVE 5, respectively, were upregulated in the roots of WT cotton plants in responses to rac-GR24 and V. dahliae infection but downregulated in the roots of both GbCCD7- and GbCCD8b-silenced cotton plants. Furthermore, GbMYC2 suppressed the expression of GbCCD7 and GbCCD8b by binding to their promoters, which might regulate the homeostasis of SLs in cotton through a negative feedback loop. We also found that GbCCD7- and GbCCD8b-silenced cotton plants were impaired in V. dahliae-induced reactive oxygen species (ROS) accumulation. Taken together, our results suggest that SLs positively regulate cotton resistance to Verticillium wilt through crosstalk with the JA- and ABA-signaling pathways and by inducing ROS accumulation.

Nanocarrier spray: a nontransgenic approach for crop engineering.

Genetic modification allows engineering of important traits in crops through expensive and tedious procedures to alter their genetic background. Recently, Thagun et al. developed a nanocarrier-based foliar spray method to translocate bioactive molecules of interest into plant cells to engineer important traits without introducing a transgene.

Enhanced phenylpropanoid metabolism underlies resistance to Fusarium oxysporum f. sp. vasinfectum race 4 infection in the cotton cultivar Pima-S6 (Gossypium barbadense L.).

Introduction: Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) is a highly pathogenic soil-borne fungus responsible for Fusarium wilt in cotton (Gossypium spp.) and represents a continuing threat to cotton production in the southwest states of the United States, including California, New Mexico, and Texas. Pima (G. barbadense L.) cotton, which is highly valued for its fiber quality, has been shown to be more susceptible to this pathogen than Upland (G. hirsutum L.) cotton. Still, some Pima cultivars present resistance to FOV4 infection. Methods: To gain insights into the FOV4-resistance mechanism, we performed comparative transcriptional and metabolomic analyses between FOV4-susceptible and FOV4-resistant Pima cotton entries. FOV4-resistant Pima-S6 and FOV4-susceptible Pima S-7 and Pima 3-79 cotton plants were infected with FOV4 in the greenhouse, and the roots harvested 11 days post-infection for further analysis. Results: We found that an enhanced root phenylpropanoid metabolism in the resistant Pima-S6 cultivar determines FOV4-resistance. Gene-ontology enrichment of phenylpropanoid biosynthesis and metabolism categories correlated with the accumulation of secondary metabolites in Pima-S6 roots. Specifically, we found esculetin, a coumarin, an inhibitor of Fusarium's growth, accumulated in the roots of Pima-S6 even under non-infected conditions. Genes related to the phenylpropanoid biosynthesis and metabolism, including phenylalanine ammonia-lyase 2 (PAL2) and pleiotropic drug resistance 12 (PDR12) transporter, were found to be upregulated in Pima-S6 roots. Discussion: Our results highlight an essential role for the phenylpropanoid synthesis pathway in FOV4 resistance in Pima-S6 cotton. These genes represent attractive research prospects for FOV4-disease resistance and breeding approaches of other cotton cultivars of economic relevance.

Twisting development, the birth of a potential new gene.

Evolution has long been considered to be a conservative process in which new genes arise from pre-existing genes through gene duplication, domain shuffling, horizontal transfer, overprinting, retrotransposition, etc. However, this view is changing as new genes originating from non-genic sequences are discovered in different organisms. Still, rather limited functional information is available. Here, we have identified TWISTED1 (TWT1), a possible de novo-originated protein-coding gene that modifies microtubule arrangement and causes helicoidal growth in Arabidopsis thaliana when its expression is increased. Interestingly, even though TWT1 is a likely recent gene, the lack of TWT1 function affects A. thaliana development. TWT1 seems to have originated from a non-genic sequence. If so, it would be one of the few examples to date of how during evolution de novo genes are integrated into developmental cellular and organismal processes.

High-temperature adaptation of an OsNRT2.3 allele is thermoregulated by small RNAs.

Climate change negatively affects crop yield, which hinders efforts to reach agricultural sustainability and food security. Here, we show that a previously unidentified allele of the nitrate transporter gene OsNRT2.3 is required to maintain high yield and high nitrogen use efficiency under high temperatures. We demonstrate that this tolerance to high temperatures in rice accessions harboring the HTNE-2 (high temperature resistant and nitrogen efficient-2) alleles from enhanced translation of the OsNRT2.3b mRNA isoform and the decreased abundance of a unique small RNA (sNRT2.3-1) derived from the 5' untranslated region of OsNRT2.3. sNRT2.3-1 binds to the OsNRT2.3a mRNA in a temperature-dependent manner. Our findings reveal that allelic variation in the 5' untranslated region of OsNRT2.3 leads to an increase in OsNRT2.3b protein levels and higher yield during high-temperature stress. Our results also provide a breeding strategy to produce rice varieties with higher grain yield and lower N fertilizer input suitable for a sustainable agriculture that is resilient against climate change.

Titanium nanoparticles activate a transcriptional response in Arabidopsis that enhances tolerance to low phosphate, osmotic stress and pathogen infection.

Titanium is a ubiquitous element with a wide variety of beneficial effects in plants, including enhanced nutrient uptake and resistance to pathogens and abiotic stresses. While there is numerous evidence supporting the beneficial effects that Ti fertilization give to plants, there is little information on which genetic signaling pathways the Ti application activate in plant tissues. In this study, we utilize RNA-seq and ionomics technologies to unravel the molecular signals that Arabidopsis plants unleash when treated with Ti. RNA-seq analysis showed that Ti activates abscisic acid and salicylic acid signaling pathways and the expression of NUCLEOTIDE BINDING SITE-LEUCINE RICH REPEAT receptors likely by acting as a chemical priming molecule. This activation results in enhanced resistance to drought, high salinity, and infection with Botrytis cinerea in Arabidopsis. Ti also grants an enhanced nutritional state, even at suboptimal phosphate concentrations by upregulating the expression of multiple nutrient and membrane transporters and by modifying or increasing the production root exudates. Our results suggest that Ti might act similarly to the beneficial element Silicon in other plant species.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton.

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a "genome shock", leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton (Gossypium hirsutum and Gossypium barbadense, AADD) and its extant AA (Gossypium arboreum) and DD (Gossypium raimondii) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii. In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense, but absent from a resynthesized hybrid of G. arboreum and G. raimondii, suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis-regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Transcriptional analysis of Ceratopteris richardii young sporophyte reveals conservation of stem cell factors in the root apical meristem.

Gene expression in roots has been assessed in different plant species in studies ranging from complete organs to specific cell layers, and more recently at the single cell level. While certain genes or functional categories are expressed in the root of all or most plant species, lineage-specific genes have also been discovered. An increasing amount of transcriptomic data is available for angiosperms, while a limited amount of data is available for ferns, and few studies have focused on fern roots. Here, we present a de novo transcriptome assembly from three different parts of the Ceratopteris richardii young sporophyte. Differential gene expression analysis of the root tip transcriptional program showed an enrichment of functional categories related to histogenesis and cell division, indicating an active apical meristem. Analysis of a diverse set of orthologous genes revealed conserved expression in the root meristem, suggesting a preserved role for different developmental roles in this tissue, including stem cell maintenance. The reconstruction of evolutionary trajectories for ground tissue specification genes suggests a high degree of conservation in vascular plants, but not for genes involved in root cap development, showing that certain genes are absent in Ceratopteris or have intricate evolutionary paths difficult to track. Overall, our results suggest different processes of conservation and divergence of genes involved in root development.