System approaches to study root hairs as a single cell plant model: current status and future perspectives

Modulating root system architecture: cross-talk between auxin and phytohormones

Root architecture is an important agronomic trait that plays an essential role in water uptake, soil compactions, nutrient recycling, plant-microbe interactions, and hormone-mediated signaling pathways. Recently, significant advancements have been made in understanding how the complex interactions of phytohormones regulate the dynamic organization of root architecture in crops. Moreover, phytohormones, particularly auxin, act as internal regulators of root development in soil, starting from the early organogenesis to the formation of root hair (RH) through diverse signaling mechanisms. However, a considerable gap remains in understanding the hormonal cross-talk during various developmental stages of roots. This review examines the dynamic aspects of phytohormone signaling, cross-talk mechanisms, and the activation of transcription factors (TFs) throughout various developmental stages of the root life cycle. Understanding these developmental processes, together with hormonal signaling and molecular engineering in crops, can improve our knowledge of root development under various environmental conditions.

Strigolactones act downstream of gibberellins to regulate fiber cell elongation and cell wall thickness in cotton (Gossypium hirsutum)

Strigolactones (SLs) are a class of phytohormones that regulate plant shoot branching and adventitious root development. However, little is known regarding the role of SLs in controlling the behavior of the smallest unit of the organism, the single cell. Here, taking advantage of a classic single-cell model offered by the cotton (Gossypium hirsutum) fiber cell, we show that SLs, whose biosynthesis is fine-tuned by gibberellins (GAs), positively regulate cell elongation and cell wall thickness by promoting the biosynthesis of very long-chain fatty acids (VLCFAs) and cellulose, respectively. Furthermore, we identified two layers of transcription factors (TFs) involved in the hierarchical regulation of this GA-SL crosstalk. The top-layer TF GROWTH-REGULATING FACTOR 4 (GhGRF4) directly activates expression of the SL biosynthetic gene DWARF27 (D27) to increase SL accumulation in fiber cells and GAs induce GhGRF4 expression. SLs induce the expression of four second-layer TF genes (GhNAC100-2, GhBLH51, GhGT2, and GhB9SHZ1), which transmit SL signals downstream to two ketoacyl-CoA synthase genes (KCS) and three cellulose synthase (CesA) genes by directly activating their transcription. Finally, the KCS and CesA enzymes catalyze the biosynthesis of VLCFAs and cellulose, respectively, to regulate development of high-grade cotton fibers. In addition to providing a theoretical basis for cotton fiber improvement, our results shed light on SL signaling in plant development at the single-cell level.

High-Resolution Translatome Analysis Reveals Cortical Cell Programs During Early Soybean Nodulation

Nodule organogenesis in legumes is regulated temporally and spatially through gene networks. Genome-wide transcriptome, proteomic, and metabolomic analyses have been used previously to define the functional role of various plant genes in the nodulation process. However, while significant progress has been made, most of these studies have suffered from tissue dilution since only a few cells/root regions respond to rhizobial infection, with much of the root non-responsive. To partially overcome this issue, we adopted translating ribosome affinity purification (TRAP) to specifically monitor the response of the root cortex to rhizobial inoculation using a cortex-specific promoter. While previous studies have largely focused on the plant response within the root epidermis (e.g., root hairs) or within developing nodules, much less is known about the early responses within the root cortex, such as in relation to the development of the nodule primordium or growth of the infection thread. We focused on identifying genes specifically regulated during early nodule organogenesis using roots inoculated with Bradyrhizobium japonicum. A number of novel nodulation gene candidates were discovered, as well as soybean orthologs of nodulation genes previously reported in other legumes. The differential cortex expression of several genes was confirmed using a promoter-GUS analysis, and RNAi was used to investigate gene function. Notably, a number of differentially regulated genes involved in phytohormone signaling, including auxin, cytokinin, and gibberellic acid (GA), were also discovered, providing deep insight into phytohormone signaling during early nodule development.

A Sight on Single-Cell Transcriptomics in Plants Through the Prism of Cell-Based Computational Modeling Approaches: Benefits and Challenges for Data Analysis

