NIN Acts as a Network Hub Controlling a Growth Module Required for Rhizobial Infection.

The symbiotic infection of root cells by nitrogen-fixing rhizobia during nodulation requires the transcription factor Nodule Inception (NIN). Our root hair transcriptomic study extends NIN's regulon to include Rhizobium Polar Growth and genes involved in cell wall modification, gibberellin biosynthesis, and a comprehensive group of nutrient (N, P, and S) uptake and assimilation genes, suggesting that NIN's recruitment to nodulation was based on its role as a growth module, a role shared with other NIN-Like Proteins. The expression of jasmonic acid genes in nin suggests the involvement of NIN in the resolution of growth versus defense outcomes. We find that the regulation of the growth module component Nodulation Pectate Lyase by NIN, and its function in rhizobial infection, are conserved in hologalegina legumes, highlighting its recruitment as a major event in the evolution of nodulation. We find that Nodulation Pectate Lyase is secreted to the infection chamber and the lumen of the infection thread. Gene network analysis using the transcription factor mutants for ERF Required for Nodulation1 and Nuclear Factor-Y Subunit A1 confirms hierarchical control of NIN over Nuclear Factor-Y Subunit A1 and shows that ERF Required for Nodulation1 acts independently to control infection. We conclude that while NIN shares functions with other NIN-Like Proteins, the conscription of key infection genes to NIN's control has made it a central regulatory hub for rhizobial infection.

SlHY5 Integrates Temperature, Light, and Hormone Signaling to Balance Plant Growth and Cold Tolerance.

During the transition from warm to cool seasons, plants experience decreased temperatures, shortened days, and decreased red/far-red (R/FR) ratios of light. The mechanism by which plants integrate these environmental cues to maintain plant growth and adaptation remains poorly understood. Here, we report that low temperature induced the transcription of PHYTOCHROME A and accumulation of LONG HYPOCOTYL5 (SlHY5, a basic Leu zipper transcription factor) in tomato (Solanum lycopersicum) plants, especially under short day conditions with low R/FR light ratios. Reverse genetic approaches and physiological analyses revealed that silencing of SlHY5 increased cold susceptibility in tomato plants, whereas overexpression of SlHY5 enhanced cold tolerance. SlHY5 directly bound to and activated the transcription of genes encoding a gibberellin-inactivation enzyme, namely GIBBERELLIN2-OXIDASE4, and an abscisic acid biosynthetic enzyme, namely 9-CIS-EPOXYCAROTENOID DIOXYGENASE6 (SlNCED6). Thus, phytochrome A-dependent SlHY5 accumulation resulted in an increased abscisic acid/gibberellin ratio, which was accompanied by growth cessation and induction of cold response. Furthermore, silencing of SlNCED6 compromises short day- and low R/FR-induced tomato resistance to cold stress. These findings provide insight into the molecular genetic mechanisms by which plants integrate environmental stimuli with hormones to coordinate their growth with impending cold temperatures. Moreover, this work reveals a molecular mechanism that plants have evolved for growth and survival in response to seasonal changes.

Cytokinin-induced parthenocarpy of San Pedro type fig (Ficus carica L.) main crop: explained by phytohormone assay and transcriptomic network comparison.

KEY MESSAGE: CPPU-induced San Pedro type fig main crop parthenocarpy exhibited constantly increasing IAA content and more significantly enriched KEGG pathways in the receptacle than in female flowers. N-(2-chloro-4-pyridyl)-N-phenylurea (CPPU) was applied to San Pedro fig (Ficus carica L.) main crop to induce parthenocarpy; the optimal effect was obtained with 25 mg L-1 application to syconia when female flowers were at anthesis. To elucidate the key expression changes in parthenocarpy conversion, significant changes in phytohormone level and transcriptome of fig female flowers and receptacles were monitored. HPLC-MS revealed increased IAA content in female flowers and receptacle 2, 4 and 10 days after treatment (DAT), decreased zeatin level in the receptacle 2, 4 and 10 DAT, decreased GA3 content 2 and 4 DAT, and increased GA3 content 10 DAT. ABA level increased 2 and 4 DAT, and decreased 10 DAT. CPPU-treated syconia released more ethylene than the control except 2 DAT. RNA-Seq and bioinformatics analysis revealed notably more differentially expressed KEGG pathways in the receptacle than in female flowers. In the phytohormone gene network, GA-biosynthesis genes GA20ox and GA3ox were upregulated, along with GA signal-transduction genes GID1 and GID2, and IAA-signaling genes AUX/IAA and GH3. ABA-biosynthesis gene NCED and signaling genes PP2C and ABF were downregulated 10 DAT. One ACO gene showed consistent upregulation in both female flowers and receptacle after CPPU treatment, and more than a dozen of ERFs demonstrated opposing changes in expression. Our results revealed early-stage spatiotemporal phytohormone and transcriptomic responses in CPPU-induced San Pedro fig main crop parthenocarpy, which could be valuable for further understanding the nature of the parthenocarpy of different fig types.

