Genome wide screens in yeast to identify potential binding sites and target genes of DNA-binding proteins

A deterministic code for transcription factor-DNA recognition through computation of binding interfaces

The recognition code between transcription factor (TF) amino acids and DNA bases remains poorly understood. Here, the determinants of TF amino acid-DNA base binding selectivity were identified through the analysis of crystals of TF-DNA complexes. Selective, high-frequency interactions were identified for the vast majority of amino acid side chains (‘structural code’). DNA binding specificities were then independently assessed by meta-analysis of random-mutagenesis studies of Zn finger-target DNA sequences. Selective, high-frequency interactions were identified for the majority of mutagenized residues (‘mutagenesis code’). The structural code and the mutagenesis code were shown to match to a striking level of accuracy (P = 3.1 × 10−33), suggesting the identification of fundamental rules of TF binding to DNA bases. Additional insight was gained by showing a geometry-dictated choice among DNA-binding TF residues with overlapping specificity. These findings indicate the existence of a DNA recognition mode whereby the physical-chemical characteristics of the interacting residues play a deterministic role. The discovery of this DNA recognition code advances our knowledge on fundamental features of regulation of gene expression and is expected to pave the way for integration with higher-order complexity approaches.

Impact of Phenol-Enriched Virgin Olive Oils on the Postprandial Levels of Circulating microRNAs Related to Cardiovascular Disease

SCOPE

We investigate the postprandial modulation of cardiovascular-related microRNAs elicited by extra virgin olive oil (EVOOs) containing different levels of their own polyphenols.

METHODS AND RESULTS

It is randomized, postprandial, parallel, double-blind study. Twelve healthy participants consumed 30 mL of EVOO containing low (L-EVOO; 250 mg total phenols kg-1 of oil), medium (M-EVOO; 500 mg total phenols kg-1 of oil), and high (H-EVOO; 750 mg total phenols kg-1 of oil) enriched EVOOs. Postprandial plasma microRNAs levels are analyzed by real-time quantitative PCR. The results show that L-EVOO intake is associated with decreased let-7e-5p and miR-328a-3p levels and increased miR-17-5p and miR-20a-5p, concentrations. M-EVOO decreases plasma let-7e-5p and increases miR-17-5p, miR-20a-5p, and miR-192-5p levels. Finally, H-EVOO decreases let-7e-5p, miR-10a-5p, miR-21-5p, and miR-26b-5p levels.

CONCLUSION

During the postprandial state, the levels of let-7e-5p decrease with EVOO regardless of polyphenol content suggesting a general response to the fatty acid composition of EVOO or/and the presence of at least 250 mg polyphenol kg-1 olive oil. Moreover, the miR-17-92 cluster increases by low and medium polyphenol content suggesting a role in fatty acid metabolism and nutrient sensing. Thus, postprandial modulation of circulating microRNAs levels could be a potential mechanism for the cardiovascular benefits associated with EVOO intake.

Foxi3 transcription factor activity is mediated by a C-terminal transactivation domain and regulated by the Protein Phosphatase 2A (PP2A) complex

The Forkhead box (FOX) family consists of at least 19 subgroups of transcription factors which are characterized by the presence of an evolutionary conserved ‘forkhead’ or ‘winged-helix’ DNA-binding domain. Despite having a conserved core DNA binding domain, FOX proteins display remarkable functional diversity and are involved in many developmental and cell specific processes. In the present study, we focus on a poorly characterized member of the Forkhead family, Foxi3, which plays a critical role in the development of the inner ear and jaw. We show that Foxi3 contains at least two important functional domains, a nuclear localization sequence (NLS) and a C-terminal transactivation domain (TAD), and that it directly binds its targets in a sequence specific manner. We also show that the transcriptional activity of Foxi3 is regulated by phosphorylation, and that the activity of Foxi3 can be attenuated by its physical interaction with the protein phosphatase 2A (PP2A) complex.

A yeast selection system for the detection of proteasomal activation

