Protocol to target a promoter region in human embryonic kidney cells using the CRISPR-dCas9 system for single-locus proteomics

The unbiased identification of less-abundant transcription factors, which direct the expression of a target gene, is technically challenging. Here, we present a protocol to analyze the locus-specific chromatin-regulating proteome using in situ capture of chromatin interactions by an inactive Cas9 (dCas9). We describe steps for designing guide RNAs and transfection, followed by precipitation of chromatin and associated proteins. In the last step, we describe the elution of DNA and proteins for PCR and mass spectrometric analysis, respectively. For complete details on the use and execution of this protocol, please refer to Alkhayer et al.1.

Combined transcriptomic and ChIPseq analyses of the Bordetella pertussis RisA regulon

The regulation of Bordetella pertussis virulence is mediated by the two-component system BvgA/S, which activates the transcription of virulence-activated genes (vags). In the avirulent phase, the vags are not expressed, but instead, virulence-repressed genes (vrgs) are expressed, under the control of another two-component system, RisA/K. Here, we combined transcriptomic and chromatin immunoprecipitation sequencing (ChIPseq) data to examine the RisA/K regulon. We performed RNAseq analyses of RisA-deficient and RisA-phosphoablative B. pertussis mutants cultivated in virulent and avirulent conditions. We confirmed that the expression of most vrgs is regulated by phosphorylated RisA. However, the expression of some, including those involved in flagellum biosynthesis and chemotaxis, requires RisA independently of phosphorylation. Many RisA-regulated genes encode proteins with regulatory functions, suggesting multiple RisA regulation cascades. By ChIPseq analyses, we identified 430 RisA-binding sites, 208 within promoter regions, 201 within open reading frames, and 21 in non-coding regions. RisA binding was demonstrated in the promoter regions of most vrgs and, surprisingly, of some vags, as well as for other genes not identified as vags or vrgs. Unexpectedly, many genes, including some vags, like prn, brpL, bipA, and cyaA, contain a BvgA-binding site and a RisA-binding site, which increases the complexity of the RisAK/BvgAS network in B. pertussis virulence regulation.IMPORTANCEThe expression of virulence-activated genes (vags) of Bordetella pertussis, the etiological agent of whooping cough, is under the transcriptional control of the two-component system BvgA/S, which allows the bacterium to switch between virulent and avirulent phases. In addition, the more recently identified two-component system RisA/K is required for the expression of B. pertussis genes, collectively named vrgs, that are repressed during the virulent phase but activated during the avirulent phase. We have characterized the RisA/K regulon by combined transcriptomic and chromatin immunoprecipitation sequencing analyses. We identified more than 400 RisA-binding sites. Many of them are localized in promoter regions, especially vrgs, but some were found within open reading frames and in non-coding regions. Surprisingly, RisA-binding sites were also found in promoter regions of some vags, illustrating the previously underappreciated complexity of virulence regulation in B. pertussis.

Contrasting effects of whole-body and hepatocyte-specific deletion of the RNA polymerase III repressor Maf1 in the mouse

MAF1 is a nutrient-sensitive, TORC1-regulated repressor of RNA polymerase III (Pol III). MAF1 downregulation leads to increased lipogenesis in Drosophila melanogaster, Caenorhabditis elegans, and mice. However, Maf1 -/- mice are lean as increased lipogenesis is counterbalanced by futile pre-tRNA synthesis and degradation, resulting in increased energy expenditure. We compared Chow-fed Maf1 -/- mice with Chow- or High Fat (HF)-fed Maf1 hep-/- mice that lack MAF1 specifically in hepatocytes. Unlike Maf1 -/- mice, Maf1 hep-/- mice become heavier and fattier than control mice with old age and much earlier under a HF diet. Liver ChIPseq, RNAseq and proteomics analyses indicate increased Pol III occupancy at Pol III genes, very few differences in mRNA accumulation, and protein accumulation changes consistent with increased lipogenesis. Futile pre-tRNA synthesis and degradation in the liver, as likely occurs in Maf1 hep-/- mice, thus seems insufficient to counteract increased lipogenesis. Indeed, RNAseq and metabolite profiling indicate that liver phenotypes of Maf1 -/- mice are strongly influenced by systemic inter-organ communication. Among common changes in the three phenotypically distinct cohorts, Angiogenin downregulation is likely linked to increased Pol III occupancy of tRNA genes in the Angiogenin promoter.

