Canalization and control in automata networks: body segmentation in Drosophila melanogaster

We present schema redescription as a methodology to characterize canalization in automata networks used to model biochemical regulation and signalling. In our formulation, canalization becomes synonymous with redundancy present in the logic of automata. This results in straightforward measures to quantify canalization in an automaton (micro-level), which is in turn integrated into a highly scalable framework to characterize the collective dynamics of large-scale automata networks (macro-level). This way, our approach provides a method to link micro- to macro-level dynamics -- a crux of complexity. Several new results ensue from this methodology: uncovering of dynamical modularity (modules in the dynamics rather than in the structure of networks), identification of minimal conditions and critical nodes to control the convergence to attractors, simulation of dynamical behaviour from incomplete information about initial conditions, and measures of macro-level canalization and robustness to perturbations. We exemplify our methodology with a well-known model of the intra- and inter cellular genetic regulation of body segmentation in Drosophila melanogaster. We use this model to show that our analysis does not contradict any previous findings. But we also obtain new knowledge about its behaviour: a better understanding of the size of its wild-type attractor basin (larger than previously thought), the identification of novel minimal conditions and critical nodes that control wild-type behaviour, and the resilience of these to stochastic interventions. Our methodology is applicable to any complex network that can be modelled using automata, but we focus on biochemical regulation and signalling, towards a better understanding of the (decentralized) control that orchestrates cellular activity -- with the ultimate goal of explaining how do cells and tissues 'compute'

Role of transcriptional regulation in the evolution of plant phenotype: A dynamic systems approach

© 2015 Wiley Periodicals, Inc. A growing body of evidence suggests that alterations in transcriptional regulation of genes involved in modulating development are an important part of phenotypic evolution, and this can be documented among species and within populations. While the effects of differential transcriptional regulation in organismal development have been preferentially studied in animal systems, this phenomenon has also been addressed in plants. In this review, we summarize evidence for cis-regulatory mutations, trans-regulatory changes and epigenetic modifications as molecular events underlying important phenotypic alterations, and thus shaping the evolution of plant development. We postulate that a mechanistic understanding of why such molecular alterations have a key role in development, morphology and evolution will have to rely on dynamic models of complex regulatory networks that consider the concerted action of genetic and nongenetic components, and that also incorporate the restrictions underlying the genotype to phenotype mapping process.CONACyT 180098, 180380, 167705, 152649 and PAPIIT UNAM IN203214-3, IN203113-3, IN203814-3. BFU2012–34821 (MINECO) to C.G. and an institutional grant from Fundación Ramón Aceres to CBMSOPeer Reviewe

Analysis of craniofacial defects in Six1/Eya1-associated Branchio-Oto-Renal Syndrome

Poster Session I - Morphogenesis: 205/B10117th ISDB 2013 cum 72nd Annual Meeting of the Society for Developmental Biology, 7th Latin American Society of Developmental Biology Meeting and 11th Congreso de la Sociedad Mexicana de Biologia del Desarrollo.Branchio-Oto-Renal (BOR) syndrome patients exhibit craniofacial and renal anomalies as well as deafness. BOR syndrome is caused by mutations in Six1 or Eya1, both of which regulate cell proliferation and differentiation. The molecular mechanism underlying the craniofacial and branchial arch (BA) defects in BOR syndrome is unclear. We have found that Hoxb3 is up-regulated in the second branchial arch (BA2) of Six1-/- mutants. Moreover, Hoxb3 over-expression in transgenic mice leads to BA abnormalities which are similar to the BA defects in Six1-/- or Eya1-/- mutants, suggesting a regulatory relationship among Six1, Eya1 and Hoxb3 genes. The aim of this study is to investigate the molecular mechanism underlying abnormal BA development in BOR syndrome using Six1 and Eya1 mutant mice. Two potential Six1 binding sites were identified on the Hoxb3 gene. In vitro and in vivo Chromatin IP assays showed that Six1 could directly bind to one of the sites specifically. Furthermore, using a chick in ovo luciferase assay we showed that Six1 could suppress gene expression through one of the specific binding sites. On the other hand, in Six1-/- mutants, we found that the Notch ligand Jag1 was up-regulated in BA2. Similarly, in Hoxb3 transgenic mice, ectopic expression of Jag1 could be also detected in BA2. To investigate the activation of Notch signaling pathway, we found that Notch intracellular domain (NICD), a direct indicator of Notch pathway activation, was up-regulated in BAs of Six1-/-; Eya1-/- double mutants. Our results indicate that Hoxb3 and Notch signaling pathway are involved in mediating the craniofacial defects of Six1/Eya1-associated Branchio-Oto-Renal Syndrome.postprin