Single-cell technology is a relatively new and promising way to obtain high-resolution transcriptomic data mostly used for animals during the last decade. However, several scientific groups developed and applied the protocols for some plant tissues. Together with deeply-developed cell-resolution imaging techniques, this achievement opens up new horizons for studying the complex mechanisms of plant tissue architecture formation. While the opportunities for integrating data from transcriptomic to morphogenetic levels in a unified system still present several difficulties, plant tissues have some additional peculiarities. One of the plants' features is that cell-to-cell communication topology through plasmodesmata forms during tissue growth and morphogenesis and results in mutual regulation of expression between neighboring cells affecting internal processes and cell domain development. Undoubtedly, we must take this fact into account when analyzing single-cell transcriptomic data. Cell-based computational modeling approaches successfully used in plant morphogenesis studies promise to be an efficient way to summarize such novel multiscale data. The inverse problem's solutions for these models computed on the real tissue templates can shed light on the restoration of individual cells' spatial localization in the initial plant organ-one of the most ambiguous and challenging stages in single-cell transcriptomic data analysis. This review summarizes new opportunities for advanced plant morphogenesis models, which become possible thanks to single-cell transcriptome data. Besides, we show the prospects of microscopy and cell-resolution imaging techniques to solve several spatial problems in single-cell transcriptomic data analysis and enhance the hybrid modeling framework opportunities.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

The role of microRNAs in the legume-Rhizobium nitrogen-fixing symbiosis

Under nitrogen starvation, most legume plants form a nitrogen-fixing symbiosis with Rhizobium bacteria. The bacteria induce the formation of a novel organ called the nodule in which rhizobia reside as intracellular symbionts and convert atmospheric nitrogen into ammonia. During this symbiosis, miRNAs are essential for coordinating the various plant processes required for nodule formation and function. miRNAs are non-coding, endogenous RNA molecules, typically 20-24 nucleotides long, that negatively regulate the expression of their target mRNAs. Some miRNAs can move systemically within plant tissues through the vascular system, which mediates, for example, communication between the stem/leaf tissues and the roots. In this review, we summarize the growing number of miRNAs that function during legume nodulation focusing on two model legumes, Lotus japonicus and Medicago truncatula, and two important legume crops, soybean (Glycine max) and common bean (Phaseolus vulgaris). This regulation impacts a variety of physiological processes including hormone signaling and spatial regulation of gene expression. The role of mobile miRNAs in regulating legume nodule number is also highlighted.

Lateral Root and Nodule Transcriptomes of Soybean

Symbiotic legume nodules and lateral roots arise away from the root meristem via dedifferentiation events. While these organs share some morphological and developmental similarities, whether legume nodules are modified lateral roots is an open question. We dissected emerging nodules, mature nodules, emerging lateral roots and young lateral roots, and constructed strand-specific RNA sequencing (RNAseq) libraries using polyA-enriched RNA preparations. Root sections above and below these organs, devoid of any lateral organs, were used to construct respective control tissue libraries. High sequence quality, predominant mapping to coding sequences, and consistency between replicates indicated that the RNAseq libraries were of a very high quality. We identified genes enriched in emerging nodules, mature nodules, emerging lateral roots and young lateral roots in soybean by comparing global gene expression profiles between each of these organs and adjacent root segments. Potential uses for this high quality transcriptome data set include generation of global gene regulatory networks to identify key regulators; metabolic pathway analyses and comparative analyses of key gene families to discover organ-specific biological processes; and identification of organ-specific alternate spliced transcripts. When combined with other similar datasets, especially from leguminous plants, these analyses can help answer questions on the evolutionary origins of root nodules and relationships between the development of different plant lateral organs.

Single-Cell Genomic Analysis in Plants

Individual cells in an organism are variable, which strongly impacts cellular processes. Advances in sequencing technologies have enabled single-cell genomic analysis to become widespread, addressing shortcomings of analyses conducted on populations of bulk cells. While the field of single-cell plant genomics is in its infancy, there is great potential to gain insights into cell lineage and functional cell types to help understand complex cellular interactions in plants. In this review, we discuss current approaches for single-cell plant genomic analysis, with a focus on single-cell isolation, DNA amplification, next-generation sequencing, and bioinformatics analysis. We outline the technical challenges of analysing material from a single plant cell, and then examine applications of single-cell genomics and the integration of this approach with genome editing. Finally, we indicate future directions we expect in the rapidly developing field of plant single-cell genomic analysis.

A comparative genomic and transcriptomic analysis at the level of isolated root hair cells reveals new conserved root hair regulatory elements

KEY MESSAGE

