Stabilizing patterning in the Drosophila segment polarity network by selecting models in silico

The segmentation of Drosophila is a prime model to study spatial patterning during embryogenesis. The spatial expression of segment polarity genes results from a complex network of interacting proteins whose expression products are maintained after successful segmentation. This prompted us to investigate the stability and robustness of this process using a dynamical model for the segmentation network based on Boolean states. The model consists of intra-cellular as well as inter-cellular interactions between adjacent cells in one spatial dimension. We quantify the robustness of the dynamical segmentation process by a systematic analysis of mutations. Our starting point consists in a previous Boolean model for Drosophila segmentation. We define mathematically the notion of dynamical robustness and show that the proposed model exhibits limited robustness in gene expression under perturbations. We applied in silico evolution (mutation and selection) and discover two classes of modified gene networks that have a more robust spatial expression pattern. We verified that the enhanced robustness of the two new models is maintained in differential equations models. By comparing the predicted model with experiments on mutated flies, we then discuss the two types of enhanced models. Drosophila patterning can be explained by modelling the underlying network of interacting genes. Here we demonstrate that simple dynamical considerations and in silico evolution can enhance the model to robustly express the expected pattern, helping to elucidate the role of further interactions.     

An analysis of the class of gene regulatory functions implied by a biochemical model

Understanding the integrated behavior of genetic regulatory networks, in which genes regulate one another's activities via RNA and protein products, is emerging as a dominant problem in systems biology. One widely studied class of models of such networks includes genes whose expression values assume Boolean values (i.e., on or off). Design decisions in the development of Boolean network models of gene regulatory systems include the topology of the network (including the distribution of input- and output-connectivity) and the class of Boolean functions used by each gene (e.g., canalizing functions, post functions, etc.). For example, evidence from simulations suggests that biologically realistic dynamics can be produced by scale-free network topologies with canalizing Boolean functions. This work seeks further insights into the design of Boolean network models through the construction and analysis of a class of models that include more concrete biochemical mechanisms than the usual abstract model, including genes and gene products, dimerization, cis-binding sites, promoters and repressors. In this model, it is assumed that the system consists of N genes, with each gene producing one protein product. Proteins may form complexes such as dimers, trimers, etc. The model also includes cis-binding sites to which proteins may bind to form activators or repressors. Binding affinities are based on structural complementarity between proteins and binding sites, with molecular binding sites modeled by bit-strings. Biochemically plausible gene expression rules are used to derive a Boolean regulatory function for each gene in the system. The result is a network model in which both topological features and Boolean functions arise as emergent properties of the interactions of components at the biochemical level. A highly biased set of Boolean functions is observed in simulations of networks of various sizes, suggesting a new characterization of the subset of Boolean functions that are likely to appear in gene regulatory networks.     

Robustness and state-space structure of Boolean gene regulatory models

Robustness to perturbation is an important characteristic of genetic regulatory systems, but the relationship between robustness and model dynamics has not been clearly quantified. We propose a method for quantifying both robustness and dynamics in terms of state-space structures, for Boolean models of genetic regulatory systems. By investigating existing models of the Drosophila melanogaster segment polarity network and the Saccharomyces cerevisiae cell-cycle network, we show that the structure of attractor basins can yield insight into the underlying decision making required of the system, and also the way in which the system maximises its robustness. In particular, gene networks implementing decisions based on a few genes have simple state-space structures, and their attractors are robust by virtue of their simplicity. Gene networks with decisions that involve many interacting genes have correspondingly more complicated state-space structures, and robustness cannot be achieved through the structure of the attractor basins, but is achieved by larger attractor basins that dominate the state space. These different types of robustness are demonstrated by the two models: the D. melanogaster segment polarity network is robust due to simple attractor basins that implement decisions based on spatial signals; the S. cerevisiae cell-cycle network has a complicated state-space structure, and is robust only due to a giant attractor basin that dominates the state space.     

Modeling segmental patterning in Drosophila: Maternal and gap genes

We propose a new mathematical model describing the establishment of maternal and gap proteins segmental patterning along the antero-posterior axis of the Drosophila early embryo. This model is based on experimental data and, without recurring to pre-defined activation thresholds, predicts qualitatively and quantitatively the expression patterns of the maternal and gap proteins, as well as the expression patterns of proteins resulting from mRNA ectopic expression and from some loss-of-function mutations. We conclude that the gap genes segmental patterning and consequent spatial organization of the embryo is determined by three main factors: (1) the initial positioning of the maternal bicoid and torso mRNA inside the egg, and subsequent diffusion of the corresponding proteins; (2) the structure of the genetic regulatory network; (3) the role of conservation laws in the establishment of steady and non-uniform spatial distributions of non-diffusing proteins.     

