Intercellular protein movement: deciphering the language of development.

Development in multicellular organisms requires the coordinated production of a large number of specialized cell types through sophisticated signaling mechanisms. Non-cell-autonomous signals are one of the key mechanisms by which organisms coordinate development. In plants, intercellular movement of transcription factors and other mobile signals, such as hormones and peptides, is essential for normal development. Through a combination of different approaches, a large number of non-cell-autonomous signals that control plant development have been identified. We review some of the transcriptional regulators that traffic between cells, as well as how changes in symplasmic continuity affect and are affected by development. We also review current models for how mobile signals move via plasmodesmata and how movement is inhibited. Finally, we consider challenges in and new tools for studying protein movement.

A bistable circuit involving SCARECROW-RETINOBLASTOMA integrates cues to inform asymmetric stem cell division.

In plants, where cells cannot migrate, asymmetric cell divisions (ACDs) must be confined to the appropriate spatial context. We investigate tissue-generating asymmetric divisions in a stem cell daughter within the Arabidopsis root. Spatial restriction of these divisions requires physical binding of the stem cell regulator SCARECROW (SCR) by the RETINOBLASTOMA-RELATED (RBR) protein. In the stem cell niche, SCR activity is counteracted by phosphorylation of RBR through a cyclinD6;1-CDK complex. This cyclin is itself under transcriptional control of SCR and its partner SHORT ROOT (SHR), creating a robust bistable circuit with either high or low SHR-SCR complex activity. Auxin biases this circuit by promoting CYCD6;1 transcription. Mathematical modeling shows that ACDs are only switched on after integration of radial and longitudinal information, determined by SHR and auxin distribution, respectively. Coupling of cell-cycle progression to protein degradation resets the circuit, resulting in a "flip flop" that constrains asymmetric cell division to the stem cell region.

Nucleo-cytoplasmic Partitioning of ARF Proteins Controls Auxin Responses in Arabidopsis thaliana.

The phytohormone auxin plays crucial roles in nearly every aspect of plant growth and development. The auxin response factor (ARF) transcription factor family regulates auxin-responsive gene expression and exhibits nuclear localization in regions of high auxin responsiveness. Here we show that the ARF7 and ARF19 proteins accumulate in micron-sized assemblies within the cytoplasm of tissues with attenuated auxin responsiveness. We found that the intrinsically disordered middle region and the folded PB1 interaction domain of ARFs drive protein assembly formation. Mutation of a single lysine within the PB1 domain abrogates cytoplasmic assemblies, promotes ARF nuclear localization, and results in an altered transcriptome and morphological defects. Our data suggest a model in which ARF nucleo-cytoplasmic partitioning regulates auxin responsiveness, providing a mechanism for cellular competence for auxin signaling.

Gene regulatory networks for compatible versus incompatible grafts identify a role for SlWOX4 during junction formation.

Grafting has been adopted for a wide range of crops to enhance productivity and resilience; for example, grafting of Solanaceous crops couples disease-resistant rootstocks with scions that produce high-quality fruit. However, incompatibility severely limits the application of grafting and graft incompatibility remains poorly understood. In grafts, immediate incompatibility results in rapid death, but delayed incompatibility can take months or even years to manifest, creating a significant economic burden for perennial crop production. To gain insight into the genetic mechanisms underlying this phenomenon, we developed a model system using heterografting of tomato (Solanum lycopersicum) and pepper (Capsicum annuum). These grafted plants express signs of anatomical junction failure within the first week of grafting. By generating a detailed timeline for junction formation, we were able to pinpoint the cellular basis for this delayed incompatibility. Furthermore, we inferred gene regulatory networks for compatible self-grafts and incompatible heterografts based on these key anatomical events, which predict core regulators for grafting. Finally, we examined the role of vascular development in graft formation and uncovered SlWOX4 as a potential regulator of graft compatibility. Following this predicted regulator up with functional analysis, we show that Slwox4 homografts fail to form xylem bridges across the junction, demonstrating that indeed, SlWOX4 is essential for vascular reconnection during grafting, and may function as an early indicator of graft failure.

Mobile PEAR transcription factors integrate positional cues to prime cambial growth.

