Gene regulation by MAPK substrate competition

Fly Me to the Micron: Microtechnologies for Drosophila Research

Multicellular model organisms, such as Drosophila melanogaster (fruit fly), are frequently used in a myriad of biological research studies due to their biological significance and global standardization. However, traditional tools used in these studies generally require manual handling, subjective phenotyping, and bulk treatment of the organisms, resulting in laborious experimental protocols with limited accuracy. Advancements in microtechnology over the course of the last two decades have allowed researchers to develop automated, high-throughput, and multifunctional experimental tools that enable novel experimental paradigms that would not be possible otherwise. We discuss recent advances in microtechnological systems developed for small model organisms using D. melanogaster as an example. We critically analyze the state of the field by comparing the systems produced for different applications. Additionally, we suggest design guidelines, operational tips, and new research directions based on the technical and knowledge gaps in the literature. This review aims to foster interdisciplinary work by helping engineers to familiarize themselves with model organisms while presenting the most recent advances in microengineering strategies to biologists.

Preformation and epigenesis converge to specify primordial germ cell fate in the early Drosophila embryo

A critical step in animal development is the specification of primordial germ cells (PGCs), the precursors of the germline. Two seemingly mutually exclusive mechanisms are implemented across the animal kingdom: epigenesis and preformation. In epigenesis, PGC specification is non-autonomous and depends on extrinsic signaling pathways. The BMP pathway provides the key PGC specification signals in mammals. Preformation is autonomous and mediated by determinants localized within PGCs. In Drosophila, a classic example of preformation, constituents of the germ plasm localized at the embryonic posterior are thought to be both necessary and sufficient for proper determination of PGCs. Contrary to this longstanding model, here we show that these localized determinants are insufficient by themselves to direct PGC specification in blastoderm stage embryos. Instead, we find that the BMP signaling pathway is required at multiple steps during the specification process and functions in conjunction with components of the germ plasm to orchestrate PGC fate.

Proper specification of primordial germ cells (PGCs) is crucial as PGCs serve as the precursors of germline stem cells. To specify PGC fate, invertebrates rely upon cell autonomous preformation involving maternally deposited germ plasm. In Drosophila melanogaster, to insulate newly formed PGCs from the adverse effects of the cell-cell signaling pathways, germ plasm determinants silence transcription and attenuate the cell cycle. However, our data on the BMP signaling pathway challenge this long-held view of PGC specification and suggest that appropriate specification of embryonic PGCs is sensitive to the BMP ligand, decapentaplegic (dpp), and its cognate receptor, thickveins. We find that PGCs are not only capable of responding to BMP signals from the soma, but also that these signals impact the proper determination of the germ cells. Based on these unanticipated similarities between mammals and flies, we propose a model integrating contribution of both the cell-autonomous (preformation) and non-autonomous (epigenesis) pathways during PGC determination. Consistent with the model, we have observed dominant genetic interactions between, oskar, the maternal determinant of PGC fate, and the BMP pathway ligand dpp.

Alternative Splicing of MAPKs in the Regulation of Signaling Specificity

The mitogen-activated protein kinase (MAPK) cascades transmit signals from extracellular stimuli to a variety of distinct cellular processes. The MAPKKs in each cascade specifically phosphorylate and activate their cognate MAPKs, indicating that this step funnels various signals into a seemingly linear pathway. Still, the effects of these cascades vary significantly, depending on the identity of the extracellular signals, which gives rise to proper outcomes. Therefore, it is clear that the specificity of the signals transmitted through the cascades is tightly regulated in order to secure the desired cell fate. Indeed, many regulatory components or processes that extend the specificity of the cascades have been identified. Here, we focus on a less discussed mechanism, that is, the role of distinct components in each tier of the cascade in extending the signaling specificity. We cover the role of distinct genes, and the alternatively spliced isoforms of MAPKKs and MAPKs, in the signaling specificity. The alternatively spliced MEK1b and ERK1c, which form an independent signaling route, are used as the main example. Unlike MEK1/2 and ERK1/2, this route’s functions are limited, including mainly the regulation of mitotic Golgi fragmentation. The unique roles of the alternatively spliced isoforms indicate that these components play an essential role in determining the proper cell fate in response to distinct stimulations.

Competitive effects in bacterial mRNA decay

In living organisms, the same enzyme catalyses the degradation of thousands of different mRNAs, but the possible influence of competing substrates has been largely ignored so far. We develop a simple mechanistic model of the coupled degradation of all cell mRNAs using the total quasi-steady-state approximation of the Michaelis-Menten framework. Numerical simulations of the model using carefully chosen parameters and analyses of rate sensitivity coefficients show how substrate competition alters mRNA decay. The model predictions reproduce and explain a number of experimental observations on mRNA decay following transcription arrest, such as delays before the onset of degradation, the occurrence of variable degradation profiles with increased non linearities and the negative correlation between mRNA half-life and concentration. The competition acts at different levels, through the initial concentration of cell mRNAs and by modifying the enzyme affinity for its targets. The consequence is a global slow down of mRNA decay due to enzyme titration and the amplification of its apparent affinity. Competition happens to stabilize weakly affine mRNAs and to destabilize the most affine ones. We believe that this mechanistic model is an interesting alternative to the exponential models commonly used for the determination of mRNA half-lives. It allows analysing regulatory mechanisms of mRNA degradation and its predictions are directly comparable to experimental data.

Kinetic analysis of substrate competition in enzymatic reactions with β-D-galactosidase by capillary electrophoresis / dynamic frontal analysis

Competitive inhibition between two substrates with an enzyme is investigated by capillary electrophoresis/dynamic frontal analysis (CE/DFA). Enzymatic hydrolyses of o-nitrophenyl β-D-galactopyranoside and p-nitrophenyl β-D-galactopyranoside with β-D-galactosidase were examined as a model competitive reaction. A sample solution containing the two substrates was injected into a capillary filled with a separation buffer containing an enzyme. Enzymatic hydrolysis occurred during the electrophoresis, and the products of o-nitrophenol and p-nitrophenol were continuously formed and resolved from the sample zone. Two-steps plateau signal was detected with the two-substrate solutions based on the difference in the effective electrophoretic mobility of o-nitrophenol and p-nitrophenol. Michaelis-Menten constants and inhibition constants were determined with the plateau heights. Usefulness of CE/DFA on competitive inhibition analysis is demonstrated in this study.

