Differential regulation of horizontally acquired and core genome genes by the bacterial modulator H-NS

Horizontal acquisition of DNA by bacteria dramatically increases genetic diversity and hence successful bacterial colonization of several niches, including the human host. A relevant issue is how this newly acquired DNA interacts and integrates in the regulatory networks of the bacterial cell. The global modulator H-NS targets both core genome and HGT genes and silences gene expression in response to external stimuli such as osmolarity and temperature. Here we provide evidence that H-NS discriminates and differentially modulates core and HGT DNA. As an example of this, plasmid R27-encoded H-NS protein has evolved to selectively silence HGT genes and does not interfere with core genome regulation. In turn, differential regulation of both gene lineages by resident chromosomal H-NS requires a helper protein: the Hha protein. Tight silencing of HGT DNA is accomplished by H-NS-Hha complexes. In contrast, core genes are modulated by H-NS homoligomers. Remarkably, the presence of Hha-like proteins is restricted to the Enterobacteriaceae. In addition, conjugative plasmids encoding H-NS variants have hitherto been isolated only from members of the family. Thus, the H-NS system in enteric bacteria presents unique evolutionary features. The capacity to selectively discriminate between core and HGT DNA may help to maintain horizontally transmitted DNA in silent form and may give these bacteria a competitive advantage in adapting to new environments, including host colonization.

Gene Regulation in The Human Immune System. Gene Expression Signatures of Th2 Cell Differentiation, Type 1 Diabetes, and Intrauterine Immune Adaptation

The human immune system is constantly interacting with the surrounding stimuli and microorganisms. However, when directed against self or harmless antigens, these vital defense mechanisms can cause great damage. In addition, the understanding the underlying mechanism of several human diseases caused by aberrant immune cell functions, for instance type 1 diabetes and allergies, remains far from being complete. In this Ph.D. study these questions were addressed using genome-wide transcriptomic analyses. Asthma and allergies are characterized by a hyperactive response of the T helper 2 (Th2) immune cells. In this study, the target genes of the STAT6 transcription factor in naïve human T cells were identified with RNAi for the first time. STAT6 was shown to act as a central activator of the genes expression upon IL-4 signaling, with both direct and indirect effects on Th2 cell transcriptome. The core transcription factor network induced by IL-4 was identified from a kinetic analysis of the transcriptome. Type 1 diabetes is an autoimmune disease influenced by both the genetic susceptibility of an individual and the disease-triggering environmental factors. To improve understanding of the autoimmune processes driving pathogenesis in the prediabetic phase in humans, a unique series of prospective whole-blood RNA samples collected from HLA-susceptible children in the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) study was studied. Changes in different timewindows of the pathogenesis process were identified, and especially the type 1 interferon response was activated early and throughout the preclinical T1D. The hygiene hypothesis states that allergic diseases, and lately also autoimmune diseases, could be prevented by infections and other microbial contacts acquired in early childhood, or even prenatally. To study the effects of the standard of hygiene on the development of neonatal immune system, cord blood samples from children born in Finland (high standard of living), Estonia (rapid economic growth) and Russian Karelia (low standard of living) were compared. Children born in Russian Karelia deviated from Finnish and Estonian children in many aspects of the neonatal immune system, which was developmentally more mature in Karelia, resembling that of older infants. The results of this thesis offer significant new information on the regulatory networks associated with immune-mediated diseases in human. The results will facilitate understanding and further research on the role of the identified target genes and mechanisms driving the allergic inflammation and type 1 diabetes, hopefully leading to a new era of drug development.

Diversité fonctionnelle du facteur de transcription Tbx5 dans le coeur

