Ligands for pheromone-sensing neurons are not conformationally activated odorant binding proteins.

Pheromones form an essential chemical language of intraspecific communication in many animals. How olfactory systems recognize pheromonal signals with both sensitivity and specificity is not well understood. An important in vivo paradigm for this process is the detection mechanism of the sex pheromone (Z)-11-octadecenyl acetate (cis-vaccenyl acetate [cVA]) in Drosophila melanogaster. cVA-evoked neuronal activation requires a secreted odorant binding protein, LUSH, the CD36-related transmembrane protein SNMP, and the odorant receptor OR67d. Crystallographic analysis has revealed that cVA-bound LUSH is conformationally distinct from apo (unliganded) LUSH. Recombinantly expressed mutant versions of LUSH predicted to enhance or diminish these structural changes produce corresponding alterations in spontaneous and/or cVA-evoked activity when infused into olfactory sensilla, leading to a model in which the ligand for pheromone receptors is not free cVA, but LUSH that is "conformationally activated" upon cVA binding. Here we present evidence that contradicts this model. First, we demonstrate that the same LUSH mutants expressed transgenically affect neither basal nor pheromone-evoked activity. Second, we compare the structures of apo LUSH, cVA/LUSH, and complexes of LUSH with non-pheromonal ligands and find no conformational property of cVA/LUSH that can explain its proposed unique activated state. Finally, we show that high concentrations of cVA can induce neuronal activity in the absence of LUSH, but not SNMP or OR67d. Our findings are not consistent with the model that the cVA/LUSH complex acts as the pheromone ligand, and suggest that pheromone molecules alone directly activate neuronal receptors.

The role of Receptor Activator of NFKB Ligand (RANKL) in mammary gland development