A comparative transcriptomic and genomic analysis between Arabidopsis thaliana and Glycine max root hair genes reveals the evolution of the expression of plant genes after speciation and whole genome duplication. Our understanding of the conservation and divergence of the expression patterns of genes between plant species is limited by the quality of the genomic and transcriptomic resources available. Specifically, the transcriptomes generated from plant organs are the reflection of the contribution of the different cell types composing the samples weighted by their relative abundances in the sample. These contributions can vary between plant species leading to the generation of datasets which are difficult to compare. To gain a deeper understanding of the evolution of gene transcription in and between plant species, we performed a comparative transcriptomic and genomic analysis at the level of one single plant cell type, the root hair cell, and between two model plants: Arabidopsis (Arabidopsis thaliana) and soybean (Glycine max). These two species, which diverged 90 million years ago, were selected as models based on the large amount of genomic and root hair transcriptomic information currently available. Our analysis revealed in detail the transcriptional divergence and conservation between soybean paralogs (i.e., the soybean genome is the product of two successive whole genome duplications) and between Arabidopsis and soybean orthologs in this single plant cell type. Taking advantage of this evolutionary study, we combined bioinformatics, molecular, cellular and microscopic tools to characterize plant promoter sequences and the discovery of two root hair regulatory elements (RHE1 and RHE2) consistently and specifically active in plant root hair cells.

Complexity and specificity of the maize (Zea mays L.) root hair transcriptome

Root hairs are tubular extensions of epidermis cells. Transcriptome profiling demonstrated that the single cell-type root hair transcriptome was less complex than the transcriptome of multiple cell-type primary roots without root hairs. In total, 831 genes were exclusively and 5585 genes were preferentially expressed in root hairs [false discovery rate (FDR) ≤1%]. Among those, the most significantly enriched Gene Ontology (GO) functional terms were related to energy metabolism, highlighting the high energy demand for the development and function of root hairs. Subsequently, the maize homologs for 138 Arabidopsis genes known to be involved in root hair development were identified and their phylogenetic relationship and expression in root hairs were determined. This study indicated that the genetic regulation of root hair development in Arabidopsis and maize is controlled by common genes, but also shows differences which need to be dissected in future genetic experiments. Finally, a maize root view of the eFP browser was implemented including the root hair transcriptome of the present study and several previously published maize root transcriptome data sets. The eFP browser provides color-coded expression levels for these root types and tissues for any gene of interest, thus providing a novel resource to study gene expression and function in maize roots.

Divergent cytosine DNA methylation patterns in single-cell, soybean root hairs

Chromatin modifications, such as cytosine methylation of DNA, play a significant role in mediating gene expression in plants, which affects growth, development, and cell differentiation. As root hairs are single-cell extensions of the root epidermis and the primary organs for water uptake and nutrients, we sought to use root hairs as a single-cell model system to measure the impact of environmental stress. We measured changes in cytosine DNA methylation in single-cell root hairs as compared with multicellular stripped roots, as well as in response to heat stress. Differentially methylated regions (DMRs) in each methylation context showed very distinct methylation patterns between cell types and in response to heat stress. Intriguingly, at normal temperature, root hairs were more hypermethylated than were stripped roots. However, in response to heat stress, both root hairs and stripped roots showed hypomethylation in each context, especially in the CHH context. Moreover, expression analysis of mRNA from similar tissues and treatments identified some associations between DMRs, genes and transposons. Taken together, the data indicate that changes in DNA methylation are directly or indirectly associated with expression of genes and transposons within the context of either specific tissues/cells or stress (heat).

Membrane glycerolipidome of soybean root hairs and its response to nitrogen and phosphate availability

Root hairs are tubular extensions of specific root epidermal cells important in plant nutrition and water absorption. To determine membrane glycerolipids in root hairs and roots may differ, as well as their respective response to nutrient availability, this study analyzed the membrane glycerolipid species in soybean root hairs and in roots stripped of root hairs, and their response to nitrogen (N) and phosphate (Pi) supplementation. The ratio of phospholipids to galactolipids was 1.5 fold higher in root hairs than in stripped roots. Under Pi deficiency, the ratio of phospholipids to galactolipids in stripped roots decreased with the greatest decrease found in the level of phosphatidylethanolamine (PE) in root hairs and stripped roots, and root hairs had an increased level of phosphatidic acid (PA). When seedlings were not supplied with N, the level of the N-containing lipids PE and phosphatidylserine in root hairs decreased whereas the level of non-N-containing lipids galactolipids and PA increased compared to N-supplied conditions. In stripped roots, the level of major membrane lipids was not different between N-sufficient and -deficient conditions. The results indicate that the membrane glycerolipidomes in root hairs are more responsive to nutrient availability than are the rest of roots.

