The combinatorics of transcriptional regulation

The ability to regulate gene expression is of central importance for the adaptability of living organisms to changes in their internal and external environment. At the transcriptional level, binding of transcription factors (TFs) in the vicinity of promoters can modulate the rate at which transcripts are produced, and as such play an important role in gene regulation. TFs with regulatory action at multiple promoters is the rule rather than the exception, with examples ranging from TFs like the cAMP receptor protein (CRP) in E. coli that regulates hundreds of different genes, to situations involving multiple copies of the same gene, such as on plasmids, or viral DNA. When the number of TFs heavily exceeds the number of binding sites, TF binding to each promoter can be regarded as independent. However, when the number of TF molecules is comparable to the number of binding sites, TF titration will result in coupling ("entanglement") between transcription of different genes. The last few decades have seen rapid advances in our ability to quantitatively measure such effects, which calls for biophysical models to explain these data. Here we develop a statistical mechanical model which takes the TF titration effect into account and use it to predict both the level of gene expression and the resulting correlation in transcription rates for a general set of promoters. To test these predictions experimentally, we create genetic constructs with known TF copy number, binding site affinities, and gene copy number; hence avoiding the need to use free fit parameters. Our results clearly prove the TF titration effect and that the statistical mechanical model can accurately predict the fold change in gene expression for the studied cases. We also generalize these experimental efforts to cover systems with multiple different genes, using the method of mRNA fluorescence in situ hybridization (FISH). Interestingly, we can use the TF titration affect as a tool to measure the plasmid copy number at different points in the cell cycle, as well as the plasmid copy number variance. Finally, we investigate the strategies of transcriptional regulation used in a real organism by analyzing the thousands of known regulatory interactions in E. coli. We introduce a "random promoter architecture model" to identify overrepresented regulatory strategies, such as TF pairs which coregulate the same genes more frequently than would be expected by chance, indicating a related biological function. Furthermore, we investigate whether promoter architecture has a systematic effect on gene expression by linking the regulatory data of E. coli to genome-wide expression censuses.

Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

The prevalence and regulation of antisense transcripts in Schizosaccharomyces pombe.

A strand-specific transcriptome sequencing strategy, directional ligation sequencing or DeLi-seq, was employed to profile antisense transcriptome of Schizosaccharomyces pombe. Under both normal and heat shock conditions, we found that polyadenylated antisense transcripts are broadly expressed while distinct expression patterns were observed for protein-coding and non-coding loci. Dominant antisense expression is enriched in protein-coding genes involved in meiosis or stress response pathways. Detailed analyses further suggest that antisense transcripts are independently regulated with respect to their sense transcripts, and diverse mechanisms might be potentially involved in the biogenesis and degradation of antisense RNAs. Taken together, antisense transcription may have profound impacts on global gene regulation in S. pombe.

Quantitative regulation of FLC via coordinated transcriptional initiation and elongation

The basis of quantitative regulation of gene expression is still poorly understood. In Arabidopsis thaliana, quantitative variation in expression of FLOWERING LOCUS C (FLC) influences the timing of flowering. In ambient temperatures, FLC expression is quantitatively modulated by a chromatin silencing mechanism involving alternative polyadenylation of antisense transcripts. Investigation of this mechanism unexpectedly showed that RNA polymerase II (Pol II) occupancy changes at FLC did not reflect RNA fold changes. Mathematical modeling of these transcriptional dynamics predicted a tight coordination of transcriptional initiation and elongation. This prediction was validated by detailed measurements of total and chromatin-bound FLC intronic RNA, a methodology appropriate for analyzing elongation rate changes in a range of organisms. Transcription initiation was found to vary ∼ 25-fold with elongation rate varying ∼ 8- to 12-fold. Premature sense transcript termination contributed very little to expression differences. This quantitative variation in transcription was coincident with variation in H3K36me3 and H3K4me2 over the FLC gene body. We propose different chromatin states coordinately influence transcriptional initiation and elongation rates and that this coordination is likely to be a general feature of quantitative gene regulation in a chromatin context.

Hydrogen peroxide sensing, signaling and regulation of transcription factors