Calculating the spatial density of regulatory chromatin interactions using multi-modal datasets from the same cell line

Here, we present a protocol for calculating the spatial density of regulatory chromatin interactions (SD-RCI) using Hi-C, ATAC-seq, and ChIP-seq datasets from the same cell line. We describe steps for selecting and preprocessing datasets, training and predicting a model to obtain regulatory chromatin interactions, and evaluating model performance. We then detail calculation of SD-RCI and visualization of the correlation between SD-RCI and gene expression. This protocol is applicable to Hi-C, ATAC-seq, and ChIP-seq data from the human cell line. For complete details on the use and execution of this protocol, please refer to Gong et al. (2023).¹.

Purification of immune-active macrophage super enhancers by chemical cross-linked chromatin immune precipitation

Isolation of extraordinarily long-length super-enhancers (SEs) using typical chromatin immune precipitation (ChIP) techniques can lead to DNA breakage due to uncontrolled cross-linking. We present a redefined ChIP technique for SE purification. After controlled paraformaldehyde-based cross-linking, glycine was used to quench the cross-linker followed by mild sonication. The sonication produced ideal fragment length of long-length SE chromatin. Presently, miR146a-5p SE of macrophages was pulled using BRD4 protein. Our protocol can reproducibly simplify the SE element isolation issues, in a quality-controlled manner. For complete details on the use and execution of this protocol, please refer to Das et al. (2021).1.

High-resolution mapping of mitotic DNA synthesis under conditions of replication stress in cultured cells

Cells experiencing DNA replication stress enter mitosis with under-replicated DNA, which activates a repair mechanism known as mitotic DNA synthesis (MiDAS). Here we describe a protocol to identify at genome wide and at high resolution the genomic sites where MiDAS occurs in cells exposed to aphidicolin. We use EdU incorporation to label nascent DNA in mitotic cells, followed by isolation of the EdU-labeled DNA and next-generation sequencing. For complete details on the use and execution of this protocol, please refer to Groelly et al. (2022)1 and Macheret et al. (2020).2.

Chromatin-bound protein colocalization analysis using bedGraph2Cluster and PanChIP

Computational pipelines for chromatin immunoprecipitation sequencing analysis can neglect colocalization events that occur in a mere subset of the genome. Here, we detail a streamlined approach for assessing colocalization of chromatin-bound proteins using the bedGraph2Cluster and PanChIP algorithms. Using histone modifications as an example, bedGraph2Cluster performs clustering analysis on chromatin binding patterns of target proteins. PanChIP then compares these clusters with a reference library of chromatin binding patterns and measures the overlap in peaks, capturing the heterogeneity in chromatin binding and colocalization patterns. For complete details on the use and execution of this protocol, please refer to Sanidas et al. (2022).1.

bulkAnalyseR: an accessible, interactive pipeline for analysing and sharing bulk multi-modal sequencing data

Bulk sequencing experiments (single- and multi-omics) are essential for exploring wide-ranging biological questions. To facilitate interactive, exploratory tasks, coupled with the sharing of easily accessible information, we present bulkAnalyseR, a package integrating state-of-the-art approaches using an expression matrix as the starting point (pre-processing functions are available as part of the package). Static summary images are replaced with interactive panels illustrating quality-checking, differential expression analysis (with noise detection) and biological interpretation (enrichment analyses, identification of expression patterns, followed by inference and comparison of regulatory interactions). bulkAnalyseR can handle different modalities, facilitating robust integration and comparison of cis-, trans- and customised regulatory networks.