Aberrant Stamen Development is Associated with Parthenocarpic Fruit Set Through Up-Regulation of Gibberellin Biosynthesis in Tomato.

Parthenocarpy, a process in which fruit set occurs without fertilization, leads to the production of seedless fruit. A number of floral homeotic mutants with abnormal stamen development exhibit parthenocarpic fruit set. Flower development is thought to repress ovary growth before anthesis. However, the mechanism of parthenocarpic fruit development caused by aberrant flower formation is poorly understood. To investigate the molecular mechanism of parthenocarpic fruit development in floral homeotic mutants, we performed functional analysis of Tomato APETALA3 (TAP3) by loss-of-function approaches. Organ-specific promoter was used to induce organ-specific loss of function in stamen and ovary/fruit. We observed increased cell expansion in tap3 mutants and TAP3-RNAi lines during parthenocarpic fruit growth. These were predominantly accompanied by the up-regulation of GA biosynthesis genes, including SlGA20ox1, SlGA20ox2, and SlGA20ox3, as well as reduced expression of the GA-inactivating gene SlGA2ox1 and the auxin signaling gene SlARF7 involved in a crosstalk between GA and auxin. These transcriptional profiles are in agreement with the GA levels in these lines. These results suggest that stamen development negatively regulates fruit set by repressing the GA biosynthesis.

NO FLOWERING IN SHORT DAY (NFL) is a bHLH transcription factor that promotes flowering specifically under short-day conditions in Arabidopsis.

Flowering in plants is a dynamic and synchronized process where various cues including age, day length, temperature and endogenous hormones fine-tune the timing of flowering for reproductive success. Arabidopsis thaliana is a facultative long day (LD) plant where LD photoperiod promotes flowering. Arabidopsis still flowers under short-day (SD) conditions, albeit much later than in LD conditions. Although factors regulating the inductive LD pathway have been extensively investigated, the non-inductive SD pathway is much less understood. Here, we identified a key basic helix-loop-helix transcription factor called NFL (NO FLOWERING IN SHORT DAY) that is essential to induce flowering specifically under SD conditions in Arabidopsis. nfl mutants do not flower under SD conditions, but flower similar to the wild type under LD conditions. The no-flowering phenotype in SD is rescued either by exogenous application of gibberellin (GA) or by introducing della quadruple mutants in the nfl background, suggesting that NFL acts upstream of GA to promote flowering. NFL is expressed at the meristematic regions and NFL is localized to the nucleus. Quantitative RT-PCR assays using apical tissues showed that GA biosynthetic genes are downregulated and the GA catabolic and receptor genes are upregulated in the nfl mutant compared with the wild type, consistent with the perturbation of the endogenous GA biosynthetic and catabolic intermediates in the mutant. Taken together, these data suggest that NFL is a key transcription factor necessary for promotion of flowering under non-inductive SD conditions through the GA signaling pathway.

The high-affinity transporter BnPHT1;4 is involved in phosphorus acquisition and mobilization for facilitating seed germination and early seedling growth of Brassica napus.