Abstract : The female reproductive hormones estrogen, progesterone and prolactin control postnatal breast development and are important to breast carcinogenesis. The mechanisms by which they elicit proliferation and morphogenesis remain poorly understood. Using the mouse as a model to study the molecular mechanisms through which hormones elicit morphogenetic changes in the mammary gland in vivo, we found the Receptor Activator of NFκB Ligand, a Tumor Necrosis Factor family member, to be strongly induced by progesterone. Recent publications suggested that hormone dependant RANKURANK signals are involved in the terminal differentiation of mammary gland alveolar buds into lobulo-alveolar structures competent for lactation. I show that in the absence of epithelial RANKL a distinct earlier stage of mammary gland development, side branch formation, is blocked. RANKL acts as a major mediator downstream of progesterone; it is required for progesterone-induced paracrine proliferation and completely rescues the mutant phenotype when ectopically expressed in progesterone receptor (PR) KO mammary epithelia. RANKL is not required for cell autonomous division of estrogen receptor alpha (ERa) /PR positive cells. Cyclin D1, previously implicated as a mediator of RANKL, is not affected by ablation of RANKL and is not required for RANKL-induced paracrine proliferation but for the cell autonomous proliferation. Gene expression arrays to find specific RANKL downstream targets have identified Id4, ElfS and one secreted metalloprotease (Adamtsl8) as potential candidates validated by Q-RT-PCR. Interestingly, Id4 and Adamtsl8 are expressed by the myoepithelial cells. Their expression additionally coincides with RANKL mRNA expression at mid pregnancy, possibly implying a functional contribution of both genes to RANKL mediated sidebranch formation. ElfS in contrast, is found to be strongly expressed by the end of pregnancy supporting recent findings of a prolactin mediated regulation. As for RANKL, this gene was in particular induced in luminal cells. Taken together, I report that progesterone is the major proliferative stimulus in the adult mammary gland eliciting proliferation of ERaJPR positive cells by a cell autonomous, cyclin D1-dependent and a paracrine RANKL-dependent mechanism. My work moreover suggests, that RANKL acts as a major orchestrator affecting different downstream mediators, through which progesterone exerts its effects concomitantly on different cellular compartments. Résumé : Les hormones sexuelles telles que l'oestrogène, la progestérone et la prolactine contrôlent le développement postnatal du sein et sont impliquées dans la cazcinogenèse. Les mécanismes par lesquels elles induisent la prolifération et la morphogénèse demeurent incompris. En utilisant la souris comme modèle, J'ai trouvé que le ligand activateur du récepteur de NFκB, une protéine de la famille du facteur de nécrose des tumeurs, peut être fortement induit par la progestérone. Les publications récentes ont suggéré que cette protéine est nécessaire à la fin de la grossesse, quand les cellules sécrétrices du lait apparaissent. Par des techniques de transplantation d'épithélium, je montre contrairement aux études précédentes, qu'en l'absence de RANKL dans l'épithélium une partie distincte du développement mammaire, la formation de branches latérales, est bloquée. La progestérone agit de manière pazacrine par l'intermédiaire de 12ANKL pour induire la prolifération tandis que la mort cellulaire n'est pas affectée. De plus, l'injection d'une protéine recombinante RANKL dans une souris mutante pour le récepteur à la progestérone induit la prolifération des cellules épithéliales en l'absence de grossesse ; la surexpression de RANKL dans ces mêmes mutants mène à une réversion complète du phénotype. Mes expériences démontrent que la progestérone induit deux types distincts de prolifération. Un premier type direct dans laquelle les cellules positives au récepteur à la progestérone prolifèrent. Cette division cellulaire est alors indépendante de RANKL mais dépendante de la cycline D1. Le second type de prolifération est induit par un mécanisme pazacrine et dépend de RANKL mais pas de la cycline D1. Ici, les cellules négatives au récepteur à la progestérone prolifèrent. Pour détecter des gènes cibles de la voie de signalisation du RANKL, un profil d'expression des gènes a été généré. Les facteurs de transcription Id4, EIf5 et une métalloprotéase sécrétée (Adamtsl8) ont été identifiés en tant que cibles potentielles. D'autres analyses de validation démontrent qu'Id4, Adamtsl8, RANKL mais pas E1f5 sont fortement exprimés au cours de la grossesse, coïncidant avec la formation de branchements latéraux induit par progestérone. EIf5 s'est avéré être exprimé vers la fin de la grossesse appuyant des résultats récents proposant une régulation par la prolactine. Le système canalaire mammaire se compose de couches cellulaires: une couche interne de cellules luminales et une externe de cellules myoépithéliale. Les expériences génétiques d'expression ont révélé que RANKL. et E1f5 sont exprimés dans la partie luminale tandis qu'Id4 et Adamtsl8 sont dans les cellules myoépithéliales. En conclusion, je prouve que la progestérone est le stimulus principal induisant la prolifération dans la glande mammaire d'adulte. Deux mécanismes de prolifération sont impliqués: l'un direct dépendant de la cycline Dl et l'autre paracrine dépendant de RANKI.. Mon travail suggère par ailleurs que RANKL agit en tant que médiateur important, par lequel la progestérone exerce ses effets sur différents compartiments cellulaires tels que la coordination de la prolifération des cellules épithéliales avec la réorganisation de la matrice extracellulaire et de la membrane basale exigées pour la morphogénèse du système canalaire latéral.

Bacterial variations on the methionine salvage pathway.