Machine learning reveals STAT motifs as predictors for GR-mediated gene repression

Glucocorticoids are potent immunosuppressive drugs, but long-term treatment leads to severe side-effects. While there is a commonly accepted model for GR-mediated gene activation, the mechanism behind repression remains elusive. Understanding the molecular action of the glucocorticoid receptor (GR) mediated gene repression is the first step towards developing novel therapies. We devised an approach that combines multiple epigenetic assays with 3D chromatin data to find sequence patterns predicting gene expression change. We systematically tested> 100 models to evaluate the best way to integrate the data types and found that GR-bound regions hold most of the information needed to predict the polarity of Dex-induced transcriptional changes. We confirmed NF-κB motif family members as predictors for gene repression and identified STAT motifs as additional negative predictors.

Computational workflow for integrative analyses of DNA replication timing, epigenomic, and transcriptomic data

Temporal profiling of DNA replication timing (RT) in combination with chromatin modifications, chromatin accessibility, and gene expression provides new insights into the causal relationships between chromatin and RT during cell cycle. Here, we describe a protocol for in-depth integrative computational analyses of Repli-seq, ATAC-seq, RNA-seq, and ChIP-seq or CUT&RUN data for multiple marks at various time points across cell cycle and changes in their interrelationships upon an experimental perturbation (e.g., knockdown or overexpression of a regulatory protein). For complete details on the use and execution of this protocol, please refer to Van Rechem et al. (2021).

The Sum of Two Halves May Be Different from the Whole-Effects of Splitting Sequencing Samples Across Lanes

The advances in high-throughput sequencing (HTS) have enabled the characterisation of biological processes at an unprecedented level of detail; most hypotheses in molecular biology rely on analyses of HTS data. However, achieving increased robustness and reproducibility of results remains a main challenge. Although variability in results may be introduced at various stages, e.g., alignment, summarisation or detection of differential expression, one source of variability was systematically omitted: the sequencing design, which propagates through analyses and may introduce an additional layer of technical variation. We illustrate qualitative and quantitative differences arising from splitting samples across lanes on bulk and single-cell sequencing. For bulk mRNAseq data, we focus on differential expression and enrichment analyses; for bulk ChIPseq data, we investigate the effect on peak calling and the peaks' properties. At the single-cell level, we concentrate on identifying cell subpopulations. We rely on markers used for assigning cell identities; both smartSeq and 10× data are presented. The observed reduction in the number of unique sequenced fragments limits the level of detail on which the different prediction approaches depend. Furthermore, the sequencing stochasticity adds in a weighting bias corroborated with variable sequencing depths and (yet unexplained) sequencing bias. Subsequently, we observe an overall reduction in sequencing complexity and a distortion in the biological signal across technologies, experimental contexts, organisms and tissues.

A cell-sorting-based protocol for cell cycle small-scale ChIP sequencing

Classic approaches to characterizing cell cycle leverage chemicals or altered nucleotide pools, which could impact chromatin states at specific phases of the cell cycle. Such approaches could induce metabolic alterations and/or DNA damage, which could reshape protein recruitment and histone modifications. In this protocol, we describe ways to fix and sort cells across the cell cycle based on their DNA content. We further detail immunoprecipitation and library preparation, allowing analysis of the epigenome by chromatin immunoprecipitation sequencing (ChIP-seq) for small numbers of cells.

For complete details on the use and execution of this protocol, please refer to Van Rechem et al. (2021).

•

Drugs are not used to arrest cells at different cell cycle phases

•

This protocol is adapted to small amounts of chromatin for ChIP-seq

•

This protocol can be adapted to separate additional cell cycle phases

•

This protocol can be adapted to any cycling cells and other histone modifications