BACKGROUND: Seed germination and seedling establishment are two of the most critical phases in plant development. However, the molecular mechanisms underlying the effect of phosphorus on seed germination and post-germinated growth of oilseed rape are unclear so far. Here, we report the role of BnPHT1;4 in seed germination and early seedling development of Brassica napus. RESULTS: Our results show that BnPHT1;4 is preferentially expressed in cotyledons of early developing seedlings. Overexpression of BnPHT1;4 in oilseed rape promoted seed germination and seedling growth. Expression levels of the genes related to ABA and GA biosynthesis and signaling were significantly altered in BnPHT1;4 transgenic seedlings. Consequently, active GA level was up-regulated, whereas ABA content was down-regulated in BnPHT1;4 transgenic seedlings. Furthermore, exogenous GA could promote seed germination of wild type, while exogenous ABA could partially recover the advanced-germination phenotype of BnPHT1;4 transgenic seeds. Total phosphorus content in cotyledons of the transgenic seedlings was decreased more rapidly than that in wild type when Pi was supplied or deficient, and Pi contents in shoots and roots of the BnPHT1;4 transgenic plants were higher than those in wild type under high and low Pi conditions. CONCLUSIONS: Our data suggest that the high-affinity transporter BnPHT1;4 is involved in phosphorus acquisition and mobilization for facilitating seed germination and seedling growth of Brassica napus by modulating ABA and GA biosynthesis.

Time-Course Transcriptomics Analysis Reveals Key Responses of Submerged Deepwater Rice to Flooding.

Water submergence is an environmental factor that limits plant growth and survival. Deepwater rice (Oryza sativa) adapts to submergence by rapidly elongating its internodes and thereby maintaining its leaves above the water surface. We performed a comparative RNA sequencing transcriptome analysis of the shoot base region, including basal nodes, internodes, and shoot apices of seedlings at two developmental stages from two varieties with contrasting deepwater growth responses. A transcriptomic comparison between deepwater rice cv C9285 and nondeepwater rice cv Taichung 65 revealed both similar and differential expression patterns between the two genotypes during submergence. The expression of genes related to gibberellin biosynthesis, trehalose biosynthesis, anaerobic fermentation, cell wall modification, and transcription factors that include ethylene-responsive factors was significantly different between the varieties. Interestingly, in both varieties, the jasmonic acid content at the shoot base decreased during submergence, while exogenous jasmonic acid inhibited submergence-induced internode elongation in cv C9285, suggesting that jasmonic acid plays a role in the submergence response of rice. Furthermore, a targeted de novo transcript assembly revealed transcripts that were specific to cv C9285, including submergence-induced biotic stress-related genes. Our multifaceted transcriptome approach using the rice shoot base region illustrates a differential response to submergence between deepwater and nondeepwater rice. Jasmonic acid metabolism appears to participate in the submergence-mediated internode elongation response of deepwater rice.

Elucidation of gibberellin biosynthesis in bacteria reveals convergent evolution.

Gibberellins (GAs) are crucial phytohormones involved in many aspects of plant growth and development, including plant-microbe interactions, which has led to GA production by plant-associated fungi and bacteria as well. While the GA biosynthetic pathways in plants and fungi have been elucidated and found to have arisen independently through convergent evolution, little has been uncovered about GA biosynthesis in bacteria. Some nitrogen-fixing, symbiotic, legume-associated rhizobia, including Bradyrhizobium japonicum-the symbiont of soybean-and Sinorhizobium fredii-a broad-host-nodulating species-contain a putative GA biosynthetic operon, or gene cluster. Through functional characterization of five unknown genes, we demonstrate that this operon encodes the enzymes necessary to produce GA9, thereby elucidating bacterial GA biosynthesis. The distinct nature of these enzymes indicates that bacteria have independently evolved a third biosynthetic pathway for GA production. Furthermore, our results also reveal a central biochemical logic that is followed in all three convergently evolved GA biosynthetic pathways.

The role of gibberellin synthase gene GhGA2ox1 in upland cotton (Gossypium hirsutum L.) responses to drought and salt stress.

Gibberellins (GAs) is one kind of important endogenous hormone in plants that regulates vegetative and reproductive growth of plants. GA2ox is a class of oxidase that plays a regulatory role in the third stage of GAs synthesis. In this paper, we cloned the GhGA2ox1 gene from an upland cotton (Gossypium hirsutum L. var. CCRI49). The results showed that the CDS of GhGA2ox1 is 996 bp, which encode 331 amino acids, which has high homology with GhGA2ox2 and NtGA2ox1. The quantitative real-time PCR showed that under the conditions of salt and drought stress, the expression of GhGA2ox1 had a higher upregulation in root, stem, and leaf of transgenic plant, compared with non-transgenic plant. Cotton plant that overexpressed the GhGA2ox1 gene showed higher drought and salt tolerance than non-transgenic cotton plant, and these results were supported by data of higher free proline, chlorophyll, and relative water content in transgenic plant compared with control plant. The expression level of antiretroviral genes, including GhEREB2, GhDREB1, GhWRKY5, and GhP5CS, was upregulated to varying degrees in transgenic plant. The above results indicate that overexpressed GhGA2ox1 gene can increases the tolerance to respond to drought and salt stress in upland cotton.

