Reduction of ETV1 is Identified as a Prominent Feature of Age-Related Cataract

PURPOSE

To identify the inactive genes in cataract lenses and explore their function in lens epithelial cells (LECs).

METHODS

Lens epithelium samples obtained from both age-related cataract (ARC) patients and normal donors were subjected to two forms of histone H3 immunoprecipitation: H3K9ac and H3K27me3 chromatin immunoprecipitation (ChIP), followed by ChIP-seq. The intersection set of "active genes in normal controls" and "repressed genes in cataract lenses" was identified. To validate the role of a specific gene, ETV1, within this set, quantitative polymerase chain reaction (qPCR), western blot, and immunofluorescence were performed using clinical lens epithelium samples. Small interference RNA (siRNA) was utilized to reduce the mRNA level of ETV1 in cultured LECs. Following this, transwell assay and western blot was performed to examine the migration ability of the cells. Furthermore, RNA-seq analysis was conducted on both cell samples with ETV1 knockdown and control cells. Additionally, the expression level of ETV1 in LECs was examined using qPCR under H2O2 treatment.

RESULTS

Six genes were identified in the intersection set of "active genes in normal controls" and "repressed genes in ARC lenses". Among these genes, ETV1 showed the most significant fold-change decrease in the cataract samples compared to the control samples. After ETV1 knockdown by siRNA in cultured LECs, reduced cell migration was observed, along with a decrease in the expression of β-Catenin and Vimentin, two specific genes associated with cell migration. In addition, under the oxidative stress induced by H2O2 treatment, the expression level of ETV1 in LECs significantly decreased.

CONCLUSIONS

Based on the findings of this study, it can be concluded that ETV1 is significantly reduced in human ARC lenses. The repression of ETV1 in ARC lenses appears to contribute to the disrupted differentiation of lens epithelium, which is likely caused by the inhibition of both cell differentiation and migration processes.

Skeletal muscle differentiation induces wide-ranging nucleosome repositioning in muscle gene promoters

In a previous report, we demonstrated that Cbx1, PurB and Sp3 inhibited cardiac muscle differentiation by increasing nucleosome density around cardiac muscle gene promoters. Since cardiac and skeletal muscle express many of the same proteins, we asked if Cbx1, PurB and Sp3 similarly regulated skeletal muscle differentiation. In a C2C12 model of skeletal muscle differentiation, Cbx1 and PurB knockdown increased myotube formation. In contrast, Sp3 knockdown inhibited myotube formation, suggesting that Sp3 played opposing roles in cardiac muscle and skeletal muscle differentiation. Consistent with this finding, Sp3 knockdown also inhibited various muscle-specific genes. The Cbx1, PurB and Sp3 proteins are believed to influence gene-expression in part by altering nucleosome position. Importantly, we developed a statistical approach to determine if changes in nucleosome positioning were significant and applied it to understanding the architecture of muscle-specific genes. Through this novel statistical approach, we found that during myogenic differentiation, skeletal muscle-specific genes undergo a set of unique nucleosome changes which differ significantly from those shown in commonly expressed muscle genes. While Sp3 binding was associated with nucleosome loss, there appeared no correlation with the aforementioned nucleosome changes. In summary, we have identified a novel role for Sp3 in skeletal muscle differentiation and through the application of quantifiable MNase-seq have discovered unique fingerprints of nucleosome changes for various classes of muscle genes during myogenic differentiation.

Principles of chromosome organization for meiotic recombination

In meiotic cells, chromosomes are organized as chromatin loop arrays anchored to a protein axis. This organization is essential to regulate meiotic recombination, from DNA double-strand break (DSB) formation to their repair. In mammals, it is unknown how chromatin loops are organized along the genome and how proteins participating in DSB formation are tethered to the chromosome axes. Here, we identify three categories of axis-associated genomic sites: PRDM9 binding sites, where DSBs form; binding sites of the insulator protein CTCF; and H3K4me3-enriched sites. We demonstrate that PRDM9 promotes the recruitment of MEI4 and IHO1, two proteins essential for DSB formation. In turn, IHO1 anchors DSB sites to the axis components HORMAD1 and SYCP3. We discovered that IHO1, HORMAD1, and SYCP3 are associated at the DSB ends during DSB repair. Our results highlight how interactions of proteins with specific genomic elements shape the meiotic chromosome organization for recombination.

Epigenetic regulation of H3K27me3 in laying hens with fatty liver hemorrhagic syndrome induced by high-energy and low-protein diets

BACKGROUND

Fatty liver hemorrhagic syndrome (FLHS) in the modern poultry industry is primarily caused by nutrition. Despite encouraging progress on FLHS, the mechanism through which nutrition influences susceptibility to FLHS is still lacking in terms of epigenetics.

RESULTS

In this study, we analyzed the genome-wide patterns of trimethylated lysine residue 27 of histone H3 (H3K27me3) enrichment by chromatin immunoprecipitation-sequencing (ChIP-seq), and examined its association with transcriptomes in healthy and FLHS hens. The study results indicated that H3K27me3 levels were increased in the FLHS hens on a genome-wide scale. Additionally, H3K27me3 was found to occupy the entire gene and the distant intergenic region, which may function as silencer-like regulatory elements. The analysis of transcription factor (TF) motifs in hypermethylated peaks has demonstrated that 23 TFs are involved in the regulation of liver metabolism and development. Transcriptomic analysis indicated that differentially expressed genes (DEGs) were enriched in fatty acid metabolism, amino acid, and carbohydrate metabolism. The hub gene identified from PPI network is fatty acid synthase (FASN). Combined ChIP-seq and transcriptome analysis revealed that the increased H3K27me3 and down-regulated genes have significant enrichment in the ECM-receptor interaction, tight junction, cell adhesion molecules, adherens junction, and TGF-beta signaling pathways.

CONCLUSIONS

Overall, the trimethylation modification of H3K27 has been shown to have significant regulatory function in FLHS, mediating the expression of crucial genes associated with the ECM-receptor interaction pathway. This highlights the epigenetic mechanisms of H3K27me3 and provides insights into exploring core regulatory targets and nutritional regulation strategies in FLHS.

Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

ChIP-seq analysis found IL21R, a target gene of GTF2I-the susceptibility gene for primary biliary cholangitis in Chinese Han

AIMS

Aimed to identify a new susceptibility gene associated with primary biliary cholangitis (PBC) in Chinese Han and investigate the possible mechanism of that gene in PBC.

METHODS

A total of 466 PBC and 694 healthy controls (HC) were included in our study, and genotyping GTF2I gene variants by Sequenom. CD19 + B cells were isolated for Chromatin immunoprecipitation sequencing (ChIP-seq). Additionally, MEME-ChIP was utilized to perform searches for known motifs and de novo motif discovery. The GTF2I ChIP-seq of hematopoietic cell line (K562) results were obtained from ENCODE (GSE176987, GSE177691). The Genomic HyperBrowser was used to determine overlap and hierarchal clustering between ours and ENCODE datasets.

RESULTS

The frequency of the rs117026326 variant T allele was significantly higher in PBC patients than that in HC (20.26% compared with 13.89%, Pc = 1.09E-04). Furthermore, we observed an elevated proportion of GTF2I binding site located in the upstream and 5' UTR of genes in PBC in comparison with HC. Additionally, an in-depth analysis of IL21R region revealed that GTF2I might bind to the IL21R promoter to regulate the expression of the IL21R, with four peaks of GTF2I binding sites, including three increased binding sites in upstream, one increased binding site in 5' UTR. Motif analysis by MEME-ChIP uncovered five significant motifs. A significant overlap between our ChIP and GSE176987, GSE17769 were found by the Genomic HyperBroswer.

CONCLUSIONS

Our study confirmed that GTF2I was associated with PBC in Chinese Han. Furthermore, our gene function analysis indicated that IL21R may be the target gene regulated by GTF2I.

Aberrant H3K4me3 modification of immune response genes in CD4+ T cells of patients with systemic lupus erythematosus

BACKGROUND

Increasing evidence has highlighted the significant role of histone modifications in pathogenesis of systemic lupus erythematosus (SLE). However, few studies have comprehensively analyzed trimethylation of histone H3 lysine 4 (H3K4me3) features at specific immune gene loci in SLE patients.

METHODS

