Phenotyping and modeling of root hydraulic architecture reveal critical determinants of axial water transport.

Water uptake by roots is a key adaptation of plants to aerial life. Water uptake depends on root system architecture (RSA) and tissue hydraulic properties that, together, shape the root hydraulic architecture. This work investigates how the interplay between conductivities along radial (e.g. aquaporins) and axial (e.g. xylem vessels) pathways determines the water transport properties of highly branched RSAs as found in adult Arabidopsis (Arabidopsis thaliana) plants. A hydraulic model named HydroRoot was developed, based on multi-scale tree graph representations of RSAs. Root water flow was measured by the pressure chamber technique after successive cuts of a same root system from the tip toward the base. HydroRoot model inversion in corresponding RSAs allowed us to concomitantly determine radial and axial conductivities, providing evidence that the latter is often overestimated by classical evaluation based on the Hagen-Poiseuille law. Organizing principles of Arabidopsis primary and lateral root growth and branching were determined and used to apply the HydroRoot model to an extended set of simulated RSAs. Sensitivity analyses revealed that water transport can be co-limited by radial and axial conductances throughout the whole RSA. The number of roots that can be sectioned (intercepted) at a given distance from the base was defined as an accessible and informative indicator of RSA. The overall set of experimental and theoretical procedures was applied to plants mutated in ESKIMO1 and previously shown to have xylem collapse. This approach will be instrumental to dissect the root water transport phenotype of plants with intricate alterations in root growth or transport functions.

PUCHI represses early meristem formation in developing lateral roots of Arabidopsis thaliana.

Lateral root organogenesis is a key process in the development of a plant's root system and its adaptation to the environment. During lateral root formation, an early phase of cell proliferation first produces a four-cell-layered primordium, and only from this stage onwards is a root meristem-like structure, expressing root stem cell niche marker genes, being established in the developing organ. Previous studies reported that the gene regulatory network controlling lateral root formation is organized into two subnetworks whose mutual inhibition may contribute to organ patterning. PUCHI encodes an AP2/ERF transcription factor expressed early during lateral root primordium development and required for correct lateral root formation. To dissect the molecular events occurring during this early phase, we generated time-series transcriptomic datasets profiling lateral root development in puchi-1 mutants and wild types. Transcriptomic and reporter analyses revealed that meristem-related genes were expressed ectopically at early stages of lateral root formation in puchi-1 mutants. We conclude that, consistent with the inhibition of genetic modules contributing to lateral root development, PUCHI represses ectopic establishment of meristematic cell identities at early stages of organ development. These findings shed light on gene network properties that orchestrate correct timing and patterning during lateral root formation.

PUCHI regulates very long chain fatty acid biosynthesis during lateral root and callus formation.

Lateral root organogenesis plays an essential role in elaborating plant root system architecture. In Arabidopsis, the AP2 family transcription factor PUCHI controls cell proliferation in lateral root primordia. To identify potential targets of PUCHI, we analyzed a time course transcriptomic dataset of lateral root formation. We report that multiple genes coding for very long chain fatty acid (VLCFA) biosynthesis enzymes are induced during lateral root development in a PUCHI-dependent manner. Significantly, several mutants perturbed in VLCFA biosynthesis show similar lateral root developmental defects as puchi-1 Moreover, puchi-1 roots display the same disorganized callus formation phenotype as VLCFA biosynthesis-deficient mutants when grown on auxin-rich callus-inducing medium. Lipidomic profiling of puchi-1 roots revealed reduced VLCFA content compared with WT. We conclude that PUCHI-regulated VLCFA biosynthesis is part of a pathway controlling cell proliferation during lateral root and callus formation.

Inference of the gene regulatory network acting downstream of CROWN ROOTLESS 1 in rice reveals a regulatory cascade linking genes involved in auxin signaling, crown root initiation, and root meristem specification and maintenance.

