genomepy: genes and genomes at your fingertips

Analyzing a functional genomics experiment, such as ATAC-, ChIP- or RNA-sequencing, requires reference data including a genome assembly and gene annotation. These resources can generally be retrieved from different organizations and in different versions. Most bioinformatic workflows require the user to supply this genomic data manually, which can be a tedious and error-prone process. Here we present genomepy, which can search, download, and preprocess the right genomic data for your analysis. Genomepy can search genomic data on NCBI, Ensembl, UCSC and GENCODE, and compare available gene annotations to enable an informed decision. The selected genome and gene annotation can be downloaded and preprocessed with sensible, yet controllable, defaults. Additional supporting data can be automatically generated or downloaded, such as aligner indexes, genome metadata and blacklists. Genomepy is freely available at https://github.com/vanheeringen-lab/genomepy under the MIT license and can be installed through pip or bioconda

GimmeMotifs: a de novo motif prediction pipeline for ChIP-sequencing experiments

Summary: Accurate prediction of transcription factor binding motifs that are enriched in a collection of sequences remains a computational challenge. Here we report on GimmeMotifs, a pipeline that incorporates an ensemble of computational tools to predict motifs de novo from ChIP-sequencing (ChIP-seq) data. Similar redundant motifs are compared using the weighted information content (WIC) similarity score and clustered using an iterative procedure. A comprehensive output report is generated with several different evaluation metrics to compare and evaluate the results. Benchmarks show that the method performs well on human and mouse ChIP-seq datasets. GimmeMotifs consists of a suite of command-line scripts that can be easily implemented in a ChIP-seq analysis pipeline

Role of p53 Serine 46 in p53 Target Gene Regulation

The tumor suppressor p53 plays a crucial role in cellular growth control inducing a plethora of different response pathways. The molecular mechanisms that discriminate between the distinct p53-responses have remained largely elusive. Here, we have analyzed the p53-regulated pathways induced by Actinomycin D and Etoposide treatment resulting in more growth arrested versus apoptotic cells respectively. We found that the genome-wide p53 DNA-binding patterns are almost identical upon both treatments notwithstanding transcriptional differences that we observed in global transcriptome analysis. To assess the role of post-translational modifications in target gene choice and activation we investigated the genome-wide level of phosphorylation of Serine 46 of p53 bound to DNA (p53-pS46) and of Serine 15 (p53-pS15). Interestingly, the extent of S46 phosphorylation of p53 bound to DNA is considerably higher in cells directed towards apoptosis while the degree of phosphorylation at S15 remains highly similar. Moreover, our data suggest that following different chemotherapeutical treatments, the amount of chromatin-associated p53 phosphorylated at S46 but not at pS15 is higher on certain apoptosis related target genes. Our data provide evidence that cell fate decisions are not made primarily on the level of general p53 DNA-binding and that post-translationally modified p53 can have distinct DNA-binding characteristics

ChIP-Chip Designs to Interrogate the Genome of Xenopus Embryos for Transcription Factor Binding and Epigenetic Regulation

Chromatin immunoprecipitation combined with genome tile path microarrays or deep sequencing can be used to study genome-wide epigenetic profiles and the transcription factor binding repertoire. Although well studied in a variety of cell lines, these genome-wide profiles have so far been little explored in vertebrate embryos. genome. In particular, a whole-genome microarray design was used to identify active promoters by close proximity to histone H3 lysine 4 trimethylation. A second microarray design features these experimentally derived promoter regions in addition to currently annotated 5′ ends of genes. These regions truly represent promoters as shown by binding of TBP, a key transcription initiation factor. embryos

Dynamics of enhancer chromatin signatures mark the transition from pluripotency to cell specification during embryogenesis