We conducted H3K4me3 chromatin immunoprecipitation sequencing (ChIP-seq) on CD4+ T cells from SLE patients and healthy controls (HC). Differential H3K4me3 peaks were identified, followed by enrichment analysis. We integrated online RNA-seq and DNA methylation datasets to explore the relationship between H3K4me3 modification, DNA methylation and gene expression. We validated several upregulated peak regions by ChIP-qPCR and confirmed their impact on gene expression using RT-qPCR. Finally, we investigated the impact of H3K4 methyltransferases KMT2A on the expression of immune response genes.

RESULTS

we identified 147 downregulated and 2701 upregulated H3K4me3 peaks in CD4+ T cells of SLE. The upregulated peaks primarily classified as gained peaks and enriched in immune response genes such as FCGR2A, C5AR1, SERPING1 and OASL. Genes with upregulated H3K4me3 and downregulated DNA methylations in the promoter were highly expressed in SLE patients. These genes, including OAS1, IFI27 and IFI44L, were enriched in immune response pathways. The IFI44L locus also showed increased H3K27ac modification, chromatin accessibility and chromatin interactions in SLE. Moreover, knockdown of KMT2A can downregulate the expression of immune response genes in T cells.

CONCLUSION

Our study uncovers dysregulated H3K4me3 modification patterns in immune response genes loci, which also exhibit downregulated DNA methylation and higher mRNA expression in CD4+ T cells of SLE patients.

The apparent loss of PRC2 chromatin occupancy as an artifact of RNA depletion

RNA has been implicated in the recruitment of chromatin modifiers, and previous studies have provided evidence in favor and against this idea. RNase treatment of chromatin is commonly used to study RNA-mediated regulation of chromatin modifiers, but the limitations of this approach remain unclear. RNase A treatment during chromatin immunoprecipitation (ChIP) reduces chromatin occupancy of the H3K27me3 methyltransferase Polycomb repressive complex 2 (PRC2). This led to suggestions of an "RNA bridge" between PRC2 and chromatin. Here, we show that RNase A treatment during ChIP causes the apparent loss of all facultative heterochromatin, including both PRC2 and H3K27me3 genome-wide. We track this observation to a gain of DNA from non-targeted chromatin, sequenced at the expense of DNA from facultative heterochromatin, which reduces ChIP signals. Our results emphasize substantial limitations in using RNase A treatment for mapping RNA-dependent chromatin occupancy and invalidate conclusions that were previously established for PRC2 based on this assay.

Apparent RNA bridging between PRC2 and chromatin is an artifact of non-specific chromatin precipitation upon RNA degradation

Polycomb repressive complex 2 (PRC2) modifies chromatin to maintain repression of genes specific for other cell lineages. In vitro, RNA inhibits PRC2 activity, but the effect of RNA on PRC2 in cells is less clear, with studies concluding that RNA either antagonizes or promotes PRC2 chromatin association. The addition of RNase A to chromatin immunoprecipitation reactions has been reported to reduce detection of PRC2 target sites, suggesting the existence of RNA bridges connecting PRC2 to chromatin. Here, we show that the apparent loss of PRC2 chromatin association after RNase A treatment is due to non-specific chromatin precipitation. RNA degradation precipitates chromatin out of solution, thereby masking enrichment of specific DNA sequences in chromatin immunoprecipitation reactions. Maintaining chromatin solubility by the addition of poly-L-glutamic acid rescues detection of PRC2 chromatin occupancy upon RNA degradation. These findings undermine support for the model that RNA bridges PRC2 and chromatin in cells.

Protocol for generating customizable and reproducible plots of sequencing coverage data using the seqNdisplayR package

The widespread usage of next-generation sequencing methods for functional genomics studies requires standardized tools for consistent visualization of the associated data. Here, we present seqNdisplayR, an R package for plotting standard sequencing data coverage within a genomic region of interest in a customizable and reproducible manner. We describe steps for installing software, preparing data files, choosing options, and plotting data. This tool is readily available for users with no prior experience with R through the "Shiny app" interface. For complete details on the use and execution of this protocol, please refer to Lykke-Andersen et al.,1 Gockert et al.,2 and Rouviere et al.3.

Epigenetic control of microglial immune responses

Microglia, the major population of brain-resident macrophages, are now recognized as a heterogeneous population comprising several cell subtypes with different (so far mostly supposed) functions in health and disease. A number of studies have performed molecular characterization of these different microglial activation states over the last years making use of "omics" technologies, that is transcriptomics, proteomics and, less frequently, epigenomics profiling. These approaches offer the possibility to identify disease mechanisms, discover novel diagnostic biomarkers, and develop new therapeutic strategies. Here, we focus on epigenetic profiling as a means to understand microglial immune responses beyond what other omics methods can offer, that is, revealing past and present molecular responses, gene regulatory networks and potential future response trajectories, and defining cell subtype-specific disease relevance through mapping non-coding genetic variants. We review the current knowledge in the field regarding epigenetic regulation of microglial identity and function, provide an exemplary analysis that demonstrates the advantages of performing joint transcriptomic and epigenomic profiling of single microglial cells and discuss how comprehensive epigenetic analyses may enhance our understanding of microglial pathophysiology.

A fluorescence-based protocol to quantitatively titrate CUT&RUN buffer components

Cleavage under targets & release using nuclease (CUT&RUN) is a technique for identifying genomic sites where proteins or histone modifications are present in chromatin in permeabilized cells. Here, we present a fluorescence-based protocol to quantitatively titrate CUT&RUN buffer components, for efficient cell permeabilization and retention of target epitopes on chromatin. We describe steps for capturing cells on concanavalin A beads and using a fluorescently labeled secondary antibody to titrate concentrations of digitonin and NaCl in CUT&RUN buffers. We then detail procedures for fluorescence imaging to identify optimal conditions. For complete details on the use and execution of this protocol, please refer to Lerner et al.1.

Shared and distinct interactions of type 1 and type 2 Epstein-Barr Nuclear Antigen 2 with the human genome

BACKGROUND

There are two major genetic types of Epstein-Barr Virus (EBV): type 1 (EBV-1) and type 2 (EBV-2). EBV functions by manipulating gene expression in host B cells, using virus-encoded gene regulatory proteins including Epstein-Barr Nuclear Antigen 2 (EBNA2). While type 1 EBNA2 is known to interact with human transcription factors (hTFs) such as RBPJ, EBF1, and SPI1 (PU.1), type 2 EBNA2 shares only ~ 50% amino acid identity with type 1 and thus may have distinct binding partners, human genome binding locations, and functions.

RESULTS

In this study, we examined genome-wide EBNA2 binding in EBV-1 and EBV-2 transformed human B cells to identify shared and unique EBNA2 interactions with the human genome, revealing thousands of type-specific EBNA2 ChIP-seq peaks. Computational predictions based on hTF motifs and subsequent ChIP-seq experiments revealed that both type 1 and 2 EBNA2 co-occupy the genome with SPI1 and AP-1 (BATF and JUNB) hTFs. However, type 1 EBNA2 showed preferential co-occupancy with EBF1, and type 2 EBNA2 preferred RBPJ. These differences in hTF co-occupancy revealed possible mechanisms underlying type-specific gene expression of known EBNA2 human target genes: MYC (shared), CXCR7 (type 1 specific), and CD21 (type 2 specific). Both type 1 and 2 EBNA2 binding events were enriched at systemic lupus erythematosus (SLE) and multiple sclerosis (MS) risk loci, while primary biliary cholangitis (PBC) risk loci were specifically enriched for type 2 peaks.

CONCLUSIONS

This study reveals extensive type-specific EBNA2 interactions with the human genome, possible differences in EBNA2 interaction partners, and a possible new role for type 2 EBNA2 in autoimmune disorders. Our results highlight the importance of considering EBV type in the control of human gene expression and disease-related investigations.

Loop-extruding Smc5/6 organizes transcription-induced positive DNA supercoils

