Vitaly V. Gursky

Canalization of gene expression in the Drosophila blastoderm by gap gene cross regulation.

Developing embryos exhibit a robust capability to reduce phenotypic variations that occur naturally or as a result of experimental manipulation. This reduction in variation occurs by an epigenetic mechanism called canalization, a phenomenon which has resisted understanding because of a lack of necessary molecular data and of appropriate gene regulation models. In recent years, quantitative gene expression data have become available for the segment determination process in the Drosophila blastoderm, revealing a specific instance of canalization. These data show that the variation of the zygotic segmentation gene expression patterns is markedly reduced compared to earlier levels by the time gastrulation begins, and this variation is significantly lower than the variation of the maternal protein gradient Bicoid. We used a predictive dynamical model of gene regulation to study the effect of Bicoid variation on the downstream gap genes. The model correctly predicts the reduced variation of the gap gene expression patterns and allows the characterization of the canalizing mechanism. We show that the canalization is the result of specific regulatory interactions among the zygotic gap genes. We demonstrate the validity of this explanation by showing that variation is increased in embryos mutant for two gap genes, Krüppel and knirps, disproving competing proposals that canalization is due to an undiscovered morphogen, or that it does not take place at all. In an accompanying article in PLoS Computational Biology (doi:10.1371/journal.pcbi.1000303), we show that cross regulation between the gap genes causes their expression to approach dynamical attractors, reducing initial variation and providing a robust output. These results demonstrate that the Bicoid gradient is not sufficient to produce gap gene borders having the low variance observed, and instead this low variance is generated by gap gene cross regulation. More generally, we show that the complex multigenic phenomenon of canalization can be understood at a quantitative and predictive level by the application of a precise dynamical model.

Canalization of gene expression and domain shifts in the Drosophila blastoderm by dynamical attractors.

The variation in the expression patterns of the gap genes in the blastoderm of the fruit fly Drosophila melanogaster reduces over time as a result of cross regulation between these genes, a fact that we have demonstrated in an accompanying article in PLoS Biology (see Manu et al., doi:10.1371/journal.pbio.1000049). This biologically essential process is an example of the phenomenon known as canalization. It has been suggested that the developmental trajectory of a wild-type organism is inherently stable, and that canalization is a manifestation of this property. Although the role of gap genes in the canalization process was established by correctly predicting the response of the system to particular perturbations, the stability of the developmental trajectory remains to be investigated. For many years, it has been speculated that stability against perturbations during development can be described by dynamical systems having attracting sets that drive reductions of volume in phase space. In this paper, we show that both the reduction in variability of gap gene expression as well as shifts in the position of posterior gap gene domains are the result of the actions of attractors in the gap gene dynamical system. Two biologically distinct dynamical regions exist in the early embryo, separated by a bifurcation at 53% egg length. In the anterior region, reduction in variation occurs because of stability induced by point attractors, while in the posterior, the stability of the developmental trajectory arises from a one-dimensional attracting manifold. This manifold also controls a previously characterized anterior shift of posterior region gap domains. Our analysis shows that the complex phenomena of canalization and pattern formation in the Drosophila blastoderm can be understood in terms of the qualitative features of the dynamical system. The result confirms the idea that attractors are important for developmental stability and shows a richer variety of dynamical attractors in developmental systems than has been previously recognized.

Mechanisms of gap gene expression canalization in the Drosophila blastoderm