Global shape of Toll activation is determined by wntD enhancer properties

Buffering variability in morphogen distribution is essential for reproducible patterning. A theoretically proposed class of mechanisms, termed "distal pinning," achieves robustness by combining local sensing of morphogen levels with global modulation of gradient spread. Here, we demonstrate a critical role for morphogen sensing by a gene enhancer, which ultimately determines the final global distribution of the morphogen and enables reproducible patterning. Specifically, we show that, while the pattern of Toll activation in the early Drosophila embryo is robust to gene dosage of its locally produced regulator, WntD, it is sensitive to a single-nucleotide change in the wntD enhancer. Thus, enhancer properties of locally produced WntD directly impinge on the global morphogen profile.

Functional divergence caused by mutations in an energetic hotspot in ERK2

The most frequent extracellular signal-regulated kinase 2 (ERK2) mutation occurring in cancers is E322K (E-K). ERK2 E-K reverses a buried charge in the ERK2 common docking (CD) site, a region that binds activators, inhibitors, and substrates. Little is known about the cellular consequences associated with this mutation, other than apparent increases in tumor resistance to pathway inhibitors. ERK2 E-K, like the mutation of the preceding aspartate (ERK2 D321N [D-N]) known as the sevenmaker mutation, causes increased activity in cells and evades inactivation by dual-specificity phosphatases. As opposed to findings in cancer cells, in developmental assays in Drosophila, only ERK2 D-N displays a significant gain of function, revealing mutation-specific phenotypes. The crystal structure of ERK2 D-N is indistinguishable from that of wild-type protein, yet this mutant displays increased thermal stability. In contrast, the crystal structure of ERK2 E-K reveals profound structural changes, including disorder in the CD site and exposure of the activation loop phosphorylation sites, which likely account for the decreased thermal stability of the protein. These contiguous mutations in the CD site of ERK2 are both required for docking interactions but lead to unpredictably different functional outcomes. Our results suggest that the CD site is in an energetically strained configuration, and this helps drive conformational changes at distal sites on ERK2 during docking interactions.

The Nuclear Translocation of Mitogen-Activated Protein Kinases: Molecular Mechanisms and Use as Novel Therapeutic Target

The mitogen-activated protein kinase (MAPK) cascades are central signaling pathways that play a central role in the regulation of most stimulated cellular processes including proliferation, differentiation, stress response and apoptosis. Currently 4 such cascades are known, each termed by its downstream MAPK components: the extracellular signal-regulated kinase 1/2 (ERK1/2), cJun-N-terminal kinase (JNK), p38 and ERK5. One of the hallmarks of these cascades is the stimulated nuclear translocation of their MAPK components using distinct mechanisms. ERK1/2 are shuttled into the nucleus by importin7, JNK and p38 by a dimer of importin3 with either importin9 or importin7, and ERK5 by importin-α/β. Dysregulation of these cascades often results in diseases, including cancer and inflammation, as well as developmental and neurological disorders. Much effort has been invested over the years in developing inhibitors to the MAPK cascades to combat these diseases. Although some inhibitors are already in clinical use or clinical trials, their effects are hampered by development of resistance or adverse side-effects. Recently, our group developed 2 myristoylated peptides: EPE peptide, which inhibits the interaction of ERK1/2 with importin7, and PERY peptide, which prevents JNK/p38 interaction with either importin7 or importin9. These peptides block the nuclear translocation of their corresponding kinases, resulting in prevention of several cancers, while the PERY peptide also inhibits inflammation-induced diseases. These peptides provide a proof of concept for the use of the nuclear translocation of MAPKs as therapeutic targets for cancer and/or inflammation.

Two Drosophilids exhibit distinct EGF pathway patterns in oogenesis

Deciphering the evolution of morphological structures is a remaining challenge in the field of developmental biology. The respiratory structures of insect eggshells, called the dorsal appendages, provide an outstanding system for exploring these processes since considerable information is known about their patterning and morphogenesis in Drosophila melanogaster and dorsal appendage number and morphology vary widely across Drosophilid species. We investigated the patterning differences that might facilitate morphogenetic differences between D. melanogaster, which produces two oar-like structures first by wrapping and then elongating the tubes via cell intercalation and cell crawling, and Scaptodrosophila lebanonensis, which produces a variable number of appendages simply by cell intercalation and crawling. Analyses of BMP pathway components thickveins and P-Mad demonstrate that anterior patterning is conserved between these species. In contrast, EGF signaling exhibits significant differences. Transcripts for the ligand encoded by gurken localize similarly in the two species, but this morphogen creates a single dorsolateral primordium in S. lebanonensis as defined by activated MAP kinase and the downstream marker broad. Expression patterns of pointed, argos, and Capicua, early steps in the EGF pathway, exhibit a heterochronic shift in S. lebanonensis relative to those seen in D. melanogaster. We demonstrate that the S. lebanonensis Gurken homolog is active in D. melanogaster but is insufficient to alter downstream patterning responses, indicating that Gurken-EGF receptor interactions do not distinguish the two species' patterning. Altogether, these results differentiate EGF signaling patterns between species and shed light on how changes to the regulation of patterning genes may contribute to different tube-forming mechanisms.

Analysis of network motifs in cellular regulation: Structural similarities, input-output relations and signal integration

Much of the complexity of regulatory networks derives from the necessity to integrate multiple signals and to avoid malfunction due to cross-talk or harmful perturbations. Hence, one may expect that the input-output behavior of larger networks is not necessarily more complex than that of smaller network motifs which suggests that both can, under certain conditions, be described by similar equations. In this review, we illustrate this approach by discussing the similarities that exist in the steady state descriptions of a simple bimolecular reaction, covalent modification cycles and bacterial two-component systems. Interestingly, in all three systems fundamental input-output characteristics such as thresholds, ultrasensitivity or concentration robustness are described by structurally similar equations. Depending on the system the meaning of the parameters can differ ranging from protein concentrations and affinity constants to complex parameter combinations which allows for a quantitative understanding of signal integration in these systems. We argue that this approach may also be extended to larger regulatory networks.