BACKGROUND: The thiomethyl group of S-adenosylmethionine is often recycled as methionine from methylthioadenosine. The corresponding pathway has been unravelled in Bacillus subtilis. However methylthioadenosine is subjected to alternative degradative pathways depending on the organism. RESULTS: This work uses genome in silico analysis to propose methionine salvage pathways for Klebsiella pneumoniae, Leptospira interrogans, Thermoanaerobacter tengcongensis and Xylella fastidiosa. Experiments performed with mutants of B. subtilis and Pseudomonas aeruginosa substantiate the hypotheses proposed. The enzymes that catalyze the reactions are recruited from a variety of origins. The first, ubiquitous, enzyme of the pathway, MtnA (methylthioribose-1-phosphate isomerase), belongs to a family of proteins related to eukaryotic intiation factor 2B alpha. mtnB codes for a methylthioribulose-1-phosphate dehydratase. Two reactions follow, that of an enolase and that of a phosphatase. While in B. subtilis this is performed by two distinct polypeptides, in the other organisms analyzed here an enolase-phosphatase yields 1,2-dihydroxy-3-keto-5-methylthiopentene. In the presence of dioxygen an aci-reductone dioxygenase yields the immediate precursor of methionine, ketomethylthiobutyrate. Under some conditions this enzyme produces carbon monoxide in B. subtilis, suggesting a route for a new gaseous mediator in bacteria. Ketomethylthiobutyrate is finally transaminated by an aminotransferase that exists usually as a broad specificity enzyme (often able to transaminate aromatic aminoacid keto-acid precursors or histidinol-phosphate). CONCLUSION: A functional methionine salvage pathway was experimentally demonstrated, for the first time, in P. aeruginosa. Apparently, methionine salvage pathways are frequent in Bacteria (and in Eukarya), with recruitment of different polypeptides to perform the needed reactions (an ancestor of a translation initiation factor and RuBisCO, as an enolase, in some Firmicutes). Many are highly dependent on the presence of oxygen, suggesting that the ecological niche may play an important role for the existence and/or metabolic steps of the pathway, even in phylogenetically related bacteria. Further work is needed to uncover the corresponding steps when dioxygen is scarce or absent (this is important to explore the presence of the pathway in Archaea). The thermophile T. tengcongensis, that thrives in the absence of oxygen, appears to possess the pathway. It will be an interesting link to uncover the missing reactions in anaerobic environments.

Robustness of Drosophila early embryo and wing imaginal disc development

Within a developing organism, cells require information on where they are in order to differentiate into the correct cell-type. Pattern formation is the process by which cells acquire and process positional cues and thus determine their fate. This can be achieved by the production and release of a diffusible signaling molecule, called a morphogen, which forms a concentration gradient: exposure to different morphogen levels leads to the activation of specific signaling pathways. Thus, in response to the morphogen gradient, cells start to express different sets of genes, forming domains characterized by a unique combination of differentially expressed genes. As a result, a pattern of cell fates and specification emerges.Though morphogens have been known for decades, it is not yet clear how these gradients form and are interpreted in order to yield highly robust patterns of gene expression. During my PhD thesis, I investigated the properties of Bicoid (Bcd) and Decapentaplegic (Dpp), two morphogens involved in the patterning of the anterior-posterior axis of Drosophila embryo and wing primordium, respectively. In particular, I have been interested in understanding how the pattern proportions are maintained across embryos of different sizes or within a growing tissue. This property is commonly referred to as scaling and is essential for yielding functional organs or organisms. In order to tackle these questions, I analysed fluorescence images showing the pattern of gene expression domains in the early embryo and wing imaginal disc. After characterizing the extent of these domains in a quantitative and systematic manner, I introduced and applied a new scaling measure in order to assess how well proportions are maintained. I found that scaling emerged as a universal property both in early embryos (at least far away from the Bcd source) and in wing imaginal discs (across different developmental stages). Since we were also interested in understanding the mechanisms underlying scaling and how it is transmitted from the morphogen to the target genes down in the signaling cascade, I also quantified scaling in mutant flies where this property could be disrupted. While scaling is largely conserved in embryos with altered bcd dosage, my modeling suggests that Bcd trapping by the nuclei as well as pre-steady state decoding of the morphogen gradient are essential to ensure precise and scaled patterning of the Bcd signaling cascade. In the wing imaginal disc, it appears that as the disc grows, the Dpp response expands and scales with the tissue size. Interestingly, scaling is not perfect at all positions in the field. The scaling of the target gene domains is best where they have a function; Spalt, for example, scales best at the position in the anterior compartment where it helps to form one of the anterior veins of the wing. Analysis of mutants for pentagone, a transcriptional target of Dpp that encodes a secreted feedback regulator of the pathway, indicates that Pentagone plays a key role in scaling the Dpp gradient activity.

Molecular analysis of a novel gene cluster encoding an insect toxin in plant-associated strains of Pseudomonas fluorescens

