Multi-level models and infrastructures for simulating biological system development

The hierarchical organisation of biological systems plays a crucial role in the pattern formation of gene expression resulting from the morphogenetic processes, where autonomous internal dynamics of cells, as well as cell-to-cell interactions through membranes, are responsible for the emergent peculiar structures of the individual phenotype. Being able to reproduce the systems dynamics at different levels of such a hierarchy might be very useful for studying such a complex phenomenon of self-organisation. The idea is to model the phenomenon in terms of a large and dynamic network of compartments, where the interplay between inter-compartment and intra-compartment events determines the emergent behaviour resulting in the formation of spatial patterns. According to these premises the thesis proposes a review of the different approaches already developed in modelling developmental biology problems, as well as the main models and infrastructures available in literature for modelling biological systems, analysing their capabilities in tackling multi-compartment / multi-level models. The thesis then introduces a practical framework, MS-BioNET, for modelling and simulating these scenarios exploiting the potential of multi-level dynamics. This is based on (i) a computational model featuring networks of compartments and an enhanced model of chemical reaction addressing molecule transfer, (ii) a logic-oriented language to flexibly specify complex simulation scenarios, and (iii) a simulation engine based on the many-species/many-channels optimised version of Gillespie’s direct method. The thesis finally proposes the adoption of the agent-based model as an approach capable of capture multi-level dynamics. To overcome the problem of parameter tuning in the model, the simulators are supplied with a module for parameter optimisation. The task is defined as an optimisation problem over the parameter space in which the objective function to be minimised is the distance between the output of the simulator and a target one. The problem is tackled with a metaheuristic algorithm. As an example of application of the MS-BioNET framework and of the agent-based model, a model of the first stages of Drosophila Melanogaster development is realised. The model goal is to generate the early spatial pattern of gap gene expression. The correctness of the models is shown comparing the simulation results with real data of gene expression with spatial and temporal resolution, acquired in free on-line sources.

Genetic Manipulation, Gene Silencing and Biomarker Development in Multiple Experimental Models.

The thesis is set in three different parts, according to the relative experimental models. First, the domestic pig (Sus scrofa) is part of the study on reproductive biotechnologies: the transgenesis technique of Sperm Mediated Gene Transfer is widely studied starting from the quality of the semen, through the study of multiple uptakes of exogenous DNA and lastly used in the production of multi-transgenic blastocysts. Finally we managed to couple the transgenesis pipeline with sperm sorting and therefore produced transgenic embryos of predetermined sex. In the second part of the thesis the attention is on the fruit fly (Drosophila melanogaster) and on its derived cell line: the S2 cells. The in vitro and in vivo models are used to develop and validate an efficient way to knock down the myc gene. First an efficient in vitro protocol is described, than we demonstrate how the decrease in myc transcript remarkably affects the ribosome biogenesis through the study of Polysome gradients, rRNA content and qPCR. In vivo we identified two optimal drivers for the conditional silencing of myc, once the flies are fed with RU486: the first one is throughout the whole body (Tubulin), while the second is a head fat body driver (S32). With these results we present a very efficient model to study the role of myc in multiple aspects of translation. In the third and last part, the focus is on human derived lung fibroblasts (hLF-1), mouse tail fibroblasts and mouse tissues. We developed an efficient assay to quantify the total protein content of the nucleus on a single cell level via fluorescence. We coupled the protocol with classical immunofluorescence so to have at the same time general and particular information, demonstrating that during senescence nuclear proteins increase by 1.8 fold either in human cells, mouse cells and mouse tissues.

Hox genes and the regulation of programmed cell death in the embryonic central nervous system of Drosophila melanogaster