Sox10 regulates enteric neural crest cell migration in the developing gut

Concurrent Sessions 1: 1.3 - Organs to organisms: Models of Human Diseases: abstract no. 1417th ISDB 2013 cum 72nd Annual Meeting of the Society for Developmental Biology, VII Latin American Society of Developmental Biology Meeting and XI Congreso de la Sociedad Mexicana de Biologia del Desarrollo. The Conference's web site is located at http://www.inb.unam.mx/isdb/Sox10 is a HMG-domain containing transcription factor which plays important roles in neural crest cell survival and differentiation. Mutations of Sox10 have been identified in patients with Waardenburg-Hirschsprung syndrome, who suffer from deafness, pigmentation defects and intestinal aganglionosis. Enteric neural crest cells (ENCCs) with Sox10 mutation undergo premature differentiation and fail to colonize the distal hindgut. It is unclear, however, whether Sox10 plays a role in the migration of ENCCs. To visualize the migration behaviour of mutant ENCCs, we generated a Sox10NGFP mouse model where EGFP is fused to the N-terminal domain of Sox10. Using time-lapse imaging, we found that ENCCs in Sox10NGFP/+ mutants displays lower migration speed and altered trajectories compared to normal controls. This behaviour was cell-autonomous, as shown by organotypic grafting of Sox10NGFP/+ gut segments onto control guts and vice versa. ENCCs encounter different extracellular matrix (ECM) molecules along the developing gut. We performed gut explant culture on various ECM and found that Sox10NGFP/+ ENCCs tend to form aggregates, particularly on fibronectin. Time-lapse imaging of single cells in gut explant culture indicated that the tightly-packed Sox10 mutant cells failed to exhibit contact inhibition of locomotion. We determined the expression of adhesion molecule families by qPCR analysis, and found integrin expression unaffected while L1-cam and selected cadherins were altered, suggesting that Sox10 mutation affects cell adhesion properties of ENCCs. Our findings identify a de novo role of Sox10 in regulating the migration behaviour of ENCCs, which has important implications for the treatment of Hirschsprung disease.postprin

Transcriptome-based Gene Networks for Systems-level Analysis of Plant Gene Functions

Present day genomic technologies are evolving at an unprecedented rate, allowing interrogation of cellular activities with increasing breadth and depth. However, we know very little about how the genome functions and what the identified genes do. The lack of functional annotations of genes greatly limits the post-analytical interpretation of new high throughput genomic datasets. For plant biologists, the problem is much severe. Less than 50% of all the identified genes in the model plant Arabidopsis thaliana, and only about 20% of all genes in the crop model Oryza sativa have some aspects of their functions assigned. Therefore, there is an urgent need to develop innovative methods to predict and expand on the currently available functional annotations of plant genes. With open-access catching the ‘pulse’ of modern day molecular research, an integration of the copious amount of transcriptome datasets allows rapid prediction of gene functions in specific biological contexts, which provide added evidence over traditional homology-based functional inference. The main goal of this dissertation was to develop data analysis strategies and tools broadly applicable in systems biology research. Two user friendly interactive web applications are presented: The Rice Regulatory Network (RRN) captures an abiotic-stress conditioned gene regulatory network designed to facilitate the identification of transcription factor targets during induction of various environmental stresses. The Arabidopsis Seed Active Network (SANe) is a transcriptional regulatory network that encapsulates various aspects of seed formation, including embryogenesis, endosperm development and seed-coat formation. Further, an edge-set enrichment analysis algorithm is proposed that uses network density as a parameter to estimate the gain or loss in correlation of pathways between two conditionally independent coexpression networks