Le cœur des vertébrés est un organe modulaire qui requiert le " patterning " complexe des champs morphogénétiques cardiogènes et la convergence coordonnée des diverses sous-populations de progéniteurs cardiogéniques. Au moins 7 facteurs de transcription de la famille T-box coopèrent au sein de ces nombreuses sous-populations de progéniteurs cardiogéniques afin de réguler la morphogenèse et l’agencement de multiples structures le long de l’ébauche cardiaque, ce qui explique que les mutations humaines de ces gènes engendrent diverses malformations congénitales cardiaques (MCCs). L’un de ces gènes T-box, Tbx5, dont l’haploinsuffisance génère le syndrome de Holt-Oram (SHO), intervient dans une grande variété de réseaux de régulation géniques (RRGs) qui orchestrent la morphogenèse des oreillettes, du ventricule gauche, de la valve mitrale, des septums inter-auriculaire et inter-ventriculaire, ainsi que du système de conduction cardiaque. La diversité des RRGs impliqués dans la formation de ces structures cardiaques suggère que Tbx5 détient une profusion de fonctions qui ne seront identifiables qu’en répertoriant ses activités moléculaires dans chaque lignée cardiaque examinée isolément. Afin d’aborder cette problématique, une ablation génétique de Tbx5 dans l’endocarde a été réalisée. Cette expérience a démontré le rôle crucial de Tbx5 dans la survie des cellules endocardiques bordant le septum primum et des cardiomyocytes au sein de cette structure embryonnaire qui contribuera à la morphogenèse du septum inter-auriculaire. En outre, cette étude a révélé l’existence d’une communication croisée entre la sous-population de cellules endocardiques Tbx5+ et le myocarde au niveau du septum primum, afin d’assurer la survie des cardiomyocytes, et ultimement de garantir la maturation du septum inter-auriculaire. Nos résultats confirment aussi l’importance de l’interdépendance génétique (Tbx5 et Gata4 ainsi que Tbx5 et Nos3) entre différents loci dans la morphogenèse de la cloison inter-auriculaire, et particulièrement de l’influence que peut avoir l’environnement sur la pénétrance et l’expressivité des communications inter-auriculaires (CIAs) dans le SHO. En outre, puisque les fonctions d’un gène dépendent ordinairement des différents isoformes qu’il peut générer, une deuxième étude a focalisé davantage sur l’aspect transcriptionnel de Tbx5. Cette approche a mené à la découverte de 6 transcrits alternatifs exhibant des fonctions à la fois communes et divergentes. La caractérisation de 2 de ces isoformes a révélé le rôle de l’isoforme long (Tbx5_v1) dans la régulation de la croissance des cardiomyocytes durant la cardiogénèse, tandis que l’isoforme court (Tbx5_v2), préférentiellement exprimé dans le cœur mature, réprime la croissance cellulaire. Il est donc entièrement concevable que les mutations de TBX5 entraînant une troncation de la région C-terminale accroissent la concentration d’une protéine mutée qui, à l’instar de Tbx5_v2, interfère avec la croissance de certaines structures cardiaques. En revanche, la divergence de fonctions de ces isoformes, caractérisée par les disparités de localisation subcellulaire et de d’interaction avec d’autres cofacteurs cardiaques, suggère que les mutations affectant davantage un isoforme favoriseraient l’émergence d’un type particulier de MCC. Finalement, un dernier objectif était d’identifier le ou les mécanisme(s) moléculaire(s) par le(s)quel(s) Tbx5 régule son principal gène cible, Nppa, et d’en extraire les indices qui éclairciraient sa fonction transcriptionnelle. Cet objectif nécessitait dans un premier lieu d’identifier les différents modules cis-régulateurs (MCRs) coordonnant la régulation transcriptionnelle de Nppa et Nppb, deux gènes natriurétiques dont l’organisation en tandem et le profil d’expression durant la cardiogénèse sont conservés dans la majorité des vertébrés. L’approche d’empreinte phylogénétique employée pour scanner le locus Nppb/Nppa a permis d’identifier trois MCRs conservés entre diverses espèces de mammifères, dont un (US3) est spécifique aux euthériens. Cette étude a corroboré que la régulation de l’expression du tandem génique Nppb/Nppa requérait l’activité transcriptionnelle d’enhancers en complément aux promoteurs de Nppa et Nppb. La concordance quasiment parfaite entre les profils d’expression de Tbx5 et de ces deux gènes natriurétiques chez les mammifères, suggère que le gradient d’expression ventriculaire de Tbx5 est interprété par le recrutement de ce facteur au niveau des différents enhancers identifiés. En somme, les études présentées dans cette thèse ont permis de clarifier la profusion de fonctions cardiaques que possède Tbx5. Certaines de ces fonctions émanent de l’épissage alternatif de Tbx5, qui favorise la synthèse d’isoformes dotés de propriétés spécifiques. Les diverses interactions combinatoires entre ces isoformes et d’autres facteurs cardiaques au sein des diverses sous-populations de progéniteurs cardiogènes contribuent à l’émergence de RRGs cardiaques divergents.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

A mathematical model of the sterol regulatory element binding protein 2 cholesterol biosynthesis pathway

