Processing and Visualization of Protec-Seq Data for the Generation of Calibrated, Genome-Wide Double Double-Strand Break (dDSB) Maps.

This chapter describes the computational pipeline for the processing and visualization of Protec-Seq data, a method for purification and genome-wide mapping of double-stranded DNA protected by a specific protein at both ends. In the published case, the protein of choice was Saccharomyces cerevisiae Spo11, a conserved topoisomerase-like enzyme that makes meiotic double-strand breaks (DSBs) to initiate homologous recombination, ensuring proper segregation of homologous chromosomes and fertility. The isolated DNA molecules were thus termed double DSB (dDSB) fragments and were found to represent 34 to several hundred base-pair long segments that are generated by Spo11 and are enriched at DSB hotspots, which are sites of topological stress. In order to allow quantitative comparisons between dDSB profiles across experiments, we implemented calibrated chromatin immunoprecipitation sequencing (ChIP-Seq) using the meiosis-competent yeast species Saccharomyces kudriavzevii as calibration strain. Here, we provide a detailed description of the computational methods for processing, analyzing, and visualizing Protec-Seq data, comprising the download of the raw data, the calibrated genome-wide alignments, and the scripted creation of either arc plots or Hi-C-style heatmaps for the illustration of chromosomal regions of interest. The workflow is based on Linux shell scripts (including wrappers for publicly available, open-source software) as well as R scripts and is highly customizable through its modular structure.

Extracting Chromosome Structural Information as One-Dimensional Metrics and Integrating Them with Epigenomics.

Hi-C is a powerful method for obtaining genome-wide chromosomal structural information. The typical Hi-C analysis utilizes a two-dimensional (2D) contact matrix, which poses challenges for quantitative comparisons, visualizations, and integrations across multiple datasets. Here, we present a protocol for extracting one-dimensional (1D) features from chromosome structure data by HiC1Dmetrics. Leveraging these 1D features enables integrated analysis of Hi-C and epigenomic data.

HNF1ß bookmarking involves Topoisomerase 1 activation and DNA topology relaxation in mitotic chromatin.

HNF1ß (HNF1B) is a transcription factor frequently mutated in patients with developmental renal disease. It binds to mitotic chromatin and reactivates gene expression after mitosis, a phenomenon referred to as bookmarking. Using a crosslinking method that circumvents the artifacts of formaldehyde, we demonstrate that HNF1ß remains associated with chromatin in a sequence-specific way in both interphase and mitosis. We identify an HNF1ß-interacting protein, BTBD2, that enables the interaction and activation of Topoisomerase 1 (TOP1) exclusively during mitosis. Our study identifies a shared microhomology domain between HNF1ß and TOP1, where a mutation, found in "maturity onset diabetes of the young" patients, disrupts their interaction. Importantly, HNF1ß recruits TOP1 and induces DNA relaxation around HNF1ß mitotic chromatin sites, elucidating its crucial role in chromatin remodeling and gene reactivation after mitotic exit. These findings shed light on how HNF1ß reactivates target gene expression after mitosis, providing insights into its crucial role in maintenance of cellular identity.

Redirecting the pioneering function of FOXA1 with covalent small molecules.

Pioneer transcription factors (TFs) bind to and open closed chromatin, facilitating engagement by other regulatory factors involved in gene activation or repression. Chemical probes are lacking for pioneer TFs, which has hindered their mechanistic investigation in cells. Here, we report the chemical proteomic discovery of electrophilic compounds that stereoselectively and site-specifically bind the pioneer TF forkhead box protein A1 (FOXA1) at a cysteine (C258) within the forkhead DNA-binding domain. We show that these covalent ligands react with FOXA1 in a DNA-dependent manner and rapidly remodel its pioneer activity in prostate cancer cells reflected in redistribution of FOXA1 binding across the genome and directionally correlated changes in chromatin accessibility. Motif analysis supports a mechanism where the ligands relax the canonical DNA-binding preference of FOXA1 by strengthening interactions with suboptimal sequences in predicted proximity to C258. Our findings reveal a striking plasticity underpinning the pioneering function of FOXA1 that can be controlled by small molecules.

Histone-like transcription factor Hfl1p in Candida albicans harmonizes nuclear and mitochondrial genomic network in regulation of energy metabolism and filamentation development.