A double-edged sword: The world according to Capicua in cancer

CIC/Capicua is an HMG‐box transcription factor that is well conserved during evolution. CIC recognizes the T(G/C)AATG(A/G)A sequence and represses its target genes, such as PEA3 family genes. The receptor tyrosine kinase/RAS/MAPK signals downregulate CIC and relieves CIC's target genes from the transrepressional activity; CIC thus acts as an important downstream molecule of the pathway and as a tumor suppressor. CIC loss‐of‐function mutations are frequently observed in several human neoplasms such as oligodendroglioma, and lung and gastric carcinoma. CIC is also involved in chromosomal translocation‐associated gene fusions in highly aggressive small round cell sarcoma that is biologically and clinically distinct from Ewing sarcoma. In these mutations, PEA3 family genes and other important target genes are upregulated, inducing malignant phenotypes. Downregulation of CIC abrogates the effect of MAPK inhibitors, suggesting its potential role as an important modifier of molecular target therapies for cancer. These data reveal the importance of CIC as a key molecule in signal transduction, carcinogenesis, and developing novel therapies.

Two faces of competition: target-mediated reverse signalling in microRNA and mitogen-activated protein kinase regulatory networks

Biomolecular regulatory networks are organised around hubs, which can interact with a large number of targets. These targets compete with each other for access to their common hubs, but whether the effect of this competition would be significant in magnitude and in function is not clear. In this review, the authors discuss recent in vivo studies that analysed the system level retroactive effects induced by target competition in microRNA and mitogen-activated protein kinase regulatory networks. The results of these studies suggest that downstream targets can regulate the overall state of their upstream regulators, and thus cannot be ignored in analysing biomolecular networks.

Mitogen-activated protein kinases participate in determination of apical-basal symmetry in Pisum sativum

Mitogen-activated protein kinases (MAPKs) are implicated in various processes in plants. Apart from response to biotic and abiotic stresses they are involved in regulation of embryo development. Although MAPKs were found to be indispensable during embryo development it has never been established at which stages of embryogenesis and in which regions of a plant embryo activated MAPKs can be observed. Here, we show that apical and basal regions display activation of the same types of MAPKs and the only difference concerns the level of their phosphorylation and cellular localization. Dually-phosphorylated MAPKs were found in nuclei of the apical region of an embryo both at the early and late cotyledonary stage while no immunofluorescence signals were detected in nuclei of the basal region. However, in this case phosphorylated MAPKs were immunodetected in cytoplasm in the apical domain of cortex cells, indicating their role in auxin transport from the basal to apical region of an embryo. Furthermore, obtained data indicate that nuclear localization of activated MAPKs may result from epigenetic modifications and polar auxin transport. The presented data and previous studies lead to the conclusion that activated MAPKs and their cellular localization define apical and basal regions during formation of an apical-basal axis.

Transcriptome sequencing reveals maelstrom as a novel target gene of the terminal-system in the red flour beetle Tribolium castaneum

Terminal regions of the Drosophila embryo are patterned by the localized activation of the Torso-RTK pathway, which promotes the down-regulation of Capicua. In the short-germ beetle Tribolium, the function of the terminal system appears to be rather different, as the pathway promotes axis elongation and in addition, is required for patterning the extraembryonic serosa at the anterior. Here we show that Torso signalling induces gene expression by relieving CAPICUA-mediated repression also in Tribolium. Given that the majority of Torso target genes remain to be identified, we established a differential gene-expression screen. A subset of 50 putative terminal target genes was screened for functions in early embryonic patterning. Of those, 13 genes show early terminal expression domains and also phenotypes were related to terminal patterning. Among others, we found the PIWI-interacting RNA factor Maelstrom to be crucial for early embryonic polarization. Tc-mael is required for proper serosal size regulation and head morphogenesis. Moreover, Tc-mael promotes growth-zone formation and axis elongation. Our results suggest that posterior patterning by Torso may be realized through Maelstrom depended activation of posterior wnt-domains.

Cell signaling regulation by protein phosphorylation: a multivariate, heterogeneous, and context-dependent process

Proper spatiotemporal regulation of protein phosphorylation in cells and tissues is required for normal development and homeostasis, but aberrant protein phosphorylation regulation leads to various diseases. The study of signaling regulation by protein phosphorylation is complicated in part by the sheer scope of the kinome and phosphoproteome, dependence of signaling protein functionality on cellular localization, and the complex multivariate relationships that exist between protein phosphorylation dynamics and the cellular phenotypes they control. Additional complexities arise from the ability of microenvironmental factors to influence phosphorylation-dependent signaling and from the tendency for some signaling processes to occur heterogeneously among cells. These considerations should be taken into account when measuring cell signaling regulation by protein phosphorylation.

A WntD-Dependent Integral Feedback Loop Attenuates Variability in Drosophila Toll Signaling

Patterning by morphogen gradients relies on the capacity to generate reproducible distribution profiles. Morphogen spread depends on kinetic parameters, including diffusion and degradation rates, which vary between embryos, raising the question of how variability is controlled. We examined this in the context of Toll-dependent dorsoventral (DV) patterning of the Drosophila embryo. We find that low embryo-to-embryo variability in DV patterning relies on wntD, a Toll-target gene expressed initially at the posterior pole. WntD protein is secreted and disperses in the extracellular milieu, associates with its receptor Frizzled4, and inhibits the Toll pathway by blocking the Toll extracellular domain. Mathematical modeling predicts that WntD accumulates until the Toll gradient narrows to its desired spread, and we support this feedback experimentally. This circuit exemplifies a broadly applicable induction-contraction mechanism, which reduces patterning variability through a restricted morphogen-dependent expression of a secreted diffusible inhibitor.

Sjögren's syndrome-associated microRNAs in CD14(+) monocytes unveils targeted TGFβ signaling

Background

