Optimal Performance Objectives in the Highly Conserved Bone Morphogenetic Protein Signaling Pathway.

Throughout development, complex networks of cell signaling pathways drive cellular decision-making across different tissues and contexts. The transforming growth factor ß (TGF-ß) pathways, including the BMP/Smad pathway, play crucial roles in these cellular responses. However, as the Smad pathway is used reiteratively throughout the life cycle of all animals, its systems-level behavior varies from one context to another, despite the pathway connectivity remaining nearly constant. For instance, some cellular systems require a rapid response, while others require high noise filtering. In this paper, we examine how the BMP- Smad pathway balances trade-offs among three such systems-level behaviors, or "Performance Objectives (POs)": response speed, noise amplification, and the sensitivity of pathway output to receptor input. Using a Smad pathway model fit to human cell data, we show that varying non-conserved parameters (NCPs) such as protein concentrations, the Smad pathway can be tuned to emphasize any of the three POs and that the concentration of nuclear phosphatase has the greatest effect on tuning the POs. However, due to competition among the POs, the pathway cannot simultaneously optimize all three, but at best must balance trade-offs among the POs. We applied the multi-objective optimization concept of the Pareto Front, a widely used concept in economics to identify optimal trade-offs among various requirements. We show that the BMP pathway efficiently balances competing POs across species and is largely Pareto optimal. Our findings reveal that varying the concentration of NCPs allows the Smad signaling pathway to generate a diverse range of POs. This insight identifies how signaling pathways can be optimally tuned for each context.

ISRES+: an improved evolutionary strategy for function minimization to estimate the free parameters of systems biology models.

MOTIVATION: Mathematical models in systems biology help generate hypotheses, guide experimental design, and infer the dynamics of gene regulatory networks. These models are characterized by phenomenological or mechanistic parameters, which are typically hard to measure. Therefore, efficient parameter estimation is central to model development. Global optimization techniques, such as evolutionary algorithms (EAs), are applied to estimate model parameters by inverse modeling, i.e. calibrating models by minimizing a function that evaluates a measure of the error between model predictions and experimental data. EAs estimate model parameters "fittest individuals" by generating a large population of individuals using strategies like recombination and mutation over multiple "generations." Typically, only a few individuals from each generation are used to create new individuals in the next generation. Improved Evolutionary Strategy by Stochastic Ranking (ISRES), proposed by Runnarson and Yao, is one such EA that is widely used in systems biology to estimate parameters. ISRES uses information at most from a pair of individuals in any generation to create a new population to minimize the error. In this article, we propose an efficient evolutionary strategy, ISRES+, which builds on ISRES by combining information from all individuals across the population and across all generations to develop a better understanding of the fitness landscape. RESULTS: ISRES+ uses the additional information generated by the algorithm during evolution to approximate the local neighborhood around the best-fit individual using linear least squares fits in one and two dimensions, enabling efficient parameter estimation. ISRES+ outperforms ISRES and results in fitter individuals with a tighter distribution over multiple runs, such that a typical run of ISRES+ estimates parameters with a higher goodness-of-fit compared with ISRES. AVAILABILITY AND IMPLEMENTATION: Algorithm and implementation: Github-https://github.com/gtreeves/isres-plus-bandodkar-2022.

C4 gene induction during de-etiolation evolved through changes in cis to allow integration with ancestral C3 gene regulatory networks.

C4 photosynthesis has evolved by repurposing enzymes found in C3 plants. Compared with the ancestral C3 state, accumulation of C4 cycle proteins is enhanced. We used de-etiolation of C4 Gynandropsis gynandra and C3 Arabidopsis thaliana to understand this process. C4 gene expression and chloroplast biogenesis in G. gynandra were tightly coordinated. Although C3 and C4 photosynthesis genes showed similar induction patterns, in G. gynandra, C4 genes were more strongly induced than orthologs from A. thaliana. In vivo binding of TGA and homeodomain as well as light-responsive elements such as G- and I-box motifs were associated with the rapid increase in transcripts of C4 genes. Deletion analysis confirmed that regions containing G- and I-boxes were necessary for high expression. The data support a model in which accumulation of transcripts derived from C4 photosynthesis genes in C4 leaves is enhanced because modifications in cis allowed integration into ancestral transcriptional networks.

Experimental and Analytical Framework for "Mix-and-Read" Assays Based on Split Luciferase.