Apical growth in plants initiates upon seed germination, whereas radial growth is primed only during early ontogenesis in procambium cells and activated later by the vascular cambium1. Although it is not known how radial growth is organized and regulated in plants, this system resembles the developmental competence observed in some animal systems, in which pre-existing patterns of developmental potential are established early on2,3. Here we show that in Arabidopsis the initiation of radial growth occurs around early protophloem-sieve-element cell files of the root procambial tissue. In this domain, cytokinin signalling promotes the expression of a pair of mobile transcription factors-PHLOEM EARLY DOF 1 (PEAR1) and PHLOEM EARLY DOF 2 (PEAR2)-and their four homologues (DOF6, TMO6, OBP2 and HCA2), which we collectively name PEAR proteins. The PEAR proteins form a short-range concentration gradient that peaks at protophloem sieve elements, and activates gene expression that promotes radial growth. The expression and function of PEAR proteins are antagonized by the HD-ZIP III proteins, well-known polarity transcription factors4-the expression of which is concentrated in the more-internal domain of radially non-dividing procambial cells by the function of auxin, and mobile miR165 and miR166 microRNAs. The PEAR proteins locally promote transcription of their inhibitory HD-ZIP III genes, and thereby establish a negative-feedback loop that forms a robust boundary that demarks the zone of cell division. Taken together, our data establish that during root procambial development there exists a network in which a module that links PEAR and HD-ZIP III transcription factors integrates spatial information of the hormonal domains and miRNA gradients to provide adjacent zones of dividing and more-quiescent cells, which forms a foundation for further radial growth.

The ALOG family members OsG1L1 and OsG1L2 regulate inflorescence branching in rice.

The architecture of the rice inflorescence is an important determinant of crop yield. The length of the inflorescence and the number of branches are among the key factors determining the number of spikelets, and thus grains, that a plant will develop. In particular, the timing of the identity transition from indeterminate branch meristem to determinate spikelet meristem governs the complexity of the inflorescence. In this context, the ALOG gene TAWAWA1 (TAW1) has been shown to delay the transition to determinate spikelet development in Oryza sativa (rice). Recently, by combining precise laser microdissection of inflorescence meristems with RNA-seq, we observed that two ALOG genes, OsG1-like 1 (OsG1L1) and OsG1L2, have expression profiles similar to that of TAW1. Here, we report that osg1l1 and osg1l2 loss-of-function CRISPR mutants have similar phenotypes to the phenotype of the previously published taw1 mutant, suggesting that these genes might act on related pathways during inflorescence development. Transcriptome analysis of the osg1l2 mutant suggested interactions of OsG1L2 with other known inflorescence architecture regulators and the data sets were used for the construction of a gene regulatory network (GRN), proposing interactions among genes potentially involved in controlling inflorescence development in rice. In this GRN, we selected the homeodomain-leucine zipper transcription factor encoding the gene OsHOX14 for further characterization. The spatiotemporal expression profiling and phenotypical analysis of CRISPR loss-of-function mutants of OsHOX14 suggests that the proposed GRN indeed serves as a valuable resource for the identification of new proteins involved in rice inflorescence development.

POPEYE intercellular localization mediates cell-specific iron deficiency responses.

Plants must tightly regulate iron (Fe) sensing, acquisition, transport, mobilization, and storage to ensure sufficient levels of this essential micronutrient. POPEYE (PYE) is an iron responsive transcription factor that positively regulates the iron deficiency response, while also repressing genes essential for maintaining iron homeostasis. However, little is known about how PYE plays such contradictory roles. Under iron-deficient conditions, pPYE:GFP accumulates in the root pericycle while pPYE:PYE-GFP is localized to the nucleus in all Arabidopsis (Arabidopsis thaliana) root cells, suggesting that PYE may have cell-specific dynamics and functions. Using scanning fluorescence correlation spectroscopy and cell-specific promoters, we found that PYE-GFP moves between different cells and that the tendency for movement corresponds with transcript abundance. While localization to the cortex, endodermis, and vasculature is required to manage changes in iron availability, vasculature and endodermis localization of PYE-GFP protein exacerbated pye-1 defects and elicited a host of transcriptional changes that are detrimental to iron mobilization. Our findings indicate that PYE acts as a positive regulator of iron deficiency response by regulating iron bioavailability differentially across cells, which may trigger iron uptake from the surrounding rhizosphere and impact root energy metabolism.

Phosphate starvation: response mechanisms and solutions.