An integrative and systems biology approach to Drosophila melanogaster transcriptomes

The availability of fully sequenced genomes of the model organisms including Drosophila, and their subsequent annotation has afforded seamless opportunities for reverse genetics in a complex model organism. With the advent of DNA microarrays to assay the levels of tens of thousands of genes in a single sample, functional genomics has been significantly aided to understand the functions in systems context. These microarrays have been employed predominantly on the RNA samples that are extracted from the whole animals for example at different developmental stages or in response to external stimuli. However, these approaches relied on the expression patterns that represent the sum of transcription coming from all the organs, which do not estimate the tissue-specificity of transcription. The purpose of this thesis is to provide tissue-specific transcriptomes of Drosophila melanogaster that were generated as part of the large FlyAtlas project using Affymetrix Drosophila GeneChips® (or microarrays). These chips, one at a time interrogate the levels of 18,500 transcripts (that represent all known genes) using 18,880 distinct probe sets in a single, total RNA sample. For each tissue, four biological replicates were analysed using the chips and the normalised signal intensities were obtained that represent the relative levels of mRNA expression. Using the transcriptomes, a general analysis was performed for potential novel insights into tissue-specific functions (Chintapalli et al., 2007) (Chapter 3). Then, a comparative analysis of epithelial tissues was performed to understand how the epithelia are organised in terms of their transcriptomes (Chapter 4). The Malpighian tubules are the Drosophila epithelial counterparts of the human kidney. They show asymmetric organisation in the body cavity. FlyAtlas segment-specific tubule transcriptomes allowed the comparison of their potential functional similarities and differences, thus to understand the asymmetry in function (Chapter 5)(Chintapalli, 2012). This identified a human Best vitelliform macular dystrophy (BVMD) disease homolog, Best2 in only the anterior pair of tubules that have the morphologically and functionally distinct enlarged initial (or distal) segment, a storage organ for Ca2+. Bestrophins were accordingly selected as candidate genes to analyse organismal functions, and thus to validate previous two theories that implicated bestrophins as Ca2+-activated Clˉ channels and/or Ca2+ channel regulators (Chapter 6). The confocal microscopy analysis of bestrophin YFP fusion proteins revealed interesting and novel localisations of bestrophins, in that Best1 was found in the apical plasma membranes, Best2 localised to peroxisomes, Best3 and Best4 were found intracellular. The salt survival analysis showed that Best1 is essential in regulating extra salt levels in the body. Furthermore, the fluid secretion analysis showed Best1’s potential role in Ca2+-dependent Clˉ function. Interestingly, the flies with reduced levels of Best2 expression showed increased ability to survive on extra salt food; the basis for this was investigated further in Chapter 7. Best2 was also found abundant in the eyes than anywhere else in the head. A comparative analysis of anterior tubule- and eye-specific transcriptomes revealed a potential overlap of Ca2+ signaling components, in that the PLCβ signaling was one. A neuropeptide Ca2+ agonist, capa1 evoked secondary cytosolic Ca2+ responses were found high in Best2 knockdowns. A quantitative PCR (qPCR) analysis of candidate Ca2+ signaling and homeostasis genes in Best2 mutants revealed their gene expression upregulation, under control-fed and salt-fed conditions than their wildtype controls, fed on similar diet regimes. The norpA that encodes PLCβ was found significantly enriched in the mutants. Cyp6a23 is another gene that was highly upregulated in Best2 mutants; it is a Drosophila homologue of human Cyp11b, a Ca2+-responsive gene implicated in renal salt wasting. Upon the downregulation of Cyp6a23, flies became sensitive to salt diet feeding. Other genes investigated and found to be upregulated in the mutants include transient-receptor-potential (trp) Ca2+ channel and retinal degeneration C (rdgC). Together, these results strongly suggest Best2 as a potential Ca2+ channel regulator, and provide fascinating insight into bestrophin function. Peroxisomal localisation of Best2 in line with the implication that peroxisomes act as dynamic regulators of cell Ca2+ homeostasis led to another aspect of the project (Chapter 8). This study identified two peroxins that are most abundant in the tubules and play essential roles in the novel cyclic nucleotide-regulated peroxisomal Ca2+ sequestration and transport pathway and that are detrimental for peroxisome biogenesis and proliferation