The structural maintenance of chromosomes (SMC) protein complexes-cohesin, condensin, and the Smc5/6 complex (Smc5/6)-are essential for chromosome function. At the molecular level, these complexes fold DNA by loop extrusion. Accordingly, cohesin creates chromosome loops in interphase, and condensin compacts mitotic chromosomes. However, the role of Smc5/6's recently discovered DNA loop extrusion activity is unknown. Here, we uncover that Smc5/6 associates with transcription-induced positively supercoiled DNA at cohesin-dependent loop boundaries on budding yeast (Saccharomyces cerevisiae) chromosomes. Mechanistically, single-molecule imaging reveals that dimers of Smc5/6 specifically recognize the tip of positively supercoiled DNA plectonemes and efficiently initiate loop extrusion to gather the supercoiled DNA into a large plectonemic loop. Finally, Hi-C analysis shows that Smc5/6 links chromosomal regions containing transcription-induced positive supercoiling in cis. Altogether, our findings indicate that Smc5/6 controls the three-dimensional organization of chromosomes by recognizing and initiating loop extrusion on positively supercoiled DNA.

Automated ChIPmentation procedure on limited biological material of the human blood fluke Schistosoma mansoni

In living cells, the genetic information stored in the DNA sequence is always associated with chromosomal and extra-chromosomal epigenetic information. Chromatin is formed by the DNA and associated proteins, in particular histones. Covalent histone modifications are important bearers of epigenetic information and as such have been increasingly studied since about the year 2000. One of the principal techniques to gather information about the association between DNA and modified histones is chromatin immunoprecipitation (ChIP), also combined with massive sequencing (ChIP-Seq). Automated ChIPmentation procedure is a convenient alternative to native chromatin immunoprecipitation (N-ChIP). It is now routinely used for ChIP-Seq in many model species, using in general roughly 10 6 cells per experiment. Such high cell numbers are sometimes difficult to produce. Using the human parasite Schistosoma mansoni, whose production requires sacrificing animals and should therefore be kept to a minimum, we show here that automated ChIPmentation is suitable for limited biological material. We define the operational limit as ≥20,000 Schistosoma cells with 30,000-300,000 cells as optimum. We also present a streamlined protocol for the preparation of ChIP input libraries.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization

The co-visualization of chromatin conformation with 1D 'omics data is key to the multi-omics driven data analysis of 3D genome organization. Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface (GUI) application, available in both web or containerized format to promote accessibility. HiCrayon is implemented in R, and includes a graphical user interface (GUI), as well as a slimmed-down web-based version that lets users quickly produce publication-ready images. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin vs. those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct 3rd chromatin compartment. Overall, we demonstrate the utility of co-visualizing 2D chromatin conformation with 1D genomic signals within the same matrix to reveal fundamental aspects of genome organization. Local version: https://github.com/JRowleyLab/HiCrayon Web version: https://jrowleylab.com/HiCrayon.

CIC-DUX4 Chromatin Profiling Reveals New Epigenetic Dependencies and Actionable Therapeutic Targets in CIC-Rearranged Sarcomas

CIC-DUX4-rearranged sarcoma (CDS) is a rare and aggressive soft tissue tumor that occurs most frequently in young adults. The key oncogenic driver of this disease is the expression of the CIC-DUX4 fusion protein as a result of chromosomal rearrangements. CIC-DUX4 displays chromatin binding properties, and is therefore believed to function as an aberrant transcription factor. However, the chromatin remodeling events induced by CIC-DUX4 are not well understood, limiting our ability to identify new mechanism-based therapeutic strategies for these patients. Here, we generated a genome-wide profile of CIC-DUX4 DNA occupancy and associated chromatin states in human CDS cell models and primary tumors. Combining chromatin profiling, proximity ligation assays, as well as genetic and pharmacological perturbations, we show that CIC-DUX4 operates as a potent transcriptional activator at its binding sites. This property is in contrast with the repressive function of the wild-type CIC protein, and is mainly mediated through the direct interaction of CIC-DUX4 with the acetyltransferase p300. In keeping with this, we show p300 to be essential for CDS tumor cell proliferation; additionally, we find its pharmacological inhibition to significantly impact tumor growth in vitro and in vivo. Taken together, our study elucidates the mechanisms underpinning CIC-DUX4-mediated transcriptional regulation.

Inducible Expression of the Restriction Enzyme Uncovered Genome-Wide Distribution and Dynamic Behavior of Histones H4K16ac and H2A.Z at DNA Double-Strand Breaks in Arabidopsis

DNA double-strand breaks (DSBs) are among the most serious types of DNA damage, causing mutations and chromosomal rearrangements. In eukaryotes, DSBs are immediately repaired in coordination with chromatin remodeling for the deposition of DSB-related histone modifications and variants. To elucidate the details of DSB-dependent chromatin remodeling throughout the genome, artificial DSBs need to be reproducibly induced at various genomic loci. Recently, a comprehensive method for elucidating chromatin remodeling at multiple DSB loci via chemically induced expression of a restriction enzyme was developed in mammals. However, this DSB induction system is unsuitable for investigating chromatin remodeling during and after DSB repair, and such an approach has not been performed in plants. Here, we established a transgenic Arabidopsis plant harboring a restriction enzyme gene Sbf I driven by a heat-inducible promoter. Using this transgenic line, we performed chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) of histones H4K16ac and H2A.Z and investigated the dynamics of these histone marks around the endogenous 623 Sbf I recognition sites. We also precisely quantified DSB efficiency at all cleavage sites using the DNA resequencing data obtained by the ChIP-seq procedure. From the results, Sbf I-induced DSBs were detected at 360 loci, which induced the transient deposition of H4K16ac and H2A.Z around these regions. Interestingly, we also observed the co-localization of H4K16ac and H2A.Z at some DSB loci. Overall, DSB-dependent chromatin remodeling was found to be highly conserved between plants and animals. These findings provide new insights into chromatin remodeling that occurs in response to DSBs in Arabidopsis.

Ascertainment Bias in the Genomic Test of Positive Selection on Regulatory Sequences

Evolution of gene expression mediated by cis-regulatory changes is thought to be an important contributor to organismal adaptation, but identifying adaptive cis-regulatory changes is challenging due to the difficulty in knowing the expectation under no positive selection. A new approach for detecting positive selection on transcription factor binding sites (TFBSs) was recently developed, thanks to the application of machine learning in predicting transcription factor (TF) binding affinities of DNA sequences. Given a TFBS sequence from a focal species and the corresponding inferred ancestral sequence that differs from the former at n sites, one can predict the TF-binding affinities of many n-step mutational neighbors of the ancestral sequence and obtain a null distribution of the derived binding affinity, which allows testing whether the binding affinity of the real derived sequence deviates significantly from the null distribution. Applying this test genomically to all experimentally identified binding sites of 3 TFs in humans, a recent study reported positive selection for elevated binding affinities of TFBSs. Here, we show that this genomic test suffers from an ascertainment bias because, even in the absence of positive selection for strengthened binding, the binding affinities of known human TFBSs are more likely to have increased than decreased in evolution. We demonstrate by computer simulation that this bias inflates the false positive rate of the selection test. We propose several methods to mitigate the ascertainment bias and show that almost all previously reported positive selection signals disappear when these methods are applied.

HBP1 promotes chicken preadipocyte proliferation via directly repressing SOCS3 transcription

Preadipocyte proliferation is an essential process in adipose development. During proliferation of preadipocytes, transcription factors play crucial roles. HMG-box protein 1 (HBP1) is an important transcription factor of cellular proliferation. However, the function and underlying mechanisms of HBP1 in the proliferation of preadipocytes remain unclear. Here, we found that the expression level of HBP1 decreased first and then increased during the proliferation of chicken preadipocytes. Knockout of HBP1 could inhibit the proliferation of preadipocytes, while overexpression of HBP1 could promote the proliferation of preadipocytes. ChIP-seq data showed that HBP1 had the unique DNA binding motif in chicken preadipocytes. By integrating ChIP-Seq and RNA-Seq, we revealed a total of 3 candidate target genes of HBP1. Furthermore, the results of ChIP-qPCR, RT-qPCR, luciferase reporter assay and EMSA showed that HBP1 could inhibit the transcription of suppressor of cytokine signaling 3 (SOCS3) by binding to its promoter. Moreover, we confirmed that SOCS3 can mediate the regulation of HBP1 on the proliferation of preadipocytes through RNAi and rescue experiments. Altogether, these data demonstrated that HBP1 directly targets SOCS3 to regulate chicken preadipocyte proliferation. Our findings expand the transcriptional regulatory network of preadipocyte proliferation, and they will be helpful in formulating a molecular breeding scheme to control excessive abdominal fat deposition and to improve meat quality in chickens.

