Predicting Ancestral Segmentation Phenotypes from Drosophila to Anopheles Using In Silico Evolution

Deriving a genetic regulatory network from an optimization principle

Many biological systems operate near the physical limits to their performance, suggesting that aspects of their behavior and underlying mechanisms could be derived from optimization principles. However, such principles have often been applied only in simplified models. Here, we explore a detailed mechanistic model of the gap gene network in the Drosophila embryo, optimizing its 50+ parameters to maximize the information that gene expression levels provide about nuclear positions. This optimization is conducted under realistic constraints, such as limits on the number of available molecules. Remarkably, the optimal networks we derive closely match the architecture and spatial gene expression profiles observed in the real organism. Our framework quantifies the tradeoffs involved in maximizing functional performance and allows for the exploration of alternative network configurations, addressing the question of which features are necessary and which are contingent. Our results suggest that multiple solutions to the optimization problem might exist across closely related organisms, offering insights into the evolution of gene regulatory networks.

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are associated with insecticide resistance in the major malaria vectors Anopheles gambiae s.l. and Anopheles funestus

Malaria remains one of the highest causes of morbidity and mortality, with 249 million cases and over 608,000 deaths in 2022. Insecticides, which target the Anopheles mosquito vector, are the primary method to control malaria. The widespread nature of resistance to the most important insecticide class, the pyrethroids, threatens the control of this disease. To reverse the stall in malaria control there is urgent need for new vector control tools, which necessitates understanding the molecular basis of pyrethroid resistance. In this study we utilised multi-omics data to identify uridine-diphosphate (UDP)-glycosyltransferases (UGTs) potentially involved in resistance across multiple Anopheles species. Phylogenetic analysis identifies sequence similarities between Anopheline UGTs and those involved in agricultural pesticide resistance to pyrethroids, pyrroles and spinosyns. Expression of five UGTs was characterised in An. gambiae and An. coluzzii to determine constitutive over-expression, induction, and tissue specificity. Furthermore, a UGT inhibitor, sulfinpyrazone, restored susceptibility to pyrethroids and DDT in An. gambiae, An. coluzzii, An. arabiensis and An. funestus, the major African malaria vectors. Taken together, this study provides clear association of UGTs with pyrethroid resistance as well as highlighting the potential use of sulfinpyrazone as a novel synergist for vector control.

Body plan evolvability: The role of variability in gene regulatory networks

Recent computational modeling of early fruit fly (Drosophila) development has characterized the degree to which gene regulation networks can be robust to natural variability. In the first few hours of development, broad spatial gradients of maternally derived transcription factors activate embryonic gap genes. These gap patterns determine the subsequent segmented insect body plan through pair-rule gene expression. Gap genes are expressed with greater spatial precision than the maternal patterns. Computational modeling of the gap-gap regulatory interactions provides a mechanistic understanding for this robustness to maternal variability in wild-type (WT) patterning. A long-standing question in evolutionary biology has been how a system which is robust, such as the developmental program creating any particular species' body plan, is also evolvable, i.e. how can a system evolve or speciate, if the WT form is strongly buffered and protected? In the present work, we use the WT model to explore the breakdown of such Waddington-type 'canalization'. What levels of variability will push the system out of the WT form; are there particular pathways in the gene regulatory mechanism which are more susceptible to losing the WT form; and when robustness is lost, what types of forms are most likely to occur (i.e. what forms lie near the WT)? Manipulating maternal effects in several different pathways, we find a common gap 'peak-to-step' pattern transition in the loss of WT. We discuss these results in terms of the evolvability of insect segmentation, and in terms of experimental perturbations and mutations which could test the model predictions. We conclude by discussing the prospects for using continuum models of pattern dynamics to investigate a wider range of evo-devo problems.

Evo-devo beyond development: Generalizing evo-devo to all levels of the phenotypic evolution

A foundational idea of evo-devo is that morphological variation is not isotropic, that is, it does not occur in all directions. Instead, some directions of morphological variation are more likely than others from DNA-level variation and these largely depend on development. We argue that this evo-devo perspective should apply not only to morphology but to evolution at all phenotypic levels. At other phenotypic levels there is no development, but there are processes that can be seen, in analogy to development, as constructing the phenotype (e.g., protein folding, learning for behavior, etc.). We argue that to explain the direction of evolution two types of arguments need to be combined: generative arguments about which phenotypic variation arises in each generation and selective arguments about which of it passes to the next generation. We explain how a full consideration of the two types of arguments improves the explanatory power of evolutionary theory. Also see the video abstract here: https://youtu.be/Egbvma_uaKc.

Dynamical systems approach to evolution-development congruence: Revisiting Haeckel's recapitulation theory

It is acknowledged that embryonic development has a tendency to proceed from common toward specific. Ernst Haeckel raised the question of why that tendency prevailed through evolution, and the question remains unsolved. Here, we revisit Haeckel's recapitulation theory, that is, the parallelism between evolution and development through numerical evolution and dynamical systems theory. By using intracellular gene expression dynamics with cell-to-cell interaction over spatially aligned cells to represent the developmental process, gene regulation networks (GRN) that govern these dynamics evolve under the selection pressure to achieve a prescribed spatial gene expression pattern. For most numerical evolutionary experiments, the evolutionary pattern changes over generations, as well as the developmental pattern changes governed by the evolved GRN exhibit remarkable similarity. Changes in both patterns consisted of several epochs where stripes are formed in a short time, whereas for other temporal regimes, the pattern hardly changes. In evolution, these quasi-stationary generations are needed to achieve relevant mutations, whereas, in development, they are due to some gene expressions that vary slowly and control the pattern change. These successive epochal changes in development and evolution are represented as common bifurcations in dynamical systems theory, regulating working network structure from feedforward subnetwork to those containing feedback loops. The congruence is the correspondence between successive acquisitions of subnetworks through evolution and changes in working subnetworks in development. Consistency of the theory with the segmentation gene-expression dynamics is discussed. Novel outlook on recapitulation and heterochrony are provided, testable experimentally by the transcriptome and network analysis.

