Search

Insects use seasonal diapause as an alternative strategy to endure adverse seasons. This developmental trajectory is induced by environmental cues like short‐day lengths in late summer and early fall, but how insects measure day length is unknown. The circadian clock has been implicated in regulating photoperiodic or seasonal responses in many insects, including the Northern house mosquito, Culex pipiens, which enters adult diapause. To investigate the potential control of diapause by circadian control, we employed ChIP‐sequencing to identify the downstream targets of a circadian transcription factor, PAR domain protein 1 (PDP1), that contribute to the hallmark features of diapause. We identified the nearest genes in a 10 kb region of the anticipated PDP1 binding sites, listed prospective targets and searched for PDP1‐specific binding sites. By examining the functional relevance to diapause‐specific behaviours and modifications such as metabolic pathways, lifespan extension, cell cycle regulation and stress tolerance, eight genes were selected as targets and validated using ChIP‐qPCR. In addition, qRT‐PCR demonstrated that the mRNA abundance of PDP1 targets increased in the heads of diapausing females during the middle of the scotophase (ZT17) compared with the early photophase (ZT1), in agreement with the peak and trough of PDP1 abundance. Thus, our investigation uncovered the mechanism by which PDP1 might generate a diapause phenotype in insects. [ABSTRACT FROM AUTHOR]

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. [Display omitted] • Activity-based protein profiling identifies covalent ligands for FOXA1 • Covalent ligands react with cysteine-258 of FOXA1 in a DNA-dependent manner • Covalent ligands remodel the pioneering activity of FOXA1 in prostate cancer cells • Covalent ligands relax the canonical DNA-binding preferences of FOXA1 Won et al. present the chemical proteomic discovery of covalent small molecules that stereoselectively and site-specifically engage the pioneer transcription factor FOXA1. These compounds rapidly remodel FOXA1 interactions with chromatin in prostate cancer cells and create corresponding changes in chromatin accessibility through relaxing the DNA-binding preferences of FOXA1. [ABSTRACT FROM AUTHOR]

Deciphering the underlying gene regulatory networks (GRNs) that govern early human embryogenesis is critical for understanding developmental mechanisms yet remains challenging due to limited sample availability and the inherent complexity of the biological processes involved. To address this, we developed InPheRNo-ChIP, a computational framework that integrates multimodal data, including RNA-seq, transcription factor (TF)–specific ChIP-seq, and phenotypic labels, to reconstruct phenotype-relevant GRNs associated with endoderm development. The core of this method is a probabilistic graphical model that models the simultaneous effect of TFs on their putative target genes to influence a particular phenotypic outcome. Unlike the majority of existing GRN inference methods that are agnostic to the phenotypic outcomes, InPheRNo-ChIP directly incorporates phenotypic information during GRN inference, enabling the distinction between lineage-specific and general regulatory interactions. We integrated data from three experimental studies and applied InPheRNo-ChIP to infer the GRN governing the differentiation of human embryonic stem cells into definitive endoderm. Benchmarking against a scRNA-seq CRISPRi study demonstrated InPheRNo-ChIP's ability to identify regulatory interactions involving endoderm markers FOXA2 , SMAD2 , and SOX17 , outperforming other methods. This highlights the importance of incorporating the phenotypic context during network inference. Furthermore, an ablation study confirms the synergistic contribution of ChIP-seq, RNA-seq, and phenotypic data, highlighting the value of multimodal integration for accurate phenotype-relevant GRN reconstruction. [ABSTRACT FROM AUTHOR]

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. [Display omitted] • 148 of the 150 examined C2H2-ZNFs bind directly to RNA in cells • Denaturing IP conditions do not abolish C2H2-ZNFs' RNA binding • iCLIP-seq for 94 RBPs, including 51 C2H2-ZNFs, reveals distinct binding preferences • C2H2-ZNFs function in RNA-related processes ranging from splicing to m6A regulation Nabeel-Shah et al. report that transcription factors of the C2H2-ZNF family also function as RNA-binding proteins that participate in multiple post-transcriptional processes, including alternative splicing, alternative polyadenylation, and m6A regulation. Integration of C2H2-ZNF iCLIP-seq data with other omics data reveals multiple modes of RNA binding by C2H2-ZNFs and uncovers mechanistic details. [ABSTRACT FROM AUTHOR]

