Analysis of Network Topologies Underlying Ethylene Growth Response Kinetics

The physiological and transcriptomic study of secondary growth in Neolamarckia cadamba stimulated by the ethylene precursor ACC

Though many biological roles of ethylene have been investigated intensively, the molecular mechanism of ethylene's action in woody plants remains unclear. In this study, we investigated the effects of exogenous 1-aminocyclopropane-1-carboxylic acid (ACC), the precursor of ethylene, on the growth of Neolamarckia cadamba seedlings, a fast-growing tropical tree. After 14 days of ACC treatment, the plants showed a reduced physiological morphology while stem diameter increased; however, this did not occur after the addition of 1-MCP. Meanwhile, the lignin content of N. cadamba also increased. Transcriptome analysis revealed that the expression of the ethylene biosynthesis and signaling genes ACC oxidase (ACO) and ethylene insensitive 3 (EIN3) were up-regulated mainly at the 6th hour and the 3rd day of the ACC treatment, respectively. The transcription levels of transcription factors, mainly in the basic helix-loop-helix (bHLH), ethylene response factor (ERF), WRKY and v-myb avian myeloblastosis viral oncogene homolog (MYB) families, involved in the ethylene signaling and secondary growth also increased significantly. Furthermore, in accordance to the increased lignification of the stem, the transcriptional level of key enzymes in the phenylalanine pathway were elevated after the ACC treatment. Our results revealed the physiological and molecular mechanisms underlying the secondary growth stimulated by exogenous ACC treatment on N. cadamba seedlings.

Ethylene signaling in plants

Ethylene is a gaseous phytohormone and the first of this hormone class to be discovered. It is the simplest olefin gas and is biosynthesized by plants to regulate plant development, growth, and stress responses via a well-studied signaling pathway. One of the earliest reported responses to ethylene is the triple response. This response is common in eudicot seedlings grown in the dark and is characterized by reduced growth of the root and hypocotyl, an exaggerated apical hook, and a thickening of the hypocotyl. This proved a useful assay for genetic screens and enabled the identification of many components of the ethylene-signaling pathway. These components include a family of ethylene receptors in the membrane of the endoplasmic reticulum (ER); a protein kinase, called constitutive triple response 1 (CTR1); an ER-localized transmembrane protein of unknown biochemical activity, called ethylene-insensitive 2 (EIN2); and transcription factors such as EIN3, EIN3-like (EIL), and ethylene response factors (ERFs). These studies led to a linear model, according to which in the absence of ethylene, its cognate receptors signal to CTR1, which inhibits EIN2 and prevents downstream signaling. Ethylene acts as an inverse agonist by inhibiting its receptors, resulting in lower CTR1 activity, which releases EIN2 inhibition. EIN2 alters transcription and translation, leading to most ethylene responses. Although this canonical pathway is the predominant signaling cascade, alternative pathways also affect ethylene responses. This review summarizes our current understanding of ethylene signaling, including these alternative pathways, and discusses how ethylene signaling has been manipulated for agricultural and horticultural applications.

Comparative RNA-Seq analysis on the regulation of cucumber sex differentiation under different ratios of blue and red light

Cucumber (Cucumis sativus L.) is a typical monoecism vegetable with individual male and female flowers, which has been used as a plant model for sex determination. It is well known that light is one of the most important environmental stimuli, which control the timing of the transition from vegetative growth to reproductive development. However, whether light controls sex determination remains elusive. To unravel this problem, we performed high-throughput RNA-Seq analyses, which compared the transcriptomes of shoot apices between R2B1(Red light:Blue light = 2:1)-treated and R4B1(Red light:Blue light = 4:1)-treated cucumber seedlings. Results showed that the higher proportion of blue light in the R2B1 treatment significantly induced the formation of female flowers and accelerated female flowering time in this whole study. The genes related to flowering time, such as flowering locus T (FT) and SUPPRESSOR OF OVEREXPRESSION OF CO1 (SOC1), were up-regulated after R2B1 treatment. Furthermore, the transcriptome analysis showed that up-regulation and down-regulation of specific DEGs (the differentially expressed genes) were primarily the result of plant hormone signal transduction after treatments. The specific DEGs related with auxin formed the highest percentage of DEGs in the plant hormone signal transduction. In addition, the expression levels of transcription factors also changed after R2B1 treatment. Thus, sex differentiation affected by light quality might be induced by plant hormone signal transduction and transcription factors. These results provide a theoretical basis for further investigation of the regulatory mechanism of female flower formation under different light qualities in cucumber seedlings.