Genome-Wide ChIPseq Analysis of AhR, COUP-TF, and HNF4 Enrichment in TCDD-Treated Mouse Liver

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor known for mediating the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds. Although the canonical mechanism of AhR activation involves heterodimerization with the aryl hydrocarbon receptor nuclear translocator, other transcriptional regulators that interact with AhR have been identified. Enrichment analysis of motifs in AhR-bound genomic regions implicated co-operation with COUP transcription factor (COUP-TF) and hepatocyte nuclear factor 4 (HNF4). The present study investigated AhR, HNF4α and COUP-TFII genomic binding and effects on gene expression associated with liver-specific function and cell differentiation in response to TCDD. Hepatic ChIPseq data from male C57BL/6 mice at 2 h after oral gavage with 30 µg/kg TCDD were integrated with bulk RNA-sequencing (RNAseq) time-course (2-72 h) and dose-response (0.01-30 µg/kg) datasets to assess putative AhR, HNF4α and COUP-TFII interactions associated with differential gene expression. Functional enrichment analysis of differentially expressed genes (DEGs) identified differential binding enrichment for AhR, COUP-TFII, and HNF4α to regions within liver-specific genes, suggesting intersections associated with the loss of liver-specific functions and hepatocyte differentiation. Analysis found that the repression of liver-specific, HNF4α target and hepatocyte differentiation genes, involved increased AhR and HNF4α binding with decreased COUP-TFII binding. Collectively, these results suggested TCDD-elicited loss of liver-specific functions and markers of hepatocyte differentiation involved interactions between AhR, COUP-TFII and HNF4α.

Chromatin Immunoprecipitation dataset of H3ac and H3K27me3 histone marks followed by DNA sequencing of Medicago truncatula embryos during control and heat stress conditions to decipher epigenetic regulation of desiccation tolerance acquisition

Desiccation tolerance (DT) is one of the most important processes that seeds need to acquire during seed maturation because it will ensure survival until seeds have favourable conditions for germinating. Moreover, in the current climate warming context, heat stress and its impact on seed maturation and quality has been increasingly studied by the scientific community. Even if the transcriptomic changes enrolled in DT acquisition and seed heat stress response are fairly known, its epigenetic control has not yet been investigated. Medicago truncatula is a model legume for studying seed molecular mechanisms, which is known to display a delay in the acquisition of seed maturation mechanisms under heat conditions, except for desiccation acquisition. Our aim was to evaluate the role of two histone marks during embryo development under control and heat stress conditions on seed maturation processes, including the DT acquisition. These histone marks have either repressive (H3K27me3) or inducible (H3ac) effects on gene transcription, respectively corresponding to markers of packed and accessible chromatins. We identified all genomic regions bound to the H3K27me3 histones at four developmental stages and to the H3ac histones at the two earlier developmental stages during seed maturation, from seed filling to mature dry seeds, collected under optimal and heat stress conditions in the model legume, Medicago truncatula (reference genotype A17). A list of genes and promoters potentially linked to these two histone marks is reported and could provide clues about the epigenetic regulation of seed maturation between control and heat stress conditions, including the desiccation tolerance acquisition.

Optimized protocols for ChIP-seq and deletion mutant construction in Pseudomonas syringae

Chromatin immunoprecipitation sequencing (ChIP-seq) is an efficient technique to identify the binding sites of transcription factors (TFs) in both eukaryotes and prokaryotes. However, its application in bacteria is very heterogeneous. In this protocol, we optimized the methods of ChIP-seq that can be widely applied to plant pathogens. We used homologous recombination to construct pK18mobsacB-Psph plasmid instead of restriction site ligation and replaced transconjugation with electroporation transformation in Pseudomonas syringae deletion mutant construction, which is more efficient and faster than previous methods.

For complete details on the use and execution of this protocol, please refer to Shao et al. (2021).

•

This approach can be used to construct TF-overexpressed Pseudomonas syringae strain

•

Protocols for ChIP-seq library construction of Pseudomonas syringae

•

A simplified procedure to construct a deletion mutant of Pseudomonas syringae

Chromatin and transcription factor profiling in rare stem cell populations using CUT&Tag