The dual role of GoPGF reveals that the pigment glands are synthetic sites of gossypol in aerial parts of cotton

Gossypol and the related terpenoids are stored in the pigment gland to protect cotton plants from biotic stresses, but little is known about the synthetic sites of these metabolites. Here, we showed that GoPGF, a key gene regulating gland formation, was expressed in gland cells and roots. The chromatin immunoprecipitation sequencing (ChIP-seq) analysis demonstrated that GoPGF targets GhJUB1 to regulate gland morphogenesis. RNA-sequencing (RNA-seq) showed high accumulation of gossypol biosynthetic genes in gland cells. Moreover, integrated analysis of the ChIP-seq and RNA-seq data revealed that GoPGF binds to the promoter of several gossypol biosynthetic genes. The cotton callus overexpressing GoPGF had dramatically increased the gossypol levels, indicating that GoPGF can directly activate the biosynthesis of gossypol. In addition, the gopgf mutant analysis revealed the existence of both GoPGF-dependent and -independent regulation of gossypol production in cotton roots. Our study revealed that the pigment glands are synthetic sites of gossypol in aerial parts of cotton and that GoPGF plays a dual role in regulating gland morphogenesis and gossypol biosynthesis. The study provides new insights for exploring the complex relationship between glands and the metabolites they store in cotton and other plant species.

NRF2 interacts with distal enhancer and inhibits nitric oxide synthase 2 expression in KRAS-driven pancreatic cancer cells

Nitric oxide is a pleiotropic free radical produced by three nitric oxide synthases (NOS1-3), of which inducible NOS2 is involved in tumor initiation and progression. In this study, RNA-seq, ChIP-seq and qRT-PCR experiments combined with bioinformatic analyses showed that NRF2 is a repressor of NOS2 gene by maintaining a distal enhancer located 22 kb downstream of TSS in an inactive state. Deletion of NRF2 leads to activation of the enhancer, which exerts a pioneering function before it is fully activated. Specifically, NRF2 controls the expression of NOS2 in response to intracellular oxidative stress and extracellular oxygen pressure. We found that abrogation of NOS2 expression by siRNAs partially reduced the ability of WT Panc-1 cells to form 3D spheroids, but strongly reduced the formation of 3D spheroids by NRF2-depleted Panc-1 cells. Mechanistically, this effect correlates with the finding that NOS2 and nitric oxide stimulate epithelial-to-mesenchymal transition in NRF2-depleted Panc-1 and MIA PaCa-2 cells. We also found that knockdown of NOS2 leads to blockade of 3D matrigel invasion of NRF2-depleted PDAC cells, demonstrating that a short-circuit in the reciprocal regulation of NOS2 and NRF2 attenuates the malignancy of PDAC cells. In summary, we show for the first time that: (i) NRF2 is a suppressor of NOS2 in pancreatic cancer cells; (ii) NRF2 binds to and inactivates an enhancer located 22 kb downstream of TSS of the NOS2 gene; (iii) activation of NOS2 requires suppression of NRF2; (iv) NOS2 is required for NRF2-depleted Panc-1 cells to maintain their malignancy and invasiveness.

Ethylene inhibits photosynthesis via temporally distinct responses in tomato plants

Ethylene is a volatile plant hormone that regulates many developmental processes and responses towards (a)biotic stress. Studies have shown that high levels of ethylene repress vegetative growth in many important crops, including tomato (Solanum lycopersicum), possibly by inhibiting photosynthesis. We investigated the temporal effects of ethylene on young tomato plants using an automated ethylene gassing system to monitor the physiological, biochemical, and molecular responses through time-course RNA-seq of a photosynthetically active source leaf. We found that ethylene evokes a dose-dependent inhibition of photosynthesis, which can be characterized by three temporally distinct phases. The earliest ethylene responses that marked the first phase and occurred a few hours after the start of the treatment were leaf epinasty and a decline in stomatal conductance, which led to lower light perception and CO2 uptake, respectively, resulting in a rapid decline of soluble sugar levels (glucose, fructose). The second phase of the ethylene effect was marked by low carbohydrate availability, which modulated plant energy metabolism to adapt by using alternative substrates (lipids and proteins) to fuel the TCA cycle. Long-term continuous exposure to ethylene led to the third phase, characterized by starch and chlorophyll breakdown, which further inhibited photosynthesis, leading to premature leaf senescence. To reveal early (3 h) ethylene-dependent regulators of photosynthesis, we performed a ChIP-seq experiment using anti-ETHYLENE INSENSITIVE 3-like 1 (EIL1) antibodies and found several candidate transcriptional regulators. Collectively, our study revealed a temporal sequence of events that led to the inhibition of photosynthesis by ethylene and identified potential transcriptional regulators responsible for this regulation.

PAPerFly: Partial Assembly-based Peak Finder for ab initio binding site reconstruction

BACKGROUND

The specific recognition of a DNA locus by a given transcription factor is a widely studied issue. It is generally agreed that the recognition can be influenced not only by the binding motif but by the larger context of the binding site. In this work, we present a novel heuristic algorithm that can reconstruct the unique binding sites captured in a sequencing experiment without using the reference genome.

RESULTS

We present PAPerFly, the Partial Assembly-based Peak Finder, a tool for the binding site and binding context reconstruction from the sequencing data without any prior knowledge. This tool operates without the need to know the reference genome of the respective organism. We employ algorithmic approaches that are used during genome assembly. The proposed algorithm constructs a de Bruijn graph from the sequencing data. Based on this graph, sequences and their enrichment are reconstructed using a novel heuristic algorithm. The reconstructed sequences are aligned and the peaks in the sequence enrichment are identified. Our approach was tested by processing several ChIP-seq experiments available in the ENCODE database and comparing the results of Paperfly and standard methods.

CONCLUSIONS

We show that PAPerFly, an algorithm tailored for experiment analysis without the reference genome, yields better results than an aggregation of ChIP-seq agnostic tools. Our tool is freely available at https://github.com/Caeph/paperfly/ or on Zenodo ( https://doi.org/10.5281/zenodo.7116424 ).

Tho2 is critical for the recruitment of Rrp6 to chromatin in response to perturbed mRNP biogenesis

The eukaryotic THO complex coordinates the assembly of so-called messenger RNA-ribonucleoprotein particles (mRNPs), a process that involves cotranscriptional coating of nascent mRNAs with proteins. Once formed, mRNPs undergo a quality control step that marks them either for active transport to the cytoplasm, or Rrp6/RNA exosome-mediated degradation in the nucleus. However, the mechanism behind the quality control of nascent mRNPs is still unclear. We investigated the cotranscriptional quality control of mRNPs in budding yeast by expressing the bacterial Rho helicase, which globally perturbs yeast mRNP formation. We examined the genome-wide binding profiles of the THO complex subunits Tho2, Thp2, Hpr1, and Mft1 upon perturbation of the mRNP biogenesis, and found that Tho2 plays two roles. In addition to its function as a subunit of the THO complex, upon perturbation of mRNP biogenesis Tho2 targets Rrp6 to chromatin via its carboxy-terminal domain. Interestingly, other THO subunits are not enriched on chromatin upon perturbation of mRNP biogenesis and are not necessary for localizing Rrp6 at its target loci. Our study highlights the potential role of Tho2 in cotranscriptional mRNP quality control, which is independent of other THO subunits. Considering that both the THO complex and the RNA exosome are evolutionarily highly conserved, our findings are likely relevant for mRNP surveillance in mammals.

Benchmarking DNA binding affinity models using allele-specific transcription factor binding data

Transcription factors (TFs) bind to DNA in a highly sequence-specific manner. This specificity can manifest itself in vivo at heterozygous loci as a difference in TF occupancy between the two alleles. When applied on a genomic scale, functional genomic assays such as ChIP-seq typically lack the statistical power to detect allele-specific binding (ASB) at the level of individual variants. To address this, we propose a framework for benchmarking sequence-to-affinity models for TF binding in terms of their ability to predict allelic imbalances in ChIP-seq counts. We show that a likelihood function based on an over-dispersed binomial distribution can aggregate evidence for allelic preference across the genome without requiring statistical significance for individual variants. This allows us to systematically compare predictive performance when multiple binding models for the same TF are available. We introduce PyProBound, an easily extensible reimplementation of the ProBound biophysically interpretable machine learning framework. Configuring PyProBound to explicitly account for a confounding sequence-specific bias in DNA fragmentation rate yields improved TF binding models when training on ChIP-seq data. We also show how our likelihood function can be leveraged to perform de novo motif discovery on the raw allele-aware ChIP-seq counts.