Notch signaling relieves the joint-suppressive activity of Defective proventriculus in the Drosophila leg

Segmentation plays crucial roles during morphogenesis. Drosophila legs are divided into segments along the proximal-distal axis by flexible structures called joints. Notch signaling is necessary and sufficient to promote leg growth and joint formation, and is activated in distal cells of each segment in everting prepupal leg discs. The homeobox gene defective proventriculus (dve) is expressed in regions both proximal and distal to the intersegmental folds at 4 h after puparium formation (APF). Dve-expressing region partly overlaps with the Notch-activated region, and they become a complementary pattern at 6 h APF. Interestingly, dve mutant legs resulted in extra joint formation at the center of each tarsal segment, and the forced expression of dve caused a jointless phenotype. We present evidence that Dve suppresses the potential joint-forming activity, and that Notch signaling represses Dve expression to form joints.     

Genetic analysis of mitochondrial protein misfolding in Drosophila melanogaster

Protein misfolding has a key role in several neurological disorders including Parkinson's disease. Although a clear mechanism for such proteinopathic diseases is well established when aggregated proteins accumulate in the cytosol, cell nucleus, endoplasmic reticulum and extracellular space, little is known about the role of protein aggregation in the mitochondria. Here we show that mutations in both human and fly PINK1 result in higher levels of misfolded components of respiratory complexes and increase in markers of the mitochondrial unfolded protein response. Through the development of a genetic model of mitochondrial protein misfolding employing Drosophila melanogaster, we show that the in vivo accumulation of an unfolded protein in mitochondria results in the activation of AMP-activated protein kinase-dependent autophagy and phenocopies of pink1 and parkin mutants. Parkin expression acts to clear mitochondria with enhanced levels of misfolded proteins by promoting their autophagic degradation in vivo, and refractory to Sigma P (ref(2)P), the Drosophila orthologue of mammalian p62, is a critical downstream effector of this quality control pathway. We show that in flies, a pathway involving pink1, parkin and ref(2)P has a role in the maintenance of a viable pool of cellular mitochondria by promoting organellar quality control.     

Expression of gooseberry-proximal in the Drosophila developing nervous system responds to cues provided by segment polarity genes

Summary Segment polarity genes define the cell states that are required for proper organization of each metameric unit of the Drosophila embryo. Among these, the gooseberry locus has been shown to be composed of two closely related genes which are expressed in an overlapping single-segment periodicity. We have used specific antibodies raised against the protein product of the gooseberry proximal (gsb-p) gene to determine the spatial distribution of this antigen in wild type embryos, and to monitor the effects of segment polarity mutants on the pattern of the gsb-p protein distribution. We find that the gsb-p protein accumulates beneath each posterior axonal commissure in the progeny of neuroblasts deriving from the epidermal compartments of wingless (wg) and engrailed (en) expression. The results of this analysis support the idea that gsb-p has a specific role in the control of cell fates during neurogenesis, and indicate that en and wg provide critical positional cues to define the domain in which gsbp will be activated. Furthermore, these data suggest that, in order to be expressed in the embryonic CNS, gsb-p may preliminarily require activity of the gooseberry-distal gene in the epidermis. Offprint requests to: S. Côté     

Genetic regulation of engrailed and wingless in Tribolium segmentation and the evolution of pair-rule segmentation

In Drosophila, primary pair-rule genes establish the parasegmental boundaries and indirectly control the periodic expression of the segment polarity genes engrailed (en) and wingless (wg) via regulation of secondary pair-rule genes. Although orthologs of some Drosophila pair-rule genes are not required for proper segmentation in Tribolium, segmental expression of Tc-en and Tc-wg is conserved. To understand how these segment polarity genes are regulated, we examined the results of expressing one or two pair-rule genes in the absence of the other known pair-rule genes. Expression of one or both of the secondary pair-rule genes, Tc-sloppy-paired (Tc-slp) and Tc-paired (Tc-prd), activated Tc-wg in the absence of the primary pair-rule genes, Tc-even-skipped (Tc-eve), Tc-runt (Tc-run) and Tc-odd-skipped (Tc-odd). Tc-eve alone failed to activate Tc-wg or Tc-en, but in combination with Tc-run or Tc-prd activated Tc-en. These results, interpreted within the pair-rule gene expression patterns, suggest separate models for the genetic regulation of the juxtaposed expression of Tc-wg and Tc-en at odd- and even-numbered parasegmental boundaries, respectively. Conserved interactions between eve and prd at the anterior boundary of odd-numbered parasegments may reflect an ancestral segmentation mechanism that functioned in every segment prior to the evolution of pair-rule segmentation.     