Phosphorus is essential to plant growth and agricultural crop yields, yet the challenges associated with phosphorus fertilization in agriculture, such as aquatic runoff pollution and poor phosphorus bioavailability, are increasingly difficult to manage. Comprehensively understanding the dynamics of phosphorus uptake and signaling mechanisms will inform the development of strategies to address these issues. This review describes regulatory mechanisms used by specific tissues in the root apical meristem to sense and take up phosphate from the rhizosphere. The major regulatory mechanisms and related hormone crosstalk underpinning phosphate starvation responses, cellular phosphate homeostasis, and plant adaptations to phosphate starvation are also discussed, along with an overview of the major mechanism of plant systemic phosphate starvation responses. Finally, this review discusses recent promising genetic engineering strategies for improving crop phosphorus use and computational approaches that may help further design strategies for improved plant phosphate acquisition. The mechanisms and approaches presented include a wide variety of species including not only Arabidopsis but also crop species such as Oryza sativa (rice), Glycine max (soybean), and Triticum aestivum (wheat) to address both general and species-specific mechanisms and strategies. The aspects of phosphorus deficiency responses and recently employed strategies of improving phosphate acquisition that are detailed in this review may provide insights into the mechanisms or phenotypes that may be targeted in efforts to improve crop phosphorus content and plant growth in low phosphorus soils.

BAM1/2 receptor kinase signaling drives CLE peptide-mediated formative cell divisions in Arabidopsis roots.

Cell division is often regulated by extracellular signaling networks to ensure correct patterning during development. In Arabidopsis, the SHORT-ROOT (SHR)/SCARECROW (SCR) transcription factor dimer activates CYCLIND6;1 (CYCD6;1) to drive formative divisions during root ground tissue development. Here, we show plasma-membrane-localized BARELY ANY MERISTEM1/2 (BAM1/2) family receptor kinases are required for SHR-dependent formative divisions and CYCD6;1 expression, but not SHR-dependent ground tissue specification. Root-enriched CLE ligands bind the BAM1 extracellular domain and are necessary and sufficient to activate SHR-mediated divisions and CYCD6;1 expression. Correspondingly, BAM-CLE signaling contributes to the restriction of formative divisions to the distal root region. Additionally, genetic analysis reveals that BAM-CLE and SHR converge to regulate additional cell divisions outside of the ground tissues. Our work identifies an extracellular signaling pathway regulating formative root divisions and provides a framework to explore this pathway in patterning and evolution.

Protein complex stoichiometry and expression dynamics of transcription factors modulate stem cell division.

Stem cells divide and differentiate to form all of the specialized cell types in a multicellular organism. In the Arabidopsis root, stem cells are maintained in an undifferentiated state by a less mitotically active population of cells called the quiescent center (QC). Determining how the QC regulates the surrounding stem cell initials, or what makes the QC fundamentally different from the actively dividing initials, is important for understanding how stem cell divisions are maintained. Here we gained insight into the differences between the QC and the cortex endodermis initials (CEI) by studying the mobile transcription factor SHORTROOT (SHR) and its binding partner SCARECROW (SCR). We constructed an ordinary differential equation model of SHR and SCR in the QC and CEI which incorporated the stoichiometry of the SHR-SCR complex as well as upstream transcriptional regulation of SHR and SCR. Our model prediction, coupled with experimental validation, showed that high levels of the SHR-SCR complex are associated with more CEI division but less QC division. Furthermore, our model prediction allowed us to propose the putative upstream SHR regulators SEUSS and WUSCHEL-RELATED HOMEOBOX 5 and to experimentally validate their roles in QC and CEI division. In addition, our model established the timing of QC and CEI division and suggests that SHR repression of QC division depends on formation of the SHR homodimer. Thus, our results support that SHR-SCR protein complex stoichiometry and regulation of SHR transcription modulate the division timing of two different specialized cell types in the root stem cell niche.

Precise transcriptional control of cellular quiescence by BRAVO/WOX5 complex in Arabidopsis roots.

Understanding stem cell regulatory circuits is the next challenge in plant biology, as these cells are essential for tissue growth and organ regeneration in response to stress. In the Arabidopsis primary root apex, stem cell-specific transcription factors BRAVO and WOX5 co-localize in the quiescent centre (QC) cells, where they commonly repress cell division so that these cells can act as a reservoir to replenish surrounding stem cells, yet their molecular connection remains unknown. Genetic and biochemical analysis indicates that BRAVO and WOX5 form a transcription factor complex that modulates gene expression in the QC cells to preserve overall root growth and architecture. Furthermore, by using mathematical modelling we establish that BRAVO uses the WOX5/BRAVO complex to promote WOX5 activity in the stem cells. Our results unveil the importance of transcriptional regulatory circuits in plant stem cell development.

tuxnet: a simple interface to process RNA sequencing data and infer gene regulatory networks.