Novel functions of the Arabidopsis transcription factor TCP5 in petal development and ethylene biosynthesis

The flowers of most dicotyledons have petals that, together with the sepals, initially protect the reproductive organs. Later during development petals are required to open the flower and to attract pollinators. This diverse set of functions demands tight temporal and spatial regulation of petal development. We studied the functioning of the Arabidopsis thaliana TCP5-like transcription factors (TFs) in petals. Overexpression of TCP5 in petal epidermal cells results in smaller petals, whereas tcp5 tcp13 tcp17 triple knockout lines have wider petals with an increased surface area. Comprehensive expression studies revealed effects of TCP5-like TFs on the expression of genes related to the cell cycle, growth regulation and organ growth. Additionally, the ethylene biosynthesis genes 1-amino-cyclopropane-1-carboxylate (ACC) synthase 2 (ACS2) and ACC oxidase 2 (ACO2) and several ETHYLENE RESPONSE FACTORS (ERFs) are found to be differentially expressed in TCP5 mutant and overexpression lines. Chromatin immunoprecipitation-quantitative PCR showed direct binding of TCP5 to the ACS2 locus in vivo. Ethylene is known to influence cell elongation, and the petal phenotype of the tcp5 tcp13 tcp17 mutant could be complemented by treatment of the plants with an ethylene pathway inhibitor. Taken together, this reveals a novel role for TCP5-like TFs in the regulation of ethylene-mediated petal development and growth.

Ethylene-Related Gene Expression Networks in Wood Formation

Thickening of tree stems is the result of secondary growth, accomplished by the meristematic activity of the vascular cambium. Secondary growth of the stem entails developmental cascades resulting in the formation of secondary phloem outwards and secondary xylem (i.e., wood) inwards of the stem. Signaling and transcriptional reprogramming by the phytohormone ethylene modifies cambial growth and cell differentiation, but the molecular link between ethylene and secondary growth remains unknown. We addressed this shortcoming by analyzing expression profiles and co-expression networks of ethylene pathway genes using the AspWood transcriptome database which covers all stages of secondary growth in aspen (Populus tremula) stems. ACC synthase expression suggests that the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) is synthesized during xylem expansion and xylem cell maturation. Ethylene-mediated transcriptional reprogramming occurs during all stages of secondary growth, as deduced from AspWood expression profiles of ethylene-responsive genes. A network centrality analysis of the AspWood dataset identified EIN3D and 11 ERFs as hubs. No overlap was found between the co-expressed genes of the EIN3 and ERF hubs, suggesting target diversification and hence independent roles for these transcription factor families during normal wood formation. The EIN3D hub was part of a large co-expression gene module, which contained 16 transcription factors, among them several new candidates that have not been earlier connected to wood formation and a VND-INTERACTING 2 (VNI2) homolog. We experimentally demonstrated Populus EIN3D function in ethylene signaling in Arabidopsis thaliana. The ERF hubs ERF118 and ERF119 were connected on the basis of their expression pattern and gene co-expression module composition to xylem cell expansion and secondary cell wall formation, respectively. We hereby establish data resources for ethylene-responsive genes and potential targets for EIN3D and ERF transcription factors in Populus stem tissues, which can help to understand the range of ethylene targeted biological processes during secondary growth.