The ubiquitin proteasome system (UPS) is a complex cellular machinery that catalyzes degradation of misfolded or damaged proteins and regulates turnover of native proteins in eukaryotic cells, thus playing a crucial role in maintaining protein homeostasis. The UPS has emerged as a drug target for a diverse range of diseases characterized by accumulation of misfolded or aggregated proteins. While enhancement of UPS activity is widely recognized as a promising strategy to prevent accumulation of aberrant, off-pathway protein conformations and ameliorate the phenotypes of a wide range of protein misfolding diseases, the molecular mechanisms underlying activation of proteasomal degradation are poorly characterized. We report the development of a yeast selection platform for genome-wide selection of UPS activators. We engineered the Saccharomyces cerevisiae selection marker orotidine-5'-phosphate decarboxylase (URA3) to function as a substrate of proteasomal degradation through fusion to UPS-sensitive tags. The resulting UPS-sensitive URA3 variant links UPS activity to cell growth. The yeast selection platform reported in this study will open the way to high-throughput, genome-wide studies aimed at identifying modulators of UPS function that might provide novel target for therapeutic applications.

The homeobox transcription factor Irxl1 negatively regulates MyoD expression and myoblast differentiation

Irxl1/Mkx (Iroquois homeobox-like 1/Mohawk) encodes a member of the TALE subfamily of homeodomain proteins. It is expressed in multiple mesoderm-derived tissues and has recently been shown to regulate tendon differentiation during mouse embryonic development. Previously we showed that knockdown of Irxl1 in zebrafish caused a deficit in neural crest cells which consequently resulted in deformation of craniofacial muscles and arch cartilages. Here, we further demonstrate that loss of Irxl1 function results in deformed somites with disordered muscle fibers and myotendinous junctions. Because expression of myoD is increased in the somites of Irxl1 knockdown morphants, we test whether Irxl1 negatively regulates myoD expression. When stable C2C12 myoblasts overexpressing Irxl1/Mkx were induced to differentiate, myotube formation was inhibited and protein levels of myoD and myosin heavy chain were decreased accordingly. A series of deletion constructs of myoD promoter fragments were tested by luciferase reporter assays, which identified a promoter fragment that is necessary and sufficient for Irxl1-mediated repression. Direct interaction of Irxl1 and myoD promoter was subsequently elucidated by yeast one-hybrid assays, electrophoretic mobility shift assays and chromatin immunoprecipitation analysis. Furthermore, mouse Mkx also binds to and represses myoD promoter. These results indicate that Irxl1/Mkx can repress myoD expression through direct binding to its promoter and may thus play a negative regulatory role in muscle differentiation.

Foxi transcription factors promote pharyngeal arch development by regulating formation of FGF signaling centers

The bones of the vertebrate face develop from transient embryonic branchial arches that are populated by cranial neural crest cells. We have characterized a mouse mutant for the Forkhead family transcription factor Foxi3, which is expressed in branchial ectoderm and endoderm. Foxi3 mutant mice are not viable and display severe branchial arch-derived facial skeleton defects, including absence of all but the most distal tip of the mandible and complete absence of the inner, middle and external ear structures. Although cranial neural crest cells of Foxi3 mutants are able to migrate, populate the branchial arches, and display some elements of correct proximo-distal patterning, they succumb to apoptosis from embryonic day 9.75 onwards. We show this cell death correlates with a delay in expression of Fgf8 in branchial arch ectoderm and a failure of neural crest cells in the arches to express FGF-responsive genes. Zebrafish foxi1 is also expressed in branchial arch ectoderm and endoderm, and morpholino knock-down of foxi1 also causes apoptosis of neural crest in the branchial arches. We show that heat shock induction of fgf3 in zebrafish arch tissue can rescue cell death in foxi1 morphants. Our results suggest that Foxi3 may play a role in the establishment of signaling centers in the branchial arches that are required for neural crest survival, patterning and the subsequent development of branchial arch derivatives.

Inferring functional transcription factor-gene binding pairs by integrating transcription factor binding data with transcription factor knockout data

Background

Chromatin immunoprecipitation (ChIP) experiments are now the most comprehensive experimental approaches for mapping the binding of transcription factors (TFs) to their target genes. However, ChIP data alone is insufficient for identifying functional binding target genes of TFs for two reasons. First, there is an inherent high false positive/negative rate in ChIP-chip or ChIP-seq experiments. Second, binding signals in the ChIP data do not necessarily imply functionality.

Methods

It is known that ChIP-chip data and TF knockout (TFKO) data reveal complementary information on gene regulation. While ChIP-chip data can provide TF-gene binding pairs, TFKO data can provide TF-gene regulation pairs. Therefore, we propose a novel network approach for identifying functional TF-gene binding pairs by integrating the ChIP-chip data with the TFKO data. In our method, a TF-gene binding pair from the ChIP-chip data is regarded to be functional if it also has high confident curated TFKO TF-gene regulatory relation or deduced hypostatic TF-gene regulatory relation.