Sjögren’s syndrome (SjS) monocytes have a pro-inflammatory phenotype, which may influence SjS pathogenesis. MicroRNAs (miRNAs) are small endogenously expressed molecules that can inhibit protein expression of their targeted genes and have important functions in regulating cell signaling responses. We profiled miRNAs in SjS monocytes to identify a SjS-specific miRNA profile and determine the potential roles of miRNAs in SjS pathogenesis.

Methods

Total RNA was extracted from healthy control (HC, n = 10), SjS (n = 18), systemic lupus erythematosus (SLE, n = 10), and rheumatoid arthritis (RA, n = 10) peripheral blood CD14⁺ monocytes for miRNA microarray analysis. To validate select miRNAs from the microarray analysis, the original cohort and a new cohort of monocyte RNA samples from HC (n = 9), SjS (n = 12), SLE (n = 8), and RA (n = 9) patients were evaluated by quantitative reverse transcription (RT)-PCR. Functional predictions of differentially expressed miRNAs were determined through miRNA target prediction database analyses. Statistical analyses performed included one-way analysis of variance with Bonferroni post tests, linear regression, and receiver operating characteristic curve analyses.

Results

MiRNAs were predominantly upregulated in SjS monocytes in comparison with controls. Quantitative RT-PCR confirmations supported co-regulation of miR-34b-3p, miR-4701-5p, miR-609, miR-300, miR-3162-3p, and miR-877-3p in SjS monocytes (13/30, 43.3 %) in comparison with SLE (1/17, 5.8 %) and RA (1/18, 5.6 %). MiRNA-target pathway predictions identified SjS-associated miRNAs appear to preferentially target the canonical TGFβ signaling pathway as opposed to pro-inflammatory interleukin-12 and Toll-like receptor/NFkB pathways.

Conclusions

Our results underscore a novel underlying molecular mechanism where SjS-associated miRNAs may collectively suppress TGFβ signaling as opposed to pro-inflammatory interleukin-12 and Toll-like receptor/NFκB pathways in SjS pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-016-0987-0) contains supplementary material, which is available to authorized users.

Microfluidics for High-Throughput Quantitative Studies of Early Development

Developmental biology has traditionally relied on qualitative analyses; recently, however, as in other fields of biology, researchers have become increasingly interested in acquiring quantitative knowledge about embryogenesis. Advances in fluorescence microscopy are enabling high-content imaging in live specimens. At the same time, microfluidics and automation technologies are increasing experimental throughput for studies of multicellular models of development. Furthermore, computer vision methods for processing and analyzing bioimage data are now leading the way toward quantitative biology. Here, we review advances in the areas of fluorescence microscopy, microfluidics, and data analysis that are instrumental to performing high-content, high-throughput studies in biology and specifically in development. We discuss a case study of how these techniques have allowed quantitative analysis and modeling of pattern formation in the Drosophila embryo.

An integrated platform for large-scale data collection and precise perturbation of live Drosophila embryos

Understanding the fundamental principles governing embryogenesis is a key goal of developmental biology. Direct observation of embryogenesis via in vivo live imaging is vital to understanding embryogenesis; yet, tedious sample preparation makes it difficult to acquire large-scale imaging data that is often required to overcome experimental and biological noises for quantitative studies. Furthermore, it is often difficult, and sometimes impossible, to incorporate environmental perturbation for understanding developmental responses to external stimuli. To address this issue, we have developed a method for high-throughput imaging of live embryos, delivering precise environmental perturbations, and unbiased data extraction. This platform includes an optimized microfluidic device specifically for live embryos and also for precise perturbations in the microenvironment of the developing embryos. In addition, we developed software for simple, yet accurate, automated segmentation of fluorescent images, and automated data extraction. Using a quantitative assessment we find that embryos develop normally within the microfluidic device. Finally, we show an application of the high-throughput assay for monitoring developmental responses to external stimuli: anoxia-induced developmental arrest in Drosophila embryos. With slight modifications, the method developed in this work can be applied to many other models of development and other stimulus-response behaviors during development.

Mapping the binding interface of ERK and transcriptional repressor Capicua using photocrosslinking

Extracellular signal-regulated kinase (ERK) coordinates cellular responses to a range of stimuli by phosphorylating its numerous substrates. One of these substrates, Capicua (Cic), is a transcriptional repressor that was first identified in Drosophila and has been implicated in a number of human diseases. Here we use a chemical biology approach to map the binding interface of ERK and Cic. The noncanonical amino acid p-azidophenylalanine (AzF) was introduced into the ERK-binding region of Drosophila Cic, and photocrosslinking and tandem mass spectrometry were used to pinpoint its binding site on ERK. We also identified the ERK-binding region of human Cic and showed that it binds to the same site on ERK despite lacking conservation with the Drosophila Cic binding region. Finally, we mapped the amino acids involved in human Cic binding to ERK using AzF-labeled ERK. These results reveal the molecular details of the ERK-Cic interaction and demonstrate that the photocrosslinking approach is complementary to existing methods for mapping kinase-substrate binding interfaces.

Imaging the dorsal-ventral axis of live and fixed Drosophila melanogaster embryos

Optimal imaging conditions are of critical importance in developmental biology, as much of the data in the discipline is acquired through microscopy. However, imaging deep sections of tissue, especially live tissue, can be a technical challenge due to light scattering and difficulties in mounting the sample. In particular, capturing high-quality images of dorsal-ventral cross sections requires "end-on" mounting to orient the anterior-posterior axis vertically. Here we present methods to mount and image dorsal-ventral cross sections of both live and fixed Drosophila melanogaster embryos. Our methods have the advantages of being rapid, allowing deep optical sections, and not requiring expensive, specialized equipment.

A load driver device for engineering modularity in biological networks

The behavior of gene modules in complex synthetic circuits is often unpredictable^1–4. Upon joining modules to create a circuit, downstream elements (such as binding sites for a regulatory protein) apply a load to upstream modules that can negatively affect circuit function^1,5. Here we devise a genetic device named a load driver that mitigates the impact of load on circuit function, and we demonstrate its behavior in Saccharomyces cerevisiae. The load driver implements the design principle of time scale separation: inclusion of the load driver’s fast phosphotransfer processes restores the capability of a slower transcriptional circuit to respond to time-varying input signals even in the presence of substantial load. Without the load driver, we observe circuit behavior that suffers from 76% delay in response time and a 25% decrease in system bandwidth due to load. With the addition of a load driver, circuit performance is almost completely restored. Load drivers will serve as fundamental building blocks in the creation of complex, higher level genetic circuits.