Crown roots (CRs) are essential components of the rice root system. Several genes involved in CR initiation or development have been identified but our knowledge about how they organize to form a gene regulatory network (GRN) is still limited. To characterize the regulatory cascades acting during CR formation, we used a systems biology approach to infer the GRN controlling CR formation downstream of CROWN ROOTLESS 1 (CRL1), coding for an ASL (asymmetric leaves-2-like)/LBD (LOB domain) transcription factor necessary for CR initiation. A time-series transcriptomic dataset was generated after synchronized induction of CR formation by dexamethasone-mediated expression of CRL1 expression in a crl1 mutant background. This time series revealed three different genome expression phases during the early steps of CR formation and was further exploited to infer a GRN using a dedicated algorithm. The predicted GRN was confronted with experimental data and 72% of the inferred links were validated. Interestingly, this network revealed a regulatory cascade linking CRL1 to other genes involved in CR initiation, root meristem specification and maintenance, such as QUIESCENT-CENTER-SPECIFIC HOMEOBOX, and in auxin signalling. This predicted regulatory cascade was validated in vivo using transient activation assays. Thus, the CRL1-dependant GRN reflects major gene regulation events at play during CR formation and constitutes a valuable source of discovery to better understand this developmental process.

A New Phenotyping Pipeline Reveals Three Types of Lateral Roots and a Random Branching Pattern in Two Cereals.

Recent progress in root phenotyping has focused mainly on increasing throughput for genetic studies, while identifying root developmental patterns has been comparatively underexplored. We introduce a new phenotyping pipeline for producing high-quality spatiotemporal root system development data and identifying developmental patterns within these data. The SmartRoot image-analysis system and temporal and spatial statistical models were applied to two cereals, pearl millet (Pennisetum glaucum) and maize (Zea mays). Semi-Markov switching linear models were used to cluster lateral roots based on their growth rate profiles. These models revealed three types of lateral roots with similar characteristics in both species. The first type corresponds to fast and accelerating roots, the second to rapidly arrested roots, and the third to an intermediate type where roots cease elongation after a few days. These types of lateral roots were retrieved in different proportions in a maize mutant affected in auxin signaling, while the first most vigorous type was absent in maize plants exposed to severe shading. Moreover, the classification of growth rate profiles was mirrored by a ranking of anatomical traits in pearl millet. Potential dependencies in the succession of lateral root types along the primary root were then analyzed using variable-order Markov chains. The lateral root type was not influenced by the shootward neighbor root type or by the distance from this root. This random branching pattern of primary roots was remarkably conserved, despite the high variability of root systems in both species. Our phenotyping pipeline opens the door to exploring the genetic variability of lateral root developmental patterns.

Inference of the Arabidopsis lateral root gene regulatory network suggests a bifurcation mechanism that defines primordia flanking and central zones.

A large number of genes involved in lateral root (LR) organogenesis have been identified over the last decade using forward and reverse genetic approaches in Arabidopsis thaliana. Nevertheless, how these genes interact to form a LR regulatory network largely remains to be elucidated. In this study, we developed a time-delay correlation algorithm (TDCor) to infer the gene regulatory network (GRN) controlling LR primordium initiation and patterning in Arabidopsis from a time-series transcriptomic data set. The predicted network topology links the very early-activated genes involved in LR initiation to later expressed cell identity markers through a multistep genetic cascade exhibiting both positive and negative feedback loops. The predictions were tested for the key transcriptional regulator AUXIN RESPONSE FACTOR7 node, and over 70% of its targets were validated experimentally. Intriguingly, the predicted GRN revealed a mutual inhibition between the ARF7 and ARF5 modules that would control an early bifurcation between two cell fates. Analyses of the expression pattern of ARF7 and ARF5 targets suggest that this patterning mechanism controls flanking and central zone specification in Arabidopsis LR primordia.

Integration of hormonal signaling networks and mobile microRNAs is required for vascular patterning in Arabidopsis roots.