Cholesterol is one of the key constituents for maintaining the cellular membrane and thus the integrity of the cell itself. In contrast high levels of cholesterol in the blood are known to be a major risk factor in the development of cardiovascular disease. We formulate a deterministic nonlinear ordinary differential equation model of the sterol regulatory element binding protein 2 (SREBP-2) cholesterol genetic regulatory pathway in an hepatocyte. The mathematical model includes a description of genetic transcription by SREBP-2 which is subsequently translated to mRNA leading to the formation of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR), a main precursor of cholesterol synthesis. Cholesterol synthesis subsequently leads to the regulation of SREBP-2 via a negative feedback formulation. Parameterised with data from the literature, the model is used to understand how SREBP-2 transcription and regulation affects cellular cholesterol concentration. Model stability analysis shows that the only positive steady-state of the system exhibits purely oscillatory, damped oscillatory or monotic behaviour under certain parameter conditions. In light of our findings we postulate how cholesterol homestasis is maintained within the cell and the advantages of our model formulation are discussed with respect to other models of genetic regulation within the literature.

GRASSIUS: A Platform for Comparative Regulatory Genomics across the Grasses

Transcription factors (TFs) are major players in gene regulatory networks and interactions between TFs and their target genes furnish spatiotemporal patterns of gene expression. Establishing the architecture of regulatory networks requires gathering information on TFs, their targets in the genome, and the corresponding binding sites. We have developed GRASSIUS (Grass Regulatory Information Services) as a knowledge-based Web resource that integrates information on TFs and gene promoters across the grasses. In its initial implementation, GRASSIUS consists of two separate, yet linked, databases. GrassTFDB holds information on TFs from maize (Zea mays), sorghum (Sorghum bicolor), sugarcane (Saccharum spp.), and rice (Oryza sativa). TFs are classified into families and phylogenetic relationships begin to uncover orthologous relationships among the participating species. This database also provides a centralized clearinghouse for TF synonyms in the grasses. GrassTFDB is linked to the grass TFome collection, which provides clones in recombination-based vectors corresponding to full-length open reading frames for a growing number of grass TFs. GrassPROMDB contains promoter and cis-regulatory element information for those grass species and genes for which enough data are available. The integration of GrassTFDB and GrassPROMDB will be accomplished through GrassRegNet as a first step in representing the architecture of grass regulatory networks. GRASSIUS can be accessed from www.grassius.org.

Reconstruction of molecular networks involved in citokine-induced myotubes atrophy integrating microRNA and mRNA expression

Pós-graduação em Ciências Biológicas (Genética) - IBB

The SUCEST-FUN regulatory network database: Designing an energy grass