ERK as a model for systems biology of enzyme kinetics in cells

A key step towards a chemical picture of enzyme catalysis was taken in 1913, when Leonor Michaelis and Maud Menten published their studies of sucrose hydrolysis by invertase. Based on a novel experimental design and a mathematical model, their work offered a quantitative view of biochemical kinetics well before the protein nature of enzymes was established and complexes with substrates could be detected. Michaelis-Menten kinetics provides a solid framework for enzyme kinetics in vitro, but what about kinetics in cells, where enzymes can be highly regulated and participate in a multitude of interactions? We discuss this question using the Extracellular Signal Regulated Kinase (ERK), which controls a myriad functions in cells, as a model of an important enzyme for which we have crystal structures, quantitative in vitro assays, and a vast list of binding partners. Despite great progress, we still cannot quantitatively predict how the rates of ERK-dependent reactions respond to genetic and pharmacological perturbations. Achieving this goal, which is important from both fundamental and practical standpoints, requires measuring the rates of enzyme reactions in their native environment and interpreting these measurements using simple but realistic mathematical models--the two elements which served as the cornerstones for Michaelis' and Menten's seminal 1913 paper.

Expression pattern and first functional characterization of riok-1 in Caenorhabditis elegans

Rio kinases are atypical serine/threonine kinases that emerge as potential cooperation partners in Ras-driven tumors. In the current study, we performed an RNAi screen in Caenorhabditis elegans to identify suppressors of oncogenic Ras signaling. Aberrant Ras/Raf signaling in C. elegans leads to the formation of a multi-vulva (Muv) phenotype. We found that depletion of riok-1, the C. elegans orthologue of the mammalian RioK1, suppressed the Muv phenotype. By using a promoter GFP construct, we could show that riok-1 is expressed in neuronal cells, the somatic gonad, the vulva, the uterus and the spermatheca. Furthermore, we observed developmental defects in the gonad upon riok-1 knockdown in a wildtype background. Our data suggest that riok-1 is a modulator of the Ras signaling pathway, suggesting implications for novel interventions in the context of Ras-driven tumors.

Downregulation of CIC does not associate with overexpression of ETV1 or MAP kinase pathway activation in gastrointestinal stromal tumors

ETV1 is a key factor in gastrointestinal stromal tumors (GIST), and is promoted by CIC downregulation in melanoma. We investigated CIC, ETV1, and the MAPK pathway in GIST. Downregulation of CIC protein levels as assessed by immunostaining was seen in 17/144 GIST, but was not associated with ETV1 or pMEK1/2 expression, KIT and PDGFRA mutations, copy number variations (CNV) of 19q13, and clinical factors. However, the data indicate that the incidence of CIC downregulation may differ for GISTs in different locations in the gastrointestinal tract, and that CNV of 19q13 is associated with shorter disease-free survival.

Competition in notch signaling with cis enriches cell fate decisions

Notch signaling is involved in cell fate choices during the embryonic development of Metazoa. Commonly, Notch signaling arises from the binding of the Notch receptor to its ligands in adjacent cells driving cell-to-cell communication. Yet, cell-autonomous control of Notch signaling through both ligand-dependent and ligand-independent mechanisms is known to occur as well. Examples include Notch signaling arising in the absence of ligand binding, and cis-inhibition of Notch signaling by titration of the Notch receptor upon binding to its ligands within a single cell. Increasing experimental evidences support that the binding of the Notch receptor with its ligands within a cell (cis-interactions) can also trigger a cell-autonomous Notch signal (cis-signaling), whose potential effects on cell fate decisions and patterning remain poorly understood. To address this question, herein we mathematically and computationally investigate the cell states arising from the combination of cis-signaling with additional Notch signaling sources, which are either cell-autonomous or involve cell-to-cell communication. Our study shows that cis-signaling can switch from driving cis-activation to effectively perform cis-inhibition and identifies under which conditions this switch occurs. This switch relies on the competition between Notch signaling sources, which share the same receptor but differ in their signaling efficiency. We propose that the role of cis-interactions and their signaling on fine-grained patterning and cell fate decisions is dependent on whether they drive cis-inhibition or cis-activation, which could be controlled during development. Specifically, cis-inhibition and not cis-activation facilitates patterning and enriches it by modulating the ratio of cells in the high-ligand expression state, by enabling additional periodic patterns like stripes and by allowing localized patterning highly sensitive to the precursor state and cell-autonomous bistability. Our study exemplifies the complexity of regulations when multiple signaling sources share the same receptor and provides the tools for their characterization.

Data-driven modeling reconciles kinetics of ERK phosphorylation, localization, and activity states

The extracellular signal‐regulated kinase (ERK) signaling pathway controls cell proliferation and differentiation in metazoans. Two hallmarks of its dynamics are adaptation of ERK phosphorylation, which has been linked to negative feedback, and nucleocytoplasmic shuttling, which allows active ERK to phosphorylate protein substrates in the nucleus and cytosol. To integrate these complex features, we acquired quantitative biochemical and live‐cell microscopy data to reconcile phosphorylation, localization, and activity states of ERK. While maximal growth factor stimulation elicits transient ERK phosphorylation and nuclear translocation responses, ERK activities available to phosphorylate substrates in the cytosol and nuclei show relatively little or no adaptation. Free ERK activity in the nucleus temporally lags the peak in nuclear translocation, indicating a slow process. Additional experiments, guided by kinetic modeling, show that this process is consistent with ERK's modification of and release from nuclear substrate anchors. Thus, adaptation of whole‐cell ERK phosphorylation is a by‐product of transient protection from phosphatases. Consistent with this interpretation, predictions concerning the dose‐dependence of the pathway response and its interruption by inhibition of MEK were experimentally confirmed.

A comparative study of Pointed and Yan expression reveals new complexity to the transcriptional networks downstream of receptor tyrosine kinase signaling

