An Integrative Systems Perspective on Plant Phosphate Research.

The case for improving crop phosphorus-use-efficiency is widely recognized. Although much is known about the molecular and regulatory mechanisms, improvements have been hampered by the extreme complexity of phosphorus (P) dynamics, which involves soil chemistry; plant-soil interactions; uptake, transport, utilization and remobilization within plants; and agricultural practices. The urgency and direction of phosphate research is also dependent upon the finite sources of P, availability of stocks to farmers and reducing environmental hazards. This work introduces integrative systems approaches as a way to represent and understand this complexity, so that meaningful links can be established between genotype, environment, crop traits and yield. It aims to provide a large set of pointers to potential genes and research practice, with a view to encouraging members of the plant-phosphate research community to adopt such approaches so that, together, we can aid efforts in global food security.

Regulatory feedback response mechanisms to phosphate starvation in rice.

Phosphorus is a growth-limiting nutrient for plants. The growing scarcity of phosphate stocks threatens global food security. Phosphate-uptake regulation is so complex and incompletely known that attempts to improve phosphorus use efficiency have had extremely limited success. This study improves our understanding of the molecular mechanisms underlying phosphate uptake by investigating the transcriptional dynamics of two regulators: the Ubiquitin ligase PHO2 and the long non-coding RNA IPS1. Temporal measurements of RNA levels have been integrated into mechanistic mathematical models using advanced statistical techniques. Models based solely on current knowledge could not adequately explain the temporal expression profiles. Further modeling and bioinformatics analysis have led to the prediction of three regulatory features: the PHO2 protein mediates the degradation of its own transcriptional activator to maintain constant PHO2 mRNA levels; the binding affinity of the transcriptional activator of PHO2 is impaired by a phosphate-sensitive transcriptional repressor/inhibitor; and the extremely high levels of IPS1 and its rapid disappearance upon Pi re-supply are best explained by Pi-sensitive RNA protection. This work offers both new opportunities for plant phosphate research that will be essential for informing the development of phosphate efficient crop varieties, and a foundation for the development of models integrating phosphate with other stress responses.

A DEMETER-like DNA demethylase governs tomato fruit ripening.

In plants, genomic DNA methylation which contributes to development and stress responses can be actively removed by DEMETER-like DNA demethylases (DMLs). Indeed, in Arabidopsis DMLs are important for maternal imprinting and endosperm demethylation, but only a few studies demonstrate the developmental roles of active DNA demethylation conclusively in this plant. Here, we show a direct cause and effect relationship between active DNA demethylation mainly mediated by the tomato DML, SlDML2, and fruit ripening- an important developmental process unique to plants. RNAi SlDML2 knockdown results in ripening inhibition via hypermethylation and repression of the expression of genes encoding ripening transcription factors and rate-limiting enzymes of key biochemical processes such as carotenoid synthesis. Our data demonstrate that active DNA demethylation is central to the control of ripening in tomato.

The circadian clock rephases during lateral root organ initiation in Arabidopsis thaliana.

The endogenous circadian clock enables organisms to adapt their growth and development to environmental changes. Here we describe how the circadian clock is employed to coordinate responses to the key signal auxin during lateral root (LR) emergence. In the model plant, Arabidopsis thaliana, LRs originate from a group of stem cells deep within the root, necessitating that new organs emerge through overlying root tissues. We report that the circadian clock is rephased during LR development. Metabolite and transcript profiling revealed that the circadian clock controls the levels of auxin and auxin-related genes including the auxin response repressor IAA14 and auxin oxidase AtDAO2. Plants lacking or overexpressing core clock components exhibit LR emergence defects. We conclude that the circadian clock acts to gate auxin signalling during LR development to facilitate organ emergence.

Multi-omics analysis identifies genes mediating the extension of cell walls in the Arabidopsis thaliana root elongation zone.