Structure-function analysis of the Drosophila retinal determination protein Dachshund

Dachshund (Dac) is a highly conserved nuclear protein that is distantly related to the Ski/Sno family of corepressor proteins. In Drosophila, Dac is necessary and sufficient for eye development and, along with Eyeless (Ey), Sine oculis (So), and Eyes absent (Eya), forms the core of the retinal determination (RD) network. In vivo and in vitro experiments suggest that members of the RD network function together in one or more complexes to regulate the expression of downstream targets. For example, Dac and Eya synergize in vivo to induce ectopic eye formation and they physically interact through conserved domains. Dac contains two highly conserved domains, named DD1 and DD2, but no function has been assigned to either of them in an in vivo context. We performed structure-function studies to understand the relationship between the conserved domains of Dac and the rest of the protein and to determine the function of each domain during development. We show that only DD1 is essential for Dac function and while DD2 facilitates DD1, it is not absolutely essential in spite of more than 500 million years of conservation. Moreover, the physical interaction between Eya and DD2 is not required for the genetic synergy between the two proteins. Finally, we show that DD1 also plays a central role for nuclear localization of Dac.     

The heterotrimeric protein Go is required for the formation of heart epithelium in Drosophila.

The gene encoding the alpha subunit of the Drosophila Go protein is expressed early in embryogenesis in the precursor cells of the heart tube, of the visceral muscles, and of the nervous system. This early expression coincides with the onset of the mesenchymal-epithelial transition to which are subjected the cardial cells and the precursor cells of the visceral musculature. This gene constitutes an appropriate marker to follow this transition. In addition, a detailed analysis of its expression suggests that the cardioblasts originate from two subpopulations of cells in each parasegment of the dorsal mesoderm that might depend on the wingless and hedgehog signaling pathways for both their determination and specification. In the nervous system, the expression of Goalpha shortly precedes the beginning of axonogenesis. Mutants produced in the Goalpha gene harbor abnormalities in the three tissues in which the gene is expressed. In particular, the heart does not form properly and interruptions in the heart epithelium are repeatedly observed, henceforth the brokenheart (bkh) name. Furthermore, in the bkh mutant embryos, the epithelial polarity of cardial cells was not acquired (or maintained) in various places of the cardiac tube. We predict that bkh might be involved in vesicular traffic of membrane proteins that is responsible for the acquisition of polarity.     

Immunohistochemical analysis of a novel dehydroepiandrosterone sulfotransferase-like protein in Drosophila neural circuits

Sulfotransferase (ST)-catalyzed sulfation plays an important role in various neuronal functions such as homeostasis of catecholamine neurotransmitters and hormones. Drosophila is a popular model for the study of memory and behavioral manifestations because it is able to mimic the intricate neuroregulation and recognition in humans. However, there has been no evidence indicating that cytosolic ST(s) is(are) present in Drosophila. The aim of this study is to investigate whether or not cytosolic ST(s) is(are) expressed in the Drosophila nervous system. Immunoblot analysis demonstrated the presence of dehydroepiandrosterone (DHEA) ST-like protein in Drosophila brain and a sensitive fluorometric assay revealed its sulfating activity toward DHEA. Immunohistochemical staining demonstrated this DHEA ST-like protein to be abundant in specific neurons as well as in several bundles of nerve fibers in Drosophila. Clarification of a possible link between ST and a neurotransmitter-mediated effect may eventually aid in designing approaches for alleviating neuronal disorders in humans.     

Functions of the segment polarity genes midline and H15 in Drosophila melanogaster neurogenesis

The Drosophila melanogaster ventral nerve cord derives from neural progenitor cells called neuroblasts. Individual neuroblasts have unique gene expression profiles and give rise to distinct clones of neurons and glia. The specification of neuroblast identity provides a cell intrinsic mechanism which ultimately results in the generation of progeny which are different from each other. Segment polarity genes have a dual function in early neurogenesis: within distinct regions of the neuroectoderm, they are required both for neuroblast formation and for the specification of neuroblast identity. Previous studies of segment polarity gene function largely focused on neuroblasts that arise within the posterior part of the segment. Here we show that the segment polarity gene midline is required for neuroblast formation in the anterior-most part of the segment. Moreover, midline contributes to the specification of anterior neuroblast identity by negatively regulating the expression of Wingless and positively regulating the expression of Mirror. In the posterior-most part of the segment, midline and its paralog, H15, have partially redundant functions in the regulation of the NB marker Eagle. Hence, the segment polarity genes midline and H15 play an important role in the development of the ventral nerve cord in the anterior- and posterior-most part of the segment.     