The regulatory mechanisms by which hydrogen peroxide (H2O2) modulates the activity of transcription factors in bacteria (OxyR and PerR), lower eukaryotes (Yap1, Maf1, Hsf1 and Msn2/4) and mammalian cells (AP-1, NRF2, CREB, HSF1, HIF-1, TP53, NF-κB, NOTCH, SP1 and SCREB-1) are reviewed. The complexity of regulatory networks increases throughout the phylogenetic tree, reaching a high level of complexity in mammalians. Multiple H2O2 sensors and pathways are triggered converging in the regulation of transcription factors at several levels: (1) synthesis of the transcription factor by upregulating transcription or increasing both mRNA stability and translation; (ii) stability of the transcription factor by decreasing its association with the ubiquitin E3 ligase complex or by inhibiting this complex; (iii) cytoplasm-nuclear traffic by exposing/masking nuclear localization signals, or by releasing the transcription factor from partners or from membrane anchors; and, (iv) DNA binding and nuclear transactivation by modulating transcription factor affinity towards DNA, co-activators or repressors, and by targeting specific regions of chromatin to activate individual genes. We also discuss how H2O2 biological specificity results from diverse thiol protein sensors, with different reactivity of their sulfhydryl groups towards H2O2, being activated by different concentrations and times of exposure to H2O2. The specific regulation of local H2O2 concentrations is also crucial and results from H2O2 localized production and removal controlled by signals. Finally, we formulate equations to extract from typical experiments quantitative data concerning H2O2 reactivity with sensor molecules. Rate constants of 140 M-1s−1 and ≥ 1.3 × 103 M-1s−1 were estimated, respectively, for the reaction of H2O2 with KEAP1 and with an unknown target that mediates NRF2 protein synthesis. In conclusion, the multitude of H2O2 targets and mechanisms provides an opportunity for highly specific effects on gene regulation that depend on the cell type and on signals received from the cellular microenvironment.

Globin Gene Expression: Role of Transcription Factors

La dérégulation de l'expression génétique est une base pathophysiologique de plusieurs maladies. On a utilisé le locus du gène β-globine humain comme modèle pour élucider le mécanisme de régulation de la transcription génétique et évaluer son expression génétique durant l'érythropoïèse. La famille des protéines 'E' est composée de facteurs de transcription qui possèdent plusieurs sites de liaison au sein de locus du gène β-globine, suggérant leur rôle potentiel dans la régulation de l’expression de ces gènes. Nous avons montré que les facteurs HEB, E2A et ETO2 interagissent d’une manière significative avec la région contrôle du Locus (LCR) et avec les promoteurs des gènes de la famille β-globine. Le recrutement de ces facteurs au locus est modifié lors de l'érythropoïèse dans les cellules souches hematopoitiques et les cellules erythroides de souris transgéniques pour le locus de la β-globine humain, ainsi que dans les cellules progénitrices hématopoïétiques humaines. De plus par cette étude, nous démontrons pour la première fois que le gène β-globine humain est dans une chromatine active et qu’il interagit avec des facteurs de transcriptions de type suppresseurs dans les cellules progénitrices lymphoïdes (voie de différentiation alternative). Cette étude a aussi été faite dans des souris ayant une génétique mutante caractérisée par l'absence des facteurs de transcription E2A ou HEB.

“Reconstruction of Gene Regulatory Network from Expression Profile of Plasma RNA Data of Colorectal Cancer Patients using Soft Computing Techniques”

Microarray data analysis is one of data mining tool which is used to extract meaningful information hidden in biological data. One of the major focuses on microarray data analysis is the reconstruction of gene regulatory network that may be used to provide a broader understanding on the functioning of complex cellular systems. Since cancer is a genetic disease arising from the abnormal gene function, the identification of cancerous genes and the regulatory pathways they control will provide a better platform for understanding the tumor formation and development. The major focus of this thesis is to understand the regulation of genes responsible for the development of cancer, particularly colorectal cancer by analyzing the microarray expression data. In this thesis, four computational algorithms namely fuzzy logic algorithm, modified genetic algorithm, dynamic neural fuzzy network and Takagi Sugeno Kang-type recurrent neural fuzzy network are used to extract cancer specific gene regulatory network from plasma RNA dataset of colorectal cancer patients. Plasma RNA is highly attractive for cancer analysis since it requires a collection of small amount of blood and it can be obtained at any time in repetitive fashion allowing the analysis of disease progression and treatment response.

