A census of human transcription factors: function, expression and evolution

Integrated approaches reveal determinants of genome-wide binding and function of the transcription factor Pho4

DNA sequences with high affinity for transcription factors occur more frequently in the genome than instances of genes bound or regulated by these factors. It is not clear what factors determine the genome-wide pattern of binding or regulation for a given transcription factor. We used an integrated approach to study how trans influences shape the binding and regulatory landscape of Pho4, a budding yeast transcription factor activated in response to phosphate limitation. We find that nucleosomes significantly restrict Pho4 binding. At nucleosome-depleted sites, competition from another transcription factor, Cbf1, determines Pho4 occupancy, raising the threshold for transcriptional activation in phosphate replete conditions and preventing Pho4 activation of genes outside the phosphate regulon during phosphate starvation. Pho4 binding is not sufficient for transcriptional activation-a cooperative interaction between Pho2 and Pho4 specifies genes that are activated. Combining these experimental observations, we are able to globally predict Pho4 binding and its functionality.

Transcriptional regulatory dynamics of the hypothalamic-pituitary-gonadal axis and its peripheral pathways as impacted by the 3-beta HSD inhibitor trilostane in zebrafish (Danio rerio)

To study mechanisms underlying generalized effects of 3β hydroxysteroid dehydrogenase (HSD3B) inhibition, reproductively mature zebrafish (Danio rerio) were exposed to trilostane at two dosages for 24, 48, or 96 h and their gonadal RNA samples profiled with Agilent zebrafish microarrays. Trilostane had substantial impact on the transcriptional dynamics of zebrafish, as reflected by a number of differentially expressed genes (DEGs) including transcription factors (TFs), altered TF networks, signaling pathways, and Gene Ontology (GO) biological processes. Changes in gene expression between a treatment and its control were mostly moderate, ranging from 1.3 to 2.0 fold. Expression of genes coding for HSD3B and many of its transcriptional regulators remained unchanged, suggesting transcriptional up-regulation is not a primary compensatory mechanism for HSD3B enzyme inhibition. While some trilostane-responsive TFs appear to share cellular functions linked to endocrine disruption, there are also many other DEGs not directly linked to steroidogenesis. Of the 65 significant TF networks, little similarity, and therefore little cross-talk, existed between them and the hypothalamic-pituitary-gonadal (HPG) axis. The most enriched GO biological processes are regulations of transcription, phosphorylation, and protein kinase activity. Most of the impacted TFs and TF networks are involved in cellular proliferation, differentiation, migration, and apoptosis. While these functions are fairly broad, their underlying TF networks may be useful to development of generalized toxicological screening methods. These findings suggest that trilostane-induced effects on fish endocrine functions are not confined to the HPG-axis alone. Its impact on corticosteroid synthesis could also have contributed to some system wide transcriptional changes in zebrafish observed in this study.

Direct inference of protein-DNA interactions using compressed sensing methods

Compressed sensing has revolutionized signal acquisition, by enabling complex signals to be measured with remarkable fidelity using a small number of so-called incoherent sensors. We show that molecular interactions, e.g., protein-DNA interactions, can be analyzed in a directly analogous manner and with similarly remarkable results. Specifically, mesoscopic molecular interactions act as incoherent sensors that measure the energies of microscopic interactions between atoms. We combine concepts from compressed sensing and statistical mechanics to determine the interatomic interaction energies of a molecular system exclusively from experimental measurements, resulting in a "de novo" energy potential. In contrast, conventional methods for estimating energy potentials are based on theoretical models premised on a priori assumptions and extensive domain knowledge. We determine the de novo energy potential for pairwise interactions between protein and DNA atoms from (i) experimental measurements of the binding affinity of protein-DNA complexes and (ii) crystal structures of the complexes. We show that the de novo energy potential can be used to predict the binding specificity of proteins to DNA with approximately 90% accuracy, compared to approximately 60% for the best performing alternative computational methods applied to this fundamental problem. This de novo potential method is directly extendable to other biomolecule interaction domains (enzymes and signaling molecule interactions) and to other classes of molecular interactions.

Genome-wide survey and expression analysis of the plant-specific NAC transcription factor family in soybean during development and dehydration stress

Plant-specific NAC transcription factors (TFs) play important roles in regulating diverse biological processes, including development, senescence, growth, cell division and responses to environmental stress stimuli. Within the soybean genome, we identified 152 full-length GmNAC TFs, including 11 membrane-bound members. In silico analysis of the GmNACs, together with their Arabidopsis and rice counterparts, revealed similar NAC architecture. Next, we explored the soybean Affymetrix array and Illumina transcriptome sequence data to analyse tissue-specific expression profiles of GmNAC genes. Phylogenetic analysis using stress-related NAC TFs from Arabidopsis and rice as seeding sequences identified 58 of the 152 GmNACs as putative stress-responsive genes, including eight previously reported dehydration-responsive GmNACs. We could design gene-specific primers for quantitative real-time PCR verification of 38 out of 50 newly predicted stress-related genes. Twenty-five and six GmNACs were found to be induced and repressed 2-fold or more, respectively, in soybean roots and/or shoots in response to dehydration. GmNAC085, whose amino acid sequence was 39%; identical to that of well-known SNAC1/ONAC2, was the most induced gene upon dehydration, showing 390-fold and 20-fold induction in shoots and roots, respectively. Our systematic analysis has identified excellent tissue-specific and/or dehydration-responsive candidate GmNAC genes for in-depth characterization and future development of improved drought-tolerant transgenic soybeans.

A gene-rich, transcriptionally active environment and the pre-deposition of repressive marks are predictive of susceptibility to KRAB/KAP1-mediated silencing

Background

KRAB-ZFPs (Krüppel-associated box domain-zinc finger proteins) are vertebrate-restricted transcriptional repressors encoded in the hundreds by the mouse and human genomes. They act via an essential cofactor, KAP1, which recruits effectors responsible for the formation of facultative heterochromatin. We have recently shown that KRAB/KAP1 can mediate long-range transcriptional repression through heterochromatin spreading, but also demonstrated that this process is at times countered by endogenous influences.

Method

To investigate this issue further we used an ectopic KRAB-based repressor. This system allowed us to tether KRAB/KAP1 to hundreds of euchromatic sites within genes, and to record its impact on gene expression. We then correlated this KRAB/KAP1-mediated transcriptional effect to pre-existing genomic and chromatin structures to identify specific characteristics making a gene susceptible to repression.

Results

We found that genes that were susceptible to KRAB/KAP1-mediated silencing carried higher levels of repressive histone marks both at the promoter and over the transcribed region than genes that were insensitive. In parallel, we found a high enrichment in euchromatic marks within both the close and more distant environment of these genes.

Conclusion

Together, these data indicate that high levels of gene activity in the genomic environment and the pre-deposition of repressive histone marks within a gene increase its susceptibility to KRAB/KAP1-mediated repression.

Endogenous tagging of the murine transcription factor Sox5 with hemaglutinin for functional studies