Modern sugarcane cultivars are complex hybrids resulting from crosses among several Saccharum species. Traditional breeding methods have been employed extensively in different countries over the past decades to develop varieties with increased sucrose yield and resistance to pests and diseases. Conventional variety improvement, however, may be limited by the narrow pool of suitable genes. Thus, molecular genetics is seen as a promising tool to assist in the process of developing improved varieties. The SUCEST-FUN Project (http://sucest-fun.org) aims to associate function with sugarcane genes using a variety of tools, in particular those that enable the study of the sugarcane transcriptome. An extensive analysis has been conducted to characterise, phenotypically, sugarcane genotypes with regard to their sucrose content, biomass and drought responses. Through the analysis of different cultivars, genes associated with sucrose content, yield, lignin and drought have been identified. Currently, tools are being developed to determine signalling and regulatory networks in grasses, and to sequence the sugarcane genome, as well as to identify sugarcane promoters. This is being implemented through the SUCEST-FUN (http://sucest-fun.org) and GRASSIUS databases (http://grassius.org), the cloning of sugarcane promoters, the identification of cis-regulatory elements (CRE) using Chromatin Immunoprecipitation-sequencing (ChIP-Seq) and the generation of a comprehensive Signal Transduction and Transcription gene catalogue (SUCAST Catalogue).

Constraint-based analysis of gene interactions using restricted boolean networks and time-series data

Abstract Background A popular model for gene regulatory networks is the Boolean network model. In this paper, we propose an algorithm to perform an analysis of gene regulatory interactions using the Boolean network model and time-series data. Actually, the Boolean network is restricted in the sense that only a subset of all possible Boolean functions are considered. We explore some mathematical properties of the restricted Boolean networks in order to avoid the full search approach. The problem is modeled as a Constraint Satisfaction Problem (CSP) and CSP techniques are used to solve it. Results We applied the proposed algorithm in two data sets. First, we used an artificial dataset obtained from a model for the budding yeast cell cycle. The second data set is derived from experiments performed using HeLa cells. The results show that some interactions can be fully or, at least, partially determined under the Boolean model considered. Conclusions The algorithm proposed can be used as a first step for detection of gene/protein interactions. It is able to infer gene relationships from time-series data of gene expression, and this inference process can be aided by a priori knowledge available.

Automatic Design of Boolean Networks for Modelling Differentiation Trees

Real living cell is a complex system governed by many process which are not yet fully understood: the process of cell differentiation is one of these. In this thesis work we make use of a cell differentiation model to develop gene regulatory networks (Boolean networks) with desired differentiation dynamics. To accomplish this task we have introduced techniques of automatic design and we have performed experiments using various differentiation trees. The results obtained have shown that the developed algorithms, except the Random algorithm, are able to generate Boolean networks with interesting differentiation dynamics. Moreover, we have presented some possible future applications and developments of the cell differentiation model in robotics and in medical research. Understanding the mechanisms involved in biological cells can gives us the possibility to explain some not yet understood dangerous disease, i.e the cancer. Le cellula è un sistema complesso governato da molti processi ancora non pienamente compresi: il differenziamento cellulare è uno di questi. In questa tesi utilizziamo un modello di differenziamento cellulare per sviluppare reti di regolazione genica (reti Booleane) con dinamiche di differenziamento desiderate. Per svolgere questo compito abbiamo introdotto tecniche di progettazione automatica e abbiamo eseguito esperimenti utilizzando vari alberi di differenziamento. I risultati ottenuti hanno mostrato che gli algoritmi sviluppati, eccetto l'algoritmo Random, sono in grado di poter generare reti Booleane con dinamiche di differenziamento interessanti. Inoltre, abbiamo presentato alcune possibili applicazioni e sviluppi futuri del modello di differenziamento in robotica e nella ricerca medica. Capire i meccanismi alla base del funzionamento cellulare può fornirci la possibilità di spiegare patologie ancora oggi non comprese, come il cancro.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

Differential dynamic properties of scleroderma fibroblasts in response to perturbation of environmental stimuli.

Diseases are believed to arise from dysregulation of biological systems (pathways) perturbed by environmental triggers. Biological systems as a whole are not just the sum of their components, rather ever-changing, complex and dynamic systems over time in response to internal and external perturbation. In the past, biologists have mainly focused on studying either functions of isolated genes or steady-states of small biological pathways. However, it is systems dynamics that play an essential role in giving rise to cellular function/dysfunction which cause diseases, such as growth, differentiation, division and apoptosis. Biological phenomena of the entire organism are not only determined by steady-state characteristics of the biological systems, but also by intrinsic dynamic properties of biological systems, including stability, transient-response, and controllability, which determine how the systems maintain their functions and performance under a broad range of random internal and external perturbations. As a proof of principle, we examine signal transduction pathways and genetic regulatory pathways as biological systems. We employ widely used state-space equations in systems science to model biological systems, and use expectation-maximization (EM) algorithms and Kalman filter to estimate the parameters in the models. We apply the developed state-space models to human fibroblasts obtained from the autoimmune fibrosing disease, scleroderma, and then perform dynamic analysis of partial TGF-beta pathway in both normal and scleroderma fibroblasts stimulated by silica. We find that TGF-beta pathway under perturbation of silica shows significant differences in dynamic properties between normal and scleroderma fibroblasts. Our findings may open a new avenue in exploring the functions of cells and mechanism operative in disease development.

The ribosome as novel target for small regulatory RNAs

Small non-protein-coding RNA (ncRNA) molecules have been recognized recently as major contributors to regulatory networks in controlling gene expression in a highly efficient manner. While the list of validated ncRNAs that regulate crucial cellular processes grows steadily, not a single ncRNA has been identified that directly interacts and regulates the ribosome during protein biosynthesis (with the notable exceptions of 7SL RNA and tmRNA). All of the recently discovered regulatory ncRNAs that act on translation (e.g. microRNAs, siRNAs or antisense RNAs) target the mRNA rather than the ribosome. This is unexpected, given the central position the ribosome plays during gene expression. Furthermore it is strongly assumed that the primordial translation system in the ‘RNA world’ most likely received direct regulatory input from ncRNA-like cofactors. The fundamental question that we would like to ask is: Does the ‘RNA world still communicate’ with the ribosome? To address this question, we have analyzed the small ncRNA interactomes of ribosomes of prokaryotic (H. volcanii, S. aureus) and unicellular eukaryotic model organisms. Deep-sequencing and subsequent bioinformatic analyses revealed thousands of putative ribosome-associated ncRNAs. For a subset of these ncRNA candidates we have gathered experimental evidence that they are expressed in a stress-dependent manner and indeed directly target the ribosome. In the archaeon H. volcanii a tRNA-derived fragment was identified to target the small ribosomal subunit upon alkaline stress in vitro and in vivo. As a consequence of ribosome binding, this tRNA-fragment reduces protein synthesis by interfering with the peptidyl transferase activity. Our data reveal the ribosome as a novel target for small regulatory ncRNAs in all domains of life. Ribosome-bound ncRNAs are capable of fine tuning translation and might represent a so far largely unexplored class of regulatory sRNAs.

The ribosome as novel target for small regulatory RNAs

Small non-protein-coding RNA (ncRNA) molecules have been recognized recently as major contributors to regulatory networks in controlling gene expression in a highly efficient manner. While the list of validated ncRNAs that regulate crucial cellular processes grows steadily, not a single ncRNA has been identified that directly interacts and regulates the ribosome during protein biosynthesis (with the notable exceptions of 7SL RNA and tmRNA). All of the recently discovered regulatory ncRNAs that act on translation (e.g. microRNAs, siRNAs or antisense RNAs) target the mRNA rather than the ribosome. This is unexpected, given the central position the ribosome plays during gene expression. Furthermore it is strongly assumed that the primordial translation system in the ‘RNA world’ most likely received direct regulatory input from ncRNA-like cofactors. The fundamental question that we would like to ask is: Does the ‘RNA world still communicate’ with the ribosome? To address this question, we have analyzed the small ncRNA interactomes of ribosomes of organisms from all three domains of life. Deep-sequencing and subsequent bioinformatic analyses revealed thousands of putative ribosome-associated ncRNAs.1,2 For a subset of these ncRNA candidates we have gathered experimental evidence that they are expressed in a stress-dependent manner and indeed directly target the ribosome. We show that some of these ribosome-bound small ncRNAs are capable of fine tuning protein synthesis in vitro and in vivo. Our data therefore reveal the ribosome as a novel target for small regulatory ncRNAs in all domains of life and suggest the existence of a so far largely unexplored mechanism of translation regulation.

An mRNA-derived ncRNA targets and regulates the ribosome

Small non-protein-coding RNA (ncRNA) molecules have been recognized recently as major contributors to regulatory networks in controlling gene expression in a highly efficient manner. While the list of validated ncRNAs that regulate crucial cellular processes grows steadily, not a single ncRNA has been identified that directly interacts and regulates the ribosome during protein biosynthesis (with the notable exceptions of 7SL RNA and tmRNA). All of the recently discovered regulatory ncRNAs that act on translation (e.g. microRNAs, siRNAs or antisense RNAs) target the mRNA rather than the ribosome. This is unexpected, given the central position the ribosome plays during gene expression. To investigate whether such a class of regulatory ncRNAs does exist we performed genomic screens for small ribosome-associated RNAs in various model organisms of all three domains [1,2]. Here we focus on the functional characterisation of an 18 nucleotide long ncRNA candidate derived from an open reading frame (ORF) of an annotated S. cerevisiae gene, which encodes a tRNA methyltransferase. Yeast cells lacking this tRNA methyltransferase showed clear growth defects in high salt containing media. Genetic analysis showed that the absence of the mRNA-derived ncRNA rather than the absence of the tRNA methyltransferase activity is responsible for the observed phenotype. Since we performed a screen for small ribosome-associated RNAs we examined the regulatory potential of the synthetic 18mer during translation in vitro and in vivo. Metabolic labeling experiments in the presence of the synthetic 18mer RNA revealed an inhibitory potential on the global protein biosynthesis rate. In vitro translation and northern blot analysis further strengthen the hypothesis, that this RNA is a ribosome-associated regulatory ncRNA. Our studies in pro- and eukaryotic model organisms reveal the ribosome as a novel target for small regulatory ncRNAs in all domains of life. Ribosome-bound ncRNAs are capable of fine tuning translation and might represent a so far largely unexplored class of regulatory ncRNAs.