This study deals with the function and regulation of programmed cell death, or apoptosis, in the development of the embryonic central nervous system of Drosophila melanogaster. The first part provides a description of apoptosis-deficient embryos, which showed that preventing apoptosis does not cause gross morphological defects in the CNS, as it appears well organized despite the presence of too many cells. An analysis of the incidence and pattern of apoptosis over the course of development discloses a partly very orderly pattern suggesting tight spatio-temporal control, but also reveals random apoptotic cells, which suggests a certain amount of plasticity in the embryo. This analysis also allowed precise identification of some of the dying neural cells in the embryo, and establishment of single cell models for studying regulation of segment-specific apoptosis in the embryonic CNS. In the second part of the work, further investigations into mechanisms controlling segment-specific apoptosis revealed the involvement of two Hox genes, Antennapedia (Antp) and Ultrabithorax (Ubx), in this process. Hox genes control the formation of segment-specific structures in their domains of expression, but also regulate organ and tissue morphogenesis. The study presented here shows that Antp and Ubx play antagonistic roles in motoneuron survival in the embryo. Ubx expression in the CNS is strongly upregulated at a late point in development, when most cells have begun to differentiate. This upregulation shortly precedes Ubx-dependent, segment-specific apoptosis of two differentiated motoneurons. It could further be demonstrated that Antp is required for proper development of the NB7-3 lineage and for survival of the NB7-3 motoneuron in the anterior thoracic segments. In segments where Antp and Ubx expression overlaps, Ubx counteracts the anti-apoptotic function of Antp, resulting in cell death. Thus, these two Hox genes play opposing roles in the survival of differentiated neurons in the late developing nervous system. They thereby contribute to establishment of correct connections between outward-projecting neurons and their targets, which is crucial for the assembly of functional neural circuits, as these have to fulfill region-specific locomotion and sensory requirements along the antero-posterior body axis.

Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

A Drosophila XPD model links cell cycle coordination with neuro-development and suggests links to cancer

XPD functions in transcription, DNA repair and in cell cycle control. Mutations in human XPD (also known as ERCC2) mainly cause three clinical phenotypes: xeroderma pigmentosum (XP), Cockayne syndrome (XP/CS) and trichothiodystrophy (TTD), and only XP patients have a high predisposition to developing cancer. Hence, we developed a fly model to obtain novel insights into the defects caused by individual hypomorphic alleles identified in human XP-D patients. This model revealed that the mutations that displayed the greatest in vivo UV sensitivity in Drosophila did not correlate with those that led to tumor formation in humans. Immunoprecipitations followed by targeted quantitative MS/MS analysis showed how different xpd mutations affected the formation or stability of different transcription factor IIH (TFIIH) subcomplexes. The XP mutants most clearly linked to high cancer risk, Xpd R683W and R601L, showed a reduced interaction with the core TFIIH and also an abnormal interaction with the Cdk-activating kinase (CAK) complex. Interestingly, these two XP alleles additionally displayed high levels of chromatin loss and free centrosomes during the rapid nuclear division phase of the Drosophila embryo. Finally, the xpd mutations showing defects in the coordination of cell cycle timing during the Drosophila embryonic divisions correlated with those human mutations that cause the neurodevelopmental abnormalities and developmental growth defects observed in XP/CS and TTD patients.

Immune-Deficient Drosophila melanogaster: A Model for the Innate Immune Response to Human Fungal Pathogens

We explored the host-pathogen interactions of the human opportunistic fungus Candida albicans using Drosophila melanogaster. We established that a Drosophila strain devoid of functional Toll receptor is highly susceptible to the human pathogen C. albicans. Using this sensitive strain, we have been able to show that a set of specific C. albicans mutants of different virulence in mammalian infection models are also impaired in virulence in Drosophila and remarkably display the same rank order of virulence. This immunodeficient insect model also revealed virulence properties undetected in an immunocompetent murine model of infection. The genetic systems available in both host and pathogen will enable the identification of host-specific components and C. albicans genes involved in the host-fungal interplay.

STUDYING AGGREGATE FORMATION BY AMYOTROPHIC LATERAL SCLEROSIS-ASSOCIATED MUTANT SOD1 PROTEIN IN DROSOPHILA MODEL