Gene expression is usually studied at the transcript level rather than at the protein level due to the lack of a specific and sensitive antibody. A way to overcome this is to fuse to the protein of interest an immunoreactive tag that has well-characterized antibodies. This epitope tagging approach is often used for in vitro experiments but for in vivo studies, the success rate of protein tagging has not been extensively analyzed and our study seeks to cover the void. A small epitope, hemaglutinin derived from the influenza virus was used to tag a transcription factor, Sox5 at the N-terminal via homologous recombination in the mouse. Sox5 is part of the Sry-related high-mobility-group box gene family and plays multiple roles in essential biological processes. Understanding of its molecular function in relation to its biological roles remains incomplete. In our study, we show that the longer isoform of Sox5 can be tagged endogenously with hemaglutinin without affecting its biological function in vivo. The tagged protein is easily and specifically detected with an anti-hemaglutinin antibody using immunohistochemistry with its expression matching the endogenous expression of Sox5. Immunoprecipitation of Sox5 was also carried out successfully using an anti-hemaglutinin antibody. The transgenic line generated from this study is predicted to be useful for future experiments such as co-immunoprecipitation and chromatin immunoprecipitation, allowing the further understanding of Sox5. Lastly, this approach can be easily employed for the investigation of other transcription factors and proteins in vivo to overcome technical limitations such as antibody cross-reactivity and to perform isoform-specific studies.

Simplified method to predict mutual interactions of human transcription factors based on their primary structure

Background

Physical interactions between transcription factors (TFs) are necessary for forming regulatory protein complexes and thus play a crucial role in gene regulation. Currently, knowledge about the mechanisms of these TF interactions is incomplete and the number of known TF interactions is limited. Computational prediction of such interactions can help identify potential new TF interactions as well as contribute to better understanding the complex machinery involved in gene regulation.

Methodology

We propose here such a method for the prediction of TF interactions. The method uses only the primary sequence information of the interacting TFs, resulting in a much greater simplicity of the prediction algorithm. Through an advanced feature selection process, we determined a subset of 97 model features that constitute the optimized model in the subset we considered. The model, based on quadratic discriminant analysis, achieves a prediction accuracy of 85.39% on a blind set of interactions. This result is achieved despite the selection for the negative data set of only those TF from the same type of proteins, i.e. TFs that function in the same cellular compartment (nucleus) and in the same type of molecular process (transcription initiation). Such selection poses significant challenges for developing models with high specificity, but at the same time better reflects real-world problems.

Conclusions

The performance of our predictor compares well to those of much more complex approaches for predicting TF and general protein-protein interactions, particularly when taking the reduced complexity of model utilisation into account.

In silico analysis of transcription factor repertoires and prediction of stress-responsive transcription factors from six major gramineae plants

The interactions between transcription factors (TFs) and cis-regulatory DNA sequences control gene expression, constituting the essential functional linkages of gene regulatory networks. The aim of this study is to identify and integrate all putative TFs from six grass species: Brachypodium distachyon, maize, rice, sorghum, barley, and wheat with significant information into an integrative database (GramineaeTFDB) for comparative genomics and functional genomics. For each TF, sequence features, promoter regions, domain alignments, GO assignment, FL-cDNA information, if available, and cross-references to various public databases and genetic resources are provided. Additionally, GramineaeTFDB possesses a tool which aids the users to search for putative cis-elements located in the promoter regions of TFs and predict the functions of the TFs using cis-element-based functional prediction approach. We also supplied hyperlinks to expression profiles of those TF genes of maize, rice, and barley, for which data are available. Furthermore, information about the availability of FOX and Ds mutant lines for rice and maize TFs, respectively, are also accessible through hyperlinks. Our study provides an important user-friendly public resource for functional analyses and comparative genomics of grass TFs, and understanding of the architecture of transcriptional regulatory networks and evolution of the TFs in agriculturally important cereal crops.

Gain, loss and divergence in primate zinc-finger genes: a rich resource for evolution of gene regulatory differences between species

The molecular changes underlying major phenotypic differences between humans and other primates are not well understood, but alterations in gene regulation are likely to play a major role. Here we performed a thorough evolutionary analysis of the largest family of primate transcription factors, the Krüppel-type zinc finger (KZNF) gene family. We identified and curated gene and pseudogene models for KZNFs in three primate species, chimpanzee, orangutan and rhesus macaque, to allow for a comparison with the curated set of human KZNFs. We show that the recent evolutionary history of primate KZNFs has been complex, including many lineage-specific duplications and deletions. We found 213 species-specific KZNFs, among them 7 human-specific and 23 chimpanzee-specific genes. Two human-specific genes were validated experimentally. Ten genes have been lost in humans and 13 in chimpanzees, either through deletion or pseudogenization. We also identified 30 KZNF orthologs with human-specific and 42 with chimpanzee-specific sequence changes that are predicted to affect DNA binding properties of the proteins. Eleven of these genes show signatures of accelerated evolution, suggesting positive selection between humans and chimpanzees. During primate evolution the most extensive re-shaping of the KZNF repertoire, including most gene additions, pseudogenizations, and structural changes occurred within the subfamily homininae. Using zinc finger (ZNF) binding predictions, we suggest potential impact these changes have had on human gene regulatory networks. The large species differences in this family of TFs stands in stark contrast to the overall high conservation of primate genomes and potentially represents a potent driver of primate evolution.

Unusually effective microRNA targeting within repeat-rich coding regions of mammalian mRNAs

MicroRNAs (miRNAs) regulate numerous biological processes by base-pairing with target messenger RNAs (mRNAs), primarily through sites in 3' untranslated regions (UTRs), to direct the repression of these targets. Although miRNAs have sometimes been observed to target genes through sites in open reading frames (ORFs), large-scale studies have shown such targeting to be generally less effective than 3' UTR targeting. Here, we show that several miRNAs each target significant groups of genes through multiple sites within their coding regions. This ORF targeting, which mediates both predictable and effective repression, arises from highly repeated sequences containing miRNA target sites. We show that such sequence repeats largely arise through evolutionary duplications and occur particularly frequently within families of paralogous C(2)H(2) zinc-finger genes, suggesting the potential for their coordinated regulation. Examples of ORFs targeted by miR-181 include both the well-known tumor suppressor RB1 and RBAK, encoding a C(2)H(2) zinc-finger protein and transcriptional binding partner of RB1. Our results indicate a function for repeat-rich coding sequences in mediating post-transcriptional regulation and reveal circumstances in which miRNA-mediated repression through ORF sites can be reliably predicted.

Interactome networks and human disease

Complex biological systems and cellular networks may underlie most genotype to phenotype relationships. Here, we review basic concepts in network biology, discussing different types of interactome networks and the insights that can come from analyzing them. We elaborate on why interactome networks are important to consider in biology, how they can be mapped and integrated with each other, what global properties are starting to emerge from interactome network models, and how these properties may relate to human disease.

Control of the embryonic stem cell state

Embryonic stem cells and induced pluripotent stem cells hold great promise for regenerative medicine. These cells can be propagated in culture in an undifferentiated state but can be induced to differentiate into specialized cell types. Moreover, these cells provide a powerful model system for studies of cellular identity and early mammalian development. Recent studies have provided insights into the transcriptional control of embryonic stem cell state, including the regulatory circuitry underlying pluripotency. These studies have, as a consequence, uncovered fundamental mechanisms that control mammalian gene expression, connect gene expression to chromosome structure, and contribute to human disease.

The characterisation of three types of genes that overlie copy number variable regions

Background

Due to the increased accuracy of Copy Number Variable region (CNV) break point mapping, it is now possible to say with a reasonable degree of confidence whether a gene (i) falls entirely within a CNV; (ii) overlaps the CNV or (iii) actually contains the CNV. We classify these as type I, II and III CNV genes respectively.

Principal Findings

Here we show that although type I genes vary in copy number along with the CNV, most of these type I genes have the same expression levels as wild type copy numbers of the gene. These genes must, therefore, be under homeostatic dosage compensation control. Looking into possible mechanisms for the regulation of gene expression we found that type I genes have a significant paucity of genes regulated by miRNAs and are not significantly enriched for monoallelically expressed genes. Type III genes, on the other hand, have a significant excess of genes regulated by miRNAs and are enriched for genes that are monoallelically expressed.