Pseudomonas fluorescens CHA0 and the related strain Pf-5 are well-characterized representatives of rhizosphere bacteria that have the capacity to protect crop plants from fungal root diseases, mainly by releasing a variety of exoproducts that are toxic to plant pathogenic fungi. Here, we report that the two plant-beneficial pseudomonads also exhibit potent insecticidal activity. Anti-insect activity is linked to a novel genomic locus encoding a large protein toxin termed Fit (for P. fluorescensinsecticidal toxin) that is related to the insect toxin Mcf (Makes caterpillars floppy) of the entomopathogen Photorhabdus luminescens, a mutualist of insect-invading nematodes. When injected into the haemocoel, even low doses of P. fluorescens CHA0 or Pf-5 killed larvae of the tobacco hornworm Manduca sexta and the greater wax moth Galleria mellonella. In contrast, mutants of CHA0 or Pf-5 with deletions in the Fit toxin gene were significantly less virulent to the larvae. When expressed from an inducible promoter in a non-toxic Escherichia coli host, the Fit toxin gene was sufficient to render the bacterium toxic to both insect hosts. Our findings establish the Fit gene products of P. fluorescens CHA0 and Pf-5 as potent insect toxins that define previously unappreciated anti-insect properties of these plant-colonizing bacteria

Comparative genome analysis of Pseudomonas knackmussii B13, the first bacterium known to degrade chloroaromatic compounds.

Pseudomonas knackmussii B13 was the first strain to be isolated in 1974 that could degrade chlorinated aromatic hydrocarbons. This discovery was the prologue for subsequent characterization of numerous bacterial metabolic pathways, for genetic and biochemical studies, and which spurred ideas for pollutant bioremediation. In this study, we determined the complete genome sequence of B13 using next generation sequencing technologies and optical mapping. Genome annotation indicated that B13 has a variety of metabolic pathways for degrading monoaromatic hydrocarbons including chlorobenzoate, aminophenol, anthranilate and hydroxyquinol, but not polyaromatic compounds. Comparative genome analysis revealed that B13 is closest to Pseudomonas denitrificans and Pseudomonas aeruginosa. The B13 genome contains at least eight genomic islands [prophages and integrative conjugative elements (ICEs)], which were absent in closely related pseudomonads. We confirm that two ICEs are identical copies of the 103 kb self-transmissible element ICEclc that carries the genes for chlorocatechol metabolism. Comparison of ICEclc showed that it is composed of a variable and a 'core' region, which is very conserved among proteobacterial genomes, suggesting a widely distributed family of so far uncharacterized ICE. Resequencing of two spontaneous B13 mutants revealed a number of single nucleotide substitutions, as well as excision of a large 220 kb region and a prophage that drastically change the host metabolic capacity and survivability.

Characterization of Fas (Apo-1, CD95)-Fas ligand interaction.

The death-inducing receptor Fas is activated when cross-linked by the type II membrane protein Fas ligand (FasL). When human soluble FasL (sFasL, containing the extracellular portion) was expressed in human embryo kidney 293 cells, the three N-linked glycans of each FasL monomer were found to be essential for efficient secretion. Based on the structure of the closely related lymphotoxin alpha-tumor necrosis factor receptor I complex, a molecular model of the FasL homotrimer bound to three Fas molecules was generated using knowledge-based protein modeling methods. Point mutations of amino acid residues predicted to affect the receptor-ligand interaction were introduced at three sites. The F275L mutant, mimicking the loss of function murine gld mutation, exhibited a high propensity for aggregation and was unable to bind to Fas. Mutants P206R, P206D, and P206F displayed reduced cytotoxicity toward Fas-positive cells with a concomitant decrease in the binding affinity for the recombinant Fas-immunoglobulin Fc fusion proteins. Although the cytotoxic activity of mutant Y218D was unaltered, mutant Y218R was inactive, correlating with the prediction that Tyr-218 of FasL interacts with a cluster of three basic amino acid side chains of Fas. Interestingly, mutant Y218F could induce apoptosis in murine, but not human cells.

Aging preferentially affects molecular pathways implicated in development and disease of myelinating glial cells