Spatiotemporal regulation of gene expression is controlled by transcription factor (TF) binding to regulatory elements, resulting in a plethora of cell types and cell states from the same genetic information. Due to the importance of regulatory elements, various sequencing methods have been developed to localise them in genomes, for example using ChIP-seq profiling of the histone mark H3K27ac that marks active regulatory regions. Moreover, multiple tools have been developed to predict TF binding to these regulatory elements based on DNA sequence. As altered gene expression is a hallmark of disease phenotypes, identifying TFs driving such gene expression programs is critical for the identification of novel drug targets.In this study, we curated 84 chromatin profiling experiments (H3K27ac ChIP-seq) where TFs were perturbed through e.g., genetic knockout or overexpression. We ran nine published tools to prioritize TFs using these real-world datasets and evaluated the performance of the methods in identifying the perturbed TFs. This allowed the nomination of three frontrunner tools, namely RcisTarget, MEIRLOP and monaLisa. Our analyses revealed opportunities and commonalities of tools that will help to guide further improvements and developments in the field.

Abstract 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. Graphical Abstract

Abstract 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.

Objective Rare study of the non-coding and regulatory regions of the genome limits our ability to decode the mechanisms of fatty liver hemorrhage syndrome (FLHS) in chickens. Methods Herein, we constructed the high-fat diet-induced FLHS chicken model to investigate the genome-wide active enhancers and transcriptome by H3K27ac target chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-Seq) profiles of normal and FLHS liver tissues. Concurrently, an integrative analysis combining ChIP-seq with RNA-Seq and a comparative analysis with chicken FLHS, rat non-alcoholic fatty liver disease (NAFLD) and human NAFLD at the transcriptome level revealed the enhancer and super enhancer target genes and conservative genes involved in metabolic processes. Results In total, 56 and 199 peak-genes were identified in upregulated peak-genes positively regulated by H3K27ac (Cor (peak-gene correlation) ≥0.5 and log2(FoldChange) ≥1) (PP) and downregulated peak-genes positively regulated by H3K27ac (Cor (peak-gene correlation) ≥0.5 and log2(FoldChange)≤−1) (PN), respectively; then we screened key regulatory targets mainly distributing in lipid metabolism (PCK1, APOA4, APOA1, INHBE) and apoptosis (KIT, NTRK2) together with MAPK and PPAR signaling pathway in FLHS. Intriguingly, PCK1 was also significantly covered in up-regulated super-enhancers (SEs), which further implied the vital role of PCK1 during the development of FLHS. Conclusion Together, our studies have identified potential therapeutic biomarkers of PCK1 and elucidated novel insights into the pathogenesis of FLHS, especially for the epigenetic perspective.