This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date; after six months, it is available under a Creative Commons License.-- et al.The generation of distinctive cell types that form different tissues and organs requires precise, temporal and spatial control of gene expression. This depends on specific cis-regulatory elements distributed in the noncoding DNA surrounding their target genes. Studies performed on mammalian embryonic stem cells and Drosophila embryos suggest that active enhancers form part of a defined chromatin landscape marked by histone H3 lysine 4 mono-methylation (H3K4me1) and histone H3 lysine 27 acetylation (H3K27ac). Nevertheless, little is known about the dynamics and the potential roles of these marks during vertebrate embryogenesis. Here, we provide genomic maps of H3K4me1/me3 and H3K27ac at four developmental time-points of zebrafish embryogenesis and analyze embryonic enhancer activity. We find that (1) changes in H3K27ac enrichment at enhancers accompany the shift from pluripotency to tissue-specific gene expression, (2) in early embryos, the peaks of H3K27ac enrichment are bound by pluripotent factors such as Nanog, and (3) the degree of evolutionary conservation is higher for enhancers that become marked by H3K27ac at the end of gastrulation, suggesting their implication in the establishment of the most conserved (phylotypic) transcriptome that is known to occur later at the pharyngula stage.We thank the Spanish and Andalusian Governments for grants (BFU2010-14839, CSD2007-00008, and Proyecto de Excelencia CVI-3488) for funding this study.Peer reviewe

WDR5, BRCA1, and BARD1 Co-regulate the DNA Damage Response and Modulate the Mesenchymal-to-Epithelial Transition during Early Reprogramming.

Differentiated cells are epigenetically stable, but can be reprogrammed to pluripotency by expression of the OSKM transcription factors. Despite significant effort, relatively little is known about the cellular requirements for reprogramming and how they affect the properties of induced pluripotent stem cells. We have performed high-content screening with small interfering RNAs targeting 300 chromatin-associated factors and extracted colony-level quantitative features. This revealed five morphological phenotypes in early reprogramming, including one displaying large round colonies exhibiting an early block of reprogramming. Using RNA sequencing, we identified transcriptional changes associated with these phenotypes. Furthermore, double knockdown epistasis experiments revealed that BRCA1, BARD1, and WDR5 functionally interact and are required for the DNA damage response. In addition, the mesenchymal-to-epithelial transition is affected in Brca1, Bard1, and Wdr5 knockdowns. Our data provide a resource of chromatin-associated factors in early reprogramming and underline colony morphology as an important high-dimensional readout for reprogramming quality.V.B. and C.B. are funded by the Stand Up to Cancer campaign for Cancer Research UK, and Cancer Research UK Program Foundation Award to C.B. (C37275/1A20146). K.M. was supported by an NWO-VIDI grant (864.12.010)

Identification of the regulatory circuit governing corneal epithelial fate determination and disease.

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity

The novel p53 target gene IRF2BP2 participates in cell survival during the p53 stress response

The tumor suppressor p53 contributes to the cellular fate after genotoxic insults, mainly through the regulation of target genes, thereby allowing e.g. repair mechanisms resulting in cell survival or inducing apoptosis. Unresolved so far is the issue, which exact mechanisms lead to one or the other cellular outcome. Here, we describe the interferon regulatory factor-2-binding protein-2 (IRF2BP2) as a new direct target gene of p53, influencing the p53-mediated cellular decision. We show that upregulation of IRF2BP2 after treatment with actinomycin D (Act.D) is dependent on functional p53 in different cell lines. This occurs in parallel with the down-regulation of the interacting partner of IRF2BP2, the interferon regulatory factor-2 (IRF2), which is known to positively influence cell growth. Analyzing the molecular functions of IRF2BP2, it appears to be able to impede on the p53-mediated transactivation of the p21- and the Bax-gene. We show here that overexpressed IRF2BP2 has an impact on the cellular stress response after Act.D treatment and that it diminishes the induction of apoptosis after doxorubicin treatment. Furthermore, the knockdown of IRF2BP2 leads to an upregulation of p21 and faster induction of apoptosis after doxorubicin as well as Act.D treatment

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