REST is a hypoxia-responsive transcriptional repressor

Cellular exposure to hypoxia results in altered gene expression in a range of physiologic and pathophysiologic states. Discrete cohorts of genes can be either up- or down-regulated in response to hypoxia. While the Hypoxia-Inducible Factor (HIF) is the primary driver of hypoxia-induced adaptive gene expression, less is known about the signalling mechanisms regulating hypoxia-dependent gene repression. Using RNA-seq, we demonstrate that equivalent numbers of genes are induced and repressed in human embryonic kidney (HEK293) cells. We demonstrate that nuclear localization of the Repressor Element 1-Silencing Transcription factor (REST) is induced in hypoxia and that REST is responsible for regulating approximately 20% of the hypoxia-repressed genes. Using chromatin immunoprecipitation assays we demonstrate that REST-dependent gene repression is at least in part mediated by direct binding to the promoters of target genes. Based on these data, we propose that REST is a key mediator of gene repression in hypoxia.

Transcriptional regulation os phenylalaline biosynthesis and utilization

Conifer trees divert large quantities of carbon into the biosynthesis of phenylpropanoids, particularly to generate lignin, an important constituent of wood. Since phenylalanine is the precursor for phenylpropanoid biosynthesis, the precise regulation of phenylalanine synthesis and utilization should occur simultaneously. This crucial pathway is finely regulated primarily at the transcriptional level. Transcriptome analyses indicate that the transcription factors (TFs) preferentially expressed during wood formation in plants belong to the MYB and NAC families. Craven-Bartle et al. (2013) have shown in conifers that Myb8 is a candidate regulator of key genes in phenylalanine biosynthesis involved in the supply of the phenylpropane carbon skeleton necessary for lignin biosynthesis. This TF is able to bind AC elements present in the promoter regions of these genes to activate transcription. Constitutive overexpression of Myb8 in white spruce increased secondary-wall thickening and led to ectopic lignin deposition (Bomal et al. 2008). In Arabidopsis, the transcriptional network controlling secondary cell wall involves NAC-domain regulators operating upstream Myb transcription factors. Functional orthologues of members of this network described have been identified in poplar and eucalyptus, but in conifers functional evidence had only been obtained for MYBs. We have identified in the P. pinaster genome 37 genes encoding NAC proteins, which 3 NAC proteins could be potential candidates to be involved in vascular development (Pascual et al. 2015). The understanding of the transcriptional regulatory network associated to phenylpropanoids and lignin biosynthesis in conifers is crucial for future applications in tree improvement and sustainable forest management. This work is supported by the projects BIO2012-33797, BIO2015-69285-R and BIO-474 References: Bomal C, et al. (2008) Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis. J Exp Bot. 59: 3925-3939. Craven-Bartle B, et al. (2013) A Myb transcription factor regulates genes of the phenylalanine pathway in maritime pine. Plant J, 74: 755-766. Pascual MB, et al. (2015) The NAC transcription factor family in maritime pine (Pinus pinaster): molecular regulation of two genes involved in stress responses. BMC Plant Biol, 15: 254.

Phase transitions in gene networks evolved under different selection rules

Mathematical models of gene regulation are a powerful tool for understanding the complex features of genetic control. While various modeling efforts have been successful at explaining gene expression dynamics, much less is known about how evolution shapes the structure of these networks. An important feature of gene regulatory networks is their stability in response to environmental perturbations. Regulatory systems are thought to have evolved to exist near the transition between stability and instability, in order to have the required stability to environmental fluctuations while also being able to achieve a wide variety of functions (corresponding to different dynamical patterns). We study a simplified model of gene network evolution in which links are added via different selection rules. These growth models are inspired by recent work on `explosive' percolation which shows that when network links are added through competitive rather than random processes, the connectivity phase transition can be significantly delayed, and when it is reached, it appears to be first order (discontinuous, e.g., going from no failure at all to large expected failure) instead of second order (continuous, e.g., going from no failure at all to very small expected failure). We find that by modifying the traditional framework for networks grown via competitive link addition to capture how gene networks evolve to avoid damage propagation, we also see significant delays in the transition that depend on the selection rules, but the transitions always appear continuous rather than `explosive'.