Divergence Involving Global Regulatory Gene Mutations in an Escherichia coli Population Evolving under Phosphate Limitation

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA-dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

MKK3/6-p38 MAPK negatively regulates murine MMP-13 gene expression induced by IL-1 beta and TNF-alpha in immortalized periodontal ligament fibroblasts

Matrix metalloprotease-13 (MMP-13) or collagenase-3 is involved in a number of pathologic processes such as tumor metastasis and angiogenesis, osteoarthritis, rheumatoid arthritis and periodontal diseases. These conditions are associated with extensive degradation of both connective tissue and bone. This report examines gene regulation mechanisms and signal transduction pathways involved in Mmp-13 expression induced by proinflammatory cytokines in periodontal ligament (PDL) fibroblasts. Mmp-13 mRNA expression was increased 10.7 and 9.5 fold after stimulation with IL-1 beta (5 ng/mL) and TNF-alpha (10 ng/mL), respectively. However, inhibition of p38 MAPKinase with SB203580 resulted in significant (p < 0.001) induction (23.2 and 18.1 fold, respectively) of Mmp-13 mRNA as assessed by real time PCR. Negative regulation of IL-1 induced Mmp-13 expression was confirmed by inhibiting p38 MAPK gene expression with siRNA. Transient transfection of dominant negative forms of MKK3 and MKK6 also resulted in increased levels of Mmp-13 mRNA after IL-1 beta stimulation. Mmp-13 mRNA expression induced by TNF-alpha was decreased by JNK and ERK inhibition. Western blot and zymogram analysis indicated that Mmp-13 protein expression induced by the proinflammatory cytokines were also upregulated by inhibition of p38 MAPK. Reporter gene experiments using stable cell lines harboring 660-bp sequence of the murine Mmp-13 proximal promoter indicated that transcriptional mechanisms were at least partially involved in this negative regulation of Mmp-13 expression by p38 MAPK and upstream MKK3/6. These results suggest a negative transcriptional regulatory mechanism mediated by p38 MAPK and upstream MKK3/6 on Mmp-13 expression induced by proinflammatory cytokines in PDL fibroblasts. (c) 2005 Elsevier B.V./International Society of Matrix Biology. All rights reserved.

Effects of BST and high energy diet on gene expression in mammary parenchyma of dairy heifers

The objective of this study was to determine the effects of dietary energy and recombinant bovine somatotropin (bST) injection to identify genes that might control mammogenesis. Total RNA was extracted from the parenchymal tissue of 32 heifers randomly assigned to one of four treatments: two diets (a standard diet and a high energy, high protein diet), each with or without bST. To perform microarray experiments, RNA samples were pooled (2 animals/pool) before reverse transcription and labeling with Cy3 or Cy5. A 4-node loop design was used to examine the differential gene expression among treatments using a bovine-specific cDNA micro array (National Bovine Functional Genomics Consortium Library, NBFGC) containing 18,263 unique expressed sequence tags (EST). Significance levels of differential gene expression among treatments were assessed using a mixed model approach. Injection of bST altered the expression of 12 % of the genes on NBFGC slide related to tissue development, whereas 6% were altered by diet. Administration of bST increases the expression of genes positively related to cell proliferation and mammary parenchyma to a greater extent than a high energy diet. © 2013 Sociedade Brasileira de Zootecnia.

Helicobacter pylori infection: Host immune response, implications on gene expression and microRNAs

Helicobacter pylori (H. pylori) infection is the most common bacterial infection worldwide. Persistent infection of the gastric mucosa leads to inflammatory processes and may remain silent for decades or progress causing more severe diseases, such as gastric adenocarcinoma. The clinical consequences of H. pylori infection are determined by multiple factors, including host genetic predisposition, gene regulation, environmental factors and heterogeneity of H. pylori virulence factors. After decades of studies of this successful relationship between pathogen and human host, various mechanisms have been elucidated. In this review, we have made an introduction on H. pylori infection and its virulence factors, and focused mainly on modulation of host immune response triggered by bacteria, changes in the pattern of gene expression in H. pylori-infected gastric mucosa, with activation of gene transcription involved in defense mechanisms, inflammatory and immunological response, cell proliferation and apoptosis. We also highlighted the role of bacteria eradication on gene expression levels. In addition, we addressed the recent involvement of different microRNAs in precancerous lesions, gastric cancer, and inflammatory processes induced by bacteria. New discoveries in this field may allow a better understanding of the role of major factors involved in the pathogenic mechanisms of H. pylori. (C) 2014 Baishideng Publishing Group Co., Limited. All rights reserved.

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

