Evaluation of the immune response to infection with metastatic and non-metastatic "Leishmania guyanensis" parasites

ABSTRACT : Les infections par le parasite Leishmania guyanensis se caractérisent par une dissémination depuis le site initial d'infection jusqu'aux tissus naso-pharyngés, responsable de la Leishmaniose à lésions secondaires muco-cutanées (LMC). Les lésions des patients atteints de LMC montrent une massive infiltration de cellules immunitaires, une réponse immunitaire élevée et la présence de parasites (bien qu'en très faible quantité). La LMC engendre une augmentation de l'expression de TNFa ainsi qu'un défaut dans le contrôle de la réponse immunitaire caractérisé par une absence de réponse à l'IL 10. La réponse immunitaire de l'hôte ainsi que la virulence du parasite sont deux facteurs reconnus pour le contrôle de l'infection. Le mécanisme de la pathogenèse de la LMC restent grandement incompris, surtout le mécanisme de dissémination de l'infection du site d'inoculation jusqu'aux sites secondaires d'infection (métastases) ainsi que les détails de la réponse de l'hôte contre le pathogène. Dans un modèle d'infection d' hamsters avec des parasites du Nouveau Monde, la classification des parasites Leishmania se fait en fonction de leur capacité à développer des métastases. Ce modéle d'infection a permis de caractériser différentes souches de parasites selon la classification de l'Organisation Mondiale de la Sante (OMS) tel que la souche de référence W>É-II/BR/78/M5313 qui est reconnue comme hautement métastatique alors que ces clones dérivés de M5313 montrent de grandes variations quand a leur capacité à créer des métastases. Les clones 13 et 21 sont métastatiques (M+) alors que les clones 3 et 17 sont nonmétastatiques (NI-). Les objectifs de cette thèse ont été d'étudier le rôle de la réponse immunitaire innée des macrophages après infection in vitro avec différents clones métastatiques et non-métastatiques du parasite L. guyanensis, ainsi que d'étudier la réponse immunitaire générée suite à une infection in vivo par les clones M+ et M- de L. guyanensis dans un modèle marin. L'analyse de la .réponse immunitaire des macrophages in vitro montrent qu'il y aune augmentation significative de leur statut d'activation après infection par des parasites M+ indiquée par la modulation des marqueurs d'activation de surface CD80, CD86 et CD40, ainsi que une augmentation significative de CXCL 10, CCLS, IL6 et TNFa au niveau transcription de l'ARNm et au niveau de la protéine. Cette phénomène d'activation a été observée chez les deux souches de souris C57BL/6 et BALB/c. L'utilisation d'un inhibiteur d'entrée des parasites (Cytochalsin D) ou d'un inhibiteur des fonctions endosomales (Chloroquine) diminue de manière significative la réponse des macrophages aux parasites M+. L'utilisation de macrophages déficients en TLR, MyD88, et TRIF a démontré que la réponse générée après infection par les parasites M+ était dépendante de la voie de signalisation de TRIF et TLR3. Lors d'infection in vivo par des parasites M5313, au moins 50% des souris BALB/c présentent un phénotype sensible caractérisé par des lésions non-nécrotiques qui ne guérissent pas, persistent plus de 13 semaines après infection et contiennent un nombre considérable de parasites. Ces souris développent une réponse immunitaire de type T helper 2 (Th2) avec un niveau élevé d'IL-4 et d'IL-10. Les autres souris ont un phénotype non-sensible, les souris développant peu ou pas de lésion, avec peu de parasites et une réponse immunitaire diminuée, caractérisée par un niveau faible d'IFNy, d'IL4 et d'IL10. De plus, les souris BALB/c infectées par un parasite L. guyanensis isolé à partir des lésions muco-cutanées d'un patient humain atteint de LMC ont démontrés un phénotype similaire aux souris infectées par la souche M5313 avec 50% des souris développant des lésions persistantes, alors qu'un parasite dérivé des lésions cutanées humains n'a montré qu'une faible sensibilité avec une lésion transitoire qui finit par guérir. Nous avons montré que la sensibilité de ces souris BALB/c dépend de l'IL-4 et de l'IL-10 car les souris IL-10-/sur fond génétique BALB/c ainsi que les souris BALB/c traitée avec de l'anti-IL4 étaient capables de contrôler l'infection par M5313. Les souris C57BL/6 sont résistantes à l'infection par le parasite M5313. Elles développent une lésion transitoire qui guérit 9 semaines après infection. Ces souris résistantes ont un très faible taux de parasites au site d'infection et développent une réponse immunitaire de type Thl avec un niveau élevé d'IFNr et peu d'IL4 et d'IL10. Les infections in vivo de souris déficientes en MyD88, TRIF, TLR3 ou TLR9 (sur fond génétique C57BL/6) ont indiqué que MyD88 et TLR9 étaient impliqués dans la résistance à l'infection par L. guyanensi, et que TRIF et TLR3 avaient un rôle important dans la sensibilité. Ce travail met en évidence le fait que la réponse immunitaire de l'hôte est modulée par le parasite selon leur caractérisation d'être soit M+ ou M-. Nous avons démontré également que plusieurs gènes et voies de signalisations étaient impliqués dans cette réponse favorisant le développement d'une LMC. ABSTRACT : Leishmania guyanensis parasites are able to disseminate from the initial site of cutaneous skin infection to the nasopharyngeal tissues resulting in destructive secondary lesions and the disease Mucocutaneous Leishmaniasis (MCL). The secondary lesions in patients have intense immune cell infiltration, elevated immune responses and the presence (albeit at low levels) of parasites. More specifically, MCL patients produce higher levels of TNFa and display impairment in their ability to control the immune response due to a defect in their ability to respond to IL10. Little is known about the pathogenesis of MCL, especially about the dissemination of the infection from the site of inoculation to secondary sites (metastasis) and the response of the host to the pathogen. The hamster model of L. guyanensis infection has previously characterized the WHO reference strain, L. guyanensis WHI/BR/78/M5313, as being highly metastatic. Clones of parasites derived from this reference strain show a differential ability to metastasize. This thesis studied the differential immune response generated by macrophages in vitro, or by mice in vivo, following infection with L. guyanensis parasites. A significant increase in the activation status of macrophages derived from C57BL/6 or BALB/c mice was observed after in vitro infection with L. guyanensis parasites when compared to non-metastatic parasites. This change in status was evidenced by the increased expression of surface activation markers, together with the chemokines, CXCL 10, CCLS, and cytokines, IL6 and TNFa. Furthermore, in vitro infection of macrophages isolated from mice deficient in either a specific Toll Like Receptor (TLR) or the adaptor molecules MyD88 or TRIF, indicated that the immune response generated following L. guyanensis metastatic parasite infection was reliant on the TRIF dependent TLR3 signalling pathway. In vivo footpad infection of BALB/c mice with the L. guyanensis M5313 parasites showed a reproducible susceptible phenotype, whereby at least 50% of infected mice developed non-healing, nonnecrosing lesions with high parasitemia that persisted over 13 weeks post infection. This phenotype was characterized by a Th2 type cytokine immune response with increased levels of IL4 and IL10 detected in the draining lymph nodes. IL 10 deficient mice on a BALB/c background, or BALB/c mice treated with anti-IL4 were able to control infection with L. guyanensis M5313 parasites, thereby proving that these cytokines were indeed implicated in the susceptibility to infection. Moreover, infection of BALB/c mice with patient isolated L. guyanensis parasites confirmed that MCL derived parasites were able to induce a susceptibility phenotype similar to that of L. guyanensis M5313. C57BL/6 mice, on the other hand, were highly resistant to infection with L. guyanensis M5313 parasites and produced transient footpad swelling that healed by week 9 post infection, together with low degrees of footpad parasitemia and a Thl polarized immune response. Infection of mice deficient in MyD88, TRIF, TLR3, and TLR9 (on a C57BL/6 background), indicated that MyD88 and TLR9 were involved in the resistance of these mice to infection, and that TRIF and TLR3 were involved in the susceptibility. This study has shown that the host response can be differentially modulated depending on the infecting parasite with several genes and pathways being identified that could be involved in promoting the development of MCL.

Role of hepatocyte wounding on " Plasmodium " liver-stage development and survival

RESUME La première étape primordiale au cycle de vie du Plasmodium dans un hôte mammifère est l'invasion des hepatocytes par des sporozoites. L'infection finale des hepatocytes est précédée de la traversée de plusieurs cellules hôtes, rompant les membranes plasmiques et ayant comme résultat la sécrétion des facteurs cytotoliques dans le micro-environnement. Ce matériel endogène libéré est fortement stimulant/immunogène et peut servir de signal de danger initiant des réponses distinctes dans diverses cellules. De nos jours, le caractère essentiel et salutaire de la migration des sporozoites comme étape d'infection du Plasmodium est vivement controversée. Ainsi, notre étude a visé à caractériser l'effet de l'interaction du parasite avec ses cellules hôtes d'un point de vue immunologique. En particulier, nous avons voulu évaluer l'effet de la perte de matériel cellulaire pendant l'infection de Plasmodium sur les hepatocytes primaires de souris et sur des cultures cellulaires HepG2. Nous avons observé que les facteurs cytotoxiques dérivés des cellules endommagés activent NF-κB - un important régulateur de réponse inflammatoires -dans des cellules voisines des cellules endommagés, qui sont des cellules hôtes potentielles pour l'infection finale du parasite. Cette activation de NF-κB s'est produite peu de temps après l'infection et a mené in vitro et in vivo à une réduction d'infection de façon dépendante du temps, un effet qui a pu être compensé par l'addition de BAY11-7082, un inhibiteur spécifique de NF-κB. De plus, aucune activation de NF-κB avec des parasites SPECT-/-, incapables de traverser les hepatocytes, n'a été observée. Nous avons montré parla suite que l'activation de NF-κB induit l'expression de l'enzyme iNOS dans les hepatocytes, qui est responsable d'une diminution des hepatocytes infectés. En outre, les hepatocytes primaires des souris MyD88-/- n'ont montré ni activation de NF-κB, ni expression d'iNOS lors de l'infection, ce qui suggère la participation des membres de famille du Toll/IL-1 récepteur dans la reconnaissance des facteurs cytosoxiques. En effet, le manque de MyD88 a augmenté significativement l'infection in vitro et in vivo. D'autre part, un rôle bénéfique pour l'activation de NF-κB a été évalué. Les cellules infectées étaient plus résistantes contre l'apoptose induite par Fas (CD95/Apo-1) que les cellules non infectées ou les cellules infectées dans lesquelles NF-κB a été bloqué par BAY11-7082 in vitro. Paradoxalement, l'expression d'iNOS contribue à la protection des cellules infectées contre l'apoptose pax Fas, puisque le traitement avec l'inhibiteur spécifique SMT (S-methylisothiourea) a rendu les cellules infectées plus susceptibles à l'apoptose. Un effet bénéfique additionnel pour le parasite est que la plupart des cellules hôtes traversées présentent des peptides du parasite aux cellules T cytotoxiques spécifiques et peuvent donc réorienter la réaction immune spécifique sur les cellules non infectées. Nous montrons que les cellules hôtes endommagés par la migration du parasite induit l'inflammation, qui limite l'ampleur de l'infection. D'autre part, nos données soutiennent que la survie du parasite Plasmodium dans le foie est assurée par une augmentation de la résistance des hepatocytes contre l'apoptose. SUMMARY The first obligatory step of the Plasmodium life cycle in the mammalian host is the invasion of hepatocytes by sporozoites. Final hepatocyte infection involves the penetration of several host cells, whose plasma membranes are ruptured in the process, resulting in the release of cytosolic factors into the microenvironment. This released endogenous material is highly stimulatory / immunogenic and can serve as a danger signal initiating distinct responses in various cells. To date, it is highly controversial whether sporozoite migration through hepatocytes is an essential and beneficial step for Plasmodium infection. Thus, our study aimed at characterizing the effect of the interaction of the parasite with its host cells from an immunological point of view In particular, we wanted to evaluate the effect of cell material leakage during Plasmodium infection on cultured mouse primary hepatocytes and HepG2 cells. We observed that wounded cell-derived cytosolic factors activate NF-κB - a main regulator of host inflammatory responses - in cells bordering wounded cells, which are potential host cells for final parasite infection. This activation of NF-κB occurred shortly after infection and led to a reduction of infection load in a time dependent manner in vitro and in viva, an effect that could be reverted by addition of the specific NF-κB inhibitor BAY11-7082. In addition, no NF-κB activation was observed when SPECT-/- parasites, which are devoid of hepatocyte traversing properties, were used. We provide further evidence that NF-κB activation causes the induction of inducible nitric oxide synthase (iNOS) expression in hepatocytes, and this is, in turn, responsible for a decrease in Plasmodium-infected hepatocytes. Furthermore, primary hepatocytes from MyD88-/- mice showed no NF-κB activation and iNOS expression upon infection, suggesting a role of the Toll/IL-1 receptor family members in sensing cytosolic factors. Indeed, lack of MyD88 significantly increased infection in vitro and in vivo. In a further complementary series of experiments, we assessed a possible beneficial role for the activation of NF-κB. Infected cells were more resistant to Fas (CD95/Apo-1)-mediated apoptosis than uninfected cells or infected cells in which NF-κB was blocked by BAYl1-7082 in vitro. Paradoxically, iNOS expression contributes to the protection of infected cells from Fas-induced apoptosis, since treatment with the specific iNOS inhibitor SMT (S-Methylisothiourea Sulfate) rendered the infected cells more susceptible to apoptosis. An additional beneficial effect of host cell traversal for the parasite is the fact that mainly traversed cells present parasite-derived peptides to specific cytotoxic T cells and therefore may redirect the specific immune response to uninfected cells. In summary, we have shown that host cells wounded by parasite migration induce inflammation, which limits the extent of parasite infection. In addition, our data support the notion that survival of Plasmodium parasites in the liver is mediated by increasing the resistance of hepatocytes to Fas-induced apoptosis.

Discovery of novel arenavirus receptors and entry factors

Arenaviruses are enveloped negative-strand RNA viruses that contain a bi-segmented genome. They are rodent-borne pathogens endemic to the Americas and Africa, with the exception of lymphocytic choriomeningitis virus (LCMV) that is world-wide distributed. The arenaviruses include numerous important human pathogens including the Old World arenavirus Lassa virus (LASV), the causative agent of a severe viral hemorrhagic fever in humans with several hundred thousand infections per year in Africa and thousands of deaths. Viruses are obligatory intracellular parasites, strictly depending on cellular processes and factors to complete their replication cycle. The binding of a virus to target cells is the first step of every viral infection, and is mainly mediated by viral proteins that can directly engage cellular receptors, providing a key determinant for viral tropism. This early step of infection represents a promising target to block the pathogen before it can take control over the host cell. Old World arenaviruses, such as LASV and LCMV, bind to host cells via attachment to their main receptor, dystroglycan (DG), an ubiquitous receptor for extracellular matrix proteins. The engagement of DG by LASV results in a fast internalization and transfer the virus to late endosomal compartment suggesting that the virus binding to DG causes marked changes in the dynamics of the receptor. These events could result in the clustering of the receptor and subsequent induction of signaling that could be modulated by the virus. Recently, numerous findings also suggest the presence of alternative receptor(s) for LASV in absence of the main DG receptor. In my first project, I was interested to investigate the effects of virus-receptor binding on the tyrosine phosphorylation of the cytoplasmic domain of DG and to test if this post-translational modification was crucial for the internalization of the LASV-receptor complex. We found that engagement of cellular DG by a recombinant LCMV expressing the envelope GP of LASV in human epithelial cells induced tyrosine phosphorylation of the cytoplasmic domain of DG. LASV GP binding to DG further resulted in dissociation of the adapter protein utrophin from virus-bound DG. Virus-induced dissociation of utrophin and consequent virus internalization were affected by the broadly specific tyrosine kinase inhibitor genistein. We speculate that the detachment of virus- bound DG from the actin-based cytoskeleton following DG phosphorylation may facilitate subsequent endocytosis of the virus-receptor complex. In the second project, I was interested to characterize the newly indentified LASV alternative receptor Axl in the context of productive arenavirus infection. In a first step, we demonstrated that Axl supports productive infection by rLCMV-LASVGP in a DG-independent manner. In line with previous studies, cell entry of rLCMV-LASVGP via Axl was less efficient when compared to functional DG. Interestingly, Axl-mediated infection showed rapid kinetics similar to DG-dependent entry. Using a panel of inhibitors, we found that Axl-mediated cell entry of rLCMV-LASVGP involved a clathrin-independent pathway that critically depended on actin and dynamin and was sensitive to EIPA but not to PAK inhibitors, compatible with a macropinocytosis-like mechanism of entry. In a next step, we aimed to investigate the molecular mechanism by which rLCMV-LASVGP recognizes Axl. Phosphatidylserine (PS) is the natural ligand of Axl via the adaptor protein Gas6. We detected the presence of PS in the envelope of Old World arenaviruses, suggesting that PS could mediate Axl-virus binding, in a mechanism of apoptotic mimicry already described for other viruses. Whether envelope PS and/or the GP of LASV plays any role in virus entry via Axl is still an open question. The molecular mechanisms underlying host cell-virus interaction are of particular interest to answer basic scientific questions as well as to apply key findings to translational research. Understanding pathogen induced-signaling and its link to invasion of the host cell is of great importance to develop drugs for therapeutic intervention against highly pathogenic viruses like LASV. - Les Arenavirus sont des virus enveloppés à ARN négatifs organisés sous forme de génome bisegmenté. Ils sont véhiculés par les rongeurs et se retrouvent de manière endémique aux Amériques et en Afrique avec l'exception du virus de la chorioméningite lymphocytaire (LCMV) qui lui est distribué mondialement. De nombreux pathogènes humains font parti de la famille des Arenavirus dont le virus de l'Ancien Monde Lassa (LASV), un agent responsable de fièvres hémorragiques sévères chez les humains. Le virus de Lassa cause plusieurs centaines de milliers d'infections par année en Afrique ainsi que des milliers de morts. De manière générale, les virus sont des parasites intracellulaires obligatoires qui dépendent strictement de processus et facteurs cellulaires pour clore leur cycle de réplication. L'attachement d'un virus à sa cellule cible représente la première étape de chaque infection virale et est principalement dirigée par des protéines virales qui interagissent directement avec leur récepteurs cellulaires respectifs fournissant ainsi un indicateur déterminant pour le tropisme d'un virus. Cette première étape de l'infection représente aussi une cible prometteuse pour bloquer le pathogène avant qu'il ne puisse prendre le contrôle de la cellule. Les Arenavirus de l'Ancien Monde comme LASV et LCMV s'attachent à la cellule hôte en se liant à leur récepteur principal, le dystroglycan (DG), un récepteur ubiquitaire pour les protéines de la matrice extracellulaire. La liaison du DG par LASV résulte en une rapide internalisation transférant le virus aux endosomes tardifs suggérant ainsi que l'attachement du virus au DG peut provoquer des changements marqués dans la dynamique moléculaire du récepteur. Ces événements sont susceptibles d'induire un regroupement du récepteur à la surface cellulaire, ainsi qu'une induction subséquente qui pourrait être, par la suite, modulée par le virus. Récemment, plusieurs découvertes suggèrent aussi la présence d'un récepteur alternatif pour LASV en l'absence du récepteur principal, le DG. Concernant mon premier projet, j'étais intéressée à étudier les effets de la liaison virus- récepteur sur la phosphorylation des acides aminés tyrosines se trouvant dans la partie cytoplasmique du DG, le but étant de tester si cette modification post-translationnelle était cruciale pour Γ internalisation du complexe LASV-DG récepteur. Nous avons découvert que l'engagement du récepteur DG par le virus recombinant LCMV, exprimant la glycoprotéine de LASV, dans des cellules épithéliales humaines induit une phosphorylation de résidu(s) tyrosine se situant dans le domaine cytoplasmique du DG. La liaison de la glycoprotéine de LASV au DG induit par la suite la dissociation de la protéine adaptatrice utrophine du complexe virus-DG récepteur. Nous avons observé que cette dissociation de l'utrophine, induite par le virus, ainsi que son internalisation, sont affectées par l'inhibiteur à large spectre des tyrosines kinases, la génistéine. Nous avons donc supposé que le détachement du virus, lié au récepteur DG, du cytosquelette d'actine suite à la phosphorylation du DG faciliterait l'endocytose subséquente du complexe virus-récepteur. Dans le second projet, j'étais intéressée à caractériser le récepteur alternatif Axl qui a été récemment identifié dans le contexte de l'infection productive des Arenavirus. Dans un premier temps, nous avons démontré que le récepteur alternatif Axl permet l'infection des cellules par le virus LCMV recombinant LASV indépendamment du récepteur DG. Conformément aux études publiées précédemment, nous avons pu observer que l'entrée du virus recombinant LASV via Axl est moins efficace que via le récepteur principal DG. De façon intéressante, nous avons aussi remarqué que l'infection autorisée par Axl manifeste une cinétique virale d'entrée similaire à celle observée avec le récepteur DG. Utilisant un éventail de différents inhibiteurs, nous avons trouvé que l'entrée du virus recombinant rLCMV-LASVGP via Axl implique une voie d'entrée indépendante de la clathrine et dépendant de manière critique de l'actine et de la dynamine. Cette nouvelle voie d'entrée est aussi sensible à l'EIPA contrairement aux inhibiteurs PAK indiquant un mécanisme d'entrée compatible avec un mécanisme de macropinocytose. L'étape suivante du projet a été d'investiguer le mécanisme moléculaire par lequel le virus recombinant rLCMV-LASVGP reconnaît le récepteur alternatif Axl. La phosphatidylsérine (PS) se trouve être un ligand naturel pour Axl via la protéine adaptatrice Gas6. Nous avons détecté la présence de PS dans l'enveloppe des Arenavirus du Vieux Monde suggérant que la PS pourrait médier la liaison du virus à Axl dans un mécanisme de mimétisme apoptotique déjà observé et décrit pour d'autres virus. Cependant, il reste encore à déterminer qui de la PS ou de la glycoprotéine de l'enveloppe virale intervient dans le processus d'entrée de LASV via le récepteur alternatif Axl. Les mécanismes moléculaires à la base de l'interaction entre virus et cellule hôte sont d'intérêts particuliers pour répondre aux questions scientifiques de base ainsi que dans l'application de découvertes clés pour la recherche translationnelle. La compréhension de la signalisation induite par les pathogènes ainsi que son lien à l'invasion de la cellule hôte est d'une importance considérable pour le développement de drogues pour l'intervention thérapeutique contre les virus hautement pathogènes comme LASV.

Cell locations AQP1, AQP4 and 9 in the non-human primate brain

Rapport de synthèse : La présence de trois canaux d'eau, appelés aquaporines AQP1, AQP4 et AQP9, a été observée dans le cerveau sain ainsi que dans plusieurs modèles des pathologies cérébrales des rongeurs. Peu est connu sur la distribution des AQP dans le cerveau des primates. Cette connaissance sera utile pour des futurs essaies médicamenteux qui visent à prévenir la formation des oedèmes cérébraux. Nous avons étudié l'expression et la distribution cellulaire des AQP1, 4 et 9 dans le cerveau primate non-humain. La distribution des AQP4 dans le cerveau primate non-humain a été observée dans des astrocytes périvasculaires, comparable à l'observation faite dans le cerveau du rongeur. Contrairement à ce qui a été décrit chez le rongeur, l'AQPI chez le primate est exprimée dans les processus et dans les prolongations périvasculaires d'un sous-type d'astrocytes, qui est avant tout localisé dans la matière blanche et dans la glia limitans et qui est peut-être impliqué dans l'homéostasie de l'eau. L'AQPI a aussi été observée dans les neurones qui innervent des vaisseaux sanguins de la pie-mère, suggérant un rôle possible dans la régularisation de la vascularisation cérébrale. Comme décrit chez le rongeur, le mRNA et les protéines de l'AQP9 ont été détectés dans des astrocytes et dans des neurones catécholaminergiques. Chez le primate, des localisations supplémentaires ont été observées dans des populations de neurones placées dans certaines zones corticales. Cet article décrit une étude détaillée sur la distribution des AQP1, 4 et 9 dans le cerveau primate non-humain. Les observations faites s'additionnent aux data déjà publié sur le cerveau du rongeur. Ces importantes différences entre les espèces doivent être considérées dans l'évaluation des médicaments qui agiront potentiellement sur des AQP des primates non-humains avant d'entrer dans la phase des essais cliniques sur des humains.

Public infrastructure and the performance of manufacturing industries: short-and long-run

We present a theoretical framework for determining the short- and long-run effects of infrastructure. While the short-run effects have been the focus of most previous studies, here we derive long-run elasticities by taking into account the adjustment of quasi-fixed inputs to their optimum levels. By considering the impact of infrastructure on private investment decisions, we observe how, apart from the direct effect on costs in the short-run, infrastructure exerts an indirect source of influence in the long-run through their effect on private capital. The model is applied to manufacturing industries in the Spanish regions

Deficit, human capital and economic growth dynamics

Long-run economic growth arouses a great interest since it can shed light on the income-path of an economy and try to explain the large differences in income we observe across countries and over time. The neoclassical model has been followed by several endogenous growth models which, contrarily to the former, seem to predict that economies with similar preferences and technological level, do not necessarily tend to converge to similar per capita income levels. This paper attempts to show a possible mechanismthrough which macroeconomic disequilibria and inefficiencies, represented by budget deficits, may hinder human capital accumulation and therefore economic growth. Using a mixed education system, deficit is characterized as a bug agent which may end up sharply reducing the resources devoted to education and training. The paper goes a step further from the literature on deficit by introducing a rich dynamic analysis of the effects of a deficit reduction on different economic aspects.Following a simple growth model and allowing for slight changes in the law of human capital accumulation, we reach a point where deficit might sharply reduce human capital accumulation. On the other hand, a deficit reduction carried on for a long time, taking that reduction as a more efficient management of the economy, may prove useful in inducing endogenous growth. Empirical evidence for a sample of countries seems to support the theoretical assumptions in the model: (1) evidence on an inverse relationship betweendeficit and human capital accumulation, (2) presence of a strongly negative associationbetween the quantity of deficit in the economy and the rate of growth. They may prove a certain role for budget deficit in economic growth

Deficit, human capital and economic growth dynamics

Long-run economic growth arouses a great interest since it can shed light on the income-path of an economy and try to explain the large differences in income we observe across countries and over time. The neoclassical model has been followed by several endogenous growth models which, contrarily to the former, seem to predict that economies with similar preferences and technological level, do not necessarily tend to converge to similar per capita income levels. This paper attempts to show a possible mechanismthrough which macroeconomic disequilibria and inefficiencies, represented by budget deficits, may hinder human capital accumulation and therefore economic growth. Using a mixed education system, deficit is characterized as a bug agent which may end up sharply reducing the resources devoted to education and training. The paper goes a step further from the literature on deficit by introducing a rich dynamic analysis of the effects of a deficit reduction on different economic aspects.Following a simple growth model and allowing for slight changes in the law of human capital accumulation, we reach a point where deficit might sharply reduce human capital accumulation. On the other hand, a deficit reduction carried on for a long time, taking that reduction as a more efficient management of the economy, may prove useful in inducing endogenous growth. Empirical evidence for a sample of countries seems to support the theoretical assumptions in the model: (1) evidence on an inverse relationship betweendeficit and human capital accumulation, (2) presence of a strongly negative associationbetween the quantity of deficit in the economy and the rate of growth. They may prove a certain role for budget deficit in economic growth

Soil shear strength under non-irrigated and irrigated short duration grazing systems

Pasture productivity can drop due to soil compaction caused by animal trampling. Physical and mechanical alterations are therefore extremely important indicators for pasture management. The objective of this research was to: draw and evaluate the Mohr failure line of a Red Yellow Latossol under different pasture cycles and natural forest; calculate apparent cohesion; observe possible physical alterations in this soil; and propose a correction factor for stocking rates based on shear strength properties. This study was conducted between March/2006 and March/2007 on the Experimental Farm of Fundação de Ensino Superior de Passos, in Passos, state of Minas Gerais. The total study area covered 6 ha, of which 2 ha were irrigated pasture, 2 ha non-irrigated pasture and 2 ha natural forest. Brachiaria brizantha cv. MG-5 Vitória was used as forage plant. The pasture area was divided into paddocks. The Mohr failure line of samples of a Red Yellow Latossol under irrigated pasture equilibrated at a tension of water content of 6 kPa indicated higher shear strength than under non-irrigated pasture. The shear strength under irrigated pasture and natural forest was higher than under non-irrigated pasture. At a tension of water content of 33 kPa no difference was found in shear strength between management and use. Possible changes in soil structure were caused by apparent cohesion. The values of the correction factor were close to 1, which may indicate a possible soil compaction in prolonged periods of management.

Public infrastructure and the performance of manufacturing industries: short-and long-run

We present a theoretical framework for determining the short- and long-run effects of infrastructure. While the short-run effects have been the focus of most previous studies, here we derive long-run elasticities by taking into account the adjustment of quasi-fixed inputs to their optimum levels. By considering the impact of infrastructure on private investment decisions, we observe how, apart from the direct effect on costs in the short-run, infrastructure exerts an indirect source of influence in the long-run through their effect on private capital. The model is applied to manufacturing industries in the Spanish regions

Force volume and stiffness tomography investigation on the dynamics of stiff material under bacterial membranes.

The determination of the characteristics of micro-organisms in clinical specimens is essential for the rapid diagnosis and treatment of infections. A thorough investigation of the nanoscale properties of bacteria can prove to be a fundamental tool. Indeed, in the latest years, the importance of high resolution analysis of the properties of microbial cell surfaces has been increasingly recognized. Among the techniques available to observe at high resolution specific properties of microscopic samples, the Atomic Force Microscope (AFM) is the most widely used instrument capable to perform morphological and mechanical characterizations of living biological systems. Indeed, AFM can routinely study single cells in physiological conditions and can determine their mechanical properties with a nanometric resolution. Such analyses, coupled with high resolution investigation of their morphological properties, are increasingly used to characterize the state of single cells. In this work, we exploit the capabilities and peculiarities of AFM to analyze the mechanical properties of Escherichia coli in order to evidence with a high spatial resolution the mechanical properties of its structure. In particular, we will show that the bacterial membrane is not mechanically uniform, but contains stiffer areas. The force volume investigations presented in this work evidence for the first time the presence and dynamics of such structures. Such information is also coupled with a novel stiffness tomography technique, suggesting the presence of stiffer structures present underneath the membrane layer that could be associated with bacterial nucleoids.

Cheminformatic studies of molecular properties and their influence on bioactivity in the context of drug discovery

Molecular shape has long been known to be an important property for the process of molecular recognition. Previous studies postulated the existence of a drug-like shape space that could be used to artificially bias the composition of screening libraries, with the aim to increase the chance of success in Hit Identification. In this work, it was analysed to which extend this assumption holds true. Normalized Principal Moments of Inertia Ratios (NPRs) have been used to describe the molecular shape of small molecules. It was investigated, whether active molecules of diverse targets are located in preferred subspaces of the NPR shape space. Results illustrated a significantly stronger clustering than could be expected by chance, with parts of the space unlikely to be occupied by active compounds. Furthermore, a strong enrichment of elongated, rather flat shapes could be observed, while globular compounds were highly underrepresented. This was confirmed for a wide range of small molecule datasets from different origins. Active compounds exhibited a high overlap in their shape distributions across different targets, making a purely shape­ based discrimination very difficult. An additional perspective was provided by comparing the shapes of protein binding pockets with those of their respective ligands. Although more globular than their ligands, it was observed that binding sites shapes exhibited a similarly skewed distribution in shape space: spherical shapes were highly underrepresented. This was different for unoccupied binding pockets of smaller size. These were on the contrary identified to possess a more globular shape. The relation between shape complementarity and exhibited bioactivity was analysed; a moderate correlation between bioactivity and parameters including pocket coverage, distance in shape space, and others could be identified, which reflects the importance of shape complementarity. However, this also suggests that other aspects are of relevance for molecular recognition. A subsequent analysis assessed if and how shape and volume information retrieved from pocket or respective reference ligands could be used as a pre-filter in a virtual screening approach. ln Lead Optimization compounds need to get optimized with respect to a variety of pararneters. Here, the availability of past success stories is very valuable, as they can guide medicinal chemists during their analogue synthesis plans. However, although of tremendous interest for the public domain, so far only large corporations had the ability to mine historical knowledge in their proprietary databases. With the aim to provide such information, the SwissBioisostere database was developed and released during this thesis. This database contains information on 21,293,355 performed substructural exchanges, corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1,948 molecular targets representing 30 target classes, and on their impact on bioactivity . A user-friendly interface was developed that provides facile access to these data and is accessible at http//www.swissbioisostere.ch. The ChEMBL database was used as primary data source of bioactivity information. Matched molecular pairs have been identified in the extracted and cleaned data. Success-based scores were developed and integrated into the database to allow re-ranking of proposed replacements by their past outcomes. It was analysed to which degree these scores correlate with chemical similarity of the underlying fragments. An unexpectedly weak relationship was detected and further investigated. Use cases of this database were envisioned, and functionalities implemented accordingly: replacement outcomes are aggregatable at the assay level, and it was shawn that an aggregation at the target or target class level could also be performed, but should be accompanied by a careful case-by-case assessment. It was furthermore observed that replacement success depends on the activity of the starting compound A within a matched molecular pair A-B. With increasing potency the probability to lose bioactivity through any substructural exchange was significantly higher than in low affine binders. A potential existence of a publication bias could be refuted. Furthermore, often performed medicinal chemistry strategies for structure-activity-relationship exploration were analysed using the acquired data. Finally, data originating from pharmaceutical companies were compared with those reported in the literature. It could be seen that industrial medicinal chemistry can access replacement information not available in the public domain. In contrast, a large amount of often-performed replacements within companies could also be identified in literature data. Preferences for particular replacements differed between these two sources. The value of combining different endpoints in an evaluation of molecular replacements was investigated. The performed studies highlighted furthermore that there seem to exist no universal substructural replacement that always retains bioactivity irrespective of the biological environment. A generalization of bioisosteric replacements seems therefore not possible. - La forme tridimensionnelle des molécules a depuis longtemps été reconnue comme une propriété importante pour le processus de reconnaissance moléculaire. Des études antérieures ont postulé que les médicaments occupent préférentiellement un sous-ensemble de l'espace des formes des molécules. Ce sous-ensemble pourrait être utilisé pour biaiser la composition de chimiothèques à cribler, dans le but d'augmenter les chances d'identifier des Hits. L'analyse et la validation de cette assertion fait l'objet de cette première partie. Les Ratios de Moments Principaux d'Inertie Normalisés (RPN) ont été utilisés pour décrire la forme tridimensionnelle de petites molécules de type médicament. Il a été étudié si les molécules actives sur des cibles différentes se co-localisaient dans des sous-espaces privilégiés de l'espace des formes. Les résultats montrent des regroupements de molécules incompatibles avec une répartition aléatoire, avec certaines parties de l'espace peu susceptibles d'être occupées par des composés actifs. Par ailleurs, un fort enrichissement en formes allongées et plutôt plates a pu être observé, tandis que les composés globulaires étaient fortement sous-représentés. Cela a été confirmé pour un large ensemble de compilations de molécules d'origines différentes. Les distributions de forme des molécules actives sur des cibles différentes se recoupent largement, rendant une discrimination fondée uniquement sur la forme très difficile. Une perspective supplémentaire a été ajoutée par la comparaison des formes des ligands avec celles de leurs sites de liaison (poches) dans leurs protéines respectives. Bien que plus globulaires que leurs ligands, il a été observé que les formes des poches présentent une distribution dans l'espace des formes avec le même type d'asymétrie que celle observée pour les ligands: les formes sphériques sont fortement sous­ représentées. Un résultat différent a été obtenu pour les poches de plus petite taille et cristallisées sans ligand: elles possédaient une forme plus globulaire. La relation entre complémentarité de forme et bioactivité a été également analysée; une corrélation modérée entre bioactivité et des paramètres tels que remplissage de poche, distance dans l'espace des formes, ainsi que d'autres, a pu être identifiée. Ceci reflète l'importance de la complémentarité des formes, mais aussi l'implication d'autres facteurs. Une analyse ultérieure a évalué si et comment la forme et le volume d'une poche ou de ses ligands de référence pouvaient être utilisés comme un pré-filtre dans une approche de criblage virtuel. Durant l'optimisation d'un Lead, de nombreux paramètres doivent être optimisés simultanément. Dans ce contexte, la disponibilité d'exemples d'optimisations réussies est précieuse, car ils peuvent orienter les chimistes médicinaux dans leurs plans de synthèse par analogie. Cependant, bien que d'un extrême intérêt pour les chercheurs dans le domaine public, seules les grandes sociétés pharmaceutiques avaient jusqu'à présent la capacité d'exploiter de telles connaissances au sein de leurs bases de données internes. Dans le but de remédier à cette limitation, la base de données SwissBioisostere a été élaborée et publiée dans le domaine public au cours de cette thèse. Cette base de données contient des informations sur 21 293 355 échanges sous-structuraux observés, correspondant à 5 586 462 remplacements uniques mesurés dans 35 039 tests contre 1948 cibles représentant 30 familles, ainsi que sur leur impact sur la bioactivité. Une interface a été développée pour permettre un accès facile à ces données, accessible à http:/ /www.swissbioisostere.ch. La base de données ChEMBL a été utilisée comme source de données de bioactivité. Une version modifiée de l'algorithme de Hussain et Rea a été implémentée pour identifier les Matched Molecular Pairs (MMP) dans les données préparées au préalable. Des scores de succès ont été développés et intégrés dans la base de données pour permettre un reclassement des remplacements proposés selon leurs résultats précédemment observés. La corrélation entre ces scores et la similarité chimique des fragments correspondants a été étudiée. Une corrélation plus faible qu'attendue a été détectée et analysée. Différents cas d'utilisation de cette base de données ont été envisagés, et les fonctionnalités correspondantes implémentées: l'agrégation des résultats de remplacement est effectuée au niveau de chaque test, et il a été montré qu'elle pourrait également être effectuée au niveau de la cible ou de la classe de cible, sous réserve d'une analyse au cas par cas. Il a en outre été constaté que le succès d'un remplacement dépend de l'activité du composé A au sein d'une paire A-B. Il a été montré que la probabilité de perdre la bioactivité à la suite d'un remplacement moléculaire quelconque est plus importante au sein des molécules les plus actives que chez les molécules de plus faible activité. L'existence potentielle d'un biais lié au processus de publication par articles a pu être réfutée. En outre, les stratégies fréquentes de chimie médicinale pour l'exploration des relations structure-activité ont été analysées à l'aide des données acquises. Enfin, les données provenant des compagnies pharmaceutiques ont été comparées à celles reportées dans la littérature. Il a pu être constaté que les chimistes médicinaux dans l'industrie peuvent accéder à des remplacements qui ne sont pas disponibles dans le domaine public. Par contre, un grand nombre de remplacements fréquemment observés dans les données de l'industrie ont également pu être identifiés dans les données de la littérature. Les préférences pour certains remplacements particuliers diffèrent entre ces deux sources. L'intérêt d'évaluer les remplacements moléculaires simultanément selon plusieurs paramètres (bioactivité et stabilité métabolique par ex.) a aussi été étudié. Les études réalisées ont souligné qu'il semble n'exister aucun remplacement sous-structural universel qui conserve toujours la bioactivité quel que soit le contexte biologique. Une généralisation des remplacements bioisostériques ne semble donc pas possible.

Do seismic waves and fluid flow sense the same permeability?

We study the discrepancy between the effective flow permeability and the effective seismic permeability, that is, the effective permeability controlling seismic attenuation due to wave-induced fluid flow, in 2D rock samples having mesoscopic heterogeneities and in the presence of strong permeability fluctuations. In order to do so, we employ a numerical oscillatory compressibility test to determine attenuation and velocity dispersion due to wave-induced fluid flow in these kinds of media and compare the responses with those obtained by replacing the heterogeneous permeability field by constant values, including the average permeability as well as the effective flow permeability of the sample. The latter is estimated in a separate upscaling procedure by solving the steady-state flow equation in the rock sample under study. Numerical experiments let us verify that attenuation levels are less significant and the attenuation peak gets broader in the presence of such strong permeability fluctuations. Moreover, we observe that for very low frequencies the effective seismic permeability is similar to the effective flow permeability, while for very high frequencies it approaches the arithmetic average of the permeability field.

Community Policing in Switzerland's Major Urban Areas: An Observational Study of the Implementation and Impact Using Geospatial Data Mining

La présente étude est à la fois une évaluation du processus de la mise en oeuvre et des impacts de la police de proximité dans les cinq plus grandes zones urbaines de Suisse - Bâle, Berne, Genève, Lausanne et Zurich. La police de proximité (community policing) est à la fois une philosophie et une stratégie organisationnelle qui favorise un partenariat renouvelé entre la police et les communautés locales dans le but de résoudre les problèmes relatifs à la sécurité et à l'ordre public. L'évaluation de processus a analysé des données relatives aux réformes internes de la police qui ont été obtenues par l'intermédiaire d'entretiens semi-structurés avec des administrateurs clés des cinq départements de police, ainsi que dans des documents écrits de la police et d'autres sources publiques. L'évaluation des impacts, quant à elle, s'est basée sur des variables contextuelles telles que des statistiques policières et des données de recensement, ainsi que sur des indicateurs d'impacts construit à partir des données du Swiss Crime Survey (SCS) relatives au sentiment d'insécurité, à la perception du désordre public et à la satisfaction de la population à l'égard de la police. Le SCS est un sondage régulier qui a permis d'interroger des habitants des cinq grandes zones urbaines à plusieurs reprises depuis le milieu des années 1980. L'évaluation de processus a abouti à un « Calendrier des activités » visant à créer des données de panel permettant de mesurer les progrès réalisés dans la mise en oeuvre de la police de proximité à l'aide d'une grille d'évaluation à six dimensions à des intervalles de cinq ans entre 1990 et 2010. L'évaluation des impacts, effectuée ex post facto, a utilisé un concept de recherche non-expérimental (observational design) dans le but d'analyser les impacts de différents modèles de police de proximité dans des zones comparables à travers les cinq villes étudiées. Les quartiers urbains, délimités par zone de code postal, ont ainsi été regroupés par l'intermédiaire d'une typologie réalisée à l'aide d'algorithmes d'apprentissage automatique (machine learning). Des algorithmes supervisés et non supervisés ont été utilisés sur les données à haute dimensionnalité relatives à la criminalité, à la structure socio-économique et démographique et au cadre bâti dans le but de regrouper les quartiers urbains les plus similaires dans des clusters. D'abord, les cartes auto-organisatrices (self-organizing maps) ont été utilisées dans le but de réduire la variance intra-cluster des variables contextuelles et de maximiser simultanément la variance inter-cluster des réponses au sondage. Ensuite, l'algorithme des forêts d'arbres décisionnels (random forests) a permis à la fois d'évaluer la pertinence de la typologie de quartier élaborée et de sélectionner les variables contextuelles clés afin de construire un modèle parcimonieux faisant un minimum d'erreurs de classification. Enfin, pour l'analyse des impacts, la méthode des appariements des coefficients de propension (propensity score matching) a été utilisée pour équilibrer les échantillons prétest-posttest en termes d'âge, de sexe et de niveau d'éducation des répondants au sein de chaque type de quartier ainsi identifié dans chacune des villes, avant d'effectuer un test statistique de la différence observée dans les indicateurs d'impacts. De plus, tous les résultats statistiquement significatifs ont été soumis à une analyse de sensibilité (sensitivity analysis) afin d'évaluer leur robustesse face à un biais potentiel dû à des covariables non observées. L'étude relève qu'au cours des quinze dernières années, les cinq services de police ont entamé des réformes majeures de leur organisation ainsi que de leurs stratégies opérationnelles et qu'ils ont noué des partenariats stratégiques afin de mettre en oeuvre la police de proximité. La typologie de quartier développée a abouti à une réduction de la variance intra-cluster des variables contextuelles et permet d'expliquer une partie significative de la variance inter-cluster des indicateurs d'impacts avant la mise en oeuvre du traitement. Ceci semble suggérer que les méthodes de géocomputation aident à équilibrer les covariables observées et donc à réduire les menaces relatives à la validité interne d'un concept de recherche non-expérimental. Enfin, l'analyse des impacts a révélé que le sentiment d'insécurité a diminué de manière significative pendant la période 2000-2005 dans les quartiers se trouvant à l'intérieur et autour des centres-villes de Berne et de Zurich. Ces améliorations sont assez robustes face à des biais dus à des covariables inobservées et covarient dans le temps et l'espace avec la mise en oeuvre de la police de proximité. L'hypothèse alternative envisageant que les diminutions observées dans le sentiment d'insécurité soient, partiellement, un résultat des interventions policières de proximité semble donc être aussi plausible que l'hypothèse nulle considérant l'absence absolue d'effet. Ceci, même si le concept de recherche non-expérimental mis en oeuvre ne peut pas complètement exclure la sélection et la régression à la moyenne comme explications alternatives. The current research project is both a process and impact evaluation of community policing in Switzerland's five major urban areas - Basel, Bern, Geneva, Lausanne, and Zurich. Community policing is both a philosophy and an organizational strategy that promotes a renewed partnership between the police and the community to solve problems of crime and disorder. The process evaluation data on police internal reforms were obtained through semi-structured interviews with key administrators from the five police departments as well as from police internal documents and additional public sources. The impact evaluation uses official crime records and census statistics as contextual variables as well as Swiss Crime Survey (SCS) data on fear of crime, perceptions of disorder, and public attitudes towards the police as outcome measures. The SCS is a standing survey instrument that has polled residents of the five urban areas repeatedly since the mid-1980s. The process evaluation produced a "Calendar of Action" to create panel data to measure community policing implementation progress over six evaluative dimensions in intervals of five years between 1990 and 2010. The impact evaluation, carried out ex post facto, uses an observational design that analyzes the impact of the different community policing models between matched comparison areas across the five cities. Using ZIP code districts as proxies for urban neighborhoods, geospatial data mining algorithms serve to develop a neighborhood typology in order to match the comparison areas. To this end, both unsupervised and supervised algorithms are used to analyze high-dimensional data on crime, the socio-economic and demographic structure, and the built environment in order to classify urban neighborhoods into clusters of similar type. In a first step, self-organizing maps serve as tools to develop a clustering algorithm that reduces the within-cluster variance in the contextual variables and simultaneously maximizes the between-cluster variance in survey responses. The random forests algorithm then serves to assess the appropriateness of the resulting neighborhood typology and to select the key contextual variables in order to build a parsimonious model that makes a minimum of classification errors. Finally, for the impact analysis, propensity score matching methods are used to match the survey respondents of the pretest and posttest samples on age, gender, and their level of education for each neighborhood type identified within each city, before conducting a statistical test of the observed difference in the outcome measures. Moreover, all significant results were subjected to a sensitivity analysis to assess the robustness of these findings in the face of potential bias due to some unobserved covariates. The study finds that over the last fifteen years, all five police departments have undertaken major reforms of their internal organization and operating strategies and forged strategic partnerships in order to implement community policing. The resulting neighborhood typology reduced the within-cluster variance of the contextual variables and accounted for a significant share of the between-cluster variance in the outcome measures prior to treatment, suggesting that geocomputational methods help to balance the observed covariates and hence to reduce threats to the internal validity of an observational design. Finally, the impact analysis revealed that fear of crime dropped significantly over the 2000-2005 period in the neighborhoods in and around the urban centers of Bern and Zurich. These improvements are fairly robust in the face of bias due to some unobserved covariate and covary temporally and spatially with the implementation of community policing. The alternative hypothesis that the observed reductions in fear of crime were at least in part a result of community policing interventions thus appears at least as plausible as the null hypothesis of absolutely no effect, even if the observational design cannot completely rule out selection and regression to the mean as alternative explanations.

Pinning of quantized vortices in helium drops by dopant atoms and molecules

Using a density functional method, we investigate the properties of liquid 4He droplets doped with atoms (Ne and Xe) and molecules ( SF6 and hydrogen cyanide). We consider the case of droplets having a quantized vortex pinned to the dopant. A liquid-drop formula is proposed that accurately describes the total energy of the complex and allows one to extrapolate the density functional results to large N. For a given impurity, we find that the formation of a dopant+vortex+4HeN complex is energetically favored below a critical size Ncr. Our results support the possibility to observe quantized vortices in helium droplets by means of spectroscopic techniques.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.