Objective To analyze the genome-wide profiles of H3K9 crotonylation (H3K9cr) in mouse embryonic stem cells (mESCs) in normal or deficient folic acid (FA) by chromatin immunoprecipitation sequencing (ChIP-seq) in order to observe the genome-wide changes in H3K9cr status and the regulated pathway genes under FA deficiency condition. Methods The mESCs were divided into normal FA group (4 mg/L, FA4) and FA deficient group (0 mg/L, FA0). Western blotting and immunohistochemical (IHC) assay were used to detect the protein level of H3K9cr in the FA deficient mESCs and the embryonic brain tissues of neural tube defects (NTDs) mice, respectively. ChIP-seq technology was carried out for sequence identification of the obtained specific binding DNA fragments. Differentially expressed genes (DEGs) were analyzed and visualized by gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathways, with the aim of analyzing the whole genome of mESCs with FA deficiency and validating the regulated genes. Results Overall H3K9cr level was significantly lower in the FA deficient mESCs and the embryonic brain tissues of NTDs mice than in the corresponding normal FA mESCs and mice (P < 0.05). Neurodevelopmental-related genes such as Bdnf, Pax6, and App were down-regulated in the FA deficient mESCs (P < 0.05). The KEGG pathway analysis indicated that the DEGs were concentrated in axon-guiding, neuroactive ligand-receptor interaction, glutaminergic synapses and other nervous system-related pathways. Conclusion FA deficiency in embryonic stem cells participates in the regulation of neurodevelopmental genes through modification of H3K9cr, suggesting that the modification of histone crotonylation may be involved in neurodevelopmental diseases induced by FA deficiency.

Summary 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Neurons are central to lifelong learning and memory, but ageing disrupts their morphology and function, leading to cognitive decline. Although epigenetic mechanisms are known to play crucial roles in learning and memory, neuron-specific genome-wide epigenetic maps into old age remain scarce, often being limited to whole-brain homogenates and confounded by glial cells. Here, we mapped H3K4me3, H3K27ac, and H3K27me3 in mouse neurons across their lifespan. This revealed stable H3K4me3 and global losses of H3K27ac and H3K27me3 into old age. We observed patterns of synaptic function gene deactivation, regulated through the loss of the active mark H3K27ac, but not H3K4me3. Alongside this, embryonic development loci lost repressive H3K27me3 in old age. This suggests a loss of a highly refined neuronal cellular identity linked to global chromatin reconfiguration. Collectively, these findings indicate a key role for epigenetic regulation in neurons that is inextricably linked with ageing. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Thiamine (vitamin B1) is essential for glucose catabolism. In the yeast species, Nakaseomyces glabratus (formerly Candida glabrata) and Saccharomyces cerevisiae , the transcription factor Pdc2 (with Thi3 and Thi2) upregulates pyruvate decarboxylase (PDC) genes and thiamine biosynthetic and acquisition (THI) genes during starvation. There have not been genome-wide analyses of Pdc2 binding. Previously, we identified small regions of Pdc2-regulated genes sufficient to confer thiamine regulation. Here, we performed deletion analyses on these regions. We observed that when the S. cerevisiae PDC5 promoter is introduced into N. glabratus , it is thiamine starvation inducible but does not require the Thi3 coregulator. The ScPDC5 promoter contains a 22-bp duplication with an AT-rich spacer between the 2 repeats, which are important for regulation. Loss of the first 22-bp element does not eliminate regulation, but the promoter becomes Thi3 dependent, suggesting cis architecture can generate a Thi3-independent, thiamine starvation inducible response. Whereas many THI promoters only have 1 copy of this element, addition of the first 22-bp element to a Thi3-dependent promoter confers Thi3 independence. Finally, we performed fluorescence anisotropy and chromatin immunoprecipitation sequencing. Pdc2 and Thi3 bind to regions that share similarity to the 22-bp element in the ScPDC5 promoter and previously identified cis elements in N. glabratus promoters. Also, while Pdc2 binds to THI and PDC promoters, neither Pdc2 nor Thi3 appears to bind the evolutionarily new NgPMU3 promoter that is regulated by Pdc2. Further study is warranted because PMU3 is required for cells to acquire thiamine from environments where thiamine is phosphorylated, such as in the human bloodstream. [ABSTRACT FROM AUTHOR]