Xylella fastidiosa inhabits the plant xylem, a nutrient-poor environment, so that mechanisms to sense and respond to adverse environmental conditions are extremely important for bacterial survival in the plant host. Although the complete genome sequences of different Xylella strains have been determined, little is known about stress responses and gene regulation in these organisms. In this work, a DNA microarray was constructed containing 2,600 ORFs identified in the genome sequencing project of Xylella fastidiosa 9a5c strain, and used to check global gene expression differences in the bacteria when it is infecting a symptomatic and a tolerant citrus tree. Different patterns of expression were found in each variety, suggesting that bacteria are responding differentially according to each plant xylem environment. The global gene expression profile was determined and several genes related to bacterial survival in stressed conditions were found to be differentially expressed between varieties, suggesting the involvement of different strategies for adaptation to the environment. The expression pattern of some genes related to the heat shock response, toxin and detoxification processes, adaptation to atypical conditions, repair systems as well as some regulatory genes are discussed in this paper. DNA microarray proved to be a powerful technique for global transcriptome analyses. This is one of the first studies of Xylella fastidiosa gene expression in vivo which helped to increase insight into stress responses and possible bacterial survival mechanisms in the nutrient-poor environment of xylem vessels.

Innate Immunity in Insects, Function and Regulation of Hemolin from Hyalophora cecropia

Insects are useful models for the study of innate immune reactions and development. The distinction between recognition mechanisms preceding the breakdown of apoptotic cells during metamorphosis, and the breakdown of cells in response to infections, is unclear. Hemolin, a Lepidopteran member of the immunoglobulin superfamily, is a candidate molecule in self/nonself recognition. This thesis investigates hemolin function and hemolin gene regulation at a molecular level. We investigated the binding and cell adhesion properties of hemolin from H. cecropia and demonstrated that the proteins could homodimerize in presence of calcium. Moreover, a higher molecular weight membrane form of hemolin was present on hemocytes. These results, taken together with an earlier finding that soluble hemolin inhibits hemocyte adhesion, indicated that the secreted hemolin could modulate hemocyte aggregation in a competitive manner in the blood. In addition, hemolin was expressed in different tissues and at different developmental stages. Since hemolin is expressed both during development and during the immune response, its different regulatory factors must act in concert. We found that the third intron contains an enhancer, through which Dif, C/EBP and HMGI synergistically activate a reporter construct in vitro. We concluded that the enhancer is used during infection, since the κB-site is crucial for an immune response. Interestingly, we also found that the active form of the steroid hormone, ecdysone, induces the hemolin gene transcription in vivo, and in addition, acts synergistically during bacterial infection. Preliminary in vivo results indicate a secondary effect of ecdysone and the importance of hormone receptor elements in the upstream promoter region of hemolin. To explore the use of Drosophila as a genetic tool for understanding hemolin function and regulation, we sought to isolate the functional homologue in this species. A fly cDNA library in yeast was screened using H. cecropia hemolin as bait. The screen was not successful. However, it did lead to the discovery of a Drosophila protein with true binding specificity for hemolin. Subsequent characterization revealed a new, highly conserved gene, which we named yippee. Yippee is distantly related to zinc finger proteins and represents a novel family of proteins present in numerous eukaryotes, including fungi, plants and humans. Notably, when the Drosophila genome sequence was revealed, no hemolin orthologue could be detected. Finally, an extensive Drosophila genome chip analysis was initiated. The goal was to investigate the Drosophila immune response, and, in contrast to earlier studies of artificially injected flies, to examine a set of natural microbes, orally and externally applied. In parallel experiments viruses, bacteria, fungi and parasites were compared to unchallenged controls. We obtained a unique set of genes that were up-regulated in the response to the parasite Octosporea muscadomesticae and to the fungus Beauveria bassiana. We expect both down-regulated and up-regulated genes to serve as a source for the discovery of new effector molecules, in particular those that are active against parasites and fungi.