Search

Background: Many mathematical models have now been employed across every area of systems biology. These models increasingly involve large numbers of unknown parameters, have complex structure which can result in substantial evaluation time relative to the needs of the analysis, and need to be compared to observed data of various forms. The correct analysis of such models usually requires a global parameter search, over a high dimensional parameter space, that incorporates and respects the most important sources of uncertainty. This can be an extremely difficult task, but it is essential for any meaningful inference or prediction to be made about any biological system. It hence represents a fundamental challenge for the whole of systems biology., Methods: Bayesian statistical methodology for the uncertainty analysis of complex models is introduced, which is designed to address the high dimensional global parameter search problem. Bayesian emulators that mimic the systems biology model but which are extremely fast to evaluate are embeded within an iterative history match: an efficient method to search high dimensional spaces within a more formal statistical setting, while incorporating major sources of uncertainty., Results: The approach is demonstrated via application to a model of hormonal crosstalk in Arabidopsis root development, which has 32 rate parameters, for which we identify the sets of rate parameter values that lead to acceptable matches between model output and observed trend data. The multiple insights into the model's structure that this analysis provides are discussed. The methodology is applied to a second related model, and the biological consequences of the resulting comparison, including the evaluation of gene functions, are described., Conclusions: Bayesian uncertainty analysis for complex models using both emulators and history matching is shown to be a powerful technique that can greatly aid the study of a large class of systems biology models. It both provides insight into model behaviour and identifies the sets of rate parameters of interest.

Root-knot nematodes (RKNs) induce giant cells (GCs) from root vascular cells inside the galls. Accompanying molecular changes as a function of infection time and across different species, and their functional impact, are still poorly understood. Thus, the transcriptomes of tomato galls and laser capture microdissected (LCM) GCs over the course of parasitism were compared with those of Arabidopsis, and functional analysis of a repressed gene was performed. Microarray hybridization with RNA from galls and LCM GCs, infection-reproduction tests and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) transcriptional profiles in susceptible and resistant (Mi-1) lines were performed in tomato. Tomato GC-induced genes include some possibly contributing to the epigenetic control of GC identity. GC-repressed genes are conserved between tomato and Arabidopsis, notably those involved in lignin deposition. However, genes related to the regulation of gene expression diverge, suggesting that diverse transcriptional regulators mediate common responses leading to GC formation in different plant species. TPX1, a cell wall peroxidase specifically involved in lignification, was strongly repressed in GCs/galls, but induced in a nearly isogenic Mi-1 resistant line on nematode infection. TPX1 overexpression in susceptible plants hindered nematode reproduction and GC expansion. Time-course and cross-species comparisons of gall and GC transcriptomes provide novel insights pointing to the relevance of gene repression during RKN establishment., (© 2013 The Authors New Phytologist © 2013 New Phytologist Trust.)

Background: The control of vascular tissue development in plants is influenced by diverse hormonal signals, but their interactions during this process are not well understood. Wild-type sterol profiles are essential for growth, tissue patterning and signalling processes in plant development, and are required for regulated vascular patterning., Methodology/principal Findings: Here we investigate the roles of sterols in vascular tissue development, through an analysis of the Arabidopsis mutants hydra1 and fackel/hydra2, which are defective in the enzymes sterol isomerase and sterol C-14 reductase respectively. We show that defective vascular patterning in the shoot is associated with ectopic cell divisions. Expression of the auxin-regulated AtHB8 homeobox gene is disrupted in mutant embryos and seedlings, associated with variably incomplete vascular strand formation and duplication of the longitudinal axis. Misexpression of the auxin reporter proIAA2ratioGUS and mislocalization of PIN proteins occurs in the mutants. Introduction of the ethylene-insensitive ein2 mutation partially rescues defective cell division, localization of PIN proteins, and vascular strand development., Conclusions: The results support a model in which sterols are required for correct auxin and ethylene crosstalk to regulate PIN localization, auxin distribution and AtHB8 expression, necessary for correct vascular development.

An important question in plant biology is how genes influence the crosstalk between hormones to regulate growth. In this study, we model POLARIS (PLS) gene function and crosstalk between auxin, ethylene and cytokinin in Arabidopsis. Experimental evidence suggests that PLS acts on or close to the ethylene receptor ETR1, and a mathematical model describing possible PLS-ethylene pathway interactions is developed, and used to make quantitative predictions about PLS-hormone interactions. Modelling correctly predicts experimental results for the effect of the pls gene mutation on endogenous cytokinin concentration. Modelling also reveals a role for PLS in auxin biosynthesis in addition to a role in auxin transport. The model reproduces available mutants, and with new experimental data provides new insights into how PLS regulates auxin concentration, by controlling the relative contribution of auxin transport and biosynthesis and by integrating auxin, ethylene and cytokinin signalling. Modelling further reveals that a bell-shaped dose-response relationship between endogenous auxin and root length is established via PLS. This combined modelling and experimental analysis provides new insights into the integration of hormonal signals in plants.