Candida albicans is an opportunistic fungal pathogen known for surviving in various nutrient-limited conditions within the host and causing infections. Our prior research revealed that Hfl1p, an archaeal histone-like or Hap5-like protein, is linked to mitochondrial ATP generation and yeast-hyphae morphogenesis. However, the specific roles of Hfl1p in these virulence behaviours, through its function in the CBF/NF-Y complex or as a DNA polymerase II subunit, remain unclear. This study explores Hfl1p's diverse functions in energy metabolism and morphogenesis. By combining proteomic analysis and phenotypic evaluations of the hfl1Δ/hfl1Δ mutant with ChIP data, we found that Hfl1p significantly impacts mitochondrial DNA-encoded CI subunits, the tricarboxylic acid (TCA) cycle, and morphogenetic pathways. This influence occurs either independently or alongside other transcription factors recognizing a conserved DNA motif (TAXXTAATTA). These findings emphasize Hfl1p's critical role in linking carbon metabolism and mitochondrial respiration to the yeast-to-filamentous form transition, enhancing our understanding of C. albicans' metabolic adaptability during morphological transition, an important pathogenic trait of this fungus. This could help identify therapeutic targets by disrupting the relationship between energy metabolism and cell morphology in C. albicans.

Deciphering genetic and nongenetic factors underlying tumour dormancy: insights from multiomics analysis of two syngeneic MRD models of melanoma and leukemia.

BACKGROUND: Tumour dormancy, a resistance mechanism employed by cancer cells, is a significant challenge in cancer treatment, contributing to minimal residual disease (MRD) and potential relapse. Despite its clinical importance, the mechanisms underlying tumour dormancy and MRD remain unclear. In this study, we employed two syngeneic murine models of myeloid leukemia and melanoma to investigate the genetic, epigenetic, transcriptomic and protein signatures associated with tumour dormancy. We used a multiomics approach to elucidate the molecular mechanisms driving MRD and identify potential therapeutic targets. RESULTS: We conducted an in-depth omics analysis encompassing whole-exome sequencing (WES), copy number variation (CNV) analysis, chromatin immunoprecipitation followed by sequencing (ChIP-seq), transcriptome and proteome investigations. WES analysis revealed a modest overlap of gene mutations between melanoma and leukemia dormancy models, with a significant number of mutated genes found exclusively in dormant cells. These exclusive genetic signatures suggest selective pressure during MRD, potentially conferring resistance to the microenvironment or therapies. CNV, histone marks and transcriptomic gene expression signatures combined with Gene Ontology (GO) enrichment analysis highlighted the potential functional roles of the mutated genes, providing insights into the pathways associated with MRD. In addition, we compared "murine MRD genes" profiles to the corresponding human disease through public datasets and highlighted common features according to disease progression. Proteomic analysis combined with multi-omics genetic investigations, revealed a dysregulated proteins signature in dormant cells with minimal genetic mechanism involvement. Pathway enrichment analysis revealed the metabolic, differentiation and cytoskeletal remodeling processes involved in MRD. Finally, we identified 11 common proteins differentially expressed in dormant cells from both pathologies. CONCLUSIONS: Our study underscores the complexity of tumour dormancy, implicating both genetic and nongenetic factors. By comparing genomic, transcriptomic, proteomic, and epigenomic datasets, our study provides a comprehensive understanding of the molecular landscape of minimal residual disease. These results provide a robust foundation for forthcoming investigations and offer potential avenues for the advancement of targeted MRD therapies in leukemia and melanoma patients, emphasizing the importance of considering both genetic and nongenetic factors in treatment strategies.

C2H2-zinc-finger transcription factors bind RNA and function in diverse post-transcriptional regulatory processes.

Cys2-His2 zinc-finger proteins (C2H2-ZNFs) constitute the largest class of DNA-binding transcription factors (TFs) yet remain largely uncharacterized. Although certain family members, e.g., GTF3A, have been shown to bind both DNA and RNA, the extent to which C2H2-ZNFs interact with-and regulate-RNA-associated processes is not known. Using UV crosslinking and immunoprecipitation (CLIP), we observe that 148 of 150 analyzed C2H2-ZNFs bind directly to RNA in human cells. By integrating CLIP sequencing (CLIP-seq) RNA-binding maps for 50 of these C2H2-ZNFs with data from chromatin immunoprecipitation sequencing (ChIP-seq), protein-protein interaction assays, and transcriptome profiling experiments, we observe that the RNA-binding profiles of C2H2-ZNFs are generally distinct from their DNA-binding preferences and that they regulate a variety of post-transcriptional processes, including pre-mRNA splicing, cleavage and polyadenylation, and m6A modification of mRNA. Our results thus define a substantially expanded repertoire of C2H2-ZNFs that bind RNA and provide an important resource for elucidating post-transcriptional regulatory programs.