The cell-adhesion molecule NEPH1 is required for maintaining the structural integrity and function of the glomerulus in the kidneys. In the nervous system of Drosophila and C. elegans, it is involved in synaptogenesis and axon branching, which are essential for establishing functional circuits. In the mammalian nervous system, the expression regulation and function of Neph1 has barely been explored. In this study, we provide a spatiotemporal characterization of Neph1 expression in mouse dorsal root ganglia (DRGs) and spinal cord. After the neurogenic phase, Neph1 is broadly expressed in the DRGs and in their putative targets at the dorsal horn of the spinal cord, comprising both GABAergic and glutamatergic neurons. Interestingly, we found that PRRXL1, a homeodomain transcription factor that is required for proper establishment of the DRG-spinal cord circuit, prevents a premature expression of Neph1 in the superficial laminae of the dorsal spinal cord at E14.5, but has no regulatory effect on the DRGs or on either structure at E16.5. By chromatin immunoprecipitation analysis of the dorsal spinal cord, we identified four PRRXL1-bound regions within the Neph1 introns, suggesting that PRRXL1 directly regulates Neph1 transcription. We also showed that Neph1 is required for branching, especially at distal neurites. Together, our work showed that Prrxl1 prevents the early expression of Neph1 in the superficial dorsal horn, suggesting that Neph1 might function as a downstream effector gene for proper assembly of the DRG-spinal nociceptive circuit. [ABSTRACT FROM AUTHOR]

Simple Summary: The emergence of next-generation sequencing (NGS) signified a revolution in biology research by capturing the significance of DNA loci and RNA molecules on a genome-wide scale. The complexity and volume of NGS data highlight the need for robust and user-friendly computational tools. The framework presented here, aPEAch, is an automated pipeline for end-to-end analysis of DNA and RNA sequencing assays, including small RNA sequencing. Implemented in Python, it allows users to customize the analysis and results, handling single or multiple replicates in batches, while also automating advanced unsupervised learning analyses. aPEAch offers quality control reports, fragment size distribution plots, and intermediate files supporting reproducibility and interoperability, along with publication-ready visualizations. With the advent of next-generation sequencing (NGS), experimental techniques that capture the biological significance of DNA loci or RNA molecules have emerged as fundamental tools for studying the epigenome and transcriptional regulation on a genome-wide scale. The volume of the generated data and the underlying complexity regarding their analysis highlight the need for robust and easy-to-use computational analytic methods that can streamline the process and provide valuable biological insights. Our solution, aPEAch, is an automated pipeline that facilitates the end-to-end analysis of both DNA- and RNA-sequencing assays, including small RNA sequencing, from assessing the quality of the input sample files to answering meaningful biological questions by exploiting the rich information embedded in biological data. Our method is implemented in Python, based on a modular approach that enables users to choose the path and extent of the analysis and the representations of the results. The pipeline can process samples with single or multiple replicates in batches, allowing the ease of use and reproducibility of the analysis across all samples. aPEAch provides a variety of sample metrics such as quality control reports, fragment size distribution plots, and all intermediate output files, enabling the pipeline to be re-executed with different parameters or algorithms, along with the publication-ready visualization of the results. Furthermore, aPEAch seamlessly incorporates advanced unsupervised learning analyses by automating clustering optimization and visualization, thus providing invaluable insight into the underlying biological mechanisms. [ABSTRACT FROM AUTHOR]

The epigenetic reader SntB was identified as an important transcriptional regulator of growth, development, and secondary metabolite synthesis in Aspergillus flavus. However, the underlying molecular mechanism is still unclear. In this study, by gene deletion and complementation, we found SntB is essential for mycelia growth, conidial production, sclerotia formation, aflatoxin synthesis, and host colonization. Chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) analysis revealed that SntB played key roles in oxidative stress response of A. flavus, influencing related gene activity, especially catC encoding catalase. SntB regulated the expression activity of catC with or without oxidative stress, and was related to the expression level of the secretory lipase (G4B84_008359). The deletion of catC showed that CatC participated in the regulation of fungal morphogenesis, reactive oxygen species (ROS) level, and aflatoxin production, and that CatC significantly regulated fungal sensitive reaction and AFB1 yield under oxidative stress. Our study revealed the potential machinery that SntB regulated fungal morphogenesis, mycotoxin anabolism, and fungal virulence through the axle of from H3K36me3 modification to fungal virulence and mycotoxin biosynthesis. The results of this study shed light into the SntB-mediated transcript regulation pathways of fungal mycotoxin anabolism and virulence, which provided potential strategy to control the contamination of A. flavus and its aflatoxins.