The use of immunodetection assays including the widely used enzyme-linked immunosorbent assay (ELISA) in applications such as point-of-care detection is often limited by the need for protein immobilization and multiple binding and washing steps. Here, we describe an experimental and analytical framework for the development of simple and modular "mix-and-read" enzymatic complementation assays based on split luciferase that enable sensitive detection and quantification of analytes in solution. In this assay, two engineered protein binders targeting nonoverlapping epitopes on the target analyte were each fused to nonactive fragments of luciferase to create biosensor probes. Binding proteins to two model targets, lysozyme and Sso6904, were isolated from a combinatorial library of Sso7d mutants using yeast surface display. In the presence of the analyte, probes were brought into close proximity, reconstituting enzymatic activity of luciferase and enabling detection of low picomolar concentrations of the analyte by chemiluminescence. Subsequently, we constructed an equilibrium binding model that relates binding affinities of the binding proteins for the target, assay parameters such as the concentrations of probes used, and assay performance (limit of detection and concentration range over which the target can be quantified). Overall, our experimental and analytical framework provides the foundation for the development of split luciferase assays for detection and quantification of various targets.

Using breeding and quantitative genetics to understand the C4 pathway.

Reducing photorespiration in C3 crops could significantly increase rates of photosynthesis and yield. One method to achieve this would be to integrate C4 photosynthesis into C3 species. This objective is challenging as it involves engineering incompletely understood traits into C3 leaves, including complex changes to their biochemistry, cell biology, and anatomy. Quantitative genetics and selective breeding offer underexplored routes to identify regulators of these processes. We first review examples of natural intraspecific variation in C4 photosynthesis as well as the potential for hybridization between C3 and C4 species. We then discuss how quantitative genetic approaches including artificial selection and genome-wide association could be used to better understand the C4 syndrome and in so doing guide the engineering of the C4 pathway into C3 crops.

Monocotyledonous plants graft at the embryonic root-shoot interface.

Grafting is possible in both animals and plants. Although in animals the process requires surgery and is often associated with rejection of non-self, in plants grafting is widespread, and has been used since antiquity for crop improvement1. However, in the monocotyledons, which represent the second largest group of terrestrial plants and include many staple crops, the absence of vascular cambium is thought to preclude grafting2. Here we show that the embryonic hypocotyl allows intra- and inter-specific grafting in all three monocotyledon groups: the commelinids, lilioids and alismatids. We show functional graft unions through histology, application of exogenous fluorescent dyes, complementation assays for movement of endogenous hormones, and growth of plants to maturity. Expression profiling identifies genes that unify the molecular response associated with grafting in monocotyledons and dicotyledons, but also gene families that have not previously been associated with tissue union. Fusion of susceptible wheat scions to oat rootstocks confers resistance to the soil-borne pathogen Gaeumannomyces graminis. Collectively, these data overturn the consensus that monocotyledons cannot form graft unions, and identify the hypocotyl (mesocotyl in grasses) as a meristematic tissue that allows this process. We conclude that graft compatibility is a shared ability among seed-bearing plants.

Discovery of Membrane-Permeating Cyclic Peptides via mRNA Display.

Small synthetic peptides capable of crossing biological membranes represent valuable tools in cell biology and drug delivery. While several cell-penetrating peptides (CPPs) of natural or synthetic origin have been reported, no peptide is currently known to cross both cytoplasmic and outer embryonic membranes. Here, we describe a method to engineer membrane-permeating cyclic peptides (MPPs) with broad permeation activity by screening mRNA display libraries of cyclic peptides against embryos at different developmental stages. The proposed method was demonstrated by identifying peptides capable of permeating Drosophila melanogaster (fruit fly) embryos and mammalian cells. The selected peptide cyclo[Glut-MRKRHASRRE-K*] showed a strong permeation activity of embryos exposed to minimal permeabilization pretreatment, as well as human embryonic stem cells and a murine fibroblast cell line. Notably, in both embryos and mammalian cells, the cyclic peptide outperformed its linear counterpart and the control MPPs. Confocal microscopy and single cell flow cytometry analysis were utilized to assess the degree of permeation both qualitatively and quantitatively. These MPPs have potential application in studying and nondisruptively controlling intracellular or intraembryonic processes.

Tunable self-cleaving ribozymes for modulating gene expression in eukaryotic systems.

Advancements in the field of synthetic biology have been possible due to the development of genetic tools that are able to regulate gene expression. However, the current toolbox of gene regulatory tools for eukaryotic systems have been outpaced by those developed for simple, single-celled systems. Here, we engineered a set of gene regulatory tools by combining self-cleaving ribozymes with various upstream competing sequences that were designed to disrupt ribozyme self-cleavage. As a proof-of-concept, we were able to modulate GFP expression in mammalian cells, and then showed the feasibility of these tools in Drosophila embryos. For each system, the fold-reduction of gene expression was influenced by the location of the self-cleaving ribozyme/upstream competing sequence (i.e. 5' vs. 3' untranslated region) and the competing sequence used. Together, this work provides a set of genetic tools that can be used to tune gene expression across various eukaryotic systems.