DaPKC-dependent phosphorylation of Crumbs is required for epithelial cell polarity in Drosophila

Both in Drosophila and vertebrate epithelial cells, the establishment of apicobasal polarity requires the apically localized, membrane-associated Par-3-Par-6-aPKC protein complex. In Drosophila, this complex colocalizes with the Crumbs-Stardust (Sdt)-Pals1-associated TJ protein (Patj) complex. Genetic and molecular analyses suggest a functional relationship between them. We show, by overexpression of a kinase-dead Drosophila atypical PKC (DaPKC), the requirement for the kinase activity of DaPKC to maintain the position of apical determinants and to restrict the localization of basolateral ones. We demonstrate a novel physical interaction between the apical complexes, via direct binding of DaPKC to both Crb and Patj, and identify Crumbs as a phosphorylation target of DaPKC. This phosphorylation of Crumbs is functionally significant. Thus, a nonphosphorylatable Crumbs protein behaves in vivo as a dominant negative. Moreover, the phenotypic effect of overexpressing wild-type Crumbs is suppressed by reducing DaPKC activity. These results provide a mechanistic framework for the functional interaction between the Par-3-Par-6-aPKC and Crumbs-Sdt-Patj complexes based in the posttranslational modification of Crb by DaPKC.     

The pleiohomeotic gene is required for maintaining expression of genes functioning in ventral appendage formation in Drosophila melanogaster

Polycomb group (PcG) proteins are negative regulators that maintain the expression of homeotic genes and affect cell proliferation. Pleiohomeotic (Pho) is a unique PcG member with a DNA-binding zinc finger motif and was proposed to recruit other PcG proteins to form a complex. The pho null mutants exhibited several mutant phenotypes such as the transformation of antennae to mesothoracic legs. We examined the effects of pho on the identification of ventral appendages and proximo-distal axis formation during postembryogenesis. In the antennal disc of the pho mutant, Antennapedia (Antp), which is a selector gene in determining leg identity, was ectopically expressed. The homothorax (hth), dachshund (dac) and Distal-less (Dll) genes involved in proximo-distal axis formation were also abnormally expressed in both the antennal and leg discs of the pho mutant. The engrailed (en) gene, which affects the formation of the anterior-posterior axis, was also misexpressed in the anterior compartment of antennal and leg discs. These mutant phenotypes were enhanced in the mutant background of Posterior sex combs (Psc) and pleiohomeotic-like (phol), which are another PcG genes. These results suggest that pho functions in maintaining expression of genes involved in the formation of ventral appendages and the proximo-distal axis.     

The role of the pro-apoptotic protein Siva in the pathogenesis of Familial Mediterranean fever: A structural and functional analysis

Familial Mediterranean fever (FMF) is an autosomal, recessive disease, attributed to mutations in MEFV gene encoding pyrin, which is characterized by recurrent, acute and self-limiting attacks of fever as well as an increased neutrophil and monocyte apoptosis. Most disease-associated mutations in MEFV gene reside on the C-terminal PRYSPRY (B30.2) domain of pyrin, an area found to interact with the pro-apoptotic protein Siva. Because apoptotic events may be contributing to endogenous inflammation we hypothesized that mutations in pyrin may affect Siva-mediated apoptosis. The confirmation of this hypothesis would be of a great biological significance since it would be demonstrated a connection between apoptosis and inflammation. We used homology modeling to construct a 3-D model of Siva protein and the constructed model of Siva defined structural elements with potential of binding other proteins to induce apoptosis. Given that Siva protein binds pyrin as shown by transfection and immunoprecipitation experiments, apoptosis was assessed by FACS and Western blotting. No differences in rates of apoptosis in myeloid cells (THP-1) upon transfection with either wt pyrin or mutant forms of pyrin were found. Patients with FMF did not display any mutations in the Siva-1 (full length) gene. Siva-1 was not linked to pyrin in the major predicted FMF gene network constructed using a literature-curated gene signature for FMF. These results suggest that Siva-mediated unprovoked apoptosis is not likely to be involved in the pathogenesis of FMF.