Abstract Background Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous (“multimappers”). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. Results In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. Conclusion Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines –currently restricted to specific contexts or communities– to ensure the reliability of genomic and transcriptomic studies.

Abstract 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.

Abstract 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.

Using machine learning (ML), we interrogated the function of all human-chimpanzee variants in 2,645 human accelerated regions (HARs), finding 43% of HARs have variants with large opposing effects on chromatin state and 14% on neurodevelopmental enhancer activity. This pattern, consistent with compensatory evolution, was confirmed using massively parallel reporter assays in chimpanzee and human neural progenitor cells. The species-specific enhancer activity of HARs was accurately predicted from the presence and absence of transcription factor footprints in each species. Despite these striking cis effects, activity of a given HAR sequence was nearly identical in human and chimpanzee cells. This suggests that HARs did not evolve to compensate for changes in the trans environment but instead altered their ability to bind factors present in both species. Thus, ML prioritized variants with functional effects on human neurodevelopment and revealed an unexpected reason why HARs may have evolved so rapidly.

Summary: Rice blast, caused by Magnaporthe oryzae, significantly impacts grain yield, necessitating the identification of broad‐spectrum resistance genes and their functional mechanisms for disease‐resistant crop breeding. Here, we report that rice with knockdown OsHDAC1 gene expression displays enhanced broad‐spectrum blast resistance without effects on plant height and tiller numbers compared to wild‐type rice, while rice overexpressing OsHDAC1 is more susceptible to M. oryzae. We identify a novel blast resistance transcription factor, OsGRAS30, which genetically acts upstream of OsHDAC1 and interacts with OsHDAC1 to suppress its enzymatic activity. This inhibition increases the histone H3K27ac level, thereby boosting broad‐spectrum blast resistance. Integrating genome‐wide mapping of OsHDAC1 and H3K27ac targets with RNA sequencing analysis unveils how OsHDAC1 mediates the expression of OsSSI2, OsF3H, OsRLR1 and OsRGA5 to regulate blast resistance. Our findings reveal that the OsGRAS30–OsHDAC1 module is critical to rice blast control. Therefore, targeting either OsHDAC1 or OsGRAS30 offers a promising approach for enhancing crop blast resistance. [ABSTRACT FROM AUTHOR]

Background: Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored. Results: In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified. Conclusion: Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines –currently restricted to specific contexts or communities– to ensure the reliability of genomic and transcriptomic studies. [ABSTRACT FROM AUTHOR]

Summary: 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. [ABSTRACT FROM AUTHOR]

To identify the inactive genes in cataract lenses and explore their function in lens epithelial cells (LECs). 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. 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. 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. [ABSTRACT FROM AUTHOR]

Background: Cervical cancer is among the highest-ranking types of cancer worldwide, with human papillomavirus (HPV) as the agent driving the malignant process. One aspect of the infection's evolution is given by epigenetic modifications, mainly DNA methylation and chromatin alteration. These processes are guided by several chromatin remodeling complexes, including NuRD. The purpose of this study was to evaluate the genome-wide binding patterns of the NuRD complex components (MBD2 and MBD3) in the presence of active HPV16 E6 and E7 oncogenes and to determine the potential of identified genes through an experimental model to differentiate between cervical precursor lesions, with the aim of establishing their utility as biomarkers. Methods: The experimental model was built using the CaSki cell line and shRNA for E6 and E7 HPV16 silencing, ChIP-seq, qRT-PCR, and Western blot analyses. Selected genes' expression was also assessed in patients. Results: Several genes have been identified to exhibit altered transcriptional activity due to the influence of HPV16 E6/E7 viral oncogenes acting through the MBD2/MBD3 NuRD complex, linking them to viral infection and cervical oncogenesis. Conclusions: The impacted genes primarily play roles in governing gene transcription, mRNA processing, and regulation of translation. Understanding these mechanisms offers valuable insights into the process of HPV-induced oncogenesis. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