Significance

Many diseases and genomic disorders are associated with CNVs so a better understanding of the different ways genes are associated with normal CNVs will help focus on candidate genes in genome wide association studies.

Biological functions and transcriptional targets of CaRF in neurons

Calcium-response factor (CaRF) is a unique DNA-binding protein first recognized as a transcription factor for its role in modulating transcription of the gene encoding Brain-Derived Neurotrophic Factor (BDNF) in neurons. Here I review evidence for the biological functions and transcriptional targets of CaRF in the brain and discuss potential mechanisms by which calcium-activated signaling pathways may modulate CaRF-dependent transcription. These data paint an emerging picture of CaRF as a regulatory nexus for signaling pathways that control aspects of synaptic development and neuronal function.

DREME: motif discovery in transcription factor ChIP-seq data

MOTIVATION

Transcription factor (TF) ChIP-seq datasets have particular characteristics that provide unique challenges and opportunities for motif discovery. Most existing motif discovery algorithms do not scale well to such large datasets, or fail to report many motifs associated with cofactors of the ChIP-ed TF.

RESULTS

We present DREME, a motif discovery algorithm specifically designed to find the short, core DNA-binding motifs of eukaryotic TFs, and optimized to analyze very large ChIP-seq datasets in minutes. Using DREME, we discover the binding motifs of the the ChIP-ed TF and many cofactors in mouse ES cell (mESC), mouse erythrocyte and human cell line ChIP-seq datasets. For example, in mESC ChIP-seq data for the TF Esrrb, we discover the binding motifs for eight cofactor TFs important in the maintenance of pluripotency. Several other commonly used algorithms find at most two cofactor motifs in this same dataset. DREME can also perform discriminative motif discovery, and we use this feature to provide evidence that Sox2 and Oct4 do not bind in mES cells as an obligate heterodimer. DREME is much faster than many commonly used algorithms, scales linearly in dataset size, finds multiple, non-redundant motifs and reports a reliable measure of statistical significance for each motif found. DREME is available as part of the MEME Suite of motif-based sequence analysis tools (http://meme.nbcr.net).

Chronic exposure to anabolic steroids induces the muscle expression of oxytocin and a more than fiftyfold increase in circulating oxytocin in cattle

Molecular mechanisms in skeletal muscle associated with anabolic steroid treatment of cattle are unclear and we aimed to characterize transcriptional changes. Cattle were chronically exposed (68 ± 20 days) to a steroid hormone implant containing 200 mg trenbolone acetate and 20 mg estradiol (Revalor-H). Biopsy samples from 48 cattle (half treated) from longissimus dorsi (LD) muscle under local anesthesia were collected. Gene expression levels were profiled by microarray, covering 16,944 unique bovine genes: 121 genes were differentially expressed (DE) due to the implant (99.99% posterior probability of not being false positives). Among DE genes, a decrease in expression of a number of fat metabolism-associated genes, likely reflecting the lipid storage activity of intramuscular adipocytes, was observed. The expression of IGF1 and genes related to the extracellular matrix, slow twitch fibers, and cell cycle (including SOX8, a satellite cell marker) was increased in the treated muscle. Unexpectedly, a very large 21- (microarray) to 97 (real time quantitative PCR)-fold higher expression of the mRNA encoding the neuropeptide hormone oxytocin was observed in treated muscle. We also observed an ∼50-fold higher level of circulating oxytocin in the plasma of treated animals at the time of biopsy. Using a coexpression network strategy OXTR was identified as more likely than IGF1R to be a major mediator of the muscle response to Revalor-H. A re-investigation of in vivo cattle LD muscle samples during early to mid-fetal development identified a >128-fold increased expression of OXT, coincident with myofiber differentiation and fusion. We propose that oxytocin may be involved in mediating the anabolic effects of Revalor-H treatment.

Systematic characterization of protein-DNA interactions

Sequence-specific protein-DNA interactions (PDIs) are critical for regulating many cellular processes, including transcription, DNA replication, repair, and rearrangement. We review recent experimental advances in high-throughput technologies designed to characterize PDIs and discuss recent studies that use these tools, including ChIP-chip/seq, SELEX-based approaches, yeast one-hybrid, bacterial one-hybrid, protein binding microarray, and protein microarray. The results of these studies have challenged some long-standing concepts of PDI and provide valuable insights into the complex transcriptional regulatory networks.

Characterization of protein-DNA interactions using protein microarrays

INTRODUCTION

Protein-DNA interactions (PDIs) are critical for many cellular processes. We present here a protocol for the identification of PDIs in vitro using protein microarray technology. The procedure involves double-stranding synthesized DNA oligonucleotides with a fluorescent-labeled primer, binding the labeled double-stranded DNA directly to the protein microarray, and analyzing binding of the resulting PDIs. This approach provides simultaneous identification of PDIs for thousands of proteins, and multiple carefully designed DNA probes can be tested in parallel, which enables a rapid mapping of PDIs on a proteome-wide scale.

Unexpected complexity of the reef-building coral Acropora millepora transcription factor network

BACKGROUND

Coral reefs are disturbed on a global scale by environmental changes including rising sea surface temperatures and ocean acidification. Little is known about how corals respond or adapt to these environmental changes especially at the molecular level. This is mostly because of the paucity of genome-wide studies on corals and the application of systems approaches that incorporate the latter. Like in any other organism, the response of corals to stress is tightly controlled by the coordinated interplay of many transcription factors.

RESULTS

Here, we develop and apply a new system-wide approach in order to infer combinatorial transcription factor networks of the reef-building coral Acropora millepora. By integrating sequencing-derived transcriptome measurements, a network of physically interacting transcription factors, and phylogenetic network footprinting we were able to infer such a network. Analysis of the network across a phylogenetically broad sample of five species, including human, reveals that despite the apparent simplicity of corals, their transcription factors repertoire and interaction networks seem to be largely conserved. In addition, we were able to identify interactions among transcription factors that appear to be species-specific lending strength to the novel concept of "Taxonomically Restricted Interactions".

CONCLUSIONS

This study provides the first look at transcription factor networks in corals. We identified a transcription factor repertoire encoded by the coral genome and found consistencies of the domain architectures of transcription factors and conserved regulatory subnetworks across eumetazoan species, providing insight into how regulatory networks have evolved.

Chromatin: constructing the big picture

Chromatin is the ensemble of genomic DNA and a large number of proteins. Various genome-wide mapping techniques have begun to reveal that, despite the tremendous complexity, chromatin organization is governed by simple principles. This review discusses the principles that drive the spatial architecture of chromatin, as well as genome-wide-binding patterns of chromatin proteins.

A user's guide to the encyclopedia of DNA elements (ENCODE)

The mission of the Encyclopedia of DNA Elements (ENCODE) Project is to enable the scientific and medical communities to interpret the human genome sequence and apply it to understand human biology and improve health. The ENCODE Consortium is integrating multiple technologies and approaches in a collective effort to discover and define the functional elements encoded in the human genome, including genes, transcripts, and transcriptional regulatory regions, together with their attendant chromatin states and DNA methylation patterns. In the process, standards to ensure high-quality data have been implemented, and novel algorithms have been developed to facilitate analysis. Data and derived results are made available through a freely accessible database. Here we provide an overview of the project and the resources it is generating and illustrate the application of ENCODE data to interpret the human genome.

Sox6, jack of all trades: a versatile regulatory protein in vertebrate development

Approximately 20,000 genes are encoded in our genome, one tenth of which are thought to be transcription factors. Considering the complexity and variety of cell types generated during development, many transcription factors likely play multiple roles. Uncovering the versatile roles of Sox6 in vertebrate development sheds some light on how an organism efficiently utilizes the limited resources of transcription factors. The structure of the Sox6 gene itself may dictate its functional versatility. First, Sox6 contains no known regulatory domains; instead, it utilizes various cofactors. Second, Sox6 has a long 3'-UTR that contains multiple microRNA targets, thus its protein level is duly adjusted by cell type-specific microRNAs. Just combining these two characteristics alone makes Sox6 extremely versatile. To date, Sox6 has been reported to regulate differentiation of tissues of mesoderm, ectoderm, and endoderm origins, making Sox6 a truly multifaceted transcription factor.

What does biologically meaningful mean? A perspective on gene regulatory network validation

Gene regulatory networks (GRNs) are rapidly being delineated, but their quality and biological meaning are often questioned. Here, I argue that biological meaning is challenging to define and discuss reasons why GRN validation should be interpreted cautiously.

Advances in analysis of transcriptional regulatory networks

A transcriptional regulatory network represents a molecular framework in which developmental or environmental cues are transformed into differential expression of genes. Transcriptional regulation is mediated by the combinatorial interplay between cis-regulatory DNA elements and trans-acting transcription factors, and is perhaps the most important mechanism for controlling gene expression. Recent innovations, most notably the method for detecting protein-DNA interactions genome-wide, can help provide a comprehensive catalog of cis-regulatory elements and their interaction with given trans-acting factors in a given condition. A transcriptional regulatory network that integrates such information can lead to a systems-level understanding of regulatory mechanisms. In this review, we will highlight the key aspects of current knowledge on eukaryotic transcriptional regulation, especially on known transcription factors and their interacting regulatory elements. Then we will review some recent technical advances for genome-wide mapping of DNA-protein interactions based on high-throughput sequencing. Finally, we will discuss the types of biological insights that can be obtained from a network-level understanding of transcription regulation as well as future challenges in the field.

Alternative splicing and nonsense-mediated RNA decay contribute to the regulation of SHOX expression

The human SHOX gene is composed of seven exons and encodes a paired-related homeodomain transcription factor. SHOX mutations or deletions have been associated with different short stature syndromes implying a role in growth and bone formation. During development, SHOX is expressed in a highly specific spatiotemporal expression pattern, the underlying regulatory mechanisms of which remain largely unknown. We have analysed SHOX expression in diverse embryonic, fetal and adult human tissues and detected expression in many tissues that were not known to express SHOX before, e.g. distinct brain regions. By using RT-PCR and comparing the results with RNA-Seq data, we have identified four novel exons (exon 2a, 7-1, 7-2 and 7-3) contributing to different SHOX isoforms, and also established an expression profile for the emerging new SHOX isoforms. Interestingly, we found the exon 7 variants to be exclusively expressed in fetal neural tissues, which could argue for a specific role of these variants during brain development. A bioinformatical analysis of the three novel 3′UTR exons yielded insights into the putative role of the different 3′UTRs as targets for miRNA binding. Functional analysis revealed that inclusion of exon 2a leads to nonsense-mediated RNA decay altering SHOX expression in a tissue and time specific manner. In conclusion, SHOX expression is regulated by different mechanisms and alternative splicing coupled with nonsense-mediated RNA decay constitutes a further component that can be used to fine tune the SHOX expression level.

Motif-guided sparse decomposition of gene expression data for regulatory module identification

Background

Genes work coordinately as gene modules or gene networks. Various computational approaches have been proposed to find gene modules based on gene expression data; for example, gene clustering is a popular method for grouping genes with similar gene expression patterns. However, traditional gene clustering often yields unsatisfactory results for regulatory module identification because the resulting gene clusters are co-expressed but not necessarily co-regulated.

Results

We propose a novel approach, motif-guided sparse decomposition (mSD), to identify gene regulatory modules by integrating gene expression data and DNA sequence motif information. The mSD approach is implemented as a two-step algorithm comprising estimates of (1) transcription factor activity and (2) the strength of the predicted gene regulation event(s). Specifically, a motif-guided clustering method is first developed to estimate the transcription factor activity of a gene module; sparse component analysis is then applied to estimate the regulation strength, and so predict the target genes of the transcription factors. The mSD approach was first tested for its improved performance in finding regulatory modules using simulated and real yeast data, revealing functionally distinct gene modules enriched with biologically validated transcription factors. We then demonstrated the efficacy of the mSD approach on breast cancer cell line data and uncovered several important gene regulatory modules related to endocrine therapy of breast cancer.

Conclusion

We have developed a new integrated strategy, namely motif-guided sparse decomposition (mSD) of gene expression data, for regulatory module identification. The mSD method features a novel motif-guided clustering method for transcription factor activity estimation by finding a balance between co-regulation and co-expression. The mSD method further utilizes a sparse decomposition method for regulation strength estimation. The experimental results show that such a motif-guided strategy can provide context-specific regulatory modules in both yeast and breast cancer studies.

Identification of genes expressed in the angiosperm female gametophyte

Until recently, identification of gene regulatory networks controlling the development of the angiosperm female gametophyte has presented a significant challenge to the plant biology community. The angiosperm female gametophyte is fairly inaccessible because it is a highly reduced structure relative to the sporophyte and is embedded within multiple layers of the sporophytic tissue of the ovule. Moreover, although mutations affecting the female gametophyte can be readily isolated, their analysis can be difficult because most affect genes involved in basic cellular processes that are also required in the diploid sporophyte. In recent years, expression-based approaches in multiple species have begun to uncover gene sets expressed in specific female gametophyte cells as a means of identifying regulatory networks controlling cell differentiation in the female gametophyte. Here, recent efforts to identify and analyse gene expression programmes in the Arabidopsis female gametophyte are reviewed.

Evolution of vertebrate appendicular structures: Insight from genetic and palaeontological data

The new body of evidence from fossils and comparative-developmental analysis of subset of appendicular patterning genes has revealed that limb elements seen in tetrapods are assembled in fish fin over evolutionary time. However, despite of deep homology in basic structure and underlying developmental system, there remains a large morphological gap between distal elements of tetrapod limb and distal fin skeleton of tetrapodomorph fish. Understanding the genetic basis of major transformations in distal-limb morphology is the next challenge for evolutionary developmental biologists. Here by integrating data from fossils, comparative-developmental and genetic studies, models are proposed describing the evolution of cis-regulatory elements as a basis for diversification of appendicular architecture. Instead of emphasizing the subset of developmental genes, for instance Hoxd genes, the focus here is on the significance of elucidating cis-regulatory elements for multiple other key molecular players of limb/fin development and genetic/molecular interactions among them, for a better understanding of the developmental and genetic basis of limb evolution.

The cardiac transcription network modulated by Gata4, Mef2a, Nkx2.5, Srf, histone modifications, and microRNAs

The transcriptome, as the pool of all transcribed elements in a given cell, is regulated by the interaction between different molecular levels, involving epigenetic, transcriptional, and post-transcriptional mechanisms. However, many previous studies investigated each of these levels individually, and little is known about their interdependency. We present a systems biology study integrating mRNA profiles with DNA–binding events of key cardiac transcription factors (Gata4, Mef2a, Nkx2.5, and Srf), activating histone modifications (H3ac, H4ac, H3K4me2, and H3K4me3), and microRNA profiles obtained in wild-type and RNAi–mediated knockdown. Finally, we confirmed conclusions primarily obtained in cardiomyocyte cell culture in a time-course of cardiac maturation in mouse around birth. We provide insights into the combinatorial regulation by cardiac transcription factors and show that they can partially compensate each other's function. Genes regulated by multiple transcription factors are less likely differentially expressed in RNAi knockdown of one respective factor. In addition to the analysis of the individual transcription factors, we found that histone 3 acetylation correlates with Srf- and Gata4-dependent gene expression and is complementarily reduced in cardiac Srf knockdown. Further, we found that altered microRNA expression in Srf knockdown potentially explains up to 45% of indirect mRNA targets. Considering all three levels of regulation, we present an Srf-centered transcription network providing on a single-gene level insights into the regulatory circuits establishing respective mRNA profiles. In summary, we show the combinatorial contribution of four DNA–binding transcription factors in regulating the cardiac transcriptome and provide evidence that histone modifications and microRNAs modulate their functional consequence. This opens a new perspective to understand heart development and the complexity cardiovascular disorders.

An evolutionary conserved orchestra of transcription factors controls cardiac development and function. More recently the contributions of epigenetic and post-transcriptional mechanisms like histone modifications and microRNAs have been identified. The interplay between these regulatory mechanisms is still an open question. However, perturbations of the cardiac transcriptome, triggered by all three levels of regulation, are underlying cardiovascular disease such as congenital heart malformations. Here, we show the impact of the interdependencies of four key transcription factors (Gata4, Mef2a, Nkx2.5, and Srf) and the contribution of activating histone modifications and microRNAs on the cardiac transcriptome. We found that even these non-paralogous transcription factors can partially compensate each other's function. Our data show that histone 3 acetylation correlates with Srf- and Gata4- dependent gene activation. Moreover, we predict a large proportion of indirect Srf targets to be regulated by Srf-dependent microRNAs, which thus might represent an important intermediate layer of regulation. Taken together, we suggest that the different levels regulating cardiac mRNA profiles have a high degree of interdependency and the potential to buffer each other, which presents a starting point to understand the phenotypic variability typically seen in complex cardiovascular disorders.

Transcriptional repression: conserved and evolved features

The regulation of gene expression by transcriptional repression is an ancient and conserved mechanism that manifests itself in diverse ways. Here we summarize conserved pathways for transcriptional repression prevalent throughout all forms of life, as well as indirect mechanisms that appear to have originated in eukaryotes, consistent with the unique chromatin environment of eukaryotic genes. The direct interactions between transcriptional repressors and the core transcriptional machinery in bacteria and archaea are sufficient to generate a sophisticated suite of mechanisms that provide flexible control. These direct interactions contrast with the activity of corepressors, which provide an additional regulatory control in eukaryotes. Their modulation of chromatin structure represents an indirect pathway to downregulate transcription, and their diversity and modulation provide additional complexity suited to the requirements of elaborate eukaryotic repression patterns. New findings indicate that corepressors are not necessarily restricted to generating a single stereotypic output, but can rather exhibit diverse functional responses depending on the context in which they are recruited, providing a hitherto unsuspected additional source of diversity in transcriptional control. Mechanisms within eukaryotes appear to be highly conserved, with novel aspects chiefly represented by addition of lineage-specific corepressor scaffolds that provide additional opportunities for recruiting the same core machinery.

Constant splice-isoform ratios in human lymphoblastoid cells support the concept of a splico-stat

Splicing generates mature transcripts from genes in pieces in eukaryotic cells. Overwhelming evidence has accumulated that alternative routes in splicing are possible for most human and mammalian genes, thereby allowing formation of different transcripts from one gene. No function has been assigned to the majority of identified alternative splice forms, and it has been assumed that they compose inert or tolerated waste from aberrant or noisy splicing. Here we demonstrate that five human transcription units (WT1, NOD2, GNAS, RABL2A, RABL2B) have constant splice-isoform ratios in genetically diverse lymphoblastoid cell lines independent of the type of alternative splicing (exon skipping, alternative donor/acceptor, tandem splice sites) and gene expression level. Even splice events that create premature stop codons and potentially trigger nonsense-mediated mRNA decay are found at constant fractions. The analyzed alternative splicing events were qualitatively but not quantitatively conserved in corresponding chimpanzee cell lines. Additionally, subtle splicing at tandem acceptor splice sites (GNAS, RABL2A/B) was highly constrained and strongly depends on the upstream donor sequence content. These results also demonstrate that unusual and unproductive splice variants are produced in a regulated manner.

A catalogue of eukaryotic transcription factor types, their evolutionary origin, and species distribution

Transcription factors (TFs) play key roles in the regulation of gene expression by binding in a sequence-specific manner to genomic DNA. In eukaryotes, DNA binding is achieved by a wide range of structural forms and motifs. TFs are typically classified by their DNA-binding domain (DBD) type. In this chapter, we catalogue and survey 91 different TF DBD types in metazoa, plants, fungi, and protists. We briefly discuss well-characterized TF families representing the major DBD superclasses. We also examine the species distributions and inferred evolutionary histories of the various families, and the potential roles played by TF family expansion and dimerization.

Identification of Transcription Factor-DNA Interactions In Vivo

Recent technological developments have revolutionized our understanding of transcriptional regulation by providing an unprecedented ability to interrogate in vivo transcription factor binding. The combination of high-throughput sequencing with chromatin precipitation of transcription factors and specifically labeled histones has allowed direct protein-DNA contacts to be visualized across genomes as large and complex as mammals at base-pair resolution. This chapter reviews the developments that led to these insights, with particular focus on examples of early protein-DNA localization experiments using genomic microarrays in mammals and yeast. Four state-of-the-art research directions are highlighted as examples of previously unimaginable frontiers now under active investigation.

Methods for Analysis of Transcription Factor DNA-Binding Specificity In Vitro

Transcription of genes during development and in response to environmental stimuli is determined by genomic DNA sequence. The DNA sequences regulating transcription are read by sequence-specific transcription factors (TFs) that recognize relatively short sequences, generally between four and twenty base pairs in length. Transcriptional regulation generally requires binding of multiple TFs in close proximity to each other. Mechanistic understanding of transcription in an organism thus requires detailed knowledge of binding affinities of all its TFs to all possible DNA sequences, and the co-operative interactions between the TFs. However, very little is known about such co-operative binding interactions, and even the simple TF-DNA binding information exists only for a very small proportion of all TFs - for example, mammals have approximately 1,300-2,000 TFs [1, 2], yet the largest public databases for TF binding specificity, Jaspar and Uniprobe [3, 4] currently list only approximately 500 moderate to high resolution profiles for human or mouse. This lack of knowledge is in part due to the fact that analysis of TF DNA binding has been laborious and expensive. In this chapter, we review methods that can be used to determine binding specificity of TFs to DNA, mainly focusing on recently developed assays that allow high-resolution analysis of TF binding specificity in relatively high throughput.

Function and Evolution of C2H2 Zinc Finger Arrays

Krüppel-type or C2H2 zinc fingers represent a dominant DNA-binding motif in eukaryotic transcription factor (TF) proteins. In Krüppel-type (KZNF) TFs, KZNF motifs are arranged in arrays of three to as many as 40 tandem units, which cooperate to define the unique DNA recognition properties of the protein. Each finger contains four amino acids located at specific positions, which are brought into direct contact with adjacent nucleotides in the DNA sequence as the KZNF array winds around the major groove of the alpha helix. This arrangement creates an intimate and potentially predictable relationship between the amino acid sequence of KZNF arrays and the nucleotide sequence of target binding sites. The large number of possible combinations and arrangements of modular KZNF motifs, and the increasing lengths of KZNF arrays in vertebrate species, has created huge repertoires of functionally unique TF proteins. The properties of this versatile DNA-binding motif have been exploited independently many times over the course of evolution, through attachment to effector motifs that confer activating, repressing or other activities to the proteins. Once created, some of these novel inventions have expanded in specific evolutionary clades, creating large families of TFs that are lineage- or species-unique. This chapter reviews the properties and their remarkable evolutionary history of eukaryotic KZNF TF proteins, with special focus on large families that dominate the TF landscapes in different metazoan species.

Transcription factor effector domains

The last decade has seen an incredible breakthrough in technologies that allow histones, transcription factors (TFs), and RNA polymerases to be precisely mapped throughout the genome. From this research, it is clear that there is a complex interaction between the chromatin landscape and the general transcriptional machinery and that the dynamic control of this interface is central to gene regulation. However, the chromatin remodeling enzymes and general TFs cannot, on their own, recognize and stably bind to promoter or enhancer regions. Rather, they are recruited to cis regulatory regions through interaction with site-specific DNA binding TFs and/or proteins that recognize epigenetic marks such as methylated cytosines or specifically modified amino acids in histones. These "recruitment" factors are modular in structure, reflecting their ability to interact with the genome via one region of the protein and to simultaneously bind to other regulatory proteins via "effector" domains. In this chapter, we provide examples of common effector domains that can function in transcriptional regulation via their ability to (a) interact with the basal transcriptional machinery and general co-activators, (b) interact with other TFs to allow cooperative binding, and (c) directly or indirectly recruit histone and chromatin modifying enzymes.

ZNF274 recruits the histone methyltransferase SETDB1 to the 3' ends of ZNF genes

Only a small percentage of human transcription factors (e.g. those associated with a specific differentiation program) are expressed in a given cell type. Thus, cell fate is mainly determined by cell type-specific silencing of transcription factors that drive different cellular lineages. Several histone modifications have been associated with gene silencing, including H3K27me3 and H3K9me3. We have previously shown that genes for the two largest classes of mammalian transcription factors are marked by distinct histone modifications; homeobox genes are marked by H3K27me3 and zinc finger genes are marked by H3K9me3. Several histone methyltransferases (e.g. G9a and SETDB1) may be involved in mediating the H3K9me3 silencing mark. We have used ChIP-chip and ChIP-seq to demonstrate that SETDB1, but not G9a, is associated with regions of the genome enriched for H3K9me3. One current model is that SETDB1 is recruited to specific genomic locations via interaction with the corepressor TRIM28 (KAP1), which is in turn recruited to the genome via interaction with zinc finger transcription factors that contain a Kruppel-associated box (KRAB) domain. However, specific KRAB-ZNFs that recruit TRIM28 (KAP1) and SETDB1 to the genome have not been identified. We now show that ZNF274 (a KRAB-ZNF that contains 5 C2H2 zinc finger domains), can interact with KAP1 both in vivo and in vitro and, using ChIP-seq, we show that ZNF274 binding sites co-localize with SETDB1, KAP1, and H3K9me3 at the 3′ ends of zinc finger genes. Knockdown of ZNF274 with siRNAs reduced the levels of KAP1 and SETDB1 recruitment to the binding sites. These studies provide the first identification of a KRAB domain-containing ZNF that is involved in recruitment of the KAP1 and SETDB1 to specific regions of the human genome.

Mammal cells double their total RNAs against diabetes, ischemia reperfusion and malaria-induced oxidative stress

Total cellular RNA level is stable usually, although it may increase gradually during growth or decrease gradually under certain stressors. However, we found that mammal cell RNAs could be doubled within 24 h in response to free heme accumulation (ischemia reperfusion and malaria infection) or a high level of glucose treatment (diabetes). Clinical investigations in rats showed that pretreatment with heme (24 h for doubling total RNAs) alleviated oxidative damages caused by diabetes, and pretreatment with glucose (24 h for trebling total RNAs) alleviated oxidative damages caused by ischemia reperfusion or malaria infection. Therefore, this rapid RNA amplification may play an important role in mammal adaptation to diabetes, ischemia reperfusion and malaria infection-derived oxidative stress. This rapid RNA amplification is derived from glucose and heme, but not from their accompanying reactive oxygen species. Hexokinases endure glucose-derived reactive oxygen species accumulation but are not related glucose-derived RNA amplification. In contrast, the TATA box-binding protein (TBP) mediates all glucose- and heme-induced RNA amplification in mammal cells.

Genome-wide binding of the orphan nuclear receptor TR4 suggests its general role in fundamental biological processes

Background

The orphan nuclear receptor TR4 (human testicular receptor 4 or NR2C2) plays a pivotal role in a variety of biological and metabolic processes. With no known ligand and few known target genes, the mode of TR4 function was unclear.

Results

We report the first genome-wide identification and characterization of TR4 in vivo binding. Using chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq), we identified TR4 binding sites in 4 different human cell types and found that the majority of target genes were shared among different cells. TR4 target genes are involved in fundamental biological processes such as RNA metabolism and protein translation. In addition, we found that a subset of TR4 target genes exerts cell-type specific functions. Analysis of the TR4 binding sites revealed that less than 30% of the peaks from any of the cell types contained the DR1 motif previously derived from in vitro studies, suggesting that TR4 may be recruited to the genome via interaction with other proteins. A bioinformatics analysis of the TR4 binding sites predicted a cis regulatory module involving TR4 and ETS transcription factors. To test this prediction, we performed ChIP-seq for the ETS factor ELK4 and found that 30% of TR4 binding sites were also bound by ELK4. Motif analysis of the sites bound by both factors revealed a lack of the DR1 element, suggesting that TR4 binding at a subset of sites is facilitated through the ETS transcription factor ELK4. Further studies will be required to investigate the functional interdependence of these two factors.