Muscle stem cells (MuSCs) are a rare stem cell population that provides myofibers with a remarkable capacity to regenerate after tissue injury. Here, we have adapted the Cleavage Under Target and Tagmentation technology to the mapping of the chromatin landscape and transcription factor binding in 50,000 activated MuSCs isolated from injured mouse hindlimb muscles. We have applied this same approach to human CD34⁺ hematopoietic stem and progenitor cells. This protocol could be adapted to any rare stem cell population.

For complete details on the use and execution of this protocol, please refer to Robinson et al. (2021).

•

Isolation of muscle stem cells from cardiotoxin-injured tissue

•

CUT&Tag-based mapping of epigenetic landscape and transcription factor enrichment

•

Adaptation of a synthetic spike-in DNA for effective normalization across samples

•

Pipeline for efficient analysis of CUT&Tag sequencing data hosted on Git-hub

Protocol for using heterologous spike-ins to normalize for technical variation in chromatin immunoprecipitation

Quantifying differential genome occupancy by chromatin immunoprecipitation (ChIP) remains challenging due to variation in chromatin fragmentation, immunoprecipitation efficiencies, and intertube variability. In this protocol, we add heterologous spike-ins from Drosophila chromatin as an internal control to the mice chromatin before immunoprecipitation to normalize for technical variation in ChIP-qPCR or ChIP-seq. The choice of spike-in depends on the evolutionary conservation of the protein of interest and the antibody used. For complete details on the use and execution of this protocol, please refer to Greulich et al. (2021).

ChIP-seq protocol for sperm cells and embryos to assess environmental impacts and epigenetic inheritance

In the field of epigenetic inheritance, delineating molecular mechanisms implicated in the transfer of paternal environmental conditions to descendants has been elusive. This protocol details how to track sperm chromatin intergenerationally. We describe mouse model design to probe chromatin states in single mouse sperm and techniques to assess pre-implantation embryo chromatin and gene expression. We place emphasis on how to obtain high-quality and quantifiable data sets in sperm and embryos, as well as highlight the limitations of working with low input.

For complete details on the use and execution of this protocol, please refer to Lismer et al. (2021).

•

Assessment of paternal environmental exposure effects on sperm chromatin by ChIP-seq

•

Tracking chromatin changes in sperm to embryos using ultra-low-input (ULI)-ChIP-seq

•

Identification of gene expression changes in embryos using low-input mRNA-seq

Genome-wide Interrogation of Protein-DNA Interactions in Mammalian Cells Using ChIPmentation

Mapping the genomic locations of chromatin-associated proteins, such as transcription factors and histone modifications, is key to understanding the mechanisms of transcriptional regulation. ChIPmentation offers a simple and robust way of investigating the genomic binding sites of a protein using relatively low-input material. Here, we present a detailed protocol for the key steps that lead to a successful ChIPmentation experiment, as well as a quick analysis pipeline to examine the data.

For complete details on the use and execution of this protocol, please refer to Schmidl et al. (2015). For example data produced by this protocol, please refer to Henriksson et al. (2019) and Zhang et al. (2019).

•

Preparation of a transposase-assisted ChIP-seq library

•

A simpler and cheaper approach than traditional ChIP-seq

•

A sensitive and reliable method for transcription factors and histone modifications

Combined RNAseq and ChIPseq Analyses of the BvgA Virulence Regulator of Bordetella pertussis

Bordetella pertussis, the etiological agent of whooping cough, remains a major global health problem. Despite the global usage of whole-cell vaccines since the 1950s and of acellular vaccines in the 1990s, it still is one of the most prevalent vaccine-preventable diseases in industrialized countries. Virulence of B. pertussis is controlled by BvgA/S, a two-component system responsible for upregulation of virulence-activated genes (vags) and downregulation of virulence-repressed genes (vrgs). By transcriptome sequencing (RNAseq) analyses, we identified more than 270 vags or vrgs, and chromatin immunoprecipitation sequencing (ChIPseq) analyses revealed 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Some vags, such as dnt and fhaL, do not contain a BvgA-binding site, suggesting indirect regulation. In contrast, several vrgs and some genes not identified by RNAseq analyses under laboratory conditions contain strong BvgA-binding sites, indicating previously unappreciated complexities of BvgA/S biology.