A common pathological hallmark of most neurodegenerative disorders is the presence of protein aggregates in the brain. Understanding the regulation of aggregate formation is thus important for elucidating disease pathogenic mechanisms and finding effective preventive avenues and cures. Amyotrophic Lateral Sclerosis (ALS), also known as Lou Gehrig’s disease, is a selective neurodegenerative disorder predominantly affecting motor neurons. The majority of ALS cases are sporadic, however, mutations in superoxide dismutase 1 (SOD1) are responsible for about 20% of familial ALS (fALS). Mutated SOD1 proteins are prone to misfold and form protein aggregates, thus representing a good candidate for studying aggregate formation. The long-term goal of this project is to identify regulators of aggregate formation by mutant SOD1 and other ALS-associated disease proteins. The specific aim of this thesis project is to assess the possibility of using the well-established Drosophila model system to study aggregation by human SOD1 (hSOD1) mutants. To this end, using wild type and the three mutant hSOD1 (A4V, G85R and G93A) most commonly found among fALS, I have generated 16 different SOD1 constructs containing either eGFP or mCherry in-frame fluorescent reporters, established and tested both cell- and animal-based Drosophila hSOD1 models. The experimental strategy allows for clear visualization of ectopic hSOD1 expression as well as versatile co-expression schemes to fully investigate protein aggregation specifically by mutant hSOD1. I have performed pilot cell-transfection experiments and verified induced expression of hSOD1 proteins. Using several tissue- or cell type-specific Gal4 lines, I have confirmed the proper expression of hSOD1 from established transgenic fly lines. Interestingly, in both Drosophila S2 cells and different fly tissues including the eye and motor neurons, robust aggregate formation by either wild type or mutant hSOD1 proteins was not observed. These preliminary observations suggest that Drosophila might not be a good experimental organism to study aggregation and toxicity of mutant hSOD1 protein. Nevertheless this preliminary conclusion implies the potential existence of a potent protective mechanism against mutant hSOD1 aggregation and toxicity in Drosophila. Thus, results from my SOD1-ALS project in Drosophila will help future studies on how to best employ this classic model organism to study ALS and other human brain degenerative diseases.

A Tropomyosin-2 Mutation Suppresses a Troponin I Myopathy in Drosophila

A suppressor mutation, D53, of the held-up2 allele of the Drosophila melanogaster Troponin I (wupA) gene is described. D53, a missense mutation, S185F, of the tropomyosin-2, Tm2, gene fully suppresses all the phenotypic effects of held-up2, including the destructive hypercontraction of the indirect flight muscles (IFMs), a lack of jumping, the progressive myopathy of the walking muscles, and reductions in larval crawling and feeding behavior. The suppressor restores normal function of the IFMs, but flight ability decreases with age and correlates with an unusual, progressive structural collapse of the myofibrillar lattice starting at the center. The S185F substitution in Tm2 is close to a troponin T binding site on tropomyosin. Models to explain suppression by D53, derived from current knowledge of the vertebrate troponin-tropomyosin complex structure and functions, are discussed. The effects of S185F are compared with those of two mutations in residues 175 and 180 of human α-tropomyosin 1 which cause familial hypertrophic cardiomyopathy (HCM).

Sexually antagonistic coevolution of a postmating-prezygotic reproductive character in desert Drosophila

Rapid divergence in postmating-prezygotic characters suggests that selection may be responsible for generating reproductive barriers between closely related species. Theoretical models indicate that this rapid divergence could be generated by a series of male adaptations and female counteradaptations by means of sexual selection or conflict, but empirical tests of particular mechanisms are generally lacking. Moreover, although a male–female genotypic interaction in mediating sperm competition attests to an active role of females, molecular or morphological evidence of the female's participation in the coevolutionary process is critically needed. Here we show that postmating-prezygotic variation among populations of cactophilic desert Drosophila reflects divergent coevolutionary trajectories between the sexes. We explicitly test the female's role in intersexual interactions by quantifying differences in a specific postmating-prezygotic reproductive character, the insemination reaction mass, in two species, Drosophila mojavensis and Drosophila arizonae. A series of interpopulation crosses confirmed that population divergence was propelled by male–female interactions, a prerequisite if the selective forces derive from sexual conflicts. An association between the reaction mass and remating and oviposition behavior argues that divergence has been propelled by sexually antagonistic coevolution, and potentially has important implications for speciation.

Subunit-destabilizing mutations in Drosophila copper/zinc superoxide dismutase: neuropathology and a model of dimer dysequilibrium.

Mutations in Cu/Zn superoxide dismutase (SOD), a hallmark of familial amyotrophic lateral sclerosis (FALS) in humans, are shown here to confer striking neuropathology in Drosophila. Heterozygotes with one wild-type and one deleted SOD allele retain the expected 50% of normal activity for this dimeric enzyme. However, heterozygotes with one wild-type and one missense SOD allele show lesser SOD activities, ranging from 37% for a heterozygote carrying a missense mutation predicted from structural models to destabilize the dimer interface, to an average of 13% for several heterozygotes carrying missense mutations predicted to destabilize the subunit fold. Genetic and biochemical evidence suggests a model of dimer dysequilibrium whereby SOD activity in missense heterozygotes is reduced through entrapment of wild-type subunits into unstable or enzymatically inactive heterodimers. This dramatic impairment of the activity of wild-type subunits in vivo has implications for our understanding of FALS and for possible therapeutic strategies.