Predicting gene regulatory networks (GRNs) from expression profiles is a common approach for identifying important biological regulators. Despite the increased use of inference methods, existing computational approaches often do not integrate RNA-sequencing data analysis, are not automated or are restricted to users with bioinformatics backgrounds. To address these limitations, we developed tuxnet, a user-friendly platform that can process raw RNA-sequencing data from any organism with an existing reference genome using a modified tuxedo pipeline (hisat 2 + cufflinks package) and infer GRNs from these processed data. tuxnet is implemented as a graphical user interface and can mine gene regulations, either by applying a dynamic Bayesian network (DBN) inference algorithm, genist, or a regression tree-based pipeline, rtp-star. We obtained time-course expression data of a PERIANTHIA (PAN) inducible line and inferred a GRN using genist to illustrate the use of tuxnet while gaining insight into the regulations downstream of the Arabidopsis root stem cell regulator PAN. Using rtp-star, we inferred the network of ATHB13, a downstream gene of PAN, for which we obtained wild-type and mutant expression profiles. Additionally, we generated two networks using temporal data from developmental leaf data and spatial data from root cell-type data to highlight the use of tuxnet to form new testable hypotheses from previously explored data. Our case studies feature the versatility of tuxnet when using different types of gene expression data to infer networks and its accessibility as a pipeline for non-bioinformaticians to analyze transcriptome data, predict causal regulations, assess network topology and identify key regulators.

GTL1 and DF1 regulate root hair growth through transcriptional repression of ROOT HAIR DEFECTIVE 6-LIKE 4 in Arabidopsis.

How plants determine the final size of growing cells is an important, yet unresolved, issue. Root hairs provide an excellent model system with which to study this as their final cell size is remarkably constant under constant environmental conditions. Previous studies have demonstrated that a basic helix-loop helix transcription factor ROOT HAIR DEFECTIVE 6-LIKE 4 (RSL4) promotes root hair growth, but how hair growth is terminated is not known. In this study, we demonstrate that a trihelix transcription factor GT-2-LIKE1 (GTL1) and its homolog DF1 repress root hair growth in Arabidopsis Our transcriptional data, combined with genome-wide chromatin-binding data, show that GTL1 and DF1 directly bind the RSL4 promoter and regulate its expression to repress root hair growth. Our data further show that GTL1 and RSL4 regulate each other, as well as a set of common downstream genes, many of which have previously been implicated in root hair growth. This study therefore uncovers a core regulatory module that fine-tunes the extent of root hair growth by the orchestrated actions of opposing transcription factors.

Characterizing the involvement of FaMADS9 in the regulation of strawberry fruit receptacle development.

FaMADS9 is the strawberry (Fragaria x ananassa) gene that exhibits the highest homology to the tomato (Solanum lycopersicum) RIN gene. Transgenic lines were obtained in which FaMADS9 was silenced. The fruits of these lines did not show differences in basic parameters, such as fruit firmness or colour, but exhibited lower Brix values in three of the four independent lines. The gene ontology MapMan category that was most enriched among the differentially expressed genes in the receptacles at the white stage corresponded to the regulation of transcription, including a high percentage of transcription factors and regulatory proteins associated with auxin action. In contrast, the most enriched categories at the red stage were transport, lipid metabolism and cell wall. Metabolomic analysis of the receptacles of the transformed fruits identified significant changes in the content of maltose, galactonic acid-1,4-lactone, proanthocyanidins and flavonols at the green/white stage, while isomaltose, anthocyanins and cuticular wax metabolism were the most affected at the red stage. Among the regulatory genes that were differentially expressed in the transgenic receptacles were several genes previously linked to flavonoid metabolism, such as MYB10, DIV, ZFN1, ZFN2, GT2, and GT5, or associated with the action of hormones, such as abscisic acid, SHP, ASR, GTE7 and SnRK2.7. The inference of a gene regulatory network, based on a dynamic Bayesian approach, among the genes differentially expressed in the transgenic receptacles at the white and red stages, identified the genes KAN1, DIV, ZFN2 and GTE7 as putative targets of FaMADS9. A MADS9-specific CArG box was identified in the promoters of these genes.

Framework for gradual progression of cell ontogeny in the Arabidopsis root meristem.

In plants, apical meristems allow continuous growth along the body axis. Within the root apical meristem, a group of slowly dividing quiescent center cells is thought to limit stem cell activity to directly neighboring cells, thus endowing them with unique properties, distinct from displaced daughters. This binary identity of the stem cells stands in apparent contradiction to the more gradual changes in cell division potential and differentiation that occur as cells move further away from the quiescent center. To address this paradox and to infer molecular organization of the root meristem, we used a whole-genome approach to determine dominant transcriptional patterns along root ontogeny zones. We found that the prevalent patterns are expressed in two opposing gradients. One is characterized by genes associated with development, the other enriched in differentiation genes. We confirmed these transcript gradients, and demonstrate that these translate to gradients in protein accumulation and gradual changes in cellular properties. We also show that gradients are genetically controlled through multiple pathways. Based on these findings, we propose that cells in the Arabidopsis root meristem gradually transition from stem cell activity toward differentiation.