Robustness of the Dorsal morphogen gradient with respect to morphogen dosage.

In multicellular organisms, the timing and placement of gene expression in a developing tissue assigns the fate of each cell in the embryo in order for a uniform field of cells to differentiate into a reproducible pattern of organs and tissues. This positional information is often achieved through the action of spatial gradients of morphogens. Spatial patterns of gene expression are paradoxically robust to variations in morphogen dosage, given that, by definition, gene expression must be sensitive to morphogen concentration. In this work we investigate the robustness of the Dorsal/NF-κB signaling module with respect to perturbations to the dosage of maternally-expressed dorsal mRNA. The Dorsal morphogen gradient patterns the dorsal-ventral axis of the early Drosophila embryo, and we found that an empirical description of the Dorsal gradient is highly sensitive to maternal dorsal dosage. In contrast, we found experimentally that gene expression patterns are highly robust. Although the components of this signaling module have been characterized in detail, how their function is integrated to produce robust gene expression patterns to variations in the dorsal maternal dosage is still unclear. Therefore, we analyzed a mechanistic model of the Dorsal signaling module and found that Cactus, a cytoplasmic inhibitor for Dorsal, must be present in the nucleus for the system to be robust. Furthermore, active Toll, the receptor that dissociates Cactus from Dorsal, must be saturated. Finally, the vast majority of robust descriptions of the system require facilitated diffusion of Dorsal by Cactus. Each of these three recently-discovered mechanisms of the Dorsal module are critical for robustness. These mechanisms synergistically contribute to changing the amplitude and shape of the active Dorsal gradient, which is required for robust gene expression. Our work highlights the need for quantitative understanding of biophysical mechanisms of morphogen gradients in order to understand emergent phenotypes, such as robustness.

Formation, interpretation, and regulation of the Drosophila Dorsal/NF-κB gradient.

The morphogen gradient of the transcription factor Dorsal in the early Drosophila embryo has become one of the most widely studied tissue patterning systems. Dorsal is a Drosophila homolog of mammalian NF-κB and patterns the dorsal-ventral axis of the blastoderm embryo into several tissue types by spatially regulating upwards of 100 zygotic genes. Recent studies using fluorescence microscopy and live imaging have quantified the Dorsal gradient and its target genes, which has paved the way for mechanistic modeling of the gradient. In this review, we describe the mechanisms behind the initiation of the Dorsal gradient and its regulation of target genes. The main focus of the review is a discussion of quantitative and computational studies of the Dl gradient system, including regulation of the Dl gradient. We conclude with a discussion of potential future directions.

Mechanism and implications of morphogen shuttling: Lessons learned from dorsal and Cactus in Drosophila.

In a developing animal, morphogen gradients act to pattern tissues into distinct domains of cell types. However, despite their prevalence in development, morphogen gradient formation is a matter of debate. In our recent publication, we showed that the Dorsal/NF-κB morphogen gradient, which patterns the DV axis of the early Drosophila embryo, is partially established by a mechanism of facilitated diffusion. This mechanism, also known as "shuttling," occurs when a binding partner of the morphogen facilitates the diffusion of the morphogen, allowing it to accumulate at a given site. In this case, the inhibitor Cactus/IκB facilitates the diffusion of Dorsal/NF-κB. In the fly embryo, we used computation and experiment to not only show that shuttling occurs in the embryo, but also that it enables the viability of embryos that inherit only one copy of dorsal maternally. In this commentary, we further discuss our evidence behind the shuttling mechanism, the previous literature data explained by the mechanism, and how it may also be critical for robustness of development. Finally, we briefly provide additional experimental data pointing toward an interaction between Dorsal and BMP signaling that is likely affected by shuttling.

Spatiotemporal control of gene expression boundaries using a feedforward loop.

BACKGROUND: A feedforward loop (FFL) is commonly observed in several biological networks. The FFL network motif has been mostly studied with respect to variation of the input signal in time, with only a few studies of FFL activity in a spatially distributed system such as morphogen-mediated tissue patterning. However, most morphogen gradients also evolve in time. RESULTS: We studied the spatiotemporal behavior of a coherent FFL in two contexts: (a) a generic, oscillating morphogen gradient and (b) the dorsal-ventral patterning of the early Drosophila embryo by a gradient of the NF-κB homolog dorsal with its early target Twist. In both models, we found features in the dynamics of the intermediate node-phase difference and noise filtering-that were largely independent of the parameterization of the models, and thus were functions of the structure of the FFL itself. In the dorsal gradient model, we also found that proper target gene expression was not possible without including the effect of maternal pioneer factor Zelda. CONCLUSIONS: An FFL buffers fluctuation to changes in the morphogen signal ensuring stable gene expression boundaries.