Probing the promiscuity of ent-kaurene oxidases via combinatorial biosynthesis.

The substrate specificity of enzymes from natural products' metabolism is a topic of considerable interest, with potential biotechnological use implicit in the discovery of promiscuous enzymes. However, such studies are often limited by the availability of substrates and authentic standards for identification of the resulting products. Here, a modular metabolic engineering system is used in a combinatorial biosynthetic approach toward alleviating this restriction. In particular, for studies of the multiply reactive cytochrome P450, ent-kaurene oxidase (KO), which is involved in production of the diterpenoid plant hormone gibberellin. Many, but not all, plants make a variety of related diterpenes, whose structural similarity to ent-kaurene makes them potential substrates for KO. Use of combinatorial biosynthesis enabled analysis of more than 20 such potential substrates, as well as structural characterization of 12 resulting unknown products, providing some insight into the underlying structure-function relationships. These results highlight the utility of this approach for investigating the substrate specificity of enzymes from complex natural products' biosynthesis.

NbGIS regulates glandular trichome initiation through GA signaling in tobacco.

KEY MESSAGE: A novel gene NbGIS positively regulates glandular trichome initiation through GA Signaling in tobacco. NbMYB123-like regulates glandular trichome initiation by acting downstream of NbGIS in tobacco. Glandular trichome is a specialized multicellular structure which has capability to synthesize and secrete secondary metabolites and protects plants from biotic and abiotic stresses. Our previous results revealed that a C2H2 zinc-finger transcription factor GIS and its sub-family genes act upstream of GL3/EGL3-GL1-TTG1 transcriptional activator complex to regulate trichome initiation in Arabidopsis. In this present study, we found that NbGIS could positively regulate glandular trichome development in Nicotiana benthamiana (tobacco). Our result demonstrated that 35S:NbGIS lines exhibited much higher densities of trichome on leaves, main stems, lateral branches and sepals than WT plants, while NbGIS:RNAi lines had the opposite phenotypes. Furthermore, our results also showed that NbGIS was required in response to GA signal to control glandular trichome initiation in Nicotiana benthamiana. In addition, our results also showed that NbGIS significantly influenced GA accumulation and expressions of marker genes of the GA biosynthesis, might result in the changes of growth and maturation in tobacco. Lastly, our results also showed that NbMYB123-like regulated glandular trichome initiation in tobacco by acting downstream of NbGIS. These findings provide new insights to discover the molecular mechanism by which C2H2 transcriptional factors regulates glandular trichome initiation through GA signaling pathway in tobacco.

The Role of Arabidopsis Inositol Polyphosphate Kinase AtIPK2ß in Glucose Suppression of Seed Germination and Seedling Development.

Seed germination and subsequent seedling development are critical phases in plants. These processes are regulated by a complex molecular network in which sugar has been reported to play an essential role. However, factors affecting sugar responses remain to be fully elucidated. In this study, we demonstrate that AtIPK2ß, known to participate in the synthesis of myo-inositol 1,2,3,4,5,6-hexakisphosphate (IP6, phytate), affects Arabidopsis responses to glucose during seed germination. The loss-of-function mutant atipk2ß showed increased sensitivity to 6% glucose and paclobutrazol (PAC). Yeast two-hybrid assay showed that AtIPK2ß interacts with sucrose non-fermenting-1-related protein kinase (SnRK1.1), and bimolecular fluorescence complementation (BiFC) and pull-down assay further confirmed this interaction. Moreover, AtIPK2ß was phosphorylated by SnRK1.1 in vitro, and the effect of restoring AtIPK2ß to yeast cells lacking IPK2 (Δipk2) was abolished by catalytically active SnRK1.1. Further analysis indicated that IP6 reduces the suppression of seed germination caused by glucose, accompanied by altered expression levels of glucose-/hormone-responsive genes. Collectively, these findings indicate that AtIPK2ß and IP6 are involved in glucose suppression of seed germination and that AtIPK2ß enzyme activity is likely to be regulated by SnRK1.1.