An optimized protocol for chromatin immunoprecipitation from murine inguinal white adipose tissue

Chromatin immunoprecipitation (ChIP) protocols have been used to reveal protein-DNA interactions of various cell types and tissues; however, optimization is required for each specific type of sample. Here, we present a ChIP protocol from murine inguinal white adipose tissue. We describe steps for tissue harvesting, crosslinking, chromatin extraction, shearing, immunoprecipitation, and purification. We then detail procedures for analysis including library preparation, sequencing, and qRT-PCR validation. For complete details on the use and execution of this protocol, please refer to Antonia Katsouda et al. (2022).1.

Understanding the role of Pax5 in development of taxane-resistant neuroendocrine like prostate cancers

Resistance to the current Androgen Receptor Signaling Inhibitor (ARSI) therapies has led to higher incidences of therapy-induced neuroendocrine-like prostate cancer (t-NEPC). This highly aggressive subtype with predominant small cell-like characteristics is resistant to taxane chemotherapies and has a dismal overall survival. t-NEPCs are mostly treated with platinum-based drugs with a combination of etoposide or taxane and have less selectivity and high systemic toxicity, which often limit their clinical potential. During t-NEPC transformation, adenocarcinomas lose their luminal features and adopt neuro-basal characteristics. Whether the adaptive neuronal characteristics of t-NEPC are responsible for such taxane resistance remains unknown. Pathway analysis from patient gene-expression databases indicates that t-NEPC upregulates various neuronal pathways associated with enhanced cellular networks. To identify transcription factor(s) (TF) that could be important for promoting the gene expression for neuronal characters in t-NEPC, we performed ATAC-Seq, acetylated-histone ChIP-seq, and RNA-seq in our NE-like cell line models and analyzed the promoters of transcriptionally active and significantly enriched neuroendocrine-like (NE-like) cancer-specific genes. Our results indicate that Pax5 could be an important transcription factor for neuronal gene expression and specific to t-NEPC. Pathway analysis revealed that Pax5 expression is involved in axonal guidance, neurotransmitter regulation, and neuronal adhesion, which are critical for strong cellular communications. Further results suggest that depletion of Pax5 disrupts cellular interaction in NE-like cells and reduces surface growth factor receptor activation, thereby, sensitizing them to taxane therapies. Moreover, t-NEPC specific hydroxymethylation of Pax5 promoter CpG islands favors Pbx1 binding to induce Pax5 expression. Based on our study, we concluded that continuous exposure to ARSI therapies leads to epigenetic modifications and Pax5 activation in t-NEPC, which promotes the expression of genes necessary to adopt taxane-resistant NE-like cancer. Thus, targeting the Pax5 axis can be beneficial for reverting their taxane sensitivity.

FoxK1 associated gene regulatory network in hepatic insulin action and its relationship to FoxO1 and insulin receptor mediated transcriptional regulation

OBJECTIVE

Insulin acts on the liver via changes in gene expression to maintain glucose and lipid homeostasis. This study aimed to the Forkhead box protein K1 (FOXK1) associated gene regulatory network as a transcriptional regulator of hepatic insulin action and to determine its role versus FoxO1 and possible actions of the insulin receptor at the DNA level.

METHODS

Genome-wide analysis of FoxK1 binding were studied by chromatin immunoprecipitation sequencing and compared to those for IR and FoxO1. These were validated by knockdown experiments and gene expression analysis.

RESULTS

Chromatin immunoprecipitation (ChIP) sequencing shows that FoxK1 binds to the proximal promoters and enhancers of over 4000 genes, and insulin enhances this interaction for about 75% of them. These include genes involved in cell cycle, senescence, steroid biosynthesis, autophagy, and metabolic regulation, including glucose metabolism and mitochondrial function and are enriched in a TGTTTAC consensus motif. Some of these genes are also bound by FoxO1. Comparing this FoxK1 ChIP-seq data to that of the insulin receptor (IR) reveals that FoxK1 may act as the transcription factor partner for some of the previously reported roles of IR in gene regulation, including for LARS1 and TIMM22, which are involved in rRNA processing and cell cycle.

CONCLUSION

These data demonstrate that FoxK1 is an important regulator of gene expression in response to insulin in liver and may act in concert with FoxO1 and IR in regulation of genes in metabolism and other important biological pathways.

A cost-free CURE: using bioinformatics to identify DNA-binding factors at a specific genomic locus

Research experiences provide diverse benefits for undergraduates. Many academic institutions have adopted course-based undergraduate research experiences (CUREs) to improve student access to research opportunities. However, potential instructors of a CURE might still face financial or practical hurdles that prevent implementation. Bioinformatics research offers an alternative that is free, safe, compatible with remote learning, and may be more accessible for students with disabilities. Here, we describe a bioinformatics CURE that leverages publicly available datasets to discover novel proteins that target an instructor-determined genomic locus of interest. We use the free, user-friendly bioinformatics platform Galaxy to map ChIP-seq datasets to a genome, which removes the computing burden from students. Both faculty and students directly benefit from this CURE, as faculty can perform candidate screens and publish CURE results. Students gain not only basic bioinformatics knowledge, but also transferable skills, including scientific communication, database navigation, and primary literature experience. The CURE is flexible and can be expanded to analyze different types of high-throughput data or to investigate different genomic loci in any species.

Ovarian ERβ cistrome and transcriptome reveal chromatin interaction with LRH-1

BACKGROUND

Estrogen receptor beta (ERβ, Esr2) plays a pivotal role in folliculogenesis and ovulation, yet its exact mechanism of action is mainly uncharacterized.

RESULTS

We here performed ERβ ChIP-sequencing of mouse ovaries followed by complementary RNA-sequencing of wild-type and ERβ knockout ovaries. By integrating the ERβ cistrome and transcriptome, we identified its direct target genes and enriched biological functions in the ovary. This demonstrated its strong impact on genes regulating organism development, cell migration, lipid metabolism, response to hypoxia, and response to estrogen. Cell-type deconvolution analysis of the bulk RNA-seq data revealed a decrease in luteal cells and an increased proportion of theca cells and a specific type of cumulus cells upon ERβ loss. Moreover, we identified a significant overlap with the gene regulatory network of liver receptor homolog 1 (LRH-1, Nr5a2) and showed that ERβ and LRH-1 extensively bound to the same chromatin locations in granulosa cells. Using ChIP-reChIP, we corroborated simultaneous ERβ and LRH-1 co-binding at the ERβ-repressed gene Greb1 but not at the ERβ-upregulated genes Cyp11a1 and Fkbp5. Transactivation assay experimentation further showed that ERβ and LRH-1 can inhibit their respective transcriptional activity at classical response elements.

CONCLUSIONS

By characterizing the genome-wide endogenous ERβ chromatin binding, gene regulations, and extensive crosstalk between ERβ and LRH-1, along with experimental corroborations, our data offer genome-wide mechanistic underpinnings of ovarian physiology and fertility.

Promoter-Specific Variants in NeuroD1 and H3K4me3 Coincident Regions and Clinical Outcomes of Small Cell Lung Cancer

BACKGROUND

Neurogenic differentiation 1 (NeuroD1) is a representative small cell lung cancer (SCLC) transcription regulator involved in the carcinogenesis and behavior of SCLC. Histone modifications play an important role in transcription, and H3 lysine 4 trimethylation (H3K4me3) is primarily associated with promoter regions.

METHODS

We investigated the association between single nucleotide polymorphisms (SNPs) in NeuroD1 and H3K4me3 coincident regions, selected using ChIP sequencing (ChIP-seq), and the clinical outcomes of 261 patients with SCLC.

RESULTS