Interrogating the topological robustness of gene regulatory circuits by randomization

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Cells are able to robustly carry out their essential biological functions, possibly because of multiple layers of tight regulation via complex, yet well-designed, gene regulatory networks involving a substantial number of genes. State-of-the-art genomics technology has enabled the mapping of these large gene networks, yet it remains a tremendous challenge to elucidate their design principles and the regulatory mechanisms underlying their biological functions such as signal processing and decision-making. One of the key barriers is the absence of accurate kinetics for the regulatory interactions, especially from in vivo experiments. To this end, we have developed a new computational modeling method, Random Circuit Perturbation (RACIPE), to explore the dynamic behaviors of gene regulatory circuits without the requirement of detailed kinetic parameters. RACIPE takes a network topology as the input, and generates an unbiased ensemble of models with varying kinetic parameters. Each model is subjected to simulation, followed by statistical analysis for the ensemble. We tested RACIPE on several gene circuits, and found that the predicted gene expression patterns from all of the models converge to experimentally observed gene state clusters. We expect RACIPE to be a powerful method to identify the role of network topology in determining network operating principles.

Combining in silico evolution and nonlinear dimensionality reduction to redesign responses of signaling networks

The rational design of network behavior is a central goal of synthetic biology. Here, we combine in silico evolution with nonlinear dimensionality reduction to redesign the responses of fixed-topology signaling networks and to characterize sets of kinetic parameters that underlie various input-output relations. We first consider the earliest part of the T cell receptor (TCR) signaling network and demonstrate that it can produce a variety of input-output relations (quantified as the level of TCR phosphorylation as a function of the characteristic TCR binding time). We utilize an evolutionary algorithm (EA) to identify sets of kinetic parameters that give rise to: (i) sigmoidal responses with the activation threshold varied over 6 orders of magnitude, (ii) a graded response, and (iii) an inverted response in which short TCR binding times lead to activation. We also consider a network with both positive and negative feedback and use the EA to evolve oscillatory responses with different periods in response to a change in input. For each targeted input-output relation, we conduct many independent runs of the EA and use nonlinear dimensionality reduction to embed the resulting data for each network in two dimensions. We then partition the results into groups and characterize constraints placed on the parameters by the different targeted response curves. Our approach provides a way (i) to guide the design of kinetic parameters of fixed-topology networks to generate novel input-output relations and (ii) to constrain ranges of biological parameters using experimental data. In the cases considered, the network topologies exhibit significant flexibility in generating alternative responses, with distinct patterns of kinetic rates emerging for different targeted responses.

Transcriptomic Analysis Implies That GA Regulates Sex Expression via Ethylene-Dependent and Ethylene-Independent Pathways in Cucumber (Cucumis sativus L.)

Sex differentiation of flower buds is an important developmental process that directly affects fruit yield of cucumber (Cucumis sativus L.). Plant hormones, such as gibberellins (GAs) and ethylene can promote development of male and female flowers, respectively, however, the regulatory mechanisms of GA-induced male flower formation and potential involvement of ethylene in this process still remain unknown. In this study, to unravel the genes and gene networks involved in GA-regulated cucumber sexual development, we performed high throughout RNA-Seq analyses that compared the transcriptomes of shoot tips between GA₃ treated and untreated gynoecious cucumber plants. Results showed that GA₃ application markedly induced male flowers but decreased ethylene production in shoot tips. Furthermore, the transcript levels of M (CsACS2) gene, ethylene receptor CsETR1 and some ethylene-responsive transcription factors were dramatically changed after GA₃ treatment, suggesting a potential involvement of ethylene in GA-regulated sex expression of cucumber. Interestingly, GA₃ down-regulated transcript of a C-class floral homeotic gene, CAG2, indicating that GA may also influence cucumber sex determination through an ethylene-independent process. These results suggest a novel model for hormone-mediated sex differentiation and provide a theoretical basis for further dissection of the regulatory mechanism of male flower formation in cucumber. Statement: We reveal that GA can regulate sex expression of cucumber via an ethylene-dependent manner, and the M (CsACS2), CsETR1, and ERFs are probably involved in this process. Moreover, CAG2, a C-class floral homeotic gene, may also participate in GA-modulated cucumber sex determination, but this pathway is ethylene-independent.