Bordetella pertussis regulates the production of its virulence factors by the two-component system BvgAS. In the virulence phase, BvgS phosphorylates BvgA, which then activates the transcription of virulence-activated genes (vags). In the avirulence phase, such as during growth in the presence of MgSO₄, BvgA is not phosphorylated and the vags are not expressed. Instead, a set of virulence-repressed genes (vrgs) is expressed. Here, we performed transcriptome sequencing (RNAseq) analyses on B. pertussis cultivated with or without MgSO₄ and on a BvgA-deficient Tohama I derivative. We observed that 146 genes were less expressed under modulating conditions or in the BvgA-deficient strain than under the nonmodulating condition, while 130 genes were more expressed. Some of the genes code for proteins with regulatory functions, suggesting a BvgA/S regulation cascade. To determine which genes are directly regulated by BvgA, we performed chromatin immunoprecipitation sequencing (ChIPseq) analyses. We identified 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Among the former, 32 are in BvgA-regulated putative promoter regions. Some vags, such as dnt and fhaL, contain no BvgA-binding site, suggesting indirect BvgA regulation. Unexpectedly, BvgA also bound to some vrg putative promoter regions. Together, these observations indicate an unrecognized complexity of BvgA/S biology.

IMPORTANCEBordetella pertussis, the etiological agent of whooping cough, remains a major global health problem. Despite the global usage of whole-cell vaccines since the 1950s and of acellular vaccines in the 1990s, it still is one of the most prevalent vaccine-preventable diseases in industrialized countries. Virulence of B. pertussis is controlled by BvgA/S, a two-component system responsible for upregulation of virulence-activated genes (vags) and downregulation of virulence-repressed genes (vrgs). By transcriptome sequencing (RNAseq) analyses, we identified more than 270 vags or vrgs, and chromatin immunoprecipitation sequencing (ChIPseq) analyses revealed 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Some vags, such as dnt and fhaL, do not contain a BvgA-binding site, suggesting indirect regulation. In contrast, several vrgs and some genes not identified by RNAseq analyses under laboratory conditions contain strong BvgA-binding sites, indicating previously unappreciated complexities of BvgA/S biology.

Context-Dependent Functions of E2F1: Cell Cycle, Cell Death, and DNA Damage Repair in Cortical Neurons

DNA damage has been reported to induce cell cycle-related neuronal death. This is significant as aberrant cell cycle re-entry of mature, post-mitotic neurons contributes to neurodegeneration. In this study, we investigate how DNA damage elicited by exposure to the topoisomerase I inhibitor camptothecin (CPT) leads to cycle-related death of cultured cortical neurons and examine the function of E2F1 in this process. CPT treatment induced cell cycle initiation of cortical neurons and elevated the expression of certain cell cycle components (e.g., cyclin D1, CDK4, E2F1) but failed to drive S phase entry or DNA synthesis. The arrest in the cell cycle is explained by the elevated expression of the CDK inhibitor p21^Cip1. Though its level was increased after CPT treatment, E2F1 did not drive treated neurons into the G1-S phase transition. E2F1 overexpression led to cell cycle activation and acute neuronal apoptosis without detectable entry of the neurons into S phase. ChIPseq analysis demonstrated that E2F1 predominantly occupies positions on or near the promoters of cell cycle related genes. Instead, in CPT-treated neurons, E2F1 preferentially regulated DNA repair related genes. Our study reveals that the functions of E2F1 in postmitotic neurons are context-dependent and offers novel insights into the role of E2F1 in DNA damage induced cycle-related neuronal death.

Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.

Genome-wide analysis of AR binding and comparison with transcript expression in primary human fetal prostate fibroblasts and cancer associated fibroblasts