The biochemical regulatory network downstream of receptor tyrosine kinase (RTK) signaling is controlled by two opposing ETS family members: the transcriptional activator Pointed (Pnt) and the transcriptional repressor Yan. A bistable switch model has been invoked to explain how pathway activation can drive differentiation by shifting the system from a high-Yan/low-Pnt activity state to a low-Yan/high-Pnt activity state. Although the model explains yan and pnt loss-of-function phenotypes in several different cell types, how Yan and Pointed protein expression dynamics contribute to these and other developmental transitions remains poorly understood. Toward this goal we have used a functional GFP-tagged Pnt transgene (Pnt-GFP) to perform a comparative study of Yan and Pnt protein expression throughout Drosophila development. Consistent with the prevailing model of the Pnt-Yan network, we found numerous instances where Pnt-GFP and Yan adopt a mutually exclusive pattern of expression. However we also observed many examples of co-expression. While some co-expression occurred in cells where RTK signaling is presumed low, other co-expression occurred in cells with high RTK signaling. The instances of co-expressed Yan and Pnt-GFP in tissues with high RTK signaling cannot be explained by the current model, and thus they provide important contexts for future investigation of how context-specific differences in RTK signaling, network topology, or responsiveness to other signaling inputs, affect the transcriptional response.

Optimized RNA ISH, RNA FISH and protein-RNA double labeling (IF/FISH) in Drosophila ovaries

In situ hybridization (ISH) is a powerful technique for detecting nucleic acids in cells and tissues. Here we describe three ISH procedures that are optimized for Drosophila ovaries: whole-mount, digoxigenin-labeled RNA ISH; RNA fluorescent ISH (FISH); and protein immunofluorescence (IF)-RNA FISH double labeling (IF/FISH). Each procedure balances conflicting requirements for permeabilization, fixation and preservation of antigenicity to detect RNA and protein expression with high resolution and sensitivity. The ISH protocol uses alkaline phosphatase-conjugated digoxigenin antibodies followed by a color reaction, whereas FISH detection involves tyramide signal amplification (TSA). To simultaneously preserve antigens for protein detection and enable RNA probe penetration for IF/FISH, we perform IF before FISH and use xylenes and detergents to permeabilize the tissue rather than proteinase K, which can damage the antigens. ISH and FISH take 3 d to perform, whereas IF/FISH takes 5 d. Probe generation takes 1 or 2 d to perform.

Calpain A modulates Toll responses by limited Cactus/IκB proteolysis

Calcium-dependent cysteine proteases of the calpain family are modulatory proteases that cleave their substrates in a limited manner. Among their substrates, calpains target vertebrate and invertebrate IκB proteins. Because proteolysis by calpains potentially generates novel protein functions, it is important to understand how this affects NFκB activity. We investigate the action of Calpain A (CalpA) on the Drosophila melanogaster IκB homologue Cactus in vivo. CalpA alters the absolute amounts of Cactus protein. Our data indicate, however, that CalpA uses additional mechanisms to regulate NFκB function. We provide evidence that CalpA interacts physically with Cactus, recognizing a Cactus pool that is not bound to Dorsal, a fly NFκB/Rel homologue. We show that proteolytic cleavage by CalpA generates Cactus fragments lacking an N-terminal region required for Toll responsiveness. These fragments are generated in vivo and display properties distinct from those of full-length Cactus. We propose that CalpA targets free Cactus, which is incorporated into and modulates Toll-responsive complexes in the embryo and immune system.

Effect of substrate competition in kinetic models of metabolic networks

Substrate competition can be found in many types of biological processes, ranging from gene expression to signal transduction and metabolic pathways. Although several experimental and in silico studies have shown the impact of substrate competition on these processes, it is still often neglected, especially in modelling approaches. Using toy models that exemplify different metabolic pathway scenarios, we show that substrate competition can influence the dynamics and the steady state concentrations of a metabolic pathway. We have additionally derived rate laws for substrate competition in reversible reactions and summarise existing rate laws for substrate competition in irreversible reactions.

Kinetics of gene derepression by ERK signaling

ERK controls gene expression in development, but mechanisms that link ERK activation to changes in transcription are not well understood. We used high-resolution analysis of signaling dynamics to study transcriptional interpretation of ERK signaling during Drosophila embryogenesis, at a stage when ERK induces transcription of intermediate neuroblasts defective (ind), a gene essential for patterning of the nerve cord. ERK induces ind by antagonizing its repression by Capicua (Cic), a transcription factor that acts as a sensor of receptor tyrosine kinases in animal development and human diseases. A recent study established that active ERK reduces the nuclear levels of Cic, but it remained unclear whether this is required for the induction of Cic target genes. We provide evidence that Cic binding sites within the regulatory DNA of ind control the spatial extent and the timing of ind expression. At the same time, we demonstrate that ERK induces ind before Cic levels in the nucleus are reduced. Based on this, we propose that ERK-dependent relief of gene repression by Cic is a two-step process, in which fast reduction of repressor activity is followed by slower changes in nuclear localization and overall protein levels. This may be a common feature of systems in which ERK induces genes by relief of transcriptional repression.

Context-dependent transcriptional interpretation of mitogen activated protein kinase signaling in the Drosophila embryo

Terminal regions of the Drosophila embryo are patterned by the localized activation of Mitogen Activated Protein Kinase (MAPK), which induces zygotic genes through relief of their repression by transcriptional repressor Capicua. The levels of MAPK activation at the anterior and posterior termini are close to each other, but the expression patterns of MAPK-target genes, such as zerknüllt (zen) and tailless (tll), display strong anterior-posterior (AP) asymmetry. This region-specific response to MAPK activation provides a clear example of context-dependent interpretation of inductive signaling, a common developmental effect that remains poorly understood. In the past, the AP asymmetry of zen expression was attributed to a mechanism that depends on MAPK substrate competition. We present data suggesting that the asymmetric expression of tll is generated by a different mechanism, based on feedforward control and multiple enhancers of the tll gene. A simple mathematical model of this mechanism correctly predicts how the wild-type expression pattern of tll changes in mutants affecting the anterior, dorsoventral, and terminal patterning systems and some of their direct targets.