The central and peripheral nervous systems are involved in multiple age-dependent neurological deficits that are often attributed to alterations in function of myelinating glial cells. However, the molecular events that underlie the age-related decline of glial cell function are unknown. We used Schwann cells as a model to study biological processes affected in glial cells by aging. We comprehensively profiled gene expression of the Schwann cellrich mouse sciatic nerve throughout life, from day of birth until senescence (840 days of age). We combined the aging data with the microarray transcriptional data obtained using nerves isolated from Schwann cell-specific neuropathy-inducing mutants MPZCre/+/Lpin1fE2−3/fE2−3 , MPZCre/+/ScapfE1/fE1 and Pmp22-null mice. The majority of age related transcripts were also affected in the analyzed mouse models of neuropathy (54.4%) and in development (59.5%) indicating a high level of overlapping in implicated molecular pathways. We observed that compared to peripheral nerve development, dynamically changing expression profiles in aging have opposite (anticorrelated) orientation while they copy the orientation of transcriptional changes observed in analyzed neuropathy models. Subsequent clustering and biological annotation of dynamically changing transcripts revealed that the processes most significantly deregulated in aging include inflammatory/immune response and lipid biosynthesis/metabolism. Importantly, the changes in these pathways were also observed in myelinated oligodendrocyte-rich optic nerves of aged mice, albeit with lower magnitude. This observation suggests that similar biological processes are affected in aging glial cells in central and peripheral nervous systems, however with different dynamics. Our data, which provide the first comprehensive comparison of molecular changes in glial cells in three distinct biological conditions comprising development, aging and disease, provide not only a new inside into the molecular alterations underlying neural system aging but also identify target pathways for potential therapeutic approaches to prevent or delay complications associated with age-related and inherited forms of neuropathies. *Current address: Department of Physiology, UCSF, San Francisco, CA, USA.

Identification of genes potentially involved in solute stress response in Sphingomonas wittichii RW1 by transposon mutant recovery.

The term water stress refers to the effects of low water availability on microbial growth and physiology. Water availability has been proposed as a major constraint for the use of microorganisms in contaminated sites with the purpose of bioremediation. Sphingomonas wittichii RW1 is a bacterium capable of degrading the xenobiotic compounds dibenzofuran and dibenzo-p-dioxin, and has potential to be used for targeted bioremediation. The aim of the current work was to identify genes implicated in water stress in RW1 by means of transposon mutagenesis and mutant growth experiments. Conditions of low water potential were mimicked by adding NaCl to the growth media. Three different mutant selection or separation method were tested which, however recovered different mutants. Recovered transposon mutants with poorer growth under salt-induced water stress carried insertions in genes involved in proline and glutamate biosynthesis, and further in a gene putatively involved in aromatic compound catabolism. Transposon mutants growing poorer on medium with lowered water potential also included ones that had insertions in genes involved in more general functions such as transcriptional regulation, elongation factor, cell division protein, RNA polymerase β or an aconitase.

Contributions of the peroxisome and β-oxidation cycle to biotin synthesis in fungi.

The first step in the synthesis of the bicyclic rings of D-biotin is mediated by 8-amino-7-oxononanoate (AON) synthase, which catalyzes the decarboxylative condensation of l-alanine and pimelate thioester. We found that the Aspergillus nidulans AON synthase, encoded by the bioF gene, is a peroxisomal enzyme with a type 1 peroxisomal targeting sequence (PTS1). Localization of AON to the peroxisome was essential for biotin synthesis because expression of a cytosolic AON variant or deletion of pexE, encoding the PTS1 receptor, rendered A. nidulans a biotin auxotroph. AON synthases with PTS1 are found throughout the fungal kingdom, in ascomycetes, basidiomycetes, and members of basal fungal lineages but not in representatives of the Saccharomyces species complex, including Saccharomyces cerevisiae. A. nidulans mutants defective in the peroxisomal acyl-CoA oxidase AoxA or the multifunctional protein FoxA showed a strong decrease in colonial growth rate in biotin-deficient medium, whereas partial growth recovery occurred with pimelic acid supplementation. These results indicate that pimeloyl-CoA is the in vivo substrate of AON synthase and that it is generated in the peroxisome via the β-oxidation cycle in A. nidulans and probably in a broad range of fungi. However, the β-oxidation cycle is not essential for biotin synthesis in S. cerevisiae or Escherichia coli. These results suggest that alternative pathways for synthesis of the pimelate intermediate exist in bacteria and eukaryotes and that Saccharomyces species use a pathway different from that used by the majority of fungi.