Fluoride Alters Gene Expression via Histone H3K27 Acetylation in Ameloblast-like LS8 Cells.

Excessive fluoride ingestion during tooth development can cause dental fluorosis. Previously, we reported that fluoride activates histone acetyltransferase (HAT) to acetylate p53, promoting fluoride toxicity in mouse ameloblast-like LS8 cells. However, the roles of HAT and histone acetylation status in fluoride-mediated gene expression remain unidentified. Here, we demonstrate that fluoride-mediated histone modification causes gene expression alterations in LS8 cells. LS8 cells were treated with or without fluoride followed by ChIP-Seq analysis of H3K27ac. Genes were identified by differential H3K27ac peaks within ±1 kb from transcription start sites. The levels of mRNA of identified genes were assessed using rea-time PCR (qPCR). Fluoride increased H3K27ac peaks associated with Bax, p21, and Mdm2 genes and upregulated their mRNA levels. Fluoride decreased H3K27ac peaks and p53, Bad, and Bcl2 had suppressed transcription. HAT inhibitors (Anacardic acid or MG149) suppressed fluoride-induced mRNA of p21 and Mdm2, while fluoride and the histone deacetylase (HDAC) inhibitor sodium butyrate increased Bad and Bcl2 expression above that of fluoride treatment alone. To our knowledge, this is the first study that demonstrates epigenetic regulation via fluoride treatment via H3 acetylation. Further investigation is required to elucidate epigenetic mechanisms of fluoride toxicity in enamel development.

A simple, robust, cost-effective, and low-input ChIP-seq method for profiling histone modifications and Pol II in plants.

Chromatin immunoprecipitation and sequencing (vs ChIP-seq) is an essential tool for epigenetic and molecular genetic studies. Although being routinely used, ChIP-seq is expensive, requires grams of plant materials, and is challenging for samples that enrich fatty acids such as seeds. Here, we developed an Ultrasensitive Plant ChIP-seq (UP-ChIP) method based on native ChIP-seq combined with Tn5 tagmentation-based library construction strategy. UP-ChIP is generally applicable for profiling both histone modification and Pol II in a wide range of plant samples, such as a single Arabidopsis seedling, a few Arabidopsis seeds, and sorted nuclei. Compared with conventional ChIP-seq, UP-ChIP is much less labor intensive and only consumes 1 µg of antibody and 10 µl of Protein-A/G conjugated beads for each IP and can work effectively with the amount of starting material down to a few milligrams. By performing UP-ChIP in various conditions and genotypes, we showed that UP-ChIP is highly reliable, sensitive, and quantitative for studying histone modifications. Detailed UP-ChIP protocol is provided. We recommend UP-ChIP as an alternative to traditional ChIP-seq for profiling histone modifications and Pol II, offering the advantages of reduced labor intensity, decreased costs, and low-sample input.

The homeodomain (PHD) protein, PdbPHD3, confers salt tolerance by regulating squamosa promoter binding protein PdbSBP3 in Populus davidiana × P. bolleana.

Plant homeodomain (PHD) proteins are a family of zinc finger transcription factors that play roles in abiotic stress tolerance. However, their mechanisms in conferring salt tolerance are largely unknown. In this study, we characterized a PHD gene, PdbPHD3, from Populus davidiana × P. bolleana (Shanxin poplar) in response to salt stress. PdbPHD3 is a nuclear protein that is strongly induced by salt and abscisic acid (ABA) treatments. Overexpression of PdbPHD3 conferred salt tolerance, while silencing of PdbPHD3 increased sensitivity to salt. PdbPHD3 could enhance the activities of superoxide dismutase and peroxidase to reduce the abundance of reactive oxygen species, and enhance the osmotic potential by increasing the proline content. Chromatin immunoprecipitation sequencing (ChIP-seq) revealed that PdbPHD3 could bind to various DNA motifs, including the G-box ("CACGTG"), PALBOXAPC ("GGACGG"), and POLLEN1LELAT52 ("TTTCTT"). ChIP-seq combined with RNA sequencing identified a transcription factor gene, squamosa promoter binding protein 3 (PdbSBP3), which is directly regulated by PdbPHD3. Overexpression and silencing of PdbSBP3 improved and decreased salt tolerance, respectively. PdbSBP3 could also regulate all the physiological changes associated with salt tolerance, similar to PdbPHD3. These results suggest that PdbPHD3 confers salt tolerance by regulating PdbSBP3 to reduce ROS accumulation and increase proline content. Therefore, the regulatory axis of PdbPHD3 and PdbSBP3 confers salt tolerance in Shanxin poplar.