Conclusions

Our data suggest that TR4 plays a pivotal role in fundamental biological processes across different cell types. In addition, the identification of cell type specific TR4 binding sites enables future studies of the pathways underlying TR4 action and its possible role in metabolic diseases.

Reprogramming cellular behavior with RNA controllers responsive to endogenous proteins

Synthetic genetic devices that interface with native cellular pathways can be used to change natural networks to implement new forms of control and behavior. The engineering of gene networks has been limited by an inability to interface with native components. We describe a class of RNA control devices that overcome these limitations by coupling increased abundance of particular proteins to targeted gene expression events through the regulation of alternative RNA splicing. We engineered RNA devices that detect signaling through the nuclear factor κB and Wnt signaling pathways in human cells and rewire these pathways to produce new behaviors, thereby linking disease markers to noninvasive sensing and reprogrammed cellular fates. Our work provides a genetic platform that can build programmable sensing-actuation devices enabling autonomous control over cellular behavior.

Accurate inference of transcription factor binding from DNA sequence and chromatin accessibility data

Accurate functional annotation of regulatory elements is essential for understanding global gene regulation. Here, we report a genome-wide map of 827,000 transcription factor binding sites in human lymphoblastoid cell lines, which is comprised of sites corresponding to 239 position weight matrices of known transcription factor binding motifs, and 49 novel sequence motifs. To generate this map, we developed a probabilistic framework that integrates cell- or tissue-specific experimental data such as histone modifications and DNase I cleavage patterns with genomic information such as gene annotation and evolutionary conservation. Comparison to empirical ChIP-seq data suggests that our method is highly accurate yet has the advantage of targeting many factors in a single assay. We anticipate that this approach will be a valuable tool for genome-wide studies of gene regulation in a wide variety of cell types or tissues under diverse conditions.