Results and conclusions

We first validated our method on a gathered ground truth set. Then we applied our method to the ChIP-chip data to identify functional TF-gene binding pairs. The biological significance of our identified functional TF-gene binding pairs was shown by assessing their functional enrichment, the prevalence of protein-protein interaction, and expression coherence. Our results outperformed the results of three existing methods across all measures. And our identified functional targets of TFs also showed statistical significance over the randomly assigned TF-gene pairs. We also showed that our method is dataset independent and can apply to ChIP-seq data and the E. coli genome. Finally, we provided an example showing the biological applicability of our notion.

Beyond the binding site: in vivo identification of tbx2, smarca5 and wnt5b as molecular targets of CNBP during embryonic development

CNBP is a nucleic acid chaperone implicated in vertebrate craniofacial development, as well as in myotonic dystrophy type 2 (DM2) and sporadic inclusion body myositis (sIBM) human muscle diseases. CNBP is highly conserved among vertebrates and has been implicated in transcriptional regulation; however, its DNA binding sites and molecular targets remain elusive. The main goal of this work was to identify CNBP DNA binding sites that might reveal target genes involved in vertebrate embryonic development. To accomplish this, we used a recently described yeast one-hybrid assay to identify DNA sequences bound in vivo by CNBP. Bioinformatic analyses revealed that these sequences are G-enriched and show high frequency of putative G-quadruplex DNA secondary structure. Moreover, an in silico approach enabled us to establish the CNBP DNA-binding site and to predict CNBP putative targets based on gene ontology terms and synexpression with CNBP. The direct interaction between CNBP and candidate genes was proved by EMSA and ChIP assays. Besides, the role of CNBP upon the identified genes was validated in loss-of-function experiments in developing zebrafish. We successfully confirmed that CNBP up-regulates tbx2b and smarca5, and down-regulates wnt5b gene expression. The highly stringent strategy used in this work allowed us to identify new CNBP target genes functionally important in different contexts of vertebrate embryonic development. Furthermore, it represents a novel approach toward understanding the biological function and regulatory networks involving CNBP in the biology of vertebrates.

Yeast-based assay identifies novel Shh/Gli target genes in vertebrate development

Background

The increasing number of developmental events and molecular mechanisms associated with the Hedgehog (Hh) pathway from Drosophila to vertebrates, suggest that gene regulation is crucial for diverse cellular responses, including target genes not yet described. Although several high-throughput, genome-wide approaches have yielded information at the genomic, transcriptional and proteomic levels, the specificity of Gli binding sites related to direct target gene activation still remain elusive. This study aims to identify novel putative targets of Gli transcription factors through a protein-DNA binding assay using yeast, and validating a subset of targets both in-vitro and in-vivo. Testing in different Hh/Gli gain- and loss-of-function scenarios we here identified known (e.g., ptc1) and novel Hh-regulated genes in zebrafish embryos.

Results

The combined yeast-based screening and MEME/MAST analysis were able to predict Gli transcription factor binding sites, and position mapping of these sequences upstream or in the first intron of promoters served to identify new putative target genes of Gli regulation. These candidates were validated by qPCR in combination with either the pharmacological Hh/Gli antagonist cyc or the agonist pur in Hh-responsive C3H10T1/2 cells. We also used small-hairpin RNAs against Gli proteins to evaluate targets and confirm specific Gli regulation their expression. Taking advantage of mutants that have been identified affecting different components of the Hh/Gli signaling system in the zebrafish model, we further analyzed specific novel candidates. Studying Hh function with pharmacological inhibition or activation complemented these genetic loss-of-function approaches. We provide evidence that in zebrafish embryos, Hh signaling regulates sfrp2, neo1, and c-myc expression in-vivo.

Conclusion

A recently described yeast-based screening allowed us to identify new Hh/Gli target genes, functionally important in different contexts of vertebrate embryonic development.

Using a yeast inverse one-hybrid system to identify functional binding sites of transcription factors