BackgroundExtensive variation in early gap gene expression in the Drosophila blastoderm is reduced over time because of gap gene cross regulation. This phenomenon is a manifestation of canalization, the ability of an organism to produce a consistent phenotype despite variations in genotype or environment. The canalization of gap gene expression can be understood as arising from the actions of attractors in the gap gene dynamical system.ResultsIn order to better understand the processes of developmental robustness and canalization in the early Drosophila embryo, we investigated the dynamical effects of varying spatial profiles of Bicoid protein concentration on the formation of the expression border of the gap gene hunchback. At several positions on the anterior-posterior axis of the embryo, we analyzed attractors and their basins of attraction in a dynamical model describing expression of four gap genes with the Bicoid concentration profile accounted as a given input in the model equations. This model was tested against a family of Bicoid gradients obtained from individual embryos. These gradients were normalized by two independent methods, which are based on distinct biological hypotheses and provide different magnitudes for Bicoid spatial variability. We showed how the border formation is dictated by the biological initial conditions (the concentration gradient of maternal Hunchback protein) being attracted to specific attracting sets in a local vicinity of the border. Different types of these attracting sets (point attractors or one dimensional attracting manifolds) define several possible mechanisms of border formation. The hunchback border formation is associated with intersection of the spatial gradient of the maternal Hunchback protein and a boundary between the attraction basins of two different point attractors. We demonstrated how the positional variability for hunchback is related to the corresponding variability of the basin boundaries. The observed reduction in variability of the hunchback gene expression can be accounted for by specific geometrical properties of the basin boundaries.ConclusionWe clarified the mechanisms of gap gene expression canalization in early Drosophila embryos. These mechanisms were specified in the case of hunchback in well defined terms of the dynamical system theory.

Mechanisms of developmental robustness

We present a review of noise buffering mechanisms responsible for developmental robustness. We focus on functions of chaperone Hsp90, miRNA, and cross-regulation of gap genes in Drosophila. The noise buffering mechanisms associated with these functions represent specific examples of the developmental canalization, reducing the phenotypical variability in presence of either genetic or environmental perturbations. We demonstrate that robustness often appears as a function of a network of interacting elements and that the system level approach is needed to understand the mechanisms of noise filtering.

Sequence-based model of gap gene regulatory network

BackgroundThe detailed analysis of transcriptional regulation is crucially important for understanding biological processes. The gap gene network in Drosophila attracts large interest among researches studying mechanisms of transcriptional regulation. It implements the most upstream regulatory layer of the segmentation gene network. The knowledge of molecular mechanisms involved in gap gene regulation is far less complete than that of genetics of the system. Mathematical modeling goes beyond insights gained by genetics and molecular approaches. It allows us to reconstruct wild-type gene expression patterns in silico, infer underlying regulatory mechanism and prove its sufficiency.ResultsWe developed a new model that provides a dynamical description of gap gene regulatory systems, using detailed DNA-based information, as well as spatial transcription factor concentration data at varying time points. We showed that this model correctly reproduces gap gene expression patterns in wild type embryos and is able to predict gap expression patterns in Kr mutants and four reporter constructs. We used four-fold cross validation test and fitting to random dataset to validate the model and proof its sufficiency in data description. The identifiability analysis showed that most model parameters are well identifiable. We reconstructed the gap gene network topology and studied the impact of individual transcription factor binding sites on the model output. We measured this impact by calculating the site regulatory weight as a normalized difference between the residual sum of squares error for the set of all annotated sites and for the set with the site of interest excluded.ConclusionsThe reconstructed topology of the gap gene network is in agreement with previous modeling results and data from literature. We showed that 1) the regulatory weights of transcription factor binding sites show very weak correlation with their PWM score; 2) sites with low regulatory weight are important for the model output; 3) functional important sites are not exclusively located in cis-regulatory elements, but are rather dispersed through regulatory region. It is of importance that some of the sites with high functional impact in hb, Kr and kni regulatory regions coincide with strong sites annotated and verified in Dnase I footprint assays.

Model with asymptotically stable dynamics for Drosophila gap gene network

We consider a model of gap gene expression during the early development of Drosophila embryo. Parameter values in the model have been obtained by fitting to experimental patterns under an additional condition that the solution be asymptotically stable at large times. The patterns at the beginning of gastrulation in such solutions are very close to an actual attractor in the model. It is shown that such solutions are more robust to perturbations of concentrations and parameter values in the model.