Among 230 SNPs, two were significantly associated with both the chemotherapy response and overall survival (OS) of patients with SCLC. RNF145 rs2043268A>G was associated with worse chemotherapy response and OS (under a recessive model, adjusted odds ratio [aOR], 0.50, 95% confidence interval [CI], 0.26-0.94, P = 0.031, and adjusted hazard ratio [aHR], 1.88, 95% CI, 1.38-2.57, P < 0.001). CINP rs762105A>G was also associated with worse chemotherapy response and OS (under a dominant model, aOR, 0.47, 95% CI, 0.23-0.99, P = 0.046, and aHR, 2.03, 95% CI, 1.47-2.82, P < 0.001). ChIP-quantitative polymerase chain reaction and luciferase assay confirmed that the two SNPs were located in the active promoter regions and influenced the promoter activity of each gene.

CONCLUSION

To summarize, among SNPs selected using ChIP-seq in promoter regions with high peaks in both NeuroD1 and H3K4me3, RNF145 rs2043268A>G and CINP rs762105A>G were associated with clinical outcomes in patients with SCLC and also affected the promoter activity of each gene.

A SNAI2/CTCF Interaction is Required for NOTCH1 Expression in Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric malignancy of the muscle with characteristics of cells blocked in differentiation. NOTCH1 is an oncogene that promotes self-renewal and blocks differentiation in the fusion negative-RMS sub-type. However, how NOTCH1 expression is transcriptionally maintained in tumors is unknown. Analyses of SNAI2 and CTCF chromatin binding and HiC analyses revealed a conserved SNAI2/CTCF overlapping peak downstream of the NOTCH1 locus marking a sub-topologically associating domain (TAD) boundary. Deletion of the SNAI2-CTCF peak showed that it is essential for NOTCH1 expression and viability of FN-RMS cells. Reintroducing constitutively activated NOTCH1-ΔE in cells with the SNAI2-CTCF peak deleted restored cell-viability. Ablation of SNAI2 using CRISPR/Cas9 reagents resulted in the loss of majority of RD and SMS-CTR FN-RMS cells. However, the few surviving clones that repopulate cultures have recovered NOTCH1. Cells that re-establish NOTCH1 expression after SNAI2 ablation are unable to differentiate robustly as SNAI2 shRNA knockdown cells; yet, SNAI2-ablated cells continued to be exquisitely sensitive to ionizing radiation. Thus, we have uncovered a novel mechanism by which SNAI2 and CTCF maintenance of a sub-TAD boundary promotes rather than represses NOTCH1 expression. Further, we demonstrate that SNAI2 suppression of apoptosis post-radiation is independent of SNAI2/NOTCH1 effects on self-renewal and differentiation.

ZNF643/ZFP69B Exerts Oncogenic Properties and Associates with Cell Adhesion and Immune Processes

The global cancer burden remains high; thus, a better understanding of the molecular mechanisms driving carcinogenesis is needed to improve current prevention and treatment options. We previously detected the ZNF643/ZFP69B gene upregulated in multiple tumors, and we speculated it may play a role in tumor biology. To test this hypothesis, we employed TCGA-centered databases to correlate ZNF643 status with various clinicopathological parameters. We also performed RNA-seq analysis and in vitro studies assessing cancer cell phenotypes, and we searched for ZNF643-bound genomic loci. Our data indicated higher levels of ZNF643 in most analyzed tumors compared to normal samples, possibly due to copy number variations. ZNF643 mRNA correlated with diverse molecular and immune subtypes and clinicopathological features (tumor stage, grade, patient survival). RNA-seq analysis revealed that ZNF643 silencing triggers the deregulation of the genes implicated in various cancer-related processes, such as growth, adhesion, and immune system. Moreover, we observed that ZNF643 positively influences cell cycle, migration, and invasion. Finally, our ChIP-seq analysis indicated that the genes associated with ZNF643 binding are linked to adhesion and immune signaling. In conclusion, our data confirm the oncogenic properties of ZNF643 and pinpoint its impact on cell adhesion and immune processes.

Seq2science: an end-to-end workflow for functional genomics analysis

Sequencing databases contain enormous amounts of functional genomics data, making them an extensive resource for genome-scale analysis. Reanalyzing publicly available data, and integrating it with new, project-specific data sets, can be invaluable. With current technologies, genomic experiments have become feasible for virtually any species of interest. However, using and integrating this data comes with its challenges, such as standardized and reproducible analysis. Seq2science is a multi-purpose workflow that covers preprocessing, quality control, visualization, and analysis of functional genomics sequencing data. It facilitates the downloading of sequencing data from all major databases, including NCBI SRA, EBI ENA, DDBJ, GSA, and ENCODE. Furthermore, it automates the retrieval of any genome assembly available from Ensembl, NCBI, and UCSC. It has been tested on a variety of species, and includes diverse workflows such as ATAC-, RNA-, and ChIP-seq. It consists of both generic as well as advanced steps, such as differential gene expression or peak accessibility analysis and differential motif analysis. Seq2science is built on the Snakemake workflow language and thus can be run on a range of computing infrastructures. It is available at https://github.com/vanheeringen-lab/seq2science.

Genes positively regulated by Mef2c in cortical neurons are enriched for common genetic variation associated with IQ and educational attainment

The myocyte enhancer factor 2 C (MEF2C) gene encodes a transcription factor important for neurogenesis and synapse development and contains common variants associated with intelligence (IQ) and educational attainment (EA). Here, we took gene expression data from the mouse cortex of a Mef2c mouse model with a heterozygous DNA binding-deficient mutation of Mef2c (Mef2c-het) and combined these data with MEF2C ChIP-seq data from cortical neurons and single-cell data from the mouse brain. This enabled us to create a set of genes that were differentially regulated in Mef2c-het mice, represented direct target genes of MEF2C and had elevated in expression in cortical neurons. We found this gene-set to be enriched for genes containing common genetic variation associated with IQ and EA. Genes within this gene-set that were down-regulated, i.e. have reduced expression in Mef2c-het mice versus controls, were specifically significantly enriched for both EA and IQ associated genes. These down-regulated genes were enriched for functionality in the adenylyl cyclase signalling system, which is known to positively regulate synaptic transmission and has been linked to learning and memory. Within the adenylyl cyclase signalling system, three genes regulated by MEF2C, CRHR1, RGS6, and GABRG3, are associated at genome-wide significant levels with IQ and/or EA. Our results indicate that genetic variation in MEF2C and its direct target genes within cortical neurons contribute to variance in cognition within the general population, and the molecular mechanisms involved include the adenylyl cyclase signalling system's role in synaptic function.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

The human adenovirus E1B-55K oncoprotein coordinates cell transformation through regulation of DNA-bound host transcription factors

The multifunctional adenovirus E1B-55K oncoprotein can induce cell transformation in conjunction with adenovirus E1A gene products. Previous data from transient expression studies and in vitro experiments suggest that these growth-promoting activities correlate with E1B-55K-mediated transcriptional repression of p53-targeted genes. Here, we analyzed genome-wide occupancies and transcriptional consequences of species C5 and A12 E1B-55Ks in transformed mammalian cells by combinatory ChIP and RNA-seq analyses. E1B-55K-mediated repression correlates with tethering of the viral oncoprotein to p53-dependent promoters via DNA-bound p53. Moreover, we found that E1B-55K also interacts with and represses transcription of numerous p53-independent genes through interactions with transcription factors that play central roles in cancer and stress signaling. Our results demonstrate that E1B-55K oncoproteins function as promiscuous transcriptional repressors of both p53-dependent and -independent genes and further support the model that manipulation of cellular transcription is central to adenovirus-induced cell transformation and oncogenesis.

A framework to validate fluorescently labeled DNA-binding proteins for single-molecule experiments

Due to the enhanced labeling capability of maleimide-based fluorescent probes, lysine-cysteine-lysine (KCK) tags are frequently added to proteins for visualization. In this study, we employed an in vitro single-molecule DNA flow-stretching assay as a sensitive way to assess the impact of the KCK tag on the property of DNA-binding proteins. Using Bacillus subtilis ParB as an example, we show that, although no noticeable changes were detected by in vivo fluorescence imaging and chromatin immunoprecipitation (ChIP) assays, the KCK tag substantially altered ParB's DNA compaction rates and its response to nucleotide binding and to the presence of the specific sequence (parS) on the DNA. While it is typically assumed that short peptide tags minimally perturb protein function, our results urge researchers to carefully validate the use of tags for protein labeling. Our comprehensive analysis can be expanded and used as a guide to assess the impacts of other tags on DNA-binding proteins in single-molecule assays.

Genomic footprinting uncovers global transcription factor responses to amino acids in Escherichia coli