As multicellular organisms grow, positional information is continually needed to regulate the pattern in which cells are arranged. In the Arabidopsis root, most cell types are organized in a radially symmetric pattern; however, a symmetry-breaking event generates bisymmetric auxin and cytokinin signaling domains in the stele. Bidirectional cross-talk between the stele and the surrounding tissues involving a mobile transcription factor, SHORT ROOT (SHR), and mobile microRNA species also determines vascular pattern, but it is currently unclear how these signals integrate. We use a multicellular model to determine a minimal set of components necessary for maintaining a stable vascular pattern. Simulations perturbing the signaling network show that, in addition to the mutually inhibitory interaction between auxin and cytokinin, signaling through SHR, microRNA165/6, and PHABULOSA is required to maintain a stable bisymmetric pattern. We have verified this prediction by observing loss of bisymmetry in shr mutants. The model reveals the importance of several features of the network, namely the mutual degradation of microRNA165/6 and PHABULOSA and the existence of an additional negative regulator of cytokinin signaling. These components form a plausible mechanism capable of patterning vascular tissues in the absence of positional inputs provided by the transport of hormones from the shoot.

Inhibition of auxin signaling in Frankia species-infected cells in Casuarina glauca nodules leads to increased nodulation.

Actinorhizal symbioses are mutualistic interactions between plants and the soil bacteria Frankia spp. that lead to the formation of nitrogen-fixing root nodules. The plant hormone auxin has been suggested to play a role in the mechanisms that control the establishment of this symbiosis in the actinorhizal tree Casuarina glauca. Here, we analyzed the role of auxin signaling in Frankia spp.-infected cells. Using a dominant-negative version of an endogenous auxin-signaling regulator, INDOLE-3-ACETIC ACID7, we established that inhibition of auxin signaling in these cells led to increased nodulation and, as a consequence, to higher nitrogen fixation per plant even if nitrogen fixation per nodule mass was similar to that in the wild type. Our results suggest that auxin signaling in Frankia spp.-infected cells is involved in the long-distance regulation of nodulation in actinorhizal symbioses.

Lateral root morphogenesis is dependent on the mechanical properties of the overlaying tissues.

In Arabidopsis, lateral root primordia (LRPs) originate from pericycle cells located deep within the parental root and have to emerge through endodermal, cortical, and epidermal tissues. These overlaying tissues place biomechanical constraints on the LRPs that are likely to impact their morphogenesis. This study probes the interplay between the patterns of cell division, organ shape, and overlaying tissues on LRP morphogenesis by exploiting recent advances in live plant cell imaging and image analysis. Our 3D/4D image analysis revealed that early stage LRPs exhibit tangential divisions that create a ring of cells corralling a population of rapidly dividing cells at its center. The patterns of division in the latter population of cells during LRP morphogenesis are not stereotypical. In contrast, statistical analysis demonstrated that the shape of new LRPs is highly conserved. We tested the relative importance of cell division pattern versus overlaying tissues on LRP morphogenesis using mutant and transgenic approaches. The double mutant aurora1 (aur1) aur2 disrupts the pattern of LRP cell divisions and impacts its growth dynamics, yet the new organ's dome shape remains normal. In contrast, manipulating the properties of overlaying tissues disrupted LRP morphogenesis. We conclude that the interaction with overlaying tissues, rather than the precise pattern of divisions, is most important for LRP morphogenesis and optimizes the process of lateral root emergence.

Short-Root regulates primary, lateral, and adventitious root development in Arabidopsis.

Short-Root (SHR) is a well-characterized regulator of radial patterning and indeterminacy of the Arabidopsis (Arabidopsis thaliana) primary root. However, its role during the elaboration of root system architecture remains unclear. We report that the indeterminate wild-type Arabidopsis root system was transformed into a determinate root system in the shr mutant when growing in soil or agar. The root growth behavior of the shr mutant results from its primary root apical meristem failing to initiate cell division following germination. The inability of shr to reactivate mitotic activity in the root apical meristem is associated with the progressive reduction in the abundance of auxin efflux carriers, PIN-FORMED1 (PIN1), PIN2, PIN3, PIN4, and PIN7. The loss of primary root growth in shr is compensated by the activation of anchor root primordia, whose tissues are radially patterned like the wild type. However, SHR function is not restricted to the primary root but is also required for the initiation and patterning of lateral root primordia. In addition, SHR is necessary to maintain the indeterminate growth of lateral and anchor roots. We conclude that SHR regulates a wide array of Arabidopsis root-related developmental processes.