Tree Peony Species Are a Novel Resource for Production of α-Linolenic Acid

Tree peony is known worldwide for its excellent ornamental and medical values, but recent reports that their seeds contain over 40% α-linolenic acid (ALA), an essential fatty acid for humans drew additional interest of biochemists. To understand the key factors that contribute to this rich accumulation of ALA, we carried out a comprehensive study of oil accumulation in developing seeds of nine wild tree peony species. The fatty acid content and composition was highly variable among the nine species; however, we selected a high- (P. rockii) and low-oil (P. lutea) accumulating species for a comparative transcriptome analysis. Similar to other oilseed transcriptomic studies, upregulation of select genes involved in plastidial fatty acid synthesis, and acyl editing, desaturation and triacylglycerol assembly in the endoplasmic reticulum was noted in seeds of P. rockii relative to P. lutea. Also, in association with the ALA content, transcript levels for fatty acid desaturases (SAD, FAD2 and FAD3), which encode for enzymes necessary for polyunsaturated fatty acid synthesis were higher in P. rockii compared to P. lutea. We further showed that the overexpression of PrFAD2 and PrFAD3 in Arabidopsis increased linoleic and α-linolenic acid content, respectively and modulated their final ratio in the seed oil. In conclusion, we identified the key steps that contribute to efficient ALA synthesis and validated the necessary desaturases in P. rockii that are responsible for not only increasing oil content but also modulating 18:2/18:3 ratio in seeds. Together, these results will aid to improve essential fatty acid content in seeds of tree peonies and other crops of agronomic interest

Intracellular transport driven by cytoskeletal motors: General mechanisms and defects

Cells are strongly out-of-equilibrium systems driven by continuous energy supply. They carry out many vital functions requiring active transport of various ingredients and organelles, some being small, others being large. The cytoskeleton, composed of three types of filaments, determines the shape of the cell and plays a role in cell motion. It also serves as a road network for the so-called cytoskeletal motors. These molecules can attach to a cytoskeletal filament, perform directed motion, possibly carrying along some cargo, and then detach. It is a central issue to understand how intracellular transport driven by molecular motors is regulated, in particular because its breakdown is one of the signatures of some neuronal diseases like the Alzheimer. We give a survey of the current knowledge on microtubule based intracellular transport. We first review some biological facts obtained from experiments, and present some modeling attempts based on cellular automata. We start with background knowledge on the original and variants of the TASEP (Totally Asymmetric Simple Exclusion Process), before turning to more application oriented models. After addressing microtubule based transport in general, with a focus on in vitro experiments, and on cooperative effects in the transportation of large cargos by multiple motors, we concentrate on axonal transport, because of its relevance for neuronal diseases. It is a challenge to understand how this transport is organized, given that it takes place in a confined environment and that several types of motors moving in opposite directions are involved. We review several features that could contribute to the efficiency of this transport, including the role of motor-motor interactions and of the dynamics of the underlying microtubule network. Finally, we discuss some still open questions.Comment: 74 pages, 43 figure