Decipher the Molecular Response of Plant Single Cell Types to Environmental Stresses

The analysis of the molecular response of entire plants or organs to environmental stresses suffers from the cellular complexity of the samples used. Specifically, this cellular complexity masks cell-specific responses to environmental stresses and logically leads to the dilution of the molecular changes occurring in each cell type composing the tissue/organ/plant in response to the stress. Therefore, to generate a more accurate picture of these responses, scientists are focusing on plant single cell type approaches. Several cell types are now considered as models such as the pollen, the trichomes, the cotton fiber, various root cell types including the root hair cell, and the guard cell of stomata. Among them, several have been used to characterize plant response to abiotic and biotic stresses. In this review, we are describing the various -omic studies performed on these different plant single cell type models to better understand plant cell response to biotic and abiotic stresses.

Nod Factor Effects on Root Hair-Specific Transcriptome of Medicago truncatula: Focus on Plasma Membrane Transport Systems and Reactive Oxygen Species Networks

Root hairs are involved in water and nutrient uptake, and thereby in plant autotrophy. In legumes, they also play a crucial role in establishment of rhizobial symbiosis. To obtain a holistic view of Medicago truncatula genes expressed in root hairs and of their regulation during the first hours of the engagement in rhizobial symbiotic interaction, a high throughput RNA sequencing on isolated root hairs from roots challenged or not with lipochitooligosaccharides Nod factors (NF) for 4 or 20 h was carried out. This provided a repertoire of genes displaying expression in root hairs, responding or not to NF, and specific or not to legumes. In analyzing the transcriptome dataset, special attention was paid to pumps, transporters, or channels active at the plasma membrane, to other proteins likely to play a role in nutrient ion uptake, NF electrical and calcium signaling, control of the redox status or the dynamic reprogramming of root hair transcriptome induced by NF treatment, and to the identification of papilionoid legume-specific genes expressed in root hairs. About 10% of the root hair expressed genes were significantly up- or down-regulated by NF treatment, suggesting their involvement in remodeling plant functions to allow establishment of the symbiotic relationship. For instance, NF-induced changes in expression of genes encoding plasma membrane transport systems or disease response proteins indicate that root hairs reduce their involvement in nutrient ion absorption and adapt their immune system in order to engage in the symbiotic interaction. It also appears that the redox status of root hair cells is tuned in response to NF perception. In addition, 1176 genes that could be considered as "papilionoid legume-specific" were identified in the M. truncatula root hair transcriptome, from which 141 were found to possess an ortholog in every of the six legume genomes that we considered, suggesting their involvement in essential functions specific to legumes. This transcriptome provides a valuable resource to investigate root hair biology in legumes and the roles that these cells play in rhizobial symbiosis establishment. These results could also contribute to the long-term objective of transferring this symbiotic capacity to non-legume plants.

Soybean Roots Grown under Heat Stress Show Global Changes in Their Transcriptional and Proteomic Profiles

Heat stress is likely to be a key factor in the negative impact of climate change on crop production. Heat stress significantly influences the functions of roots, which provide support, water, and nutrients to other plant organs. Likewise, roots play an important role in the establishment of symbiotic associations with different microorganisms. Despite the physiological relevance of roots, few studies have examined their response to heat stress. In this study, we performed genome-wide transcriptomic and proteomic analyses on isolated root hairs, which are a single, epidermal cell type, and compared their response to stripped roots. On average, we identified 1849 and 3091 genes differentially regulated in root hairs and stripped roots, respectively, in response to heat stress. Our gene regulatory module analysis identified 10 key modules that might control the majority of the transcriptional response to heat stress. We also conducted proteomic analysis on membrane fractions isolated from root hairs and compared these responses to stripped roots. These experiments identified a variety of proteins whose expression changed within 3 h of application of heat stress. Most of these proteins were predicted to play a significant role in thermo-tolerance, as well as in chromatin remodeling and post-transcriptional regulation. The data presented represent an in-depth analysis of the heat stress response of a single cell type in soybean.