High-resolution genome-wide in vivo footprinting of diverse transcription factors in human cells

Regulation of gene transcription in diverse cell types is determined largely by varied sets of cis-elements where transcription factors bind. Here we demonstrate that data from a single high-throughput DNase I hypersensitivity assay can delineate hundreds of thousands of base-pair resolution in vivo footprints in human cells that precisely mark individual transcription factor-DNA interactions. These annotations provide a unique resource for the investigation of cis-regulatory elements. We find that footprints for specific transcription factors correlate with ChIP-seq enrichment and can accurately identify functional versus nonfunctional transcription factor motifs. We also find that footprints reveal a unique evolutionary conservation pattern that differentiates functional footprinted bases from surrounding DNA. Finally, detailed analysis of CTCF footprints suggests multiple modes of binding and a novel DNA binding motif upstream of the primary binding site.

EBP1 is a novel E2F target gene regulated by transforming growth factor-β

Regulation of gene expression requires transcription factor binding to specific DNA elements, and a large body of work has focused on the identification of such sequences. However, it is becoming increasingly clear that eukaryotic transcription factors can exhibit widespread, nonfunctional binding to genomic DNA sites. Conversely, some of these proteins, such as E2F, can also modulate gene expression by binding to non-consensus elements. E2F comprises a family of transcription factors that play key roles in a wide variety of cellular functions, including survival, differentiation, activation during tissue regeneration, metabolism, and proliferation. E2F factors bind to the Erb3-binding protein 1 (EBP1) promoter in live cells. We now show that E2F binding to the EBP1 promoter occurs through two tandem DNA elements that do not conform to typical consensus E2F motifs. Exogenously expressed E2F1 activates EBP1 reporters lacking one, but not both sites, suggesting a degree of redundancy under certain conditions. E2F1 increases the levels of endogenous EBP1 mRNA in breast carcinoma and other transformed cell lines. In contrast, in non-transformed primary epidermal keratinocytes, E2F, together with the retinoblastoma family of proteins, appears to be involved in decreasing EBP1 mRNA abundance in response to growth inhibition by transforming growth factor-β1. Thus, E2F is likely a central coordinator of multiple responses that culminate in regulation of EBP1 gene expression, and which may vary depending on cell type and context.