A Rubisco-binding protein is required for normal pyrenoid number and starch sheath morphology in Chlamydomonas reinhardtii.

A phase-separated, liquid-like organelle called the pyrenoid mediates CO2 fixation in the chloroplasts of nearly all eukaryotic algae. While most algae have 1 pyrenoid per chloroplast, here we describe a mutant in the model alga Chlamydomonas that has on average 10 pyrenoids per chloroplast. Characterization of the mutant leads us to propose a model where multiple pyrenoids are favored by an increase in the surface area of the starch sheath that surrounds and binds to the liquid-like pyrenoid matrix. We find that the mutant's phenotypes are due to disruption of a gene, which we call StArch Granules Abnormal 1 (SAGA1) because starch sheath granules, or plates, in mutants lacking SAGA1 are more elongated and thinner than those of wild type. SAGA1 contains a starch binding motif, suggesting that it may directly regulate starch sheath morphology. SAGA1 localizes to multiple puncta and streaks in the pyrenoid and physically interacts with the small and large subunits of the carbon-fixing enzyme Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase), a major component of the liquid-like pyrenoid matrix. Our findings suggest a biophysical mechanism by which starch sheath morphology affects pyrenoid number and CO2-concentrating mechanism function, advancing our understanding of the structure and function of this biogeochemically important organelle. More broadly, we propose that the number of phase-separated organelles can be regulated by imposing constraints on their surface area.

The engineering principles of combining a transcriptional incoherent feedforward loop with negative feedback.

BACKGROUND: Regulation of gene expression is of paramount importance in all living systems. In the past two decades, it has been discovered that certain motifs, such as the feedforward motif, are overrepresented in gene regulatory circuits. Feedforward loops are also ubiquitous in process control engineering, and are nearly always structured so that one branch has the opposite effect of the other, which is a structure known as an "incoherent" feedforward loop in biology. In engineered systems, feedforward control loops are subject to several engineering constraints, including that (1) they are finely-tuned so that the system returns to the original steady state after a disturbance occurs (perfect adaptation), (2) they are typically only implemented in the combination with negative feedback, and (3) they can greatly improve the stability and dynamical characteristics of the conjoined negative feedback loop. On the other hand, in biology, incoherent feedforward loops can serve many purposes, one of which may be perfect adaptation. It is an open question as to whether those that achieve perfect adaptation are subject to the above engineering principles. RESULTS: We analyzed an incoherent feedforward gene regulatory motif from the standpoint of the above engineering principles. In particular, we showed that an incoherent feedforward loop Type 1 (I1-FFL), from within a gene regulatory circuit, can be finely-tuned for perfect adaptation after a stimulus, and that the robustness of this behavior is increased by the presence of moderate negative feedback. In addition, we analyzed the advantages of adding a feedforward loop to a system that already operated under negative feedback, and found that the dynamical properties of the combined feedforward/feedback system were superior. CONCLUSIONS: Our analysis shows that many of the engineering principles used in engineering design of feedforward control are also applicable to feedforward loops in biological systems. We speculate that principles found in other domains of engineering may also be applicable to analogous structures in biology.

Natural Variation within a Species for Traits Underpinning C4 Photosynthesis.

Engineering C4 photosynthesis into C3 crops could substantially increase their yield by alleviating photorespiratory losses. This objective is challenging because the C4 pathway involves complex modifications to the biochemistry, cell biology, and anatomy of leaves. Forward genetics has provided limited insight into the mechanistic basis of these properties, and there have been no reports of significant quantitative intraspecific variation of C4 attributes that would allow trait mapping. Here, we show that accessions of the C4 species Gynandropsis gynandra collected from locations across Africa and Asia exhibit natural variation in key characteristics of C4 photosynthesis. Variable traits include bundle sheath size and vein density, gas-exchange parameters, and carbon isotope discrimination associated with the C4 state. The abundance of transcripts encoding core enzymes of the C4 cycle also showed significant variation. Traits relating to water use showed more quantitative variation than those associated with carbon assimilation. We propose that variation in these traits likely adapted the hydraulic system for increased water use efficiency rather than improving carbon fixation, indicating that selection pressure may drive C4 diversity in G. gynandra by modifying water use rather than photosynthesis. The accessions analyzed can be easily crossed and produce fertile offspring. Our findings, therefore, indicate that natural variation within this C4 species is sufficiently large to allow genetic mapping of key C4 traits and regulators.