Geometric models for robust encoding of dynamical information into embryonic patterns

During development, cells gradually assume specialized fates via changes of transcriptional dynamics, sometimes even within the same developmental stage. For anterior-posterior (AP) patterning in metazoans, it has been suggested that the gradual transition from a dynamic genetic regime to a static one is encoded by different transcriptional modules. In that case, the static regime has an essential role in pattern formation in addition to its maintenance function. In this work, we introduce a geometric approach to study such transition. We exhibit two types of genetic regime transitions arising through local or global bifurcations, respectively. We find that the global bifurcation type is more generic, more robust, and better preserves dynamical information. This could parsimoniously explain common features of metazoan segmentation, such as changes of periods leading to waves of gene expressions, ‘speed/frequency-gradient’ dynamics, and changes of wave patterns. Geometric approaches appear as possible alternatives to gene regulatory networks to understand development.

Arthropod segmentation

There is now compelling evidence that many arthropods pattern their segments using a clock-and-wavefront mechanism, analogous to that operating during vertebrate somitogenesis. In this Review, we discuss how the arthropod segmentation clock generates a repeating sequence of pair-rule gene expression, and how this is converted into a segment-polarity pattern by ‘timing factor’ wavefronts associated with axial extension. We argue that the gene regulatory network that patterns segments may be relatively conserved, although the timing of segmentation varies widely, and double-segment periodicity appears to have evolved at least twice. Finally, we describe how the repeated evolution of a simultaneous (Drosophila-like) mode of segmentation within holometabolan insects can be explained by heterochronic shifts in timing factor expression plus extensive pre-patterning of the pair-rule genes.

Speeding up anterior-posterior patterning of insects by differential initialization of the gap gene cascade

Recently, it was shown that anterior-posterior patterning genes in the red flour beetle Tribolium castaneum are expressed sequentially in waves. However, in the fruit fly Drosophila melanogaster, an insect with a derived mode of embryogenesis compared to Tribolium, anterior-posterior patterning genes quickly and simultaneously arise as mature gene expression domains that, afterwards, undergo slight posterior-to-anterior shifts. This raises the question of how a fast and simultaneous mode of patterning, like that of Drosophila, could have evolved from a rather slow sequential mode of patterning, like that of Tribolium. In this paper, we propose a mechanism for this evolutionary transition based on a switch from a uniform to a gradient-mediated initialization of the gap gene cascade by maternal Hb. The model is supported by computational analyses and experiments.

φ-evo: A program to evolve phenotypic models of biological networks

Molecular networks are at the core of most cellular decisions, but are often difficult to comprehend. Reverse engineering of network architecture from their functions has proved fruitful to classify and predict the structure and function of molecular networks, suggesting new experimental tests and biological predictions. We present φ-evo, an open-source program to evolve in silico phenotypic networks performing a given biological function. We include implementations for evolution of biochemical adaptation, adaptive sorting for immune recognition, metazoan development (somitogenesis, hox patterning), as well as Pareto evolution. We detail the program architecture based on C, Python 3, and a Jupyter interface for project configuration and network analysis. We illustrate the predictive power of φ-evo by first recovering the asymmetrical structure of the lac operon regulation from an objective function with symmetrical constraints. Second, we use the problem of hox-like embryonic patterning to show how a single effective fitness can emerge from multi-objective (Pareto) evolution. φ-evo provides an efficient approach and user-friendly interface for the phenotypic prediction of networks and the numerical study of evolution itself.

Modelling Time-Dependent Acquisition of Positional Information

Theoretical and computational modelling are crucial to understand dynamics of embryonic development. In this tutorial chapter, we describe two models of gene networks performing time-dependent acquisition of positional information under control of a dynamic morphogen: a toy-model of a bistable gene under control of a morphogen, allowing for the numerical computation of a simple Waddington's epigenetic landscape, and a recently published model of gap genes in Tribolium under control of multiple enhancers. We present detailed commented implementations of the models using python and jupyter notebooks.

Dynamic patterning by the Drosophila pair-rule network reconciles long-germ and short-germ segmentation

Drosophila segmentation is a well-established paradigm for developmental pattern formation. However, the later stages of segment patterning, regulated by the "pair-rule" genes, are still not well understood at the system level. Building on established genetic interactions, I construct a logical model of the Drosophila pair-rule system that takes into account the demonstrated stage-specific architecture of the pair-rule gene network. Simulation of this model can accurately recapitulate the observed spatiotemporal expression of the pair-rule genes, but only when the system is provided with dynamic "gap" inputs. This result suggests that dynamic shifts of pair-rule stripes are essential for segment patterning in the trunk and provides a functional role for observed posterior-to-anterior gap domain shifts that occur during cellularisation. The model also suggests revised patterning mechanisms for the parasegment boundaries and explains the aetiology of the even-skipped null mutant phenotype. Strikingly, a slightly modified version of the model is able to pattern segments in either simultaneous or sequential modes, depending only on initial conditions. This suggests that fundamentally similar mechanisms may underlie segmentation in short-germ and long-germ arthropods.