Our knowledge of transcriptional responses to changes in nutrient availability comes primarily from few well-studied transcription factors (TFs), often lacking an unbiased genome-wide perspective. Leveraging recent advances allowing bacterial genomic footprinting, we comprehensively mapped the genome-wide regulatory responses of Escherichia coli to exogenous leucine, methionine, alanine, and lysine. The global TF Lrp was found to individually sense three amino acids and mount three different target gene responses. Overall, 531 genes had altered RNA polymerase occupancy, and 32 TFs responded directly or indirectly to the presence of amino acids, including regulators of membrane and osmotic pressure homeostasis. About 70% of the detected TF-DNA interactions had not been reported before. We thus identified 682 previously unknown TF-binding locations, for a subset of which the involved TFs were identified by affinity purification. This comprehensive map of amino acid regulation illustrates the incompleteness of the known transcriptional regulation network, even in E. coli.

Brucella MucR acts as an H-NS-like protein to silence virulence genes and structure the nucleoid

Histone-like nucleoid structuring (H-NS) and H-NS-like proteins serve as global gene silencers and work with antagonistic transcriptional activators (counter-silencers) to properly coordinate the expression of virulence genes in pathogenic bacteria. In Brucella, MucR has been proposed as a novel H-NS-like gene silencer, but direct experimental evidence is lacking. Here, we show that MucR serves as an H-NS-like silencer of the Brucella abortus genes encoding the polar autotransporter adhesins BtaE and BmaC, the c-di-GMP-specific phosphodiesterase BpdB, and the quorum-sensing regulator BabR. We also demonstrate that the MarR-type transcriptional activator MdrA can displace MucR from the btaE promoter, supporting the existence of MucR counter-silencers in Brucella. Moreover, our chromatin immunoprecipitation (ChIP)-seq analysis identified 546 MucR enrichment peaks along the genome, including in the promoters of the genes encoding the Type IV secretion machinery and effectors and the quorum-sensing regulator VjbR. Importantly, MucR ChIP-seq peaks overlap with the previously described binding sites for the transcriptional activators VjbR, BvrR, and CtrA suggesting that these regulators serve as MucR counter-silencers and work in concert with MucR to coordinate virulence gene expression in Brucella. In addition, using chromosome conformation capture (Hi-C), we show that like H-NS in Escherichia coli, MucR alters the global structure of the Brucella nucleoid. Finally, a copy of the E. coli hns rescues the distinctive growth defect and elevated btaE expression of a B. abortus mucR mutant. Together, these findings solidify the role of MucR as a novel type of H-NS-like protein and suggest that MucR's gene-silencing properties play a key role in virulence in Brucella. IMPORTANCE Histone-like nucleoid structuring (H-NS) and H-NS-like proteins coordinate host-associated behaviors in many pathogenic bacteria, often through forming silencer/counter-silencer pairs with signal-responsive transcriptional activators to tightly control gene expression. Brucella and related bacteria do not encode H-NS or homologs of known H-NS-like proteins, and it is unclear if they have other proteins that perform analogous functions during pathogenesis. In this work, we provide compelling evidence for the role of MucR as a novel H-NS-like protein in Brucella. We show that MucR possesses many of the known functions attributed to H-NS and H-NS-like proteins, including the formation of silencer/counter-silencer pairs to control virulence gene expression and global structuring of the nucleoid. These results uncover a new role for MucR as a nucleoid structuring protein and support the importance of temporal control of gene expression in Brucella and related bacteria.

The two-component system WalKR provides an essential link between cell wall homeostasis and DNA replication in Staphylococcus aureus

Among the 16 two-component systems in the opportunistic human pathogen Staphylococcus aureus, only WalKR is essential. Like the orthologous systems in other Bacillota, S. aureus WalKR controls autolysins involved in peptidoglycan remodeling and is therefore intimately involved in cell division. However, despite the importance of WalKR in S. aureus, the basis for its essentiality is not understood and the regulon is poorly defined. Here, we defined a consensus WalR DNA-binding motif and the direct WalKR regulon by using functional genomics, including chromatin immunoprecipitation sequencing, with a panel of isogenic walKR mutants that had a spectrum of altered activities. Consistent with prior findings, the direct regulon includes multiple autolysin genes. However, this work also revealed that WalR directly regulates at least five essential genes involved in lipoteichoic acid synthesis (ltaS): translation (rplK), DNA compaction (hup), initiation of DNA replication (dnaA, hup) and purine nucleotide metabolism (prs). Thus, WalKR in S. aureus serves as a polyfunctional regulator that contributes to fundamental control over critical cell processes by coordinately linking cell wall homeostasis with purine biosynthesis, protein biosynthesis, and DNA replication. Our findings further address the essentiality of this locus and highlight the importance of WalKR as a bona fide target for novel anti-staphylococcal therapeutics. IMPORTANCE The opportunistic human pathogen Staphylococcus aureus uses an array of protein sensing systems called two-component systems (TCS) to sense environmental signals and adapt its physiology in response by regulating different genes. This sensory network is key to S. aureus versatility and success as a pathogen. Here, we reveal for the first time the full extent of the regulatory network of WalKR, the only staphylococcal TCS that is indispensable for survival under laboratory conditions. We found that WalKR is a master regulator of cell growth, coordinating the expression of genes from multiple, fundamental S. aureus cellular processes, including those involved in maintaining cell wall metabolism, protein biosynthesis, nucleotide metabolism, and the initiation of DNA replication.

Iron Oxide Nanoparticles with and without Cobalt Functionalization Provoke Changes in the Transcription Profile via Epigenetic Modulation of Enhancer Activity

Despite the progress in the field of nanotoxicology, much about the cellular mechanisms that mediate the adverse effects of nanoparticles (NPs) and, in particular, the possible role of epigenetics in nanotoxicity, remains to be clarified. Therefore, we studied the changes occurring in the genome-wide distribution of H3K27ac, H3K4me1, H3K9me2, and H3K27me3 histone modifications and compared them with the transcriptome after exposing NIH3T3 cells to iron-based magnetic NPs (i.e., Fe2O3 and Fe2O3@Co NPs). We found that the transcription response is mainly due to changes in the genomic distribution of H3K27ac that can modulate the activity of enhancers. We propose that alteration of the epigenetic landscape is a key mechanism in defining the gene expression program changes resulting in nanotoxicity. With this approach, it is possible to construct a data set of genomic regions that could be useful for defining toxicity in a manner that is more comprehensive than what is possible with the present toxicology assays.

Transcription factor-binding k-mer analysis clarifies the cell type dependency of binding specificities and cis-regulatory SNPs in humans

BACKGROUND

Transcription factors (TFs) exhibit heterogeneous DNA-binding specificities in individual cells and whole organisms under natural conditions, and de novo motif discovery usually provides multiple motifs, even from a single chromatin immunoprecipitation-sequencing (ChIP-seq) sample. Despite the accumulation of ChIP-seq data and ChIP-seq-derived motifs, the diversity of DNA-binding specificities across different TFs and cell types remains largely unexplored.

RESULTS

Here, we applied MOCCS2, our k-mer-based motif discovery method, to a collection of human TF ChIP-seq samples across diverse TFs and cell types, and systematically computed profiles of TF-binding specificity scores for all k-mers. After quality control, we compiled a set of TF-binding specificity score profiles for 2,976 high-quality ChIP-seq samples, comprising 473 TFs and 398 cell types. Using these high-quality samples, we confirmed that the k-mer-based TF-binding specificity profiles reflected TF- or TF-family dependent DNA-binding specificities. We then compared the binding specificity scores of ChIP-seq samples with the same TFs but with different cell type classes and found that half of the analyzed TFs exhibited differences in DNA-binding specificities across cell type classes. Additionally, we devised a method to detect differentially bound k-mers between two ChIP-seq samples and detected k-mers exhibiting statistically significant differences in binding specificity scores. Moreover, we demonstrated that differences in the binding specificity scores between k-mers on the reference and alternative alleles could be used to predict the effect of variants on TF binding, as validated by in vitro and in vivo assay datasets. Finally, we demonstrated that binding specificity score differences can be used to interpret disease-associated non-coding single-nucleotide polymorphisms (SNPs) as TF-affecting SNPs and provide candidates responsible for TFs and cell types.

CONCLUSIONS