How gap genes make their domains: An analytical study based on data driven approximations.

We consider a mathematical formulation of the problem of protein production during segment determination in the Drosophila blastoderm, together with some preliminary results of its analytical study. We reformulate the spatial difference equations as a set of nonlinear partial differential equations and obtain their dimensionless form in the continuum limit. Using previous results obtained by the gene circuit method, we find an asymptotic statement of the problem with a small parameter. Some results of the comparison method applied to the model are obtained, and exact stationary upper solutions are derived. They exhibit distinctive features of localized bell-shaped structures. (c) 2001 American Institute of Physics.

A New Stochastic Model for Subgenomic Hepatitis C Virus Replication Considers Drug Resistant Mutants

As an RNA virus, hepatitis C virus (HCV) is able to rapidly acquire drug resistance, and for this reason the design of effective anti-HCV drugs is a real challenge. The HCV subgenomic replicon-containing cells are widely used for experimental studies of the HCV genome replication mechanisms, for drug testing in vitro and in studies of HCV drug resistance. The NS3/4A protease is essential for virus replication and, therefore, it is one of the most attractive targets for developing specific antiviral agents against HCV. We have developed a stochastic model of subgenomic HCV replicon replication, in which the emergence and selection of drug resistant mutant viral RNAs in replicon cells is taken into account. Incorporation into the model of key NS3 protease mutations leading to resistance to BILN-2061 (A156T, D168V, R155Q), VX-950 (A156S, A156T, T54A) and SCH 503034 (A156T, A156S, T54A) inhibitors allows us to describe the long term dynamics of the viral RNA suppression for various inhibitor concentrations. We theoretically showed that the observable difference between the viral RNA kinetics for different inhibitor concentrations can be explained by differences in the replication rate and inhibitor sensitivity of the mutant RNAs. The pre-existing mutants of the NS3 protease contribute more significantly to appearance of new resistant mutants during treatment with inhibitors than wild-type replicon. The model can be used to interpret the results of anti-HCV drug testing on replicon systems, as well as to estimate the efficacy of potential drugs and predict optimal schemes of their usage.

Modeling of the Drosophila gap-gene network with the variation of the Bcd morphogen

Expression patterns of segmentation genes form under the influence of gradients of maternal transcription factors, which initiate spatially local expression in the segmentation gene cascade. Bcd acts as one of these activators. A model of the regulation in the gap-gene network was studied by varying the Bcd concentration. The topology that is known for the gap-gene network was found to be insufficient for explaining the experimental finding that the hb anterior expression domain shifts when the Bcd concentration changes in the embryo. Modeling that was performed to comply with this experimental finding yielded a new topology, which determined the proper shifts of the hb expression domain. The result indicates that interactions of hb, Kr, and gt act as key regulatory factors to ensure the correct behavior of the hb expression pattern upon changes in Bcd concentration. This study made it possible to specify the limits of the validity of phenomenological models of gene networks.

Course 9 A survey of gene circuit approach applied to modelling of segment determination in fruit fly

Publisher Summary This chapter discusses the ability of gene circuit method to interpret and predict regulatory mechanisms using (as an example) the segment determination system in fruit fly Drosophila. In modern genetics the regulatory interactions in multicellular organisms are inferred by a comparison of mutant and wild type phenotypes. An important example of use of this method is the deduction that virtually all of the pair-rule class of segmentation genes in the fruit fly Drosophila melanogaster is regulated by members of gap gene class. This conclusion follows directly from observations of segmentation gene expression patterns at gastrulation in a variety of gap and pair-rule mutants. The gene circuit is a data driven mathematical modeling method, whose main aim is to extract information about dynamical regulatory interactions between transcription factors from given gene expression patterns. This is achieved in four steps: (1) formulation of a mathematical modeling framework, (2) collection of gene expression data, (3) fitting of the model to expression data to obtain regulatory parameters, and (4) biological analysis of the resulting gene circuits.