TcoF-DB: dragon database for human transcription co-factors and transcription factor interacting proteins

The initiation and regulation of transcription in eukaryotes is complex and involves a large number of transcription factors (TFs), which are known to bind to the regulatory regions of eukaryotic DNA. Apart from TF–DNA binding, protein–protein interaction involving TFs is an essential component of the machinery facilitating transcriptional regulation. Proteins that interact with TFs in the context of transcription regulation but do not bind to the DNA themselves, we consider transcription co-factors (TcoFs). The influence of TcoFs on transcriptional regulation and initiation, although indirect, has been shown to be significant with the functionality of TFs strongly influenced by the presence of TcoFs. While the role of TFs and their interaction with regulatory DNA regions has been well-studied, the association between TFs and TcoFs has so far been given less attention. Here, we present a resource that is comprised of a collection of human TFs and the TcoFs with which they interact. Other proteins that have a proven interaction with a TF, but are not considered TcoFs are also included. Our database contains 157 high-confidence TcoFs and additionally 379 hypothetical TcoFs. These have been identified and classified according to the type of available evidence for their involvement in transcriptional regulation and their presence in the cell nucleus. We have divided TcoFs into four groups, one of which contains high-confidence TcoFs and three others contain TcoFs which are hypothetical to different extents. We have developed the Dragon Database for Human Transcription Co-Factors and Transcription Factor Interacting Proteins (TcoF-DB). A web-based interface for this resource can be freely accessed at http://cbrc.kaust.edu.sa/tcof/ and http://apps.sanbi.ac.za/tcof/.