Examining NF-κB genomic interactions by ChIP-seq and CUT&Tag.

An understanding of the mechanisms and logic by which transcription factors coordinate gene regulation requires delineation of their genomic interactions at a genome-wide scale. Chromatin immunoprecipitation-sequencing (ChIP-seq) and more recent techniques, including CUT&Tag, typically reveal thousands of genomic interactions by transcription factors, but without insight into their functional roles. Due to cost and time considerations, optimization of ChIP experimental conditions is typically carried out only with representative interaction sites rather than through genome-wide analyses. Here, we describe insights gained from the titration of two chemical crosslinking reagents in genome-wide ChIP-seq experiments examining two members of the NF-κB family of transcription factors: RelA and c-Rel. We also describe a comparison of ChIP-seq and CUT&Tag. Our results highlight the large impact of ChIP-seq experimental conditions on the number of interactions detected, on the enrichment of consensus and non-consensus DNA motifs for the factor, and on the frequency with which the genomic interactions detected are located near potential target genes. We also found considerable consistency between ChIP-seq and CUT&Tag results, but with a substantial fraction of genomic interactions detected with only one of the two techniques. Together, the results demonstrate the dramatic impact of experimental conditions on the results obtained in a genome-wide analysis of transcription factor binding, highlighting the need for further scrutiny of the functional significance of these condition-dependent differences.

SPO-Seq: An Accessible Method for Efficient Evaluation of Spo11 Catalytic Activity and Profiling Meiotic DSB Hotspots in Saccharomyces cerevisiae.

Meiotic recombination is a key process facilitating the formation of crossovers and the exchange of genetic material between homologous chromosomes in early meiosis. This involves controlled double-strand breaks (DSBs) formation catalyzed by Spo11. DSBs exhibit a preferential location in specific genomic regions referred to as hotspots, and their variability is tied to varying Spo11 activity levels. We have refined a ChIP-Seq technique, called SPO-Seq, to map Spo11-specific DSB formation in Saccharomyces cerevisiae. The chapter describes our streamlined approach and the developed bioinformatic tools for processing data and comparing with existing DSB hotspot maps. Through this combined experimental and computational approach, we aim to enhance our understanding of meiotic recombination and genetic exchange processes in budding yeast, with the potential to expand this methodology to other organisms by applying a few modifications.

Efficient Enrichment of Synchronized Mouse Spermatocytes Suitable for Genome-Wide Analysis.

Chromatin undergoes extensive remodeling during meiosis, leading to specific patterns of gene expression and chromosome organization, which ultimately controls fundamental meiotic processes such as recombination and homologous chromosome associations. Recent game-changing advances have been made by analysis of chromatin binding sites of meiotic specific proteins genome-wide in mouse spermatocytes. However, further progress is still highly dependent on the reliable isolation of sufficient quantities of spermatocytes at specific stages of prophase I. Here, we describe a combination of methodologies we adapted for rapid and reliable isolation of synchronized fixed mouse spermatocytes. We show that chromatin isolated from these cells can be used to study chromatin-binding sites by ChIP-seq. High-quality data we obtained from INO80 ChIP-seq in zygotene cells was used for functional analysis of chromatin-binding sites.

Native ChIP: Studying the Genome-Wide Distribution of Histone Modifications in Cells and Tissue.

For the genome-wide mapping of histone modifications, chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing remains the benchmark method. While crosslinked ChIP can be used for all kinds of targets, native ChIP is predominantly used for strong and direct DNA interactors like histones and their modifications. Here we describe a native ChIP protocol that can be used for cells and tissue material.

Bioinformatics Core Workflow for ChIP-Seq Data Analysis.

Chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (-seq) has been the most common genomics method for studying DNA-protein interactions in the last decade. ChIP-seq technology became standard both experimentally and computationally. This chapter presents a core workflow that covers data processing and initial analytical steps of ChIP-seq data. We provide a step-by-step protocol of the commands as well as a fully assembled Snakemake workflow. Along the protocol, we discuss key tool parameters, quality control, output reports, and preliminary results.

Differential Analysis of Protein-DNA Binding Using ChIP-Seq Data.