Plant cell wall composition is important for regulating growth rates, especially in roots. However, neither analyses of cell wall composition nor transcriptomes on their own can comprehensively reveal which genes and processes are mediating growth and cell elongation rates. This study reveals the benefits of carrying out multiple analyses in combination. Sections of roots from five anatomically and functionally defined zones in Arabidopsis thaliana were prepared and divided into three biological replicates. We used glycan microarrays and antibodies to identify the major classes of glycans and glycoproteins present in the cell walls of these sections, and identified the expected decrease in pectin and increase in xylan from the meristematic zone (MS), through the rapid and late elongation zones (REZ, LEZ) to the maturation zone and the rest of the root, including the emerging lateral roots. Other compositional changes included extensin and xyloglucan levels peaking in the REZ and increasing levels of arabinogalactan-proteins (AGP) epitopes from the MS to the LEZ, which remained high through the subsequent mature zones. Immuno-staining using the same antibodies identified the tissue and (sub)cellular localization of many epitopes. Extensins were localized in epidermal and cortex cell walls, while AGP glycans were specific to different tissues from root-hair cells to the stele. The transcriptome analysis found several gene families peaking in the REZ. These included a large family of peroxidases (which produce the reactive oxygen species (ROS) needed for cell expansion), and three xyloglucan endo-transglycosylase/hydrolase genes (XTH17, XTH18, and XTH19). The significance of the latter may be related to a role in breaking and re-joining xyloglucan cross-bridges between cellulose microfibrils, a process which is required for wall expansion. Knockdowns of these XTHs resulted in shorter root lengths, confirming a role of the corresponding proteins in root extension growth.

Mechanical modelling quantifies the functional importance of outer tissue layers during root elongation and bending.

Root elongation and bending require the coordinated expansion of multiple cells of different types. These processes are regulated by the action of hormones that can target distinct cell layers. We use a mathematical model to characterise the influence of the biomechanical properties of individual cell walls on the properties of the whole tissue. Taking a simple constitutive model at the cell scale which characterises cell walls via yield and extensibility parameters, we derive the analogous tissue-level model to describe elongation and bending. To accurately parameterise the model, we take detailed measurements of cell turgor, cell geometries and wall thicknesses. The model demonstrates how cell properties and shapes contribute to tissue-level extensibility and yield. Exploiting the highly organised structure of the elongation zone (EZ) of the Arabidopsis root, we quantify the contributions of different cell layers, using the measured parameters. We show how distributions of material and geometric properties across the root cross-section contribute to the generation of curvature, and relate the angle of a gravitropic bend to the magnitude and duration of asymmetric wall softening. We quantify the geometric factors which lead to the predominant contribution of the outer cell files in driving root elongation and bending.

Systems analysis of auxin transport in the Arabidopsis root apex.

Auxin is a key regulator of plant growth and development. Within the root tip, auxin distribution plays a crucial role specifying developmental zones and coordinating tropic responses. Determining how the organ-scale auxin pattern is regulated at the cellular scale is essential to understanding how these processes are controlled. In this study, we developed an auxin transport model based on actual root cell geometries and carrier subcellular localizations. We tested model predictions using the DII-VENUS auxin sensor in conjunction with state-of-the-art segmentation tools. Our study revealed that auxin efflux carriers alone cannot create the pattern of auxin distribution at the root tip and that AUX1/LAX influx carriers are also required. We observed that AUX1 in lateral root cap (LRC) and elongating epidermal cells greatly enhance auxin's shootward flux, with this flux being predominantly through the LRC, entering the epidermal cells only as they enter the elongation zone. We conclude that the nonpolar AUX1/LAX influx carriers control which tissues have high auxin levels, whereas the polar PIN carriers control the direction of auxin transport within these tissues.

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.

Growth-induced hormone dilution can explain the dynamics of plant root cell elongation.