Multivariate quantitative genetics and the Lek Paradox: Genetic variance in male sexually selected traits of Drosophila serrata under field control

Single male sexually selected traits have been found to exhibit substantial genetic variance, even though natural and sexual selection are predicted to deplete genetic variance in these traits. We tested whether genetic variance in multiple male display traits of Drosophila serrata was maintained under field conditions. A breeding design involving 300 field-reared males and their laboratory-reared offspring allowed the estimation of the genetic variance-covariance matrix for six male cuticular hydrocarbons (CHCs) under field conditions. Despite individual CHCs displaying substantial genetic variance under field conditions, the vast majority of genetic variance in CHCs was not closely associated with the direction of sexual selection measured on field phenotypes. Relative concentrations of three CHCs correlated positively with body size in the field, but not under laboratory conditions, suggesting condition-dependent expression of CHCs under field conditions. Therefore condition dependence may not maintain genetic variance in preferred combinations of male CHCs under field conditions, suggesting that the large mutational target supplied by the evolution of condition dependence may not provide a solution to the lek paradox in this species. Sustained sexual selection may be adequate to deplete genetic variance in the direction of selection, perhaps as a consequence of the low rate of favorable mutations expected in multiple trait systems.

Filling in gaps of Drosophila melanogaster urate degradation metabolic pathway using metabolomics approaches: towards the core metabolome of the fruit fly

The primary goal of systems biology is to integrate complex omics data, and data obtained from traditional experimental studies in order to provide a holistic understanding of organismal function. One way of achieving this aim is to generate genome-scale metabolic models (GEMs), which contain information on all metabolites, enzyme-coding genes, and biochemical reactions in a biological system. Drosophila melanogaster GEM has not been reconstructed to date. Constraint-free genome-wide metabolic model of the fruit fly has been reconstructed in our lab, identifying gaps, where no enzyme was identified and metabolites were either only produced or consume. The main focus of the work presented in this thesis was to develop a pipeline for efficient gap filling using metabolomics approaches combined with standard reverse genetics methods, using 5-hydroxyisourate hydrolase (5-HIUH) as an example. 5-HIUH plays a role in urate degradation pathway. Inability to degrade urate can lead to inborn errors of metabolism (IEMs) in humans, including hyperuricemia. Based on sequence analysis Drosophila CG30016 gene was hypothesised to encode 5- HIUH. CG30016 knockout flies were examined to identify Malpighian tubules phenotype, and shortened lifespan might reflect kidney disorders in hyperuricemia in humans. Moreover, LC-MS analysis of mutant tubules revealed that CG30016 is involved in purine metabolism, and specifically urate degradation pathway. However, the exact role of the gene has not been identified, and the complete method for gap filling has not been developed. Nevertheless, thanks to the work presented here, we are a step closer towards the development of a gap-filling pipeline in Drosophila melanogaster GEM. Importantly, the areas that require further optimisation were identified and are the focus of future research. Moreover, LC-MS analysis confirmed that tubules rather than the whole fly were more suitable for metabolomics analysis of purine metabolism. Previously, Dow/Davies lab has generated the most complete tissue-specific transcriptomic atlas for Drosophila – FlyAtlas.org, which provides data on gene expression across multiple tissues of adult fly and larva. FlyAtlas revealed that transcripts of many genes are enriched in specific Drosophila tissues, and that it is possible to deduce the functions of individual tissues within the fly. Based on FlyAtlas data, it has become clear that the fly (like other metazoan species) must be considered as a set of tissues, each 2 with its own distinct transcriptional and functional profile. Moreover, it revealed that for about 30% of the genome, reverse genetic methods (i.e. mutation in an unknown gene followed by observation of phenotype) are only useful if specific tissues are investigated. Based on the FlyAtlas findings, we aimed to build a primary tissue-specific metabolome of the fruit fly, in order to establish whether different Drosophila tissues have different metabolomes and if they correspond to tissue-specific transcriptome of the fruit fly (FlyAtlas.org). Different fly tissues have been dissected and their metabolome elucidated using LC-MS. The results confirmed that tissue metabolomes differ significantly from each other and from the whole fly, and that some of these differences can be correlated to the tissue function. The results illustrate the need to study individual tissues as well as the whole organism. It is clear that some metabolites that play an important role in a given tissue might not be detected in the whole fly sample because their abundance is much lower in comparison to other metabolites present in all tissues, which prevent the detection of the tissue-specific compound.