Identification of a Poor-Prognosis BRAF-Mutant-Like Population of Patients With Colon Cancer.

PURPOSE Our purpose was development and assessment of a BRAF-mutant gene expression signature for colon cancer (CC) and the study of its prognostic implications. Materials and METHODS A set of 668 stage II and III CC samples from the PETACC-3 (Pan-European Trails in Alimentary Tract Cancers) clinical trial were used to assess differential gene expression between c.1799T>A (p.V600E) BRAF mutant and non-BRAF, non-KRAS mutant cancers (double wild type) and to construct a gene expression-based classifier for detecting BRAF mutant samples with high sensitivity. The classifier was validated in independent data sets, and survival rates were compared between classifier positive and negative tumors. Results A 64 gene-based classifier was developed with 96% sensitivity and 86% specificity for detecting BRAF mutant tumors in PETACC-3 and independent samples. A subpopulation of BRAF wild-type patients (30% of KRAS mutants, 13% of double wild type) showed a gene expression pattern and had poor overall survival and survival after relapse, similar to those observed in BRAF-mutant patients. Thus they form a distinct prognostic subgroup within their mutation class. CONCLUSION A characteristic pattern of gene expression is associated with and accurately predicts BRAF mutation status and, in addition, identifies a population of BRAF mutated-like KRAS mutants and double wild-type patients with similarly poor prognosis. This suggests a common biology between these tumors and provides a novel classification tool for cancers, adding prognostic and biologic information that is not captured by the mutation status alone. These results may guide therapeutic strategies for this patient segment and may help in population stratification for clinical trials.

Anchoring of both PKA and 14-3-3 inhibits the Rho-GEF activity of the AKAP-Lbc signaling complex.

A-kinase anchoring proteins (AKAPs) target the cAMP-regulated protein kinase (PKA) to its physiological substrates. We recently identified a novel anchoring protein, called AKAP-Lbc, which functions as a PKA-targeting protein as well as a guanine nucleotide exchange factor (GEF) for RhoA. We demonstrated that AKAP-Lbc Rho-GEF activity is stimulated by the alpha subunit of the heterotrimeric G protein G12. Here, we identified 14-3-3 as a novel regulatory protein interacting with AKAP-Lbc. Elevation of the cellular concentration of cAMP activates the PKA holoenzyme anchored to AKAP-Lbc, which phosphorylates the anchoring protein on the serine 1565. This phosphorylation event induces the recruitment of 14-3-3, which inhibits the Rho-GEF activity of AKAP-Lbc. AKAP-Lbc mutants that fail to interact with PKA or with 14-3-3 show a higher basal Rho-GEF activity as compared to the wild-type protein. This suggests that, under basal conditions, 14-3-3 maintains AKAP-Lbc in an inactive state. Therefore, while it is known that AKAP-Lbc activity can be stimulated by Galpha12, in this study we demonstrated that it is inhibited by the anchoring of both PKA and 14-3-3.

Aging preferentially affects molecular pathways implicated in development and disease of myelinating glial cells

The central and peripheral nervous systems are involved in multiple agedependent neurological deficits that are often attributed to alterations in function of myelinating glial cells. However, the molecular events that underlie the age-related decline of glial cell function are unknown. We used Schwann cells as a model to study biological processes affected in glial cells by aging. We comprehensively profiled gene expression of the Schwann cell-rich mouse sciatic nerve throughout life, from day of birth until senescence (840 days of age). We combined the aging data with the microarray transcriptional data obtained using nerves isolated from Schwann cell-specific neuropathy-inducing mutants MPZCre/þ/Lpin1fE2-3/fE2-3, MPZCre/þ/ScapfE1/fE1 and Pmp22-null mice. A majority of age related transcripts were also affected in the analyzed mouse models of neuropathy (54.4%) and in development (59.5%) indicating a high level of overlapping in implicated molecular pathways. We observed that compared to peripheral nerve development, dynamically changing expression profiles in aging have opposite (anticorrelated) orientation while they copy the orientation of transcriptional changes observed in analyzed neuropathy models. Subsequent clustering and biological annotation of dynamically changing transcripts revealed that the processes most significantly deregulated in aging include inflammatory/ immune response and lipid biosynthesis/metabolism. Importantly, the changes in these pathways were also observed in myelinated oligodendrocyte- rich optic nerves of aged mice, albeit with lower magnitude. This observation suggests that similar biological processes are affected in aging glial cells in central and peripheral nervous systems, however with different dynamics. Our data, which provide the first comprehensive comparison of molecular changes in glial cells in three distinct biological conditions comprising development, aging and disease, provide not only a new inside into the molecular alterations underlying neural system aging but also identify target pathways for potential therapeutical approaches to prevent or delay complications associated with age-related and inherited forms of neuropathies.