Plant Systems Biology: Lessons from Teaching.

Systems biology is the study of biological interactions. These interactions exist between biological entities at every scale, from genes to population, and create incredibly complex networks of feedbacks responsible for emerging behaviors. To study these behaviors, biologists can use models based on mathematical and computational formalisms grounded on vast existing corpus of theoretical work. This chapter develops an overview of this process of plant systems biology study from the point of view of a teaching course, and introduces the methods and studies presented in this second edition of the "Plant Systems Biology" book series.

Future Challenges in Plant Systems Biology.

Plant systems biology is currently facing several important challenges, whose nature depend on the considered frame of reference and associated scale. This review covers some of the issues associated respectively with the molecular, tissue, and whole-plant scales, as well as discusses the potential for latest advances in synthetic biology and machine-learning methods to be of use in the future of plant systems biology.

Plant Systems Biology: Further Reading and Resources.

This chapter compiles a list of useful references for aspiring plant systems biologists. It is structured in four sections focusing on systems biology books, reviews for the different types of models and resources at each plant scale, online resources, and plant systems biology communities.

How to Use the TDCor Algorithm to Infer Gene Regulatory Networks from Time Series Transcriptomic Data.

Over the last few decades, many genes have been functionally characterized and shown to be involved in various metabolic, developmental, and signaling pathways. However it still remains unclear how all these genes and pathways integrate into a unique regulatory network to coordinate the development and the growth, or the response to the environment. This is why unraveling the topology of gene regulatory networks (GRN) has become central to our understanding of all these processes. The recent advancement of high-throughput methods has provided enormous amount of -omics data. These data can now be exploited for rapid network reconstruction with statistical inference methods. We recently published a new GRN inference algorithm called TDCor which reconstructs GRN from time-series transcriptomic data. The algorithm has been released in the form of an R package. Here, I describe into details how to install and use the package.

Postembryonic Organogenesis in Plants: Experimental Induction of New Shoot and Root Organs.

Postembryonic organogenesis is a critical component in plant root and shoot development and its adaptation to the environment. Decades of scientific analyses have yielded a wealth of experimental data about the cellular and molecular processes orchestrating the postembryonic formation of new shoot and root organs. Among these, distribution and signaling of the plant hormone auxin play a prominent role. Systems biology approaches are now particularly interesting to study the emerging properties of such complex and dynamic regulatory networks. To fully explore the precise kinetics of these organogenesis processes, efficient protocols for the synchronized induction of shoot and root organogenesis are extremely valuable. Two protocols for shoot and root organ induction are detailed.

Auxin carriers localization drives auxin accumulation in plant cells infected by Frankia in Casuarina glauca actinorhizal nodules.

Actinorhizal symbioses are mutualistic interactions between plants and the soil bacteria Frankia that lead to the formation of nitrogen-fixing root nodules. Little is known about the signaling mechanisms controlling the different steps of the establishment of the symbiosis. The plant hormone auxin has been suggested to play a role. Here we report that auxin accumulates within Frankia-infected cells in actinorhizal nodules of Casuarina glauca. Using a combination of computational modeling and experimental approaches, we establish that this localized auxin accumulation is driven by the cell-specific expression of auxin transporters and by Frankia auxin biosynthesis in planta. Our results indicate that the plant actively restricts auxin accumulation to Frankia-infected cells during the symbiotic interaction.

Plant systems biology: network matters.

Systems biology is all about networks. A recent trend has been to associate systems biology exclusively with the study of gene regulatory or protein-interaction networks. However, systems biology approaches can be applied at many other scales, from the subatomic to the ecosystem scales. In this review, we describe studies at the sub-cellular, tissue, whole plant and crop scales and highlight how these studies can be related to systems biology. We discuss the properties of system approaches at each scale as well as their current limits, and pinpoint in each case advances unique to the considered scale but representing potential for the other scales. We conclude by examining plant models bridging different scales and considering the future prospects of plant systems biology.