The androgen receptor (AR) is a transcription factor, and key regulator of prostate development and cancer, which has discrete functions in stromal versus epithelial cells. AR expressed in mesenchyme is necessary and sufficient for prostate development while loss of stromal AR is predictive of prostate cancer progression. Many studies have characterized genome-wide binding of AR in prostate tumour cells but none have used primary mesenchyme or stroma. We applied ChIPseq to identify genomic AR binding sites in primary human fetal prostate fibroblasts and patient derived cancer associated fibroblasts, as well as the WPMY1 cell line overexpressing AR. We identified AR binding sites that were specific to fetal prostate fibroblasts (7534), cancer fibroblasts (629), WPMY1-AR (2561) as well as those common among all (783). Primary fibroblasts had a distinct AR binding profile versus prostate cancer cell lines and tissue, and showed a localisation to gene promoter binding sites 1 kb upstream of the transcriptional start site, as well as non-classical AR binding sequence motifs. We used RNAseq to define transcribed genes associated with AR binding sites and derived cistromes for embryonic and cancer fibroblasts as well as a cistrome common to both. These were compared to several in vivo ChIPseq and transcript expression datasets; which identified subsets of AR targets that were expressed in vivo and regulated by androgens. This analysis enabled us to deconvolute stromal AR targets active in stroma within tumour samples. Taken together, our data suggest that the AR shows significantly different genomic binding site locations in primary prostate fibroblasts compared to that observed in tumour cells. Validation of our AR binding site data with transcript expression in vitro and in vivo suggests that the AR target genes we have identified in primary fibroblasts may contribute to clinically significant and biologically important AR-regulated changes in prostate tissue.

The Transcription Factor NFATc1 Supports the Rejection of Heterotopic Heart Allografts

The immune suppressants cyclosporin A (CsA) and tacrolimus (FK506) are used worldwide in transplantation medicine to suppress graft rejection. Both CsA and FK506 inhibit the phosphatase calcineurin (CN) whose activity controls the immune receptor-mediated activation of lymphocytes. Downstream targets of CN in lymphocytes are the nuclear factors of activated T cells (NFATs). We show here that the activity of NFATc1, the most prominent NFAT factor in activated lymphocytes supports the acute rejection of heterotopic heart allografts. While ablation of NFATc1 in T cells prevented graft rejection, ectopic expression of inducible NFATc1/αA isoform led to rejection of heart allografts in recipient mice. Acceptance of transplanted hearts in mice bearing NFATc1-deficient T cells was accompanied by a reduction in number and cytotoxicity of graft infiltrating cells. In CD8⁺ T cells, NFATc1 controls numerous intracellular signaling pathways that lead to the metabolic switch to aerobic glycolysis and the expression of numerous lymphokines, chemokines, and their receptors, including Cxcr3 that supports the rejection of allogeneic heart transplants. These findings favors NFATc1 as a molecular target for the development of new strategies to control the cytotoxicity of T cells upon organ transplantation.

Regulatory integration of Hox factor activity with T-box factors in limb development

In tetrapods, Tbx4, Tbx5 and Hox cluster genes are crucial for forelimb and hindlimb development and mutations in these genes are responsible for congenital limb defects. The molecular basis of their integrated mechanisms of action in the context of limb development remains poorly understood. We studied Tbx4 and Hoxc10 owing to their overlapping loss-of-function phenotypes and colocalized expression in mouse hindlimb buds. We report an extensive overlap between Tbx4 and Hoxc10 genome occupancy and their putative target genes. Tbx4 and Hoxc10 interact directly with each other, have the ability to bind to a previously unrecognized T-box-Hox composite DNA motif and show synergistic activity when acting on reporter genes. Pitx1, the master regulator for hindlimb specification, also shows extensive genomic colocalization with Tbx4 and Hoxc10. Genome occupancy by Tbx4 in hindlimb buds is similar to Tbx5 occupancy in forelimbs. By contrast, another Hox factor, Hoxd13, also interacts with Tbx4/Tbx5 but antagonizes Tbx4/Tbx5-dependent transcriptional activity. Collectively, the modulation of Tbx-dependent activity by Hox factors acting on common DNA targets may integrate different developmental processes for the balanced formation of proportionate limbs.