Development and application of computational methodologies in qualitative modeling

La présente thèse s'intitule "Développent et Application des Méthodologies Computationnelles pour la Modélisation Qualitative". Elle comprend tous les différents projets que j'ai entrepris en tant que doctorante. Plutôt qu'une mise en oeuvre systématique d'un cadre défini a priori, cette thèse devrait être considérée comme une exploration des méthodes qui peuvent nous aider à déduire le plan de processus regulatoires et de signalisation. Cette exploration a été mue par des questions biologiques concrètes, plutôt que par des investigations théoriques. Bien que tous les projets aient inclus des systèmes divergents (réseaux régulateurs de gènes du cycle cellulaire, réseaux de signalisation de cellules pulmonaires) ainsi que des organismes (levure à fission, levure bourgeonnante, rat, humain), nos objectifs étaient complémentaires et cohérents. Le projet principal de la thèse est la modélisation du réseau de l'initiation de septation (SIN) du S.pombe. La cytokinèse dans la levure à fission est contrôlée par le SIN, un réseau signalant de protéines kinases qui utilise le corps à pôle-fuseau comme échafaudage. Afin de décrire le comportement qualitatif du système et prédire des comportements mutants inconnus, nous avons décidé d'adopter l'approche de la modélisation booléenne. Dans cette thèse, nous présentons la construction d'un modèle booléen étendu du SIN, comprenant la plupart des composantes et des régulateurs du SIN en tant que noeuds individuels et testable expérimentalement. Ce modèle utilise des niveaux d'activité du CDK comme noeuds de contrôle pour la simulation d'évènements du SIN à différents stades du cycle cellulaire. Ce modèle a été optimisé en utilisant des expériences d'un seul "knock-out" avec des effets phénotypiques connus comme set d'entraînement. Il a permis de prédire correctement un set d'évaluation de "knock-out" doubles. De plus, le modèle a fait des prédictions in silico qui ont été validées in vivo, permettant d'obtenir de nouvelles idées de la régulation et l'organisation hiérarchique du SIN. Un autre projet concernant le cycle cellulaire qui fait partie de cette thèse a été la construction d'un modèle qualitatif et minimal de la réciprocité des cyclines dans la S.cerevisiae. Les protéines Clb dans la levure bourgeonnante présentent une activation et une dégradation caractéristique et séquentielle durant le cycle cellulaire, qu'on appelle communément les vagues des Clbs. Cet évènement est coordonné avec la courbe d'activation inverse du Sic1, qui a un rôle inhibitoire dans le système. Pour l'identification des modèles qualitatifs minimaux qui peuvent expliquer ce phénomène, nous avons sélectionné des expériences bien définies et construit tous les modèles minimaux possibles qui, une fois simulés, reproduisent les résultats attendus. Les modèles ont été filtrés en utilisant des simulations ODE qualitatives et standardisées; seules celles qui reproduisaient le phénotype des vagues ont été gardées. L'ensemble des modèles minimaux peut être utilisé pour suggérer des relations regulatoires entre les molécules participant qui peuvent ensuite être testées expérimentalement. Enfin, durant mon doctorat, j'ai participé au SBV Improver Challenge. Le but était de déduire des réseaux spécifiques à des espèces (humain et rat) en utilisant des données de phosphoprotéines, d'expressions des gènes et des cytokines, ainsi qu'un réseau de référence, qui était mis à disposition comme donnée préalable. Notre solution pour ce concours a pris la troisième place. L'approche utilisée est expliquée en détail dans le dernier chapitre de la thèse. -- The present dissertation is entitled "Development and Application of Computational Methodologies in Qualitative Modeling". It encompasses the diverse projects that were undertaken during my time as a PhD student. Instead of a systematic implementation of a framework defined a priori, this thesis should be considered as an exploration of the methods that can help us infer the blueprint of regulatory and signaling processes. This exploration was driven by concrete biological questions, rather than theoretical investigation. Even though the projects involved divergent systems (gene regulatory networks of cell cycle, signaling networks in lung cells), as well as organisms (fission yeast, budding yeast, rat, human), our goals were complementary and coherent. The main project of the thesis is the modeling of the Septation Initiation Network (SIN) in S.pombe. Cytokinesis in fission yeast is controlled by the SIN, a protein kinase signaling network that uses the spindle pole body as scaffold. In order to describe the qualitative behavior of the system and predict unknown mutant behaviors we decided to adopt a Boolean modeling approach. In this thesis, we report the construction of an extended, Boolean model of the SIN, comprising most SIN components and regulators as individual, experimentally testable nodes. The model uses CDK activity levels as control nodes for the simulation of SIN related events in different stages of the cell cycle. The model was optimized using single knock-out experiments of known phenotypic effect as a training set, and was able to correctly predict a double knock-out test set. Moreover, the model has made in silico predictions that have been validated in vivo, providing new insights into the regulation and hierarchical organization of the SIN. Another cell cycle related project that is part of this thesis was to create a qualitative, minimal model of cyclin interplay in S.cerevisiae. CLB proteins in budding yeast present a characteristic, sequential activation and decay during the cell cycle, commonly referred to as Clb waves. This event is coordinated with the inverse activation curve of Sic1, which has an inhibitory role in the system. To generate minimal qualitative models that can explain this phenomenon, we selected well-defined experiments and constructed all possible minimal models that, when simulated, reproduce the expected results. The models were filtered using standardized qualitative ODE simulations; only the ones reproducing the wave-like phenotype were kept. The set of minimal models can be used to suggest regulatory relations among the participating molecules, which will subsequently be tested experimentally. Finally, during my PhD I participated in the SBV Improver Challenge. The goal was to infer species-specific (human and rat) networks, using phosphoprotein, gene expression and cytokine data and a reference network provided as prior knowledge. Our solution to the challenge was selected as in the final chapter of the thesis.