Chromatin immunoprecipitation in combination with next-generation sequencing (ChIP-Seq) allows probing of protein-DNA binding in a rapid and genome-wide fashion. Herein we describe the required steps to preprocess ChIP-Seq data and to analyze the differential binding of proteins to DNA for perturbation experiments. In these experiments, different conditions are compared to find the underlying biological mechanisms caused by the stimulus or treatment. In addition, we provide a sample analysis using the steps outlined in the chapter.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) of Epigenetic Regulators.

Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) allows for the identification of genomic targeting of DNA-binding proteins. Cleavage Under Targets and Release Using Nuclease (CUT&RUN) modifies this process by including a nuclease to digest DNA around a protein of interest. The result is a higher signal-to-noise ratio and decreased required starting material. This allows for high-fidelity sequence identification from as few as 500 cells, enabling chromatin profiling of precious tissue samples or primary cell types, as well as less abundant chromatin-binding proteins: all at significantly increased throughput.

Integrative Analysis of CUT&Tag and RNA-Seq Data Through Bioinformatics: A Unified Workflow for Enhanced Insights.

Cleavage Under Targets and Tagmentation (CUT&Tag) is a recent methodology used for robust epigenomic profiling that, unlike conventional chromatin immunoprecipitation (ChIP-Seq), requires only a limited amount of cells as starting material. RNA sequencing (RNA-Seq) reveals the presence and quantity of RNA in a biological sample, describing the continuously changing cellular transcriptome. The integrated analysis of transcriptional activity, histone modifications, and chromatin accessibility via CUT&Tag is still in its infancy compared to the well-established ChIP-Seq. This chapter describes a robust bioinformatics methodology and workflow to perform an integrative CUT&Tag/RNA-Seq analysis.

Role of Csdc2 in Regulating Secondary Hair Follicle Growth in Cashmere Goats.

Cashmere goats possess two types of hair follicles, with the secondary hair follicles producing valuable cashmere fiber used for textiles. The growth of cashmere exhibits a seasonal pattern arising from photoperiod change. Transcription factors play crucial roles during this process. The transcription factor, cold-shock domain, containing C2 (Csdc2) plays a crucial role in modulating cell proliferation and differentiation. Our preceding research indicated that the expression of Csdc2 changes periodically during anagen to telogen. However, the mechanisms of Csdc2 in regulating SHF growth remain unclear. Here, we found that the knockdown of Csdc2 inhibits the proliferation of dermal papilla cells. ChIP-Seq analysis showed that Csdc2 had a unique DNA binding motif in SHFs. Through conjoint analysis of ChIP-Seq and RNA-Seq, we revealed a total of 25 candidate target genes of Csdc2. Notably, we discovered a putative Csdc2 binding site within roundabout guidance receptor 2 (Robo2) on chromosome 1 of the goat genome. Furthermore, qRT-PCR and dual-luciferase reporter assay confirmed Csdc2's positive regulatory influence on Robo2. These findings expand the research field of hair follicle transcriptional regulatory networks, offering insights into molecular breeding strategies to enhance cashmere production in goats.

Overlooked poor-quality patient samples in sequencing data impair reproducibility of published clinically relevant datasets.

BACKGROUND: Reproducibility is a major concern in biomedical studies, and existing publication guidelines do not solve the problem. Batch effects and quality imbalances between groups of biological samples are major factors hampering reproducibility. Yet, the latter is rarely considered in the scientific literature. RESULTS: Our analysis uses 40 clinically relevant RNA-seq datasets to quantify the impact of quality imbalance between groups of samples on the reproducibility of gene expression studies. High-quality imbalance is frequent (14 datasets; 35%), and hundreds of quality markers are present in more than 50% of the datasets. Enrichment analysis suggests common stress-driven effects among the low-quality samples and highlights a complementary role of transcription factors and miRNAs to regulate stress response. Preliminary ChIP-seq results show similar trends. Quality imbalance has an impact on the number of differential genes derived by comparing control to disease samples (the higher the imbalance, the higher the number of genes), on the proportion of quality markers in top differential genes (the higher the imbalance, the higher the proportion; up to 22%) and on the proportion of known disease genes in top differential genes (the higher the imbalance, the lower the proportion). We show that removing outliers based on their quality score improves the resulting downstream analysis. CONCLUSIONS: Thanks to a stringent selection of well-designed datasets, we demonstrate that quality imbalance between groups of samples can significantly reduce the relevance of differential genes, consequently reducing reproducibility between studies. Appropriate experimental design and analysis methods can substantially reduce the problem.