Pitx1 directly modulates the core limb development program to implement hindlimb identity

Forelimbs (FLs) and hindlimbs (HLs) develop complex musculoskeletal structures that rely on the deployment of a conserved developmental program. Pitx1, a transcription factor gene with expression restricted to HL and absent from FL, plays an important role in generating HL features. The genomic mechanisms by which Pitx1 effects HL identity remain poorly understood. Here, we use expression profiling and analysis of direct Pitx1 targets to characterize the HL- and FL-restricted genetic programs in mouse and situate the Pitx1-dependent gene network within the context of limb-specific gene regulation. We show that Pitx1 is a crucial component of a narrow network of HL-restricted regulators, acting on a developmental program that is shared between FL and HL. Pitx1 targets sites that are in a similar chromatin state in FL and HL and controls expression of patterning genes as well as the chondrogenic program, consistent with impaired chondrogenesis in Pitx1^-/- HL. These findings support a model in which multifactorial actions of a limited number of HL regulators redirect the generic limb development program in order to generate the unique structural features of the limb.

Futures Challenges in Thyroid Hormone Signaling Research

The canonical pathway of thyroid hormone signaling involves its binding to nuclear receptors (TRs) acting directly on the transcription of a number of genes. Recent genome-wide studies revealed that chromatin occupancy by TR is not sufficient for transactivation of gene expression. Reciprocally, in some cases, DNA binding by TR may not be required for cellular response. This leaves many new questions to be addressed in future research.

bPeaks: a bioinformatics tool to detect transcription factor binding sites from ChIPseq data in yeasts and other organisms with small genomes

Peak calling is a critical step in ChIPseq data analysis. Choosing the correct algorithm as well as optimized parameters for a specific biological system is an essential task. In this article, we present an original peak-calling method (bPeaks) specifically designed to detect transcription factor (TF) binding sites in small eukaryotic genomes, such as in yeasts. As TF interactions with DNA are strong and generate high binding signals, bPeaks uses simple parameters to compare the sequences (reads) obtained from the immunoprecipitation (IP) with those from the control DNA (input). Because yeasts have small genomes (<20 Mb), our program has the advantage of using ChIPseq information at the single nucleotide level and can explore, in a reasonable computational time, results obtained with different sets of parameter values. Graphical outputs and text files are provided to rapidly assess the relevance of the detected peaks. Taking advantage of the simple promoter structure in yeasts, additional functions were implemented in bPeaks to automatically assign the peaks to promoter regions and retrieve peak coordinates on the DNA sequence for further predictions of regulatory motifs, enriched in the list of peaks. Applications of the bPeaks program to three different ChIPseq datasets from Saccharomyces cerevisiae, Candida albicans and Candida glabrata are presented. Each time, bPeaks allowed us to correctly predict the DNA binding sequence of the studied TF and provided relevant lists of peaks. The bioinformatics tool bPeaks is freely distributed to academic users. Supplementary data, together with detailed tutorials, are available online: http://bpeaks.gene-networks.net.

D-peaks: a visual tool to display ChIP-seq peaks along the genome

ChIP-sequencing is a method of choice to localize the positions of protein binding sites on DNA on a whole genomic scale. The deciphering of the sequencing data produced by this novel technique is challenging and it is achieved by their rigorous interpretation using dedicated tools and adapted visualization programs. Here, we present a bioinformatics tool (D-peaks) that adds several possibilities (including, user-friendliness, high-quality, relative position with respect to the genomic features) to the well-known visualization browsers or databases already existing. D-peaks is directly available through its web interface http://rsat.ulb.ac.be/dpeaks/ as well as a command line tool.