Primary root protophloem differentiation requires balanced phosphatidylinositol-4,5-biphosphate levels and systemically affects root branching.

Protophloem is a specialized vascular tissue in growing plant organs, such as root meristems. In Arabidopsis mutants with impaired primary root protophloem differentiation, brevis radix (brx) and octopus (ops), meristematic activity and consequently overall root growth are strongly reduced. Second site mutation in the protophloem-specific presumed phosphoinositide 5-phosphatase COTYLEDON VASCULAR PATTERN 2 (CVP2), but not in its homolog CVP2-LIKE 1 (CVL1), partially rescues brx defects. Consistent with this finding, CVP2 hyperactivity in a wild-type background recreates a brx phenotype. Paradoxically, however, while cvp2 or cvl1 single mutants display no apparent root defects, the root phenotype of cvp2 cvl1 double mutants is similar to brx or ops, although, as expected, cvp2 cvl1 seedlings contain more phosphatidylinositol-4,5-biphosphate. Thus, tightly balanced phosphatidylinositol-4,5-biphosphate levels appear essential for proper protophloem differentiation. Genetically, OPS acts downstream of phosphatidylinositol-4,5-biphosphate levels, as cvp2 mutation cannot rescue ops defects, whereas increased OPS dose rescues cvp2 cvl1 defects. Finally, all three mutants display higher density and accelerated emergence of lateral roots, which correlates with increased auxin response in the root differentiation zone. This phenotype is also created by application of peptides that suppress protophloem differentiation, CLAVATA3/EMBRYO SURROUNDING REGION 26 (CLE26) and CLE45. Thus, local changes in the primary root protophloem systemically shape overall root system architecture.