PIL5 represses floral transition in Arabidopsis under long day conditions.

PHYTOCHROME INTERACING FACTOR 3 LIKE 5 (PIL5), also named PHYTOCHROME INTERACTING FACTOR 1 (PIF1) is an important b-HLH transcription factor in Arabidopsis thaliana. Here we show that mutant of pil5-1 displays early flowering phenotype. We demonstrate that the expressions of the major flowering promoter genes [FLOWERING LOCUS T (FT), SUPPRESOR OF OVEREXPRESSION OF CO 1 (SOC1), and LEAFY (LFY)] are upregulated in the mutant of pil5-1. There is a significant increase of the mRNA of PIL5 in the mutants of co2-1, ft-10, soc1-2, and lfy-4. These changes provide the molecular evidence that PIL5 interacts with the flowering regulators to control flowering time. Moreover, it is shown in our results that PIL5 mutation mediates the increased contents of gibberellic acid (GA). Which is further supported by the qRT-PCR analysis, an increased transcriptome level of the GA biosynthesis genes (GA3ox1, GA3ox2, GA20ox1, GA20ox2, and GA20ox3) has been observed in the pil5-1 mutants as compared to the wild type. Collectively, PIL5 is involved in floral transition interacting with flowering integrators and GA.

Isoprenoid-derived plant signaling molecules: biosynthesis and biological importance.

MAIN CONCLUSION: The present review summarizes current knowledge of the biosynthesis and biological importance of isoprenoid-derived plant signaling compounds. Cellular organisms use chemical signals for intercellular communication to coordinate their growth, development, and responses to environmental cues. The skeletons of majority of plant signaling molecules, mediators of plant intercellular 'broadcasting', are built from C5 units of isoprene and therefore belong to a huge and diverse group of natural substances called isoprenoids (terpenoids). They fill many important roles in nature. This review summarizes current knowledge of the biosynthesis and biological importance of a group of isoprenoid-derived plant signaling compounds.

From network to phenotype: the dynamic wiring of an Arabidopsis transcriptional network induced by osmotic stress.

Plants have established different mechanisms to cope with environmental fluctuations and accordingly fine-tune their growth and development through the regulation of complex molecular networks. It is largely unknown how the network architectures change and what the key regulators in stress responses and plant growth are. Here, we investigated a complex, highly interconnected network of 20 Arabidopsis transcription factors (TFs) at the basis of leaf growth inhibition upon mild osmotic stress. We tracked the dynamic behavior of the stress-responsive TFs over time, showing the rapid induction following stress treatment, specifically in growing leaves. The connections between the TFs were uncovered using inducible overexpression lines and were validated with transient expression assays. This study resulted in the identification of a core network, composed of ERF6, ERF8, ERF9, ERF59, and ERF98, which is responsible for most transcriptional connections. The analyses highlight the biological function of this core network in environmental adaptation and its redundancy. Finally, a phenotypic analysis of loss-of-function and gain-of-function lines of the transcription factors established multiple connections between the stress-responsive network and leaf growth.

A brassinosteroid responsive miRNA-target module regulates gibberellin biosynthesis and plant development.

Plant growth and development are highly coordinated by hormones, including brassinosteroid (BR) and gibberellin (GA). Although much progress has been made in understanding the fundamental signaling transduction in BR and GA, their relationship remains elusive in rice. Here, we show that BR suppresses the level of OsmiR159d, which cleaves the target OsGAMYBL2 gene. The OsmiR159d-OsGAMYBL2 pair functions as an early BR-responsive module regulating the expression of BU1, a BR-regulated gene involved in BR signaling, and CPS1 and GA3ox2, two genes in GA biosynthesis, by binding to the promoters of these genes. Furthermore, OsGSK2, a key negative player in BR signaling, interacts with OsGAMYBL2 and prevents it from being degraded under 24-epibrassinolide treatment, whereas SLR1, a rice DELLA protein negatively regulating GA signaling, interacts with OsGAMYBL2 and prevents OsGAMYBL2 from binding to the target gene promoter. GA signaling induces degradation of OsGAMYBL2 and, consequently, enhances BR signaling. These results demonstrate that a BR-responsive module acts as a common component functioning in both BR and GA pathways, which connects BR signaling and GA biosynthesis, and thus coordinates the regulation of BR and GA in plant growth and development.