Our study provides a basis for investigating the regulation of gene expression in a TF-, TF family-, or cell-type-dependent manner. Furthermore, our differential analysis of binding-specificity scores highlights noncoding disease-associated variants in humans.

ChIP-seq identifies McSLC35E2 as a novel target gene of McNrf2 in Mytilus coruscus, highlighting its role in the regulation of oxidative stress response in marine mollusks

NF-E2-related factor 2 (Nrf2) plays a crucial role in the oxidative regulatory process, which could trigger hundreds of antioxidant elements to confront xenobiotics. In the previous study, we identified Nrf2 from the marine mussel Mytilus coruscus, and the findings demonstrated that McNrf2 effectively protected the mussels against oxidative stress induced by benzopyrene (Bap). In order to delve deeper into the underlying mechanism, we utilized Chromatin Immunoprecipitation followed by sequencing (ChIP-seq) technology to systematically identify potential novel target genes of McNrf2. A total of 3,465 potential target genes were screened, of which 219 owned binding sites located within the promoter region. During subsequent experimental verification, it was found that McSLC35E2, a candidate target gene of McNrf2, exhibited negative regulation by McNrf2, as confirmed through dual luciferase and qRT-PCR detection. Further, the enzyme activity tests demonstrated that McNrf2 could counteract Bap induced oxidative stress by inhibiting McSLC35E2. The current study provides valuable insights into the application of ChIP-seq technology in the research of marine mollusks, advancing our understanding of the key role of Nrf2 in antioxidant defense mechanisms, and highlighting the significance of SLC35E2 in the highly sophisticated regulation of oxidative stress response in marine invertebrates.

Runt-related Transcription Factors and Gene Regulatory Mechanisms in Skeletal Development and Diseases

PURPOSE OF REVIEW

Runt-related transcription factors (RUNX) play critical roles in skeletal development, metabolism, and diseases. In mammals, three RUNX members, namely RUNX1, RUNX2, and RUNX3, play distinct and redundant roles, although RUNX2 is a dominant factor in skeletal development and several skeletal diseases. This review is to provide an overview of the current understanding of RUNX-mediated transcriptional regulation in different skeletal cell types.

RECENT FINDINGS

Advances in chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) have revealed genome-wide RUNX-mediated gene regulatory mechanisms, including their association with cis-regulatory elements and putative target genes. Further studies with genome-wide analysis and biochemical assays have shed light on RUNX-mediated pioneering action and involvements of RUNX2 in lipid-lipid phase separation. Emerging multi-layered mechanisms of RUNX-mediated gene regulations help us better understanding of skeletal development and diseases, which also provides clues to think how genome-wide studies can help develop therapeutic strategies for skeletal diseases.

Ligand Activation of the Aryl Hydrocarbon Receptor Upregulates Epidermal Uridine Diphosphate Glucose Ceramide Glucosyltransferase and Glucosylceramides

Ligand activation of the aryl hydrocarbon receptor (AHR) accelerates keratinocyte differentiation and the formation of the epidermal permeability barrier. Several classes of lipids, including ceramides, are critical to the epidermal permeability barrier. In normal human epidermal keratinocytes, the AHR ligand, 2,3,7,8-tetrachlorodibenzo-p-dioxin, increased RNA levels of ceramide metabolism and transport genes: uridine diphosphate glucose ceramide glucosyltransferase (UGCG), ABCA12, GBA1, and SMPD1. Levels of abundant skin ceramides were also increased by 2,3,7,8-tetrachlorodibenzo-p-dioxin. These included the metabolites synthesized by UGCG, glucosylceramides, and acyl glucosylceramides. Chromatin immunoprecipitation-sequence analysis and luciferase reporter assays identified UGCG as a direct AHR target. The AHR antagonist, GNF351, inhibited the 2,3,7,8-tetrachlorodibenzo-p-dioxin-mediated RNA and transcriptional increases. Tapinarof, an AHR ligand approved for the treatment of psoriasis, increased UGCG RNA, protein, and its lipid metabolites hexosylceramides as well as increased the RNA expression of ABCA12, GBA1, and SMPD1. In Ahr-null mice, Ugcg RNA and hexosylceramides were lower than those in the wild type. These results indicate that the AHR regulates the expression of UGCG, a ceramide-metabolizing enzyme required for ceramide trafficking, keratinocyte differentiation, and epidermal permeability barrier formation.

Exploring the Genomic Landscape: An In-depth ChIP-seq Analysis Protocol for Uncovering Protein-DNA Interactions

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a widely employed technique for investigating protein-DNA interactions. However, the absence of a standardized and clear workflow necessitates researchers to independently assemble methodologies from diverse resources. This lack of uniformity hampers reproducibility and makes version control a complex endeavor, thereby limiting the accessibility of ChIP-seq analyses to individuals with extensive training in bioinformatics. In light of these challenges, we have developed an executable protocol that addresses these limitations. Our protocol encompasses all aspects of ChIP-seq analysis, ranging from quality control of raw reads to peak calling and downstream functional analyses. We have implemented two distinct approaches for peak calling, providing researchers with flexibility to choose the most suitable method for their specific experimental needs. This protocol will contribute to the scientific community by providing a standardized and clear resource that will enhance the reproducibility and accessibility of ChIP-seq analyses. © 2023 Wiley Periodicals LLC. Basic Protocol: ChIP-seq analysis workflow Alternative Protocol: Call differentially enriched peaks by using MACS3.

Integrative analyses of RNA-seq and ChIP-seq Reveal MITF as a Target Gene of TFPI-2 in MDA231 Cells

Breast cancer is the most prevalent cancer in female patients worldwide. Tissue factor pathway inhibitor 2 (TFPI-2) is identified as an important tumor suppressor in various cancers. Recent studies have shown that TFPI-2 translocates into the nucleus, where it modulates the transcription of the matrix metalloproteinase-2 (MMP-2) gene. However, its biological role and molecular mechanisms in the progression of breast cancer remain unclear. In this study, we identified 5125 differentially expressed genes (DEGs) from RNA sequencing (RNA-seq) in TFPI-2-overexpressing MDA231 cells compared with control cells. Gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) analysis shown that cell cycle, cell differentiation, proteoglycans in cancer, and pathways associated with cancer were highly enriched in downregulated DEGs. Integration of the RNA-seq and ChIP-sequencing (ChIP-seq) data identified 73 genes directly controlled by TFPI-2 in MDA231 cells. Among them, melanocyte inducing transcription factor (MITF) gene expression was repressed by TFPI-2, which was further verified by a luciferase reporter assay and ChIP-quantitative PCR. Our study provides evidence of a novel role of TFPI-2 in human breast cancer involving targeting of the MITF.

A Micro-evolutionary Change in Target Binding Sites as a Key Determinant of Ultrabithorax Function in Drosophila

Hox genes encode Homeodomain-containing transcription factors, which specify segmental identities along the anterior-posterior axis. Functional changes in Hox genes have been directly implicated in the evolution of body plans across the metazoan lineage. The Hox protein Ultrabithorax (Ubx) is expressed and required in developing third thoracic (T3) segments in holometabolous insects studied so far, particularly, of the order Coleoptera, Lepidoptera and Diptera. Ubx function is key to specify differential development of the second (T2) and T3 thoracic segments in these insects. While Ubx is expressed in the third thoracic segment in developing larvae of Hymenopteran Apis mellifera, the morphological differences between T2 and T3 are subtle. To identify evolutionary changes that are behind the differential function of Ubx in Drosophila and Apis, which are diverged for more than 350 million years, we performed comparative analyses of genome wide Ubx-binding sites between these two insects. Our studies reveal that a motif with a TAAAT core is a preferred binding site for Ubx in Drosophila, but not in Apis. Biochemical and transgenic assays suggest that in Drosophila, the TAAAT core sequence in the Ubx binding sites is required for Ubx-mediated regulation of two of its target genes studied here; CG13222, a gene that is normally upregulated by Ubx and vestigial (vg), whose expression is repressed by Ubx in T3. Interestingly, changing the TAAT site to a TAAAT site was sufficient to bring an otherwise unresponsive enhancer of the vg gene from Apis under the control of Ubx in a Drosophila transgenic assay. Taken together, our results suggest an evolutionary mechanism by which critical wing patterning genes might have come under the regulation of Ubx in the Dipteran lineage.