In the elongation zone of the Arabidopsis thaliana plant root, cells undergo rapid elongation, increasing their length by ∼10-fold over 5 h while maintaining a constant radius. Although progress is being made in understanding how this growth is regulated, little consideration has been given as to how cell elongation affects the distribution of the key regulating hormones. Using a multiscale mathematical model and measurements of growth dynamics, we investigate the distribution of the hormone gibberellin in the root elongation zone. The model quantifies how rapid cell expansion causes gibberellin to dilute, creating a significant gradient in gibberellin levels. By incorporating the gibberellin signaling network, we simulate how gibberellin dilution affects the downstream components, including the growth-repressing DELLA proteins. We predict a gradient in DELLA that provides an explanation of the reduction in growth exhibited as cells move toward the end of the elongation zone. These results are validated at the molecular level by comparing predicted mRNA levels with transcriptomic data. To explore the dynamics further, we simulate perturbed systems in which gibberellin levels are reduced, considering both genetically modified and chemically treated roots. By modeling these cases, we predict how these perturbations affect gibberellin and DELLA levels and thereby provide insight into their altered growth dynamics.

Root gravitropism is regulated by a transient lateral auxin gradient controlled by a tipping-point mechanism.

Gravity profoundly influences plant growth and development. Plants respond to changes in orientation by using gravitropic responses to modify their growth. Cholodny and Went hypothesized over 80 years ago that plants bend in response to a gravity stimulus by generating a lateral gradient of a growth regulator at an organ's apex, later found to be auxin. Auxin regulates root growth by targeting Aux/IAA repressor proteins for degradation. We used an Aux/IAA-based reporter, domain II (DII)-VENUS, in conjunction with a mathematical model to quantify auxin redistribution following a gravity stimulus. Our multidisciplinary approach revealed that auxin is rapidly redistributed to the lower side of the root within minutes of a 90° gravity stimulus. Unexpectedly, auxin asymmetry was rapidly lost as bending root tips reached an angle of 40° to the horizontal. We hypothesize roots use a "tipping point" mechanism that operates to reverse the asymmetric auxin flow at the midpoint of root bending. These mechanistic insights illustrate the scientific value of developing quantitative reporters such as DII-VENUS in conjunction with parameterized mathematical models to provide high-resolution kinetics of hormone redistribution.

Statistical evaluation of transcriptomic data generated using the Affymetrix one-cycle, two-cycle and IVT-Express RNA labelling protocols with the Arabidopsis ATH1 microarray.

BACKGROUND: Microarrays are a powerful tool used for the determination of global RNA expression. There is an increasing requirement to focus on profiling gene expression in tissues where it is difficult to obtain large quantities of material, for example individual tissues within organs such as the root, or individual isolated cells. From such samples, it is difficult to produce the amount of RNA required for labelling and hybridisation in microarray experiments, thus a process of amplification is usually adopted. Despite the increasing use of two-cycle amplification for transcriptomic analyses on the Affymetrix ATH1 array, there has been no report investigating any potential bias in gene representation that may occur as a result. RESULTS: Here we compare transcriptomic data generated using Affymetrix one-cycle (standard labelling protocol), two-cycle (small-sample protocol) and IVT-Express protocols with the Affymetrix ATH1 array using Arabidopsis root samples. Results obtained with each protocol are broadly similar. However, we show that there are 35 probe sets (of a total of 22810) that are misrepresented in the two-cycle data sets. Of these, 33 probe sets were classed as mis-amplified when comparisons of two independent publicly available data sets were undertaken. CONCLUSIONS: Given the unreliable nature of the highlighted probes, we caution against using data associated with the corresponding genes in analyses involving transcriptomic data generated with two-cycle amplification protocols. We have shown that the Affymetrix IVT-E labelling protocol produces data with less associated bias than the two-cycle protocol, and as such, would recommend this kit for new experiments that involve small samples.

Inference of gene regulatory networks using boolean-network inference methods.

The modeling of genetic networks especially from microarray and related data has become an important aspect of the biosciences. This review takes a fresh look at a specific family of models used for constructing genetic networks, the so-called Boolean networks. The review outlines the various different types of Boolean network developed to date, from the original Random Boolean Network to the current Probabilistic Boolean Network. In addition, some of the different inference methods available to infer these genetic networks are also examined. Where possible, particular attention is paid to input requirements as well as the efficiency, advantages and drawbacks of each method. Though the Boolean network model is one of many models available for network inference today, it is well established and remains a topic of considerable interest in the field of genetic network inference. Hybrids of Boolean networks with other approaches may well be the way forward in inferring the most informative networks.