Transcriptome profiling of laser-captured crown root primordia reveals new pathways activated during early stages of crown root formation in rice.

Crown roots constitute the main part of the rice root system. Several key genes involved in crown root initiation and development have been identified by functional genomics approaches. Nevertheless, these approaches are impaired by functional redundancy and mutant lethality. To overcome these limitations, organ targeted transcriptome analysis can help to identify genes involved in crown root formation and early development. In this study, we generated an atlas of genes expressed in developing crown root primordia in comparison with adjacent stem cortical tissue at three different developmental stages before emergence, using laser capture microdissection. We identified 3975 genes differentially expressed in crown root primordia. About 30% of them were expressed at the three developmental stages, whereas 10.5%, 19.5% and 12.8% were specifically expressed at the early, intermediate and late stages, respectively. Sorting them by functional ontology highlighted an active transcriptional switch during the process of crown root primordia formation. Cross-analysis with other rice root development-related datasets revealed genes encoding transcription factors, chromatin remodeling factors, peptide growth factors, and cell wall remodeling enzymes that are likely to play a key role during crown root primordia formation. This atlas constitutes an open primary data resource for further studies on the regulation of crown root initiation and development.

Flux-based transport enhancement as a plausible unifying mechanism for auxin transport in meristem development.

Plants continuously generate new organs through the activity of populations of stem cells called meristems. The shoot apical meristem initiates leaves, flowers, and lateral meristems in highly ordered, spiralled, or whorled patterns via a process called phyllotaxis. It is commonly accepted that the active transport of the plant hormone auxin plays a major role in this process. Current hypotheses propose that cellular hormone transporters of the PIN family would create local auxin maxima at precise positions, which in turn would lead to organ initiation. To explain how auxin transporters could create hormone fluxes to distinct regions within the plant, different concepts have been proposed. A major hypothesis, canalization, proposes that the auxin transporters act by amplifying and stabilizing existing fluxes, which could be initiated, for example, by local diffusion. This convincingly explains the organised auxin fluxes during vein formation, but for the shoot apical meristem a second hypothesis was proposed, where the hormone would be systematically transported towards the areas with the highest concentrations. This implies the coexistence of two radically different mechanisms for PIN allocation in the membrane, one based on flux sensing and the other on local concentration sensing. Because these patterning processes require the interaction of hundreds of cells, it is impossible to estimate on a purely intuitive basis if a particular scenario is plausible or not. Therefore, computational modelling provides a powerful means to test this type of complex hypothesis. Here, using a dedicated computer simulation tool, we show that a flux-based polarization hypothesis is able to explain auxin transport at the shoot meristem as well, thus providing a unifying concept for the control of auxin distribution in the plant. Further experiments are now required to distinguish between flux-based polarization and other hypotheses.

Development of a model estimating root length density from root impacts on a soil profile in pearl millet (Pennisetum glaucum (L.) R. Br). Application to measure root system response to water stress in field conditions.

Pearl millet is able to withstand dry and hot conditions and plays an important role for food security in arid and semi-arid areas of Africa and India. However, low soil fertility and drought constrain pearl millet yield. One target to address these constraints through agricultural practices or breeding is root system architecture. In this study, in order to easily phenotype the root system in field conditions, we developed a model to predict root length density (RLD) of pearl millet plants from root intersection densities (RID) counted on a trench profile in field conditions. We identified root orientation as an important parameter to improve the relationship between RID and RLD. Root orientation was notably found to depend on soil depth and to differ between thick roots (more anisotropic with depth) and fine roots (isotropic at all depths). We used our model to study pearl millet root system response to drought and showed that pearl millet reorients its root growth toward deeper soil layers that retain more water in these conditions. Overall, this model opens ways for the characterization of the impact of environmental factors and management practices on pearl millet root system development.