Microfluidic trap array for massively parallel imaging of Drosophila embryos

Here we describe a protocol for the fabrication and use of a microfluidic device to rapidly orient >700 Drosophila embryos in parallel for end-on imaging. The protocol describes master microfabrication (∼1 d), polydimethylsiloxane molding (few hours), system setup and device operation (few minutes) and imaging (depending on application). Our microfluidics-based approach described here is one of the first to facilitate rapid orientation for end-on imaging, and it is a major breakthrough for quantitative studies on Drosophila embryogenesis. The operating principle of the embryo trap is based on passive hydrodynamics, and it does not require direct manipulation of embryos by the user; biologists following the protocol should be able to repeat these procedures. The compact design and fabrication materials used allow the device to be used with traditional microscopy setups and do not require specialized fixtures. Furthermore, with slight modification, this array can be applied to the handling of other model organisms and oblong objects.

The Capicua repressor--a general sensor of RTK signaling in development and disease

Receptor tyrosine kinase (RTK) signaling pathways control multiple cellular decisions in metazoans, often by regulating the expression of downstream genes. In Drosophila melanogaster and other systems, E-twenty-six (ETS) transcription factors are considered to be the predominant nuclear effectors of RTK pathways. Here, we highlight recent progress in identifying the HMG-box protein Capicua (CIC) as a key sensor of RTK signaling in both Drosophila and mammals. Several studies have shown that CIC functions as a repressor of RTK-responsive genes, keeping them silent in the absence of signaling. Following the activation of RTK signaling, CIC repression is relieved, and this allows the expression of the targeted gene in response to local or ubiquitous activators. This regulatory switch is essential for several RTK responses in Drosophila, from the determination of cell fate to cell proliferation. Furthermore, increasing evidence supports the notion that this mechanism is conserved in mammals, where CIC has been implicated in cancer and neurodegeneration. In addition to summarizing our current knowledge on CIC, we also discuss the implications of these findings for our understanding of RTK signaling specificity in different biological processes.

Lack of tailless leads to an increase in expression variability in Drosophila embryos

Developmental processes are robust, or canalised: dynamic patterns of gene expression across space and time are regulated reliably and precisely in the presence of genetic and environmental perturbations. It remains unclear whether canalisation relies on specific regulatory factors (such as heat-shock proteins), or whether it is based on more general redundancy and distributed robustness at the network level. The latter explanation implies that mutations in many regulatory factors should exhibit loss of canalisation. Here, we present a quantitative characterisation of segmentation gene expression patterns in mutants of the terminal gap gene tailless (tll) in Drosophila melanogaster. Our analysis provides new insights into the dynamic mechanisms underlying gap gene regulation, and reveals significantly increased variability of gene expression in the mutant compared to the wild-type background. We show that both position and timing of posterior segmentation gene expression domains vary strongly from embryo-to-embryo in tll mutants. This variability must be caused by a vulnerability in the regulatory system which is hidden or buffered in the wild-type, but becomes uncovered by the deletion of tll. Our analysis provides evidence that loss of canalisation in mutants could be more widespread than previously thought.

Dissecting sources of quantitative gene expression pattern divergence between Drosophila species

Gene expression patterns can diverge between species due to changes in a gene's regulatory DNA or changes in the proteins, e.g., transcription factors (TFs), that regulate the gene. We developed a modeling framework to uncover the sources of expression differences in blastoderm embryos of three Drosophila species, focusing on the regulatory circuit controlling expression of the hunchback (hb) posterior stripe. Using this framework and cellular-resolution expression measurements of hb and its regulating TFs, we found that changes in the expression patterns of hb's TFs account for much of the expression divergence. We confirmed our predictions using transgenic D. melanogaster lines, which demonstrate that this set of orthologous cis-regulatory elements (CREs) direct similar, but not identical, expression patterns. We related expression pattern differences to sequence changes in the CRE using a calculation of the CRE's TF binding site content. By applying this calculation in both the transgenic and endogenous contexts, we found that changes in binding site content affect sensitivity to regulating TFs and that compensatory evolution may occur in circuit components other than the CRE.

Torso RTK controls Capicua degradation by changing its subcellular localization

The transcriptional repressor Capicua (Cic) controls multiple aspects of Drosophila embryogenesis and has been implicated in vertebrate development and human diseases. Receptor tyrosine kinases (RTKs) can antagonize Cic-dependent gene repression, but the mechanisms responsible for this effect are not fully understood. Based on genetic and imaging studies in the early Drosophila embryo, we found that Torso RTK signaling can increase the rate of Cic degradation by changing its subcellular localization. We propose that Cic is degraded predominantly in the cytoplasm and show that Torso reduces the stability of Cic by controlling the rates of its nucleocytoplasmic transport. This model accounts for the experimentally observed spatiotemporal dynamics of Cic in the early embryo and might explain RTK-dependent control of Cic in other developmental contexts.

Mathematical study of the role of Delta/Notch lateral inhibition during primary branching of Drosophila trachea development

A wide range of cellular developmental processes employ intercellular signaling via the Delta/Notch lateral inhibitory pathway to achieve stable spatial patterning. Recent genetic experiments have shown the importance of Delta/Notch lateral inhibition for regulating the number of tip cells in the tracheal primary branching of Drosophila. To examine the role of Delta/Notch regulation in the tip-cell selection, we analyzed a mathematical model of a simple lateral inhibitory system having input signals. Mathematical and numerical analyses revealed that the lateral inhibition did not amplify the signal difference between neighboring cells over the parameter ranges in which the spatial pattern of tip selection was realized. We also show that the number of tip cells becomes less affected by a fluctuation of the input gradient signal as the lateral inhibition becomes stronger. In addition, we demonstrate that the lateral inhibitory regulation enhances the robustness of the tip-cell selection compared with a system regulated by self-inhibition, an alternative means of inhibitory regulation. These results suggest that the lateral inhibition promotes the robustness of tip-cell selection in the tracheal development of Drosophila.

MAP kinase signalling cascades and transcriptional regulation