Ancient duons may underpin spatial patterning of gene expression in C4 leaves.

If the highly efficient C4 photosynthesis pathway could be transferred to crops with the C3 pathway there could be yield gains of up to 50%. It has been proposed that the multiple metabolic and developmental modifications associated with C4 photosynthesis are underpinned by relatively few master regulators that have allowed the evolution of C4 photosynthesis more than 60 times in flowering plants. Here we identify a component of one such regulator that consists of a pair of cis-elements located in coding sequence of multiple genes that are preferentially expressed in bundle sheath cells of C4 leaves. These motifs represent duons as they play a dual role in coding for amino acids as well as controlling the spatial patterning of gene expression associated with the C4 leaf. They act to repress transcription of C4 photosynthesis genes in mesophyll cells. These duons are also present in the C3 model Arabidopsis thaliana, and, in fact, are conserved in all land plants and even some algae that use C3 photosynthesis. C4 photosynthesis therefore appears to have coopted an ancient regulatory code to generate the spatial patterning of gene expression that is a hallmark of C4 photosynthesis. This intragenic transcriptional regulatory sequence could be exploited in the engineering of efficient photosynthesis of crops.

A facilitated diffusion mechanism establishes the Drosophila Dorsal gradient.

The transcription factor NF-κB plays an important role in the immune system, apoptosis and inflammation. Dorsal, a Drosophila homolog of NF-κB, patterns the dorsal-ventral axis in the blastoderm embryo. During this stage, Dorsal is sequestered outside the nucleus by the IκB homolog Cactus. Toll signaling on the ventral side breaks the Dorsal/Cactus complex, allowing Dorsal to enter the nucleus to regulate target genes. Fluorescent data show that Dorsal accumulates on the ventral side of the syncytial blastoderm. Here, we use modeling and experimental studies to show that this accumulation is caused by facilitated diffusion, or shuttling, of the Dorsal/Cactus complex. We also show that active Toll receptors are limiting in wild-type embryos, which is a key factor in explaining global Dorsal gradient formation. Our results suggest that shuttling is necessary for viability of embryos from mothers with compromised dorsal levels. Therefore, Cactus not only has the primary role of regulating Dorsal nuclear import, but also has a secondary role in shuttling. Given that this mechanism has been found in other, independent, systems, we suggest that it might be more prevalent than previously thought.

Regulatory gateways for cell-specific gene expression in C4 leaves with Kranz anatomy.

C4 photosynthesis is a carbon-concentrating mechanism that increases delivery of carbon dioxide to RuBisCO and as a consequence reduces photorespiration. The C4 pathway is therefore beneficial in environments that promote high photorespiration. This pathway has evolved many times, and involves restricting gene expression to either mesophyll or bundle sheath cells. Here we review the regulatory mechanisms that control cell-preferential expression of genes in the C4 cycle. From this analysis, it is clear that the C4 pathway has a complex regulatory framework, with control operating at epigenetic, transcriptional, post-transcriptional, translational, and post-translational levels. Some genes of the C4 pathway are regulated at multiple levels, and we propose that this ensures robust expression in each cell type. Accumulating evidence suggests that multiple genes of the C4 pathway may share the same regulatory mechanism. The control systems for C4 photosynthesis gene expression appear to operate in C3 plants, and so it appears that pre-existing mechanisms form the basis of C4 photosynthesis gene expression.

Analyzing negative feedback using a synthetic gene network expressed in the Drosophila melanogaster embryo.

BACKGROUND: A complex network of gene interactions controls gene regulation throughout development and the life of the organisms. Insights can be made into these processes by studying the functional interactions (or "motifs") which make up these networks. RESULTS: We sought to understand the functionality of one of these network motifs, negative feedback, in a multi-cellular system. This was accomplished using a synthetic network expressed in the Drosophila melanogaster embryo using the yeast proteins Gal4 (a transcriptional activator) and Gal80 (an inhibitor of Gal4 activity). This network is able to produce an attenuation or shuttling phenotype depending on the Gal80/Gal4 ratio. This shuttling behavior was validated by expressing Gal3, which inhibits Gal80, to produce a localized increase in free Gal4 and therefore signaling. Mathematical modeling was used to demonstrate the capacity for negative feedback to produce these varying outputs. CONCLUSIONS: The capacity of a network motif to exhibit different phenotypes due to minor changes to the network in multi-cellular systems was shown. This work demonstrates the importance of studying network motifs in multi-cellular systems.