Cux1 transcription factor is induced in inflammatory bowel disease and protects against experimental colitis

BACKGROUND

Cux1 is a ubiquitous transcriptional factor that has been associated with cell proliferation, migration, invasion, and differentiation. Cux1 is an effector of the transforming growth factor beta (TGFβ) pathway, PAR(2) receptor signaling, and cellular migration, mechanisms intimately related to inflammatory bowel diseases (IBD).

METHODS

CD1 mice treated with dextran sulfate sodium (DSS) in drinking water and cultured intestinal epithelial cells were used to determine Cux1 expression under inflammatory conditions. A commercial cDNA library was used to monitor CUX1 expression in IBD patients. The Cux1(ΔHD/ΔHD) hypomorphic mouse model (Cux1ΔHD) treated with DSS in drinking water was used and the disease severity assessed.

RESULTS

Cux1 expression increased in cultured intestinal epithelial cells stimulated with tumor necrosis factor alpha (TNFα), in the mouse intestinal epithelium during experimental colitis and in human IBD patient samples. DSS-induced colitis in Cux1ΔHD mice was more severe according to clinical observations such as weight loss, colon length, and rectal bleeding. Histological observations confirmed an increase of IBD-related morphological changes including ulceration and mucosal infiltration of leukocytes in Cux1ΔHD mice. An increased number of pSer(276)-RelA-positive cells and higher expression levels of proinflammatory cytokines were also measured in the colon of Cux1ΔHD diseased animals. Elevated levels of Cxcl1 were measured before and after DSS-treatment and a greater neutrophilic infiltration was quantified in DSS-treated Cux1ΔHD mice. Finally, mucosal healing was significantly impaired in Cux1ΔHD mice during recovery from DSS treatment.

CONCLUSIONS

CUX1 is increased in response to inflammatory stress and its nuclear expression is crucial to protect against DSS-induced colitis and subsequent mucosal healing.

Genome-wide analysis of two-component systems and prediction of stress-responsive two-component system members in soybean

In plants, the two-component systems (TCSs) play important roles in regulating diverse biological processes, including responses to environmental stress stimuli. Within the soybean genome, the TCSs consist of at least 21 histidine kinases, 13 authentic and pseudo-phosphotransfers and 18 type-A, 15 type-B, 3 type-C and 11 pseudo-response regulator proteins. Structural and phylogenetic analyses of soybean TCS members with their Arabidopsis and rice counterparts revealed similar architecture of their TCSs. We identified a large number of closely homologous soybean TCS genes, which likely resulted from genome duplication. Additionally, we analysed tissue-specific expression profiles of those TCS genes, whose data are available from public resources. To predict the putative regulatory functions of soybean TCS members, with special emphasis on stress-responsive functions, we performed comparative analyses from all the TCS members of soybean, Arabidopsis and rice and coupled these data with annotations of known abiotic stress-responsive cis-elements in the promoter region of each soybean TCS gene. Our study provides insights into the architecture and a solid foundation for further functional characterization of soybean TCS elements. In addition, we provide a new resource for studying the conservation and divergence among the TCSs within plant species and/or between plants and other organisms.

Human microRNA oncogenes and tumor suppressors show significantly different biological patterns: from functions to targets

MicroRNAs (miRNAs) are small noncoding RNAs which play essential roles in many important biological processes. Therefore, their dysfunction is associated with a variety of human diseases, including cancer. Increasing evidence shows that miRNAs can act as oncogenes or tumor suppressors, and although there is great interest in research into these cancer-associated miRNAs, little is known about them. In this study, we performed a comprehensive analysis of putative human miRNA oncogenes and tumor suppressors. We found that miRNA oncogenes and tumor suppressors clearly show different patterns in function, evolutionary rate, expression, chromosome distribution, molecule size, free energy, transcription factors, and targets. For example, miRNA oncogenes are located mainly in the amplified regions in human cancers, whereas miRNA tumor suppressors are located mainly in the deleted regions. miRNA oncogenes tend to cleave target mRNAs more frequently than miRNA tumor suppressors. These results indicate that these two types of cancer-associated miRNAs play different roles in cancer formation and development. Moreover, the patterns identified here can discriminate novel miRNA oncogenes and tumor suppressors with a high degree of accuracy. This study represents the first large-scale bioinformatic analysis of human miRNA oncogenes and tumor suppressors. Our findings provide help for not only understanding of miRNAs in cancer but also for the specific identification of novel miRNAs as miRNA oncogenes and tumor suppressors. In addition, the data presented in this study will be valuable for the study of both miRNAs and cancer.