ABSTRACT The tumor suppressor p53, primarily functioning as a transcription factor, has exhibited antiviral capabilities against various viruses in chickens, including infectious bursal disease virus (IBDV), avian leukosis virus subgroup J (ALV-J), and avian infectious laryngotracheitis virus (ILTV). Nevertheless, the existence of a universal antiviral mechanism employed by chicken p53 (chp53) against these viruses remains uncertain. This study conducted a comprehensive comparison of molecular networks involved in chp53’s antiviral function against IBDV, ALV-J, and ILTV. This was achieved through an integrated analysis of ChIP-seq data, examining chp53’s genome-wide chromatin occupancy, and RNA-seq data from chicken cells infected with these viruses. The consistent observation of chp53 target gene enrichment in metabolic pathways, confirmed via ChIP-qPCR, suggests a ubiquitous regulation of host cellular metabolism by chp53 across different viruses. Further genome binding motif conservation analysis and transcriptional co-factor prediction suggest conserved transcriptional regulation mechanism by which chp53 regulates host cellular metabolism during viral infection. These findings offer novel insights into the antiviral role of chp53 and propose that targeting the virus-host metabolic interaction through regulating p53 could serve as a universal strategy for antiviral therapies in chickens.IMPORTANCEThe current study conducted a comprehensive analysis, comparing molecular networks underlying chp53’s antiviral role against infectious bursal disease virus (IBDV), avian leukosis virus subgroup J (ALV-J), and avian infectious laryngotracheitis virus (ILTV). This was achieved through a combined assessment of ChIP-seq and RNA-seq data obtained from infected chicken cells. Notably, enrichment of chp53 target genes in metabolic pathways was consistently observed across viral infections, indicating a universal role of chp53 in regulating cellular metabolism during diverse viral infections. These findings offer novel insights into the antiviral capabilities of chicken p53, laying a foundation for the potential development of broad-spectrum antiviral therapies in chickens.

Abstract Background Discovery biological motifs plays a fundamental role in understanding regulatory mechanisms. Computationally, they can be efficiently represented as kmers, making the counting of these elements a critical aspect for ensuring not only the accuracy but also the efficiency of the analytical process. This is particularly useful in scenarios involving large data volumes, such as those generated by the ChIP-seq protocol. Against this backdrop, we introduce biomapp::chip, a tool specifically designed to optimize the discovery of biological motifs in large data volumes. Results We conducted a comprehensive set of comparative tests with state-of-the-art algorithms. Our analyses revealed that biomapp::chip outperforms existing approaches in various metrics, excelling both in terms of performance and accuracy. The tests demonstrated a higher detection rate of significant motifs and also greater agility in the execution of the algorithm. Furthermore, the smt component played a vital role in the system’s efficiency, proving to be both agile and accurate in kmer counting, which in turn improved the overall efficacy of our tool. Conclusion biomapp::chip represent real advancements in the discovery of biological motifs, particularly in large data volume scenarios, offering a relevant alternative for the analysis of ChIP-seq data and have the potential to boost future research in the field. This software can be found at the following address: (https://github.com/jadermcg/biomapp-chip).

Abstract 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.

Abstract Background Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress. Results In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P

Avian pathogenic Escherichia coli (APEC) belonging to extraintestinal pathogenic E. coli (ExPEC) can cause severe infections in extraintestinal tissues in birds and humans, such as the lungs and blood. MprA (microcin production regulation, locus A, herein renamed AbsR, a blood survival regulator), a member of the MarR (multiple antibiotic resistance regulator) transcriptional regulator family, governs the expression of capsule biosynthetic genes in human ExPEC and represents a promising druggable target for antimicrobials. However, a deep understanding of the AbsR regulatory mechanism as well as its regulon is lacking. In this study, we present a systems-level analysis of the APEC AbsR regulon using ChIP-Seq (chromatin immunoprecipitation sequencing) and RNA-Seq (RNA sequencing) methods. We found that AbsR directly regulates 99 genes and indirectly regulates 667 genes. Furthermore, we showed that: 1) AbsR contributes to antiphagocytotic effects by macrophages and virulence in a mouse model for systemic infection by directly activating the capsular gene cluster; 2) AbsR positively impacts biofilm formation via direct regulation of the T2SS (type II secretion system) but plays a marginal role in virulence; and 3) AbsR directly upregulates the acid tolerance signaling system EvgAS to withstand acid stress but is dispensable in ExPEC virulence. Finally, our data indicate that the role of AbsR in virulence gene regulation is relatively conserved in ExPEC strains. Altogether, this study provides a comprehensive analysis of the AbsR regulon and regulatory mechanism, and our data suggest that AbsR likely influences virulence primarily through the control of capsule production. Interestingly, we found that AbsR severely represses the expression of the type I-F CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR associated) systems, which could have implications in CRISPR biology and application.