The MAP kinase (MAPK) signalling pathways play fundamental roles in a wide range of cellular processes and are often deregulated in disease states. One major mode of action for these pathways is in controlling gene expression, in particular through regulating transcription. In this review, we discuss recent significant advances in this area. In particular we focus on the mechanisms by which MAPKs are targeted to the nucleus and chromatin, and once there, how they impact on chromatin structure and subsequent gene regulation. We also discuss how systems biology approaches have contributed to our understanding of MAPK signaling networks, and also how the MAPK pathways intersect with other regulatory pathways in the nucleus. Finally, we summarise progress in studying the physiological functions of key MAPK transcriptional targets.

RTK signaling modulates the Dorsal gradient

The dorsoventral (DV) axis of the Drosophila embryo is patterned by a nuclear gradient of the Rel family transcription factor, Dorsal (Dl), that activates or represses numerous target genes in a region-specific manner. Here, we demonstrate that signaling by receptor tyrosine kinases (RTK) reduces nuclear levels and transcriptional activity of Dl, both at the poles and in the mid-body of the embryo. These effects depend on wntD, which encodes a Dl antagonist belonging to the Wingless/Wnt family of secreted factors. Specifically, we show that, via relief of Groucho- and Capicua-mediated repression, the Torso and EGFR RTK pathways induce expression of WntD, which in turn limits Dl nuclear localization at the poles and along the DV axis. Furthermore, this RTK-dependent control of Dl is important for restricting expression of its targets in both contexts. Thus, our results reveal a new mechanism of crosstalk, whereby RTK signals modulate the spatial distribution and activity of a developmental morphogen in vivo.

A system of repressor gradients spatially organizes the boundaries of Bicoid-dependent target genes

The homeodomain (HD) protein Bicoid (Bcd) is thought to function as a gradient morphogen that positions boundaries of target genes via threshold-dependent activation mechanisms. Here, we analyze 66 Bcd-dependent regulatory elements and show that their boundaries are positioned primarily by repressive gradients that antagonize Bcd-mediated activation. A major repressor is the pair-rule protein Runt (Run), which is expressed in an opposing gradient and is necessary and sufficient for limiting Bcd-dependent activation. Evidence is presented that Run functions with the maternal repressor Capicua and the gap protein Kruppel as the principal components of a repression system that correctly orders boundaries throughout the anterior half of the embryo. These results put conceptual limits on the Bcd morphogen hypothesis and demonstrate how the Bcd gradient functions within the gene network that patterns the embryo.

CIC and FUBP1 mutations in oligodendrogliomas, oligoastrocytomas and astrocytomas

CIC and FUBP1 mutations have recently been detected in oligodendrogliomas but not in oligoastrocytomas. However, allelic losses in the regions on chromosomal arms 19q and 1p harboring CIC and FUBP1 are a common feature of both, oligodendrogliomas and oligoastrocytomas. To resolve this discrepancy, we analyzed CIC and FUBP1 mutations in a set of primary brain tumors including 18 oligodendrogliomas and 42 oligoastrocytomas. In addition, we analyzed 10 astrocytomas and 16 glioblastomas with allelic losses on 19q as well as a set of 12 medulloblastomas for CIC mutations. CIC mutations were found in 15/18 oligodendrogliomas, 14/42 oligoastrocytomas and 3/10 preselected astrocytomas. With the exception of a single case, all CIC mutations occurred in tumors with combined 1p/19q losses. In contrast to oligodendrogliomas where CIC mutations were always detected along with 1p/19q co-deletion, CIC mutations were only found in 52 % of the 1p/19q co-deleted oligoastrocytomas. FUBP1 mutations were detected in 7/61 tumors, all presenting with CIC mutations. FUBP1 mutations appear to cluster in the DNA binding domain spanning exons 5-14. CIC and FUBP1 mutations exclusively occurred in presence of either IDH1 or IDH2 mutations. Our data confirm CIC and FUBP1 mutations in oligodendrogliomas and demonstrate the presence of these mutations in oligoastrocytomas.

Substrate-dependent control of ERK phosphorylation can lead to oscillations

The extracellular signal-regulated kinase (ERK) controls cellular processes by phosphorylating multiple substrates. The ERK protein can use the same domains to interact with phosphatases, which dephosphorylate and deactivate ERK, and with substrates, which connect ERK to its downstream effects. As a consequence, substrates can compete with phosphatases and control the level of ERK phosphorylation. We propose that this effect can qualitatively change the dynamics of a network that controls ERK activation. On its own, this network can be bistable, but in a larger system, where ERK accelerates the degradation of a substrate competing with a phosphatase, this network can oscillate. Previous studies proposed that oscillatory ERK signaling requires a negative feedback in which active ERK reduces the rate at which it is phosphorylated by upstream kinase. We argue that oscillations can also emerge even when this rate is constant, due to substrate-dependent control of ERK phosphorylation.

Evaluating the Drosophila Bicoid morphogen gradient system through dissecting the noise in transcriptional bursts

MOTIVATION

We describe a statistical model to dissect the noise in transcriptional bursts in a developmental system.

RESULTS

We assume that, at any given moment of time, each copy of a native gene inside a cell can exist in either a bursting (active) or non-bursting (inactive) state. The experimentally measured total noise in the transcriptional states of a gene in a population of cells can be mathematically dissected into two contributing components: internal and external. While internal noise quantifies the stochastic nature of transcriptional bursts, external noise is caused by cell-to-cell differences including fluctuations in activator concentration. We use our developed methods to analyze the Drosophila Bicoid (Bcd) morphogen gradient system. For its target gene hunchback (hb), the noise properties can be recapitulated by a simplified gene regulatory model in which Bcd acts as the only input, suggesting that the external noise in hb transcription is primarily derived from fluctuations in the Bcd activator input. However, such a simplified gene regulatory model is insufficient to predict the noise properties of another Bcd target gene, orthodenticle (otd), suggesting that otd transcription is sensitive to additional external fluctuations beyond those in Bcd. Our results show that analysis of the relationship between input and output noise can reveal important insights into how a morphogen gradient system works. Our study also advances the knowledge about transcription at a fundamental level.

CONTACT

jun.ma@cchmc.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.