Determining the specificity of protein-DNA interactions

Proteins, such as many transcription factors, that bind to specific DNA sequences are essential for the proper regulation of gene expression. Identifying the specific sequences that each factor binds can help to elucidate regulatory networks within cells and how genetic variation can cause disruption of normal gene expression, which is often associated with disease. Traditional methods for determining the specificity of DNA-binding proteins are slow and laborious, but several new high-throughput methods can provide comprehensive binding information much more rapidly. Combined with in vivo determinations of transcription factor binding locations, this information provides more detailed views of the regulatory circuitry of cells and the effects of variation on gene expression.

Functional modularity of nuclear hormone receptors in a Caenorhabditis elegans metabolic gene regulatory network

We present the first gene regulatory network (GRN) that pertains to post-developmental gene expression. Specifically, we mapped a transcription regulatory network of Caenorhabditis elegans metabolic gene promoters using gene-centered yeast one-hybrid assays. We found that the metabolic GRN is enriched for nuclear hormone receptors (NHRs) compared with other gene-centered regulatory networks, and that these NHRs organize into functional network modules.

The NHR family has greatly expanded in nematodes; C. elegans has 284 NHRs, whereas humans have only 48. We show that the NHRs in the metabolic GRN have metabolic phenotypes, suggesting that they do not simply function redundantly.

The mediator subunit MDT-15 preferentially interacts with NHRs that occur in the metabolic GRN.

We describe an NHR circuit that responds to nutrient availability and propose a model for the evolution and organization of NHRs in C. elegans metabolic regulatory networks.

Physical and/or regulatory interactions between transcription factors (TFs) and their target genes are essential to establish body plans of multicellular organisms during development, and these interactions have been studied extensively in the context of GRNs. The precise control of differential gene expression is also of critical importance to maintain physiological homeostasis, and many metabolic disorders such as obesity and diabetes coincide with substantial changes in gene expression. Much work has focused on the GRNs that control metazoan development; however, the design principles and organization of the GRNs that control systems physiology remain largely unexplored.

In this study, we present the first gene-centered GRN that includes ∼70 genes involved in C. elegans metabolism and physiology, 100 TFs and more than 500 protein–DNA interactions between them. The resulting metabolic GRN is enriched for NHRs, compared with other gene-centered regulatory networks. NHRs are well-known regulators of lipid meta-qj;bolism in mammals. The transcriptional activity of NHRs can be modified by diffusible ligands, which allows these TFs to function as molecular sensors and rapidly alter the expression of their target genes. Interestingly, NHRs comprise the largest family of TFs in nematodes; the C. elegans genome encodes 284 NHRs, most of which are uncharacterized. Furthermore, their organization in GRNs has not yet been investigated. In our study, we show that the C. elegans NHRs that we retrieved in the metabolic GRN organize into network modules, and that most of these NHRs function to maintain lipid homeostasis in the nematode. Interestingly, network modularity has been proposed to facilitate rapid and robust changes in gene expression. Our results suggest that the C. elegans metabolic GRN may have evolved by combining NHR family expansion with the specific modular wiring of NHRs to enable the rapid adaptation of the animal to different environmental cues.

NHRs can interact with transcriptional cofactors such as chromatin remodeling complexes and Mediator components. For instance, the C. elegans Mediator subunit, MDT-15, can interact with NHR-49 to regulate the expression of its target genes. To find all the TFs that MDT-15 can interact with, we performed systematic yeast two-hybrid assays with MDT-15 versus 755 full-length TFs. We found that MDT-15 preferentially associates with NHRs, and specifically with those NHRs that confer a metabolic phenotype and that occur in the metabolic GRN. This illustrates the central role of MDT-15 in the regulation of metabolic gene expression.

Using a variety of genetic and biochemical approaches, we characterized NHR-86 in more detail. NHR-86 participates in one of the two NHR modules, and has a high-flux capacity; that is it has both a high incoming and a high outgoing degree. We obtained an nhr-86 mutant and generated an NHR-86 antibody, and showed that NHR-86 functions as an auto-repressor in vivo and that nhr-86 mutant animals store abnormally high levels of body fat.

Finally, we discovered a novel NHR circuit that responds to nutrient availability. In this circuit NHR-45 regulates the activity of nhr-178 promoter in two distinct physiologically important tissues: the intestine and the hypodermis. Both of these NHRs are required to maintain lipid homeostasis in C. elegans. The expression of nhr-178 is responsive to the nutritional status of the animal, which switches between ON and OFF states in the hypodermis. We found that NHR-45 activity is necessary to control this switch in the hypodermis. Interestingly, NHR-45 has opposite effects on the activity of the nhr-178 promoter in these tissues: NHR-45 activates this promoter in the intestine, but represses it in the hypodermis.

Altogether our study leads to a model in which the expansion of the NHR family, TFs that have the capacity to act as fast molecular sensors, is combined with a modular network organization to enable rapid and robust responses to various environmental cues.

Gene regulatory networks (GRNs) provide insights into the mechanisms of differential gene expression at a systems level. GRNs that relate to metazoan development have been studied extensively. However, little is still known about the design principles, organization and functionality of GRNs that control physiological processes such as metabolism, homeostasis and responses to environmental cues. In this study, we report the first experimentally mapped metazoan GRN of Caenorhabditis elegans metabolic genes. This network is enriched for nuclear hormone receptors (NHRs). The NHR family has greatly expanded in nematodes: humans have 48 NHRs, but C. elegans has 284, most of which are uncharacterized. We find that the C. elegans metabolic GRN is highly modular and that two GRN modules predominantly consist of NHRs. Network modularity has been proposed to facilitate a rapid response to different cues. As NHRs are metabolic sensors that are poised to respond to ligands, this suggests that C. elegans GRNs evolved to enable rapid and adaptive responses to different cues by a concurrence of NHR family expansion and modular GRN wiring.

PPARG: Gene Expression Regulation and Next-Generation Sequencing for Unsolved Issues

Peroxisome proliferator-activated receptor gamma (PPARγ) is one of the most extensively studied ligand-inducible transcription factors (TFs), able to modulate its transcriptional activity through conformational changes. It is of particular interest because of its pleiotropic functions: it plays a crucial role in the expression of key genes involved in adipogenesis, lipid and glucid metabolism, atherosclerosis, inflammation, and cancer. Its protein isoforms, the wide number of PPARγ target genes, ligands, and coregulators contribute to determine the complexity of its function. In addition, the presence of genetic variants is likely to affect expression levels of target genes although the impact of PPARG gene variations on the expression of target genes is not fully understood. The introduction of massively parallel sequencing platforms-in the Next Generation Sequencing (NGS) era-has revolutionized the way of investigating the genetic causes of inherited diseases. In this context, DNA-Seq for identifying-within both coding and regulatory regions of PPARG gene-novel nucleotide variations and haplotypes associated to human diseases, ChIP-Seq for defining a PPARγ binding map, and RNA-Seq for unraveling the wide and intricate gene pathways regulated by PPARG, represent incredible steps toward the understanding of PPARγ in health and disease.