Embryonic genome activation (EGA) occurs during preimplantation development and is characterized by the initiation of de novo transcription from the embryonic genome. Despite its importance, the regulation of EGA and the transcription factors involved in this process are poorly understood. Paired-like homeobox (PRDL) family proteins are implicated as potential transcriptional regulators of EGA, yet the PRDL-mediated gene regulatory networks remain uncharacterized. To investigate the function of PRDL proteins, we are identifying the molecular interactions and the functions of a subset family of the Eutherian Totipotent Cell Homeobox (ETCHbox) proteins, seven PRDL family proteins and six other transcription factors (TFs), all suggested to participate in transcriptional regulation during preimplantation. Using mass spectrometry-based interactomics methods, AP-MS and proximity-dependent biotin labeling, and chromatin immunoprecipitation sequencing we derive the comprehensive regulatory networks of these preimplantation TFs. By these interactomics tools we identify more than a thousand high-confidence interactions for the 21 studied bait proteins with more than 300 interacting proteins. We also establish that TPRX2, currently assigned as pseudogene, is a transcriptional activator. Synopsis: This study characterizes the molecular interactions and functions of human PRDL and other transcription factors in embryonic genome activation (EGA), revealing over a thousand high-confidence interactions and indicating key roles in transcription regulation. PRDL Family in EGA: The data support a crucial role for PRDL family proteins in embryonic genome activation, contributing to our understanding of the regulation of EGA. Interactomics Analysis: The study employs advanced mass spectrometry-based interactomics methods, identifying over a thousand high-confidence interactions involving PRDL proteins, thus mapping a comprehensive network of preimplantation transcription factors. This study characterizes the molecular interactions and functions of human PRDL and other transcription factors in embryonic genome activation (EGA), revealing over a thousand high-confidence interactions and indicating key roles in transcription regulation. [ABSTRACT FROM AUTHOR]

Simple Summary: CIC-DUX4-rearranged sarcoma (CDS) is a rare and aggressive soft tissue tumor predominantly affecting young adults for which the only oncogenic driver is the CIC-DUX4 fusion protein, resulting from chromosomal rearrangements. The molecular mechanisms by which CIC-DUX4 drives CDS development remain largely unknown, hindering the development of targeted therapies. This study conducted genome-wide profiling of CIC-DUX4 DNA occupancy and associated chromatin states in CDS cell models and tumors. Our findings reveal CIC-DUX4 as a potent transcriptional activator at its binding sites, facilitated by direct interaction with the acetyltransferase p300. Significantly, inhibiting p300 impedes CDS tumor cell proliferation, suggesting a potential therapeutic avenue. This research enhances understanding of CIC-DUX4-mediated transcriptional regulation and proposes targeting p300 as a promising strategy for the clinical management of CDS. 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. [ABSTRACT FROM AUTHOR]