Predicting gene regulatory networks by combining spatial and temporal gene expression data in Arabidopsis root stem cells.

Identifying the transcription factors (TFs) and associated networks involved in stem cell regulation is essential for understanding the initiation and growth of plant tissues and organs. Although many TFs have been shown to have a role in the Arabidopsis root stem cells, a comprehensive view of the transcriptional signature of the stem cells is lacking. In this work, we used spatial and temporal transcriptomic data to predict interactions among the genes involved in stem cell regulation. To accomplish this, we transcriptionally profiled several stem cell populations and developed a gene regulatory network inference algorithm that combines clustering with dynamic Bayesian network inference. We leveraged the topology of our networks to infer potential major regulators. Specifically, through mathematical modeling and experimental validation, we identified PERIANTHIA (PAN) as an important molecular regulator of quiescent center function. The results presented in this work show that our combination of molecular biology, computational biology, and mathematical modeling is an efficient approach to identify candidate factors that function in the stem cells.

Auxin minimum triggers the developmental switch from cell division to cell differentiation in the Arabidopsis root.

In multicellular organisms, a stringent control of the transition between cell division and differentiation is crucial for correct tissue and organ development. In the Arabidopsis root, the boundary between dividing and differentiating cells is positioned by the antagonistic interaction of the hormones auxin and cytokinin. Cytokinin affects polar auxin transport, but how this impacts the positional information required to establish this tissue boundary, is still unknown. By combining computational modeling with molecular genetics, we show that boundary formation is dependent on cytokinin's control on auxin polar transport and degradation. The regulation of both processes shapes the auxin profile in a well-defined auxin minimum. This auxin minimum positions the boundary between dividing and differentiating cells, acting as a trigger for this developmental transition, thus controlling meristem size.

Genetic Architecture and Molecular Networks Underlying Leaf Thickness in Desert-Adapted Tomato Solanum pennellii.

Thicker leaves allow plants to grow in water-limited conditions. However, our understanding of the genetic underpinnings of this highly functional leaf shape trait is poor. We used a custom-built confocal profilometer to directly measure leaf thickness in a set of introgression lines (ILs) derived from the desert tomato Solanum pennellii and identified quantitative trait loci. We report evidence of a complex genetic architecture of this trait and roles for both genetic and environmental factors. Several ILs with thick leaves have dramatically elongated palisade mesophyll cells and, in some cases, increased leaf ploidy. We characterized the thick IL2-5 and IL4-3 in detail and found increased mesophyll cell size and leaf ploidy levels, suggesting that endoreduplication underpins leaf thickness in tomato. Next, we queried the transcriptomes and inferred dynamic Bayesian networks of gene expression across early leaf ontogeny in these lines to compare the molecular networks that pattern leaf thickness. We show that thick ILs share S. pennellii-like expression profiles for putative regulators of cell shape and meristem determinacy as well as a general signature of cell cycle-related gene expression. However, our network data suggest that leaf thickness in these two lines is patterned at least partially by distinct mechanisms. Consistent with this hypothesis, double homozygote lines combining introgression segments from these two ILs show additive phenotypes, including thick leaves, higher ploidy levels, and larger palisade mesophyll cells. Collectively, these data establish a framework of genetic, anatomical, and molecular mechanisms that pattern leaf thickness in desert-adapted tomato.

Exposure to heavy metal stress triggers changes in plasmodesmatal permeability via deposition and breakdown of callose.

Both plants and animals must contend with changes in their environment. The ability to respond appropriately to these changes often underlies the ability of the individual to survive. In plants, an early response to environmental stress is an alteration in plasmodesmatal permeability with accompanying changes in cell to cell signaling. However, the ways in which plasmodesmata are modified, the molecular players involved in this regulation, and the biological significance of these responses are not well understood. Here, we examine the effects of nutrient scarcity and excess on plasmodesmata-mediated transport in the Arabidopsis thaliana root and identify two CALLOSE SYNTHASES and two ß-1,3-GLUCANASES as key regulators of these processes. Our results suggest that modification of plasmodesmata-mediated signaling underlies the ability of the plant to maintain root growth and properly partition nutrients when grown under conditions of excess nutrients.