Characterization Of A Protein-protein Interaction Network Of The Cbl-interacting Protein Kinase 8 From Sugarcane.

Plants are sessile organisms and have evolved to tolerate a constantly changing environment. After the onset of different stress conditions, calcineurin B-like (CBL) proteins can sense calcium signals and activate CBL-interacting protein kinase (CIPK) proteins, which can phosphorylate downstream proteins to reestablish plant homeostasis. Previous studies in the bioenergy crop sugarcane showed that the ScCIPK8 gene is induced by drought stress and is also related to sucrose content. Here, we have characterized the protein-protein interactions of ScCIPK8 with six CBL proteins (ScCBL1, ScCBL2, ScCBL3, ScCBL6, ScCBL9, and ScCBL10). Yeast two-hybrid assays showed that ScCIPK8 interacts with ScCBL1, ScCBL3, and ScCBL6. Bimolecular fluorescence complementation assays confirmed in planta the interactions that were observed in yeast cells. These findings give insights on the regulatory networks related to sugar accumulation and drought stress responses in sugarcane.

Chicken skeletal muscle-associated macroarray for gene discovery

Macro- and microarrays are well-established technologies to determine gene functions through repeated measurements of transcript abundance. We constructed a chicken skeletal muscle-associated array based on a muscle-specific EST database, which was used to generate a tissue expression dataset of similar to 4500 chicken genes across 5 adult tissues (skeletal muscle, heart, liver, brain, and skin). Only a small number of ESTs were sufficiently well characterized by BLAST searches to determine their probable cellular functions. Evidence of a particular tissue-characteristic expression can be considered an indication that the transcript is likely to be functionally significant. The skeletal muscle macroarray platform was first used to search for evidence of tissue-specific expression, focusing on the biological function of genes/transcripts, since gene expression profiles generated across tissues were found to be reliable and consistent. Hierarchical clustering analysis revealed consistent clustering among genes assigned to 'developmental growth', such as the ontology genes and germ layers. Accuracy of the expression data was supported by comparing information from known transcripts and tissue from which the transcript was derived with macroarray data. Hybridization assays resulted in consistent tissue expression profile, which will be useful to dissect tissue-regulatory networks and to predict functions of novel genes identified after extensive sequencing of the genomes of model organisms. Screening our skeletal-muscle platform using 5 chicken adult tissues allowed us identifying 43 'tissue-specific' transcripts, and 112 co-expressed uncharacterized transcripts with 62 putative motifs. This platform also represents an important tool for functional investigation of novel genes; to determine expression pattern according to developmental stages; to evaluate differences in muscular growth potential between chicken lines, and to identify tissue-specific genes.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background: Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor sigma(S) (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain. Results: Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms. Conclusions: The ppGpp and sigma(S) variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

Gene expression regulation and lineage evolution: the North and South tale of the hybrid polyploid Squalius alburnoides complex

The evolution of hybrid polyploid vertebrates, their viability and their perpetuation over evolutionary time have always been questions of great interest. However, little is known about the impact of hybridization and polyploidization on the regulatory networks that guarantee the appropriate quantitative and qualitative gene expression programme. The Squalius alburnoides complex of hybrid fish is an attractive system to address these questions, as it includes a wide variety of diploid and polyploid forms, and intricate systems of genetic exchange. Through the study of genome-specific allele expression of seven housekeeping and tissue-specific genes, we found that a gene copy silencing mechanism of dosage compensation exists throughout the distribution range of the complex. Here we show that the allele-specific patterns of silencing vary within the complex, according to the geographical origin and the type of genome involved in the hybridization process. In southern populations, triploids of S. alburnoides show an overall tendency for silencing the allele from the minority genome, while northern population polyploids exhibit preferential biallelic gene expression patterns, irrespective of genomic composition. The present findings further suggest that gene copy silencing and variable expression of specific allele combinations may be important processes in vertebrate polyploid evolution.

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.

Giving SBML support to ByoDyn

Report for the scientific stay at the California Institute of Technology during the summer of 2005. ByoDyn is a tool for simulating the dynamical expression of gene regulatory networks (GRNs) and for parameter estimation in uni- and multicellular models. A software support was carried out describing GRNs in the Systems Biology Markup Language (SBML). This one is a computer format for representing and storing computational models of biochemical pathways in software tools and databases. Supporting this format gives ByoDyn a wide range of possibilities to study the dynamical properties of multiple regulatory pathways.

Malignant transformation of DMBA/TPA-induced papillomas and nevi in the skin of mice selectively lacking retinoid-X-receptor alpha in epidermal keratinocytes.

Retinoid-X-receptor alpha (RXRalpha), a member of the nuclear receptor (NR) superfamily, is a ligand-dependent transcriptional regulatory factor. It plays a crucial role in NR signalling through heterodimerization with some 15 NRs. We investigated the role of RXRalpha and its partners on mouse skin tumor formation and malignant progression upon topical DMBA/TPA treatment. In mutants selectively ablated for RXRalpha in keratinocytes, epidermal tumors increased in size and number, and frequently progressed to carcinomas. As keratinocyte-selective peroxisome proliferator-activated receptor gamma (PPARgamma) ablation had similar effects, RXRalpha/PPARgamma heterodimers most probably mediate epidermal tumor suppression. Keratinocyte-selective RXRalpha-null and vitamin-D-receptor null mice also exhibited more numerous dermal melanocytic growths (nevi) than control mice, but only nevi from RXRalpha mutant mice progressed to invasive human-melanoma-like tumors. Distinct RXRalpha-mediated molecular events appear therefore to be involved, in keratinocytes, in cell-autonomous suppression of epidermal tumorigenesis and malignant progression, and in non-cell-autonomous suppression of nevi formation and progression. Our study emphasizes the crucial role of keratinocytes in chemically induced epidermal and melanocytic tumorigenesis, and raises the possibility that they could play a similar role in UV-induced tumorigenesis, notably in nevi formation and progression to melanoma.

The loss of circadian PAR bZip transcription factors results in epilepsy.

DBP (albumin D-site-binding protein), HLF (hepatic leukemia factor), and TEF (thyrotroph embryonic factor) are the three members of the PAR bZip (proline and acidic amino acid-rich basic leucine zipper) transcription factor family. All three of these transcriptional regulatory proteins accumulate with robust circadian rhythms in tissues with high amplitudes of clock gene expression, such as the suprachiasmatic nucleus (SCN) and the liver. However, they are expressed at nearly invariable levels in most brain regions, in which clock gene expression only cycles with low amplitude. Here we show that mice deficient for all three PAR bZip proteins are highly susceptible to generalized spontaneous and audiogenic epilepsies that frequently are lethal. Transcriptome profiling revealed pyridoxal kinase (Pdxk) as a target gene of PAR bZip proteins in both liver and brain. Pyridoxal kinase converts vitamin B6 derivatives into pyridoxal phosphate (PLP), the coenzyme of many enzymes involved in amino acid and neurotransmitter metabolism. PAR bZip-deficient mice show decreased brain levels of PLP, serotonin, and dopamine, and such changes have previously been reported to cause epilepsies in other systems. Hence, the expression of some clock-controlled genes, such as Pdxk, may have to remain within narrow limits in the brain. This could explain why the circadian oscillator has evolved to generate only low-amplitude cycles in most brain regions.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

A hierarchical approach to model parameter optimization for developmental systems.

In the context of Systems Biology, computer simulations of gene regulatory networks provide a powerful tool to validate hypotheses and to explore possible system behaviors. Nevertheless, modeling a system poses some challenges of its own: especially the step of model calibration is often difficult due to insufficient data. For example when considering developmental systems, mostly qualitative data describing the developmental trajectory is available while common calibration techniques rely on high-resolution quantitative data. Focusing on the calibration of differential equation models for developmental systems, this study investigates different approaches to utilize the available data to overcome these difficulties. More specifically, the fact that developmental processes are hierarchically organized is exploited to increase convergence rates of the calibration process as well as to save computation time. Using a gene regulatory network model for stem cell homeostasis in Arabidopsis thaliana the performance of the different investigated approaches is evaluated, documenting considerable gains provided by the proposed hierarchical approach.

Bristle pattern formation in the thorax formation of Drosophila: a systems level approach

A fundamental question in developmental biology is how tissues are patterned to give rise to differentiated body structures with distinct morphologies. The Drosophila wing disc offers an accessible model to understand epithelial spatial patterning. It has been studied extensively using genetic and molecular approaches. Bristle patterns on the thorax, which arise from the medial part of the wing disc, are a classical model of pattern formation, dependent on a pre-pattern of trans-activators and –repressors. Despite of decades of molecular studies, we still only know a subset of the factors that determine the pre-pattern. We are applying a novel and interdisciplinary approach to predict regulatory interactions in this system. It is based on the description of expression patterns by simple logical relations (addition, subtraction, intersection and union) between simple shapes (graphical primitives). Similarities and relations between primitives have been shown to be predictive of regulatory relationships between the corresponding regulatory factors in other Systems, such as the Drosophila egg. Furthermore, they provide the basis for dynamical models of the bristle-patterning network, which enable us to make even more detailed predictions on gene regulation and expression dynamics. We have obtained a data-set of wing disc expression patterns which we are now processing to obtain average expression patterns for each gene. Through triangulation of the images we can transform the expression patterns into vectors which can easily be analysed by Standard clustering methods. These analyses will allow us to identify primitives and regulatory interactions. We expect to identify new regulatory interactions and to understand the basic Dynamics of the regulatory network responsible for thorax patterning. These results will provide us with a better understanding of the rules governing gene regulatory networks in general, and provide the basis for future studies of the evolution of the thorax-patterning network in particular.

Sleep deprivation effects on circadian clock gene expression in the cerebral cortex parallel electroencephalographic differences among mouse strains.

Sleep deprivation (SD) results in increased electroencephalographic (EEG) delta power during subsequent non-rapid eye movement sleep (NREMS) and is associated with changes in the expression of circadian clock-related genes in the cerebral cortex. The increase of NREMS delta power as a function of previous wake duration varies among inbred mouse strains. We sought to determine whether SD-dependent changes in circadian clock gene expression parallel this strain difference described previously at the EEG level. The effects of enforced wakefulness of incremental durations of up to 6 h on the expression of circadian clock genes (bmal1, clock, cry1, cry2, csnk1epsilon, npas2, per1, and per2) were assessed in AKR/J, C57BL/6J, and DBA/2J mice, three strains that exhibit distinct EEG responses to SD. Cortical expression of clock genes subsequent to SD was proportional to the increase in delta power that occurs in inbred strains: the strain that exhibits the most robust EEG response to SD (AKR/J) exhibited dramatic increases in expression of bmal1, clock, cry2, csnkIepsilon, and npas2, whereas the strain with the least robust response to SD (DBA/2) exhibited either no change or a decrease in expression of these genes and cry1. The effect of SD on circadian clock gene expression was maintained in mice in which both of the cryptochrome genes were genetically inactivated. cry1 and cry2 appear to be redundant in sleep regulation as elimination of either of these genes did not result in a significant deficit in sleep homeostasis. These data demonstrate transcriptional regulatory correlates to previously described strain differences at the EEG level and raise the possibility that genetic differences underlying circadian clock gene expression may drive the EEG differences among these strains.