The hydra mutants of Arabidopsis are characterized by a pleiotropic phenotype that shows defective embryonic and seedling cell patterning, morphogenesis, and root growth. We demonstrate that the HYDRA1 gene encodes a Delta8-Delta7 sterol isomerase, whereas HYDRA2 encodes a sterol C14 reductase, previously identified as the FACKEL gene product. Seedlings mutant for each gene are similarly defective in the concentrations of the three major Arabidopsis sterols. Promoter::reporter gene analysis showed misexpression of the auxin-regulated DR5 and ACS1 promoters and of the epidermal cell file-specific GL2 promoter in the mutants. The mutants exhibit enhanced responses to auxin. The phenotypes can be rescued partially by inhibition of auxin and ethylene signaling but not by exogenous sterols or brassinosteroids. We propose a model in which correct sterol profiles are required for regulated auxin and ethylene signaling through effects on membrane function.

Development can be considered to comprise the co-ordinated regulation of patterning at different levels: patterning of cells to form tissues, patterning of tissues to form organs, and patterning of organs to generate the characteristic architecture of the organism. These processes are expected, in turn, to be mediated by the precise spatial and temporal regulation of patterns of gene expression during development, which depend on appropriate signalling mechanisms. In order to investigate molecular events of morphogenesis in plants, we have utilized a system of promoter trap insertional mutagenesis in Arabidopsis, to generate both phenotypic mutants and gene fusions that represent markers useful in studying the regulation of patterning. A screen of transgenic seedlings containing a T-DNA promoter trap has led to the identification of mutants defective in seedling shape and embryonic development, and of GUS fusion genes that are expressed in spatially restricted patterns. Mutants have been crossed with marker lines expressing cell type-specific GUS activities, to investigate their cellular organization. For example, the POLARIS marker gene is expressed in the embryonic and seedling root tip. When crossed with hydra, which lacks an embryonic root, and with emb30, which lacks both embryonic and seedling roots, it is nevertheless expressed in the correct relative position, and we hypothesize that it represents a novel marker of root positional information, independent of root morphogenesis.

Mutations in the HYDRA1 (HYD1) gene of Arabidopsis thaliana can prevent normal morphological development of embryos and seedlings. Three allelic mutants (hydra 1-1, hydra1-2 and hydra1-3) have been identified, and in each the seedling is characterized by having a variable number of cotyledons, a short and wide hypocotyl and a much reduced root system. hydra1 embryos appear to develop normally to the octant stage, but fail to establish a distinct protoderm and lack bilateral symmetry, developing multiple cotyledonary primordia of irregular size and shape. Cells of the embryo proper, but not the suspensor, exhibit abnormalities in size and shape. The hydra1 embryo fails to develop an embryonic root, but embryos and seedlings express molecular markers of apical-basal polarity. Mutant seedlings produce leaves to form a small cabbage-like habit and may occasionally produce sterile flowers, though the mutation is commonly seedling-lethal. hydra1 seedlings exhibit abnormal radial patterning, but nevertheless express at least one molecular marker of vascular cell differentiation. A model is proposed in which the HYDRA1 protein functions as an essential component of the cell expansion system.

To investigate mechanisms involved in establishing polar organization in Arabidopsis embryos and seedlings, we used promoter trapping to identify molecular markers (beta-glucuronidase fusion genes) expressed in spatially restricted patterns along the apical-basal axis. Three markers were identified that are expressed, respectively, in the embryonic and seedling root tip (POLARIS), cotyledons and shoot and root apices (EXORDIUM), and root cap (COLUMELLA). Each marker was crossed into the mutants hydra and emb30, which are defective in embryonic and seedling morphogenesis. All three markers were expressed in hydra mutants in patterns similar to those observed in phenotypically wild-type embryos and seedlings. In emb30 mutants, the EXORDIUM marker was expressed in cotyledons but not in the expected position of shoot and root meristems, and the marker COLUMELLA was not expressed at all, which is consistent with the view that the emb30 mutant, but not hydra, lacks shoot and root meristems. However, POLARIS was expressed in the basal part of hydra embryos lacking an embryonic root and in the basal parts of both hydra and emb30 seedlings. Expression of POLARIS is inducible by exogenous auxin and suppressed by cytokinin but is unaffected by inhibitors of polar auxin transport or cell division. We conclude that POLARIS differentiates positional aspects of polar development from structural aspects.

The technique of promoter trapping has been exploited to identify markers of embryogenesis in Arabidopsis thaliana. A population of transgenic A. thaliana was generated containing the promoter trap vector pdeltagusBin19, following Agrobacterium tumefaciens-mediated transformation. This vector contains at the T-DNA left border, a promoterless gusA gene, which is activated following integration downstream of a native gene promoter that directs transgene transcription. Screening of a population of 430 independent transgenic lines revealed that 74 lines (17.2%) exhibited GUS activity in siliques, as determined by fluorimetric assay. Histochemical GUS analysis was used to identify lines that expressed GUS in embryos, and three lines that were demonstrated to contain single T-DNA inserts were analysed in detail. Each line showed a distinct pattern of GUS fusion activity. Fusion transcripts were identified, demonstrating that transcription was initiated in the genomic DNA flanking the T-DNA left border. Inverse PCR was used to clone the T-DNA flanking sequences, and for one line a corresponding cDNA was identified, demonstrating that the tagged sequences are transcribed. The markers represent the earliest embryonic genes known to be expressed in plants.