Activation of Mitochondrial Protein Phosphatase SLP2 by MIA40 Regulates Seed Germination.

Reversible protein phosphorylation catalyzed by protein kinases and phosphatases represents the most prolific and well-characterized posttranslational modification known. Here, we demonstrate that Arabidopsis (Arabidopsis thaliana) Shewanella-like protein phosphatase 2 (AtSLP2) is a bona fide Ser/Thr protein phosphatase that is targeted to the mitochondrial intermembrane space (IMS) where it interacts with the mitochondrial oxidoreductase import and assembly protein 40 (AtMIA40), forming a protein complex. Interaction with AtMIA40 is necessary for the phosphatase activity of AtSLP2 and is dependent on the formation of disulfide bridges on AtSLP2. Furthermore, by utilizing atslp2 null mutant, AtSLP2 complemented and AtSLP2 overexpressing plants, we identify a function for the AtSLP2-AtMIA40 complex in negatively regulating gibberellic acid-related processes during seed germination. Results presented here characterize a mitochondrial IMS-localized protein phosphatase identified in photosynthetic eukaryotes as well as a protein phosphatase target of the highly conserved eukaryotic MIA40 IMS oxidoreductase.

Brassinosteroids Are Master Regulators of Gibberellin Biosynthesis in Arabidopsis.

Plant growth and development are highly regulated processes that are coordinated by hormones including the brassinosteroids (BRs), a group of steroids with structural similarity to steroid hormones of mammals. Although it is well understood how BRs are produced and how their signals are transduced, BR targets, which directly confer the hormone's growth-promoting effects, have remained largely elusive. Here, we show that BRs regulate the biosynthesis of gibberellins (GAs), another class of growth-promoting hormones, in Arabidopsis thaliana. We reveal that Arabidopsis mutants deficient in BR signaling are severely impaired in the production of bioactive GA, which is correlated with defective GA biosynthetic gene expression. Expression of the key GA biosynthesis gene GA20ox1 in the BR signaling mutant bri1-301 rescues many of its developmental defects. We provide evidence that supports a model in which the BR-regulated transcription factor BES1 binds to a regulatory element in promoters of GA biosynthesis genes in a BR-induced manner to control their expression. In summary, our study underscores a role of BRs as master regulators of GA biosynthesis and shows that this function is of major relevance for the growth and development of vascular plants.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis.

KEY MESSAGE: GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

Parameter inference to motivate asymptotic model reduction: An analysis of the gibberellin biosynthesis pathway.

Developing effective strategies to use models in conjunction with experimental data is essential to understand the dynamics of biological regulatory networks. In this study, we demonstrate how combining parameter estimation with asymptotic analysis can reveal the key features of a network and lead to simplified models that capture the observed network dynamics. Our approach involves fitting the model to experimental data and using the profile likelihood to identify small parameters and cases where model dynamics are insensitive to changing particular individual parameters. Such parameter diagnostics provide understanding of the dominant features of the model and motivate asymptotic model reductions to derive simpler models in terms of identifiable parameter groupings. We focus on the particular example of biosynthesis of the plant hormone gibberellin (GA), which controls plant growth and has been mutated in many current crop varieties. This pathway comprises two parallel series of enzyme-substrate reactions, which have previously been modelled using the law of mass action (Middleton et al., 2012). Considering the GA20ox-mediated steps, we analyse the identifiability of the model parameters using published experimental data; the analysis reveals the ratio between enzyme and GA levels to be small and motivates us to perform a quasi-steady state analysis to derive a reduced model. Fitting the parameters in the reduced model reveals additional features of the pathway and motivates further asymptotic analysis which produces a hierarchy of reduced models. Calculating the Akaike information criterion and parameter confidence intervals enables us to select a parsimonious model with identifiable parameters. As well as demonstrating the benefits of combining parameter estimation and asymptotic analysis, the analysis shows how GA biosynthesis is limited by the final GA20ox-mediated steps in the pathway and generates a simple mathematical description of this part of the GA biosynthesis pathway.