Background & Aims: Acinar-to-ductal metaplasia (ADM) serves as a precursor event in the development of pancreatic ductal adenocarcinoma (PDAC) upon constitutive environmental and genetical stress. While the role of ADM in PDAC progression has been established, the molecular mechanisms underlying human ADM remain elusive. We previously demonstrated the induction of ADM in human acinar cells through the transforming growth factor beta (TGFβ) signaling pathway. We aim to investigate the interaction between TGFβ and Hippo pathways in mediating ADM. Methods: RNA-sequencing was conducted on sorted normal primary human acinar, ductal, and AD (acinar cells that have undergone ADM) cells. ATAC-seq analysis was utilized to reveal the chromatin accessibility in these three cell types. ChIP-Seq of YAP1, SMAD4, and H3K27ac was performed to identify the gene targets of YAP1 and SMAD4. The role of YAP1/TAZ in ADM-driven cell proliferation, as well as in oncogenic KRAS driven proliferation, was assessed using sphere formation assay. Results: AD cells have a unique transcription profile, with upregulated genes in open chromatin states in acinar cells. YAP1 and SMAD4 co-occupy the loci of ADM-related genes, including PROM1, HES1, and MMP7, co-regulating biological functions such as cell adhesion, cell migration, and inflammation. Overexpression of YAP1/TAZ promoted acinar cell proliferation but still required the TGFβ pathway. YAP1/TAZ were also crucial for TGFβ-induced sphere formation and were necessary for KRAS-induced proliferation. Conclusions: Our study reveals the intricate transition between acinar and AD states in human pancreatic tissues. It unveils the complex interaction between the Hippo and TGF-β pathways during ADM, highlighting the pivotal role of YAP1/TAZ and SMAD4 in PDAC initiation. [ABSTRACT FROM AUTHOR]

Summary: 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Objective: To provide novel insights into the aetiology of non-syndromic cleft lip with or without cleft palate (NSCL/P) by integrating multi-omics data and exploring susceptibility genes associated with NSCL/P. Methods: A two-stage genome-wide association study (GWAS) of NSCL/P was performed, involving a total of 1,069 cases and 1,724 controls. Using promoter capture Hi-C (pCHi-C) datasets in human embryonic stem cells (hESC) and chromatin immunoprecipitation sequencing (ChIP-seq) in craniofacial tissues, we filtered out single nucleotide polymorphisms (SNPs) with active cis-regulation and their target genes. Additionally, we employed expression quantitative trait loci (eQTL) analysis to identify candidate genes. Results: Thirteen SNPs were identified as cis-regulation units associated with the risk of NSCL/P. Five of these were proven to be active in chromatin states in early human craniofacial development (rs7218002: odds ratio [OR] 1.50, P = 8.14E-08; rs835367: OR 0.78, P = 3.48E05; rs77022994: OR 0.55, P = 1.05E-04; rs961470: OR 0.73, P = 1.38E-04; rs17314727: OR 0.73, P = 1.85E-04). Additionally, pCHi-C and eQTL analysis prioritised three candidate genes (rs7218002: NTN1, rs835367: FGGY, LINC01135). NTN1 and FGGY were expressed in mouse orofacial development. Deficiencies in NTN1, FGGY and LINC01135 were associated with cleft palate and cleft lip, abnormal facial shape and bifid uvula, and abnormality of the face, respectively. Conclusion: Our study identified five SNPs (rs7218002, rs835367, rs77022994, rs961470 and rs17314727) and three susceptibility genes (NTN1, FGGY and LINC01135) associated with NSCL/P. These findings contribute to a better understanding of the genetic factors involved. [ABSTRACT FROM AUTHOR]