Binding of transcription factors to promoters is a necessary step to initiate transcription. From an evolutionary standpoint, the regulatory proteins and their binding sites are considered to have molecularly coevolved. We developed an efficient yeast strategy, an "inverse one-hybrid system", to identify binding targets of transcription factors globally in a genome of interest. The technique consists of a yeast strain expressing a -transcription factor of interest mated to yeast containing a library of random genomic fragments cloned upstream of a reporter gene (URA3). Positive growth on media without uracil denotes a fragment being bound by the transcription factor, e.g., zebrafish FoxI1. The bound fragments in hundreds of positive clones are sequenced and retested for their binding activities using a colony PCR and sequencing strategy. The resulting tools allow for rapid and genomic-wide identification of transcriptional binding targets.

MicroRNAs in Development and Disease

MicroRNAs (miRNAs) are a class of posttranscriptional regulators that have recently introduced an additional level of intricacy to our understanding of gene regulation. There are currently over 10,000 miRNAs that have been identified in a range of species including metazoa, mycetozoa, viridiplantae, and viruses, of which 940, to date, are found in humans. It is estimated that more than 60% of human protein-coding genes harbor miRNA target sites in their 3′ untranslated region and, thus, are potentially regulated by these molecules in health and disease. This review will first briefly describe the discovery, structure, and mode of function of miRNAs in mammalian cells, before elaborating on their roles and significance during development and pathogenesis in the various mammalian organs, while attempting to reconcile their functions with our existing knowledge of their targets. Finally, we will summarize some of the advances made in utilizing miRNAs in therapeutics.

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.

Lysine120 interactions with p53 response elements can allosterically direct p53 organization

p53 can serve as a paradigm in studies aiming to figure out how allosteric perturbations in transcription factors (TFs) triggered by small changes in DNA response element (RE) sequences, can spell selectivity in co-factor recruitment. p53-REs are 20-base pair (bp) DNA segments specifying diverse functions. They may be located near the transcription start sites or thousands of bps away in the genome. Their number has been estimated to be in the thousands, and they all share a common motif. A key question is then how does the p53 protein recognize a particular p53-RE sequence among all the similar ones? Here, representative p53-REs regulating diverse functions including cell cycle arrest, DNA repair, and apoptosis were simulated in explicit solvent. Among the major interactions between p53 and its REs involving Lys120, Arg280 and Arg248, the bps interacting with Lys120 vary while the interacting partners of other residues are less so. We observe that each p53-RE quarter site sequence has a unique pattern of interactions with p53 Lys120. The allosteric, DNA sequence-induced conformational and dynamic changes of the altered Lys120 interactions are amplified by the perturbation of other p53-DNA interactions. The combined subtle RE sequence-specific allosteric effects propagate in the p53 and in the DNA. The resulting amplified allosteric effects far away are reflected in changes in the overall p53 organization and in the p53 surface topology and residue fluctuations which play key roles in selective co-factor recruitment. As such, these observations suggest how similar p53-RE sequences can spell the preferred co-factor binding, which is the key to the selective gene transactivation and consequently different functional effects.

Preferred drifting along the DNA major groove and cooperative anchoring of the p53 core domain: mechanisms and scenarios

While the importance of specific p53-DNA binding is broadly accepted, the recognition process is still not fully understood. Figuring out the initial tetrameric p53-DNA association and the swift and cooperative search for specific binding sites is crucial for understanding the transactivation mechanism and selectivity. To gain insight into the p53-DNA binding process, here we have carried out explicit solvent molecular dynamic (MD) simulations of several p53 core domain-DNA conformations with the p53 and the DNA separated by varying distances. p53 approached the DNA, bound non-specifically, and quickly drifted along the DNA surface to find the major groove, cooperatively anchoring in a way similar to the specific binding observed in the crystal structure. Electrostatics was the major driving force behind the p53 movement. Mechanistically, this is a cooperative process: key residues, particularly Lys120 and Arg280 acted as sensors; upon finding their hydrogen-bonding partners, they lock in, anchoring p53 into the major groove. Concomitantly, the DNA adopted a conformation that facilitated p53 easy access. The initial non-specific core domain-DNA contacts assist in shifting the DNA and the p53 substrates toward conformations "ready" for specific major groove binding, with subsequent optimization of the interactions. This work is an invited contribution for the special issue of the Journal of Molecular Recognition dedicated to Professor Martin Karplus.

Cooperativity dominates the genomic organization of p53-response elements: a mechanistic view

p53-response elements (p53-REs) are organized as two repeats of a palindromic DNA segment spaced by 0 to 20 base pairs (bp). Several experiments indicate that in the vast majority of the human p53-REs there are no spacers between the two repeats; those with spacers, particularly with sizes beyond two nucleotides, are rare. This raises the question of what it indicates about the factors determining the p53-RE genomic organization. Clearly, given the double helical DNA conformation, the orientation of two p53 core domain dimers with respect to each other will vary depending on the spacer size: a small spacer of 0 to 2 bps will lead to the closest p53 dimer-dimer orientation; a 10-bp spacer will locate the p53 dimers on the same DNA face but necessitate DNA looping; while a 5-bp spacer will position the p53 dimers on opposite DNA faces. Here, via conformational analysis we show that when there are 0-2 bp spacers, p53-DNA binding is cooperative; however, cooperativity is greatly diminished when there are spacers with sizes beyond 2 bp. Cooperative binding is broadly recognized to be crucial for biological processes, including transcriptional regulation. Our results clearly indicate that cooperativity of the p53-DNA association dominates the genomic organization of the p53-REs, raising questions of the structural organization and functional roles of p53-REs with larger spacers. We further propose that a dynamic landscape scenario of p53 and p53-REs can better explain the selectivity of the degenerate p53-REs. Our conclusions bear on the evolutionary preference of the p53-RE organization and as such, are expected to have broad implications to other multimeric transcription factor response element organization.

Functional analysis of transcription factors in Arabidopsis

Transcription factors (TFs) regulate the expression of genes at the transcriptional level. Modification of TF activity dynamically alters the transcriptome, which leads to metabolic and phenotypic changes. Thus, functional analysis of TFs using 'omics-based' methodologies is one of the most important areas of the post-genome era. In this mini-review, we present an overview of Arabidopsis TFs and introduce strategies for the functional analysis of plant TFs, which include both traditional and recently developed technologies. These strategies can be assigned to five categories: bioinformatic analysis; analysis of molecular function; expression analysis; phenotype analysis; and network analysis for the description of entire transcriptional regulatory networks.

Age and diet affect gene expression profile in canine skeletal muscle

We evaluated gene transcription in canine skeletal muscle (biceps femoris) using microarray analysis to identify effects of age and diet on gene expression. Twelve female beagles were used (six 1-year olds and six 12-year olds) and they were fed one of two experimental diets for 12 months. One diet contained primarily plant-based protein sources (PPB), whereas the second diet contained primarily animal-based protein sources (APB). Affymetrix GeneChip Canine Genome Arrays were used to hybridize extracted RNA. Age had the greatest effect on gene transcription (262 differentially expressed genes), whereas the effect of diet was relatively small (22 differentially expressed genes). Effects of age (regardless of diet) were most notable on genes related to metabolism, cell cycle and cell development, and transcription function. All these genes were predominantly down-regulated in geriatric dogs. Age-affected genes that were differentially expressed on only one of two diets were primarily noted in the PPB diet group (144/165 genes). Again, genes related to cell cycle (22/35) and metabolism (15/19) had predominantly decreased transcription in geriatric dogs, but 6/8 genes related to muscle development had increased expression. Effects of diet on muscle gene expression were mostly noted in geriatric dogs, but no consistent patterns in transcription were observed. The insight these data provide into gene expression profiles of canine skeletal muscle as affected by age, could serve as a foundation for future research pertaining to age-related muscle diseases.

Characterization of genome-wide p53-binding sites upon stress response

The tumor suppressor p53 is a sequence-specific transcription factor, which regulates the expression of target genes involved in different stress responses. To understand p53's essential transcriptional functions, unbiased analysis of its DNA-binding repertoire is pivotal. In a genome-wide tiling ChIP-on-chip approach, we have identified and characterized 1546 binding sites of p53 upon Actinomycin D treatment. Among those binding sites were known as well as novel p53 target sites, which included regulatory regions of potentially novel transcripts. Using this collection of genome-wide binding sites, a new high-confidence algorithm was developed, p53scan, to identify the p53 consensus-binding motif. Strikingly, this motif was present in the majority of all bound sequences with 83% of all binding sites containing the motif. In the surrounding sequences of the binding sites, several motifs for potential regulatory cobinders were identified. Finally, we show that the majority of the genome-wide p53 target sites can also be bound by overexpressed p63 and p73 in vivo, suggesting that they can possibly play an important role at p53 binding sites. This emphasizes the possible interplay of p53 and its family members in the context of target gene binding. Our study greatly expands the known, experimentally validated p53 binding site repertoire and serves as a valuable knowledgebase for future research.