Search

Background: During the latter stages of their development, mammalian oocytes under dramatic chromatin reconfiguration, transitioning from a non-surrounded nucleolus (NSN) to a surrounded nucleolus (SN) stage, and concomitant transcriptional silencing. Although the NSN-SN transition is known to be essential for developmental competence of the oocyte, less is known about the accompanying molecular changes. Here we examine the changes in the transcriptome and DNA methylation during the NSN to SN transition in mouse oocytes. Results: To study the transcriptome and DNA methylation dynamics during the NSN to SN transition, we used single-cell (sc)M&T-seq to generate scRNA-seq and sc-bisulphite-seq (scBS-seq) data from GV oocytes classified as NSN or SN by Hoechst staining of their nuclei. Transcriptome analysis showed a lower number of detected transcripts in SN compared with NSN oocytes as well as downregulation of 576 genes, which were enriched for processes related to mRNA processing. We used the transcriptome data to generate a classifier that can infer chromatin stage in scRNA-seq datasets. The classifier was successfully tested in multiple published datasets of mouse models with a known skew in NSN: SN ratios. Analysis of the scBS-seq data showed increased DNA methylation in SN compared to NSN oocytes, which was most pronounced in regions with intermediate levels of methylation. Overlap with chromatin immunoprecipitation and sequencing (ChIP-seq) data for the histone modifications H3K36me3, H3K4me3 and H3K27me3 showed that regions gaining methylation in SN oocytes are enriched for overlapping H3K36me3 and H3K27me3, which is an unusual combination, as these marks do not typically coincide. Conclusions: We characterise the transcriptome and DNA methylation changes accompanying the NSN-SN transition in mouse oocytes. We develop a classifier that can be used to infer chromatin status in single-cell or bulk RNA-seq data, enabling identification of altered chromatin transition in genetic knock-outs, and a quality control to identify skewed NSN-SN proportions that could otherwise confound differential gene expression analysis. We identify late-methylating regions in SN oocytes that are associated with an unusual combination of chromatin modifications, which may be regions with high chromatin plasticity and transitioning between H3K27me3 and H3K36me3, or reflect heterogeneity on a single-cell level. [ABSTRACT FROM AUTHOR]

As a response regulator of the EsrA-EsrB two-component system, EsrB is conserved in Hafniaceae and plays a crucial role in virulence and pathogenicity. EsrB possesses DNA binding abilities, enabling it to regulate the transcription of virulence genes to confront different stresses and achieve systematic infections. Here, ChIP-seq analysis of EsrB in Dulbecco's Modified Eagle's Medium (DMEM) (mimicking in vivo environments) revealed that EsrB preferred to bind to virulence-associated promoters with a distinct 7'-4-7" pseudopalindromic DNA motif and interact with metabolic-related promoters with a high AT DNA motif. The crystal structure of the C-terminal of EsrB (EsrBC) was solved at 2.20-Å resolution. Specifically, Lys181 enabled the DNA-binding affinity of EsrB and promoted the in vitro and in vivo pathogenicity of Edwardsiella piscicida. Moreover, EsrB directly regulated the expression of genes associated with basal metabolism, including iron and tricarboxylic acid (TCA) cycles. Furthermore, EsrB enhanced iron transport capability and the enzyme activity of succinate dehydrogenase and pyruvate dehydrogenase in DMEM. Collectively, our structural and ChIP-seq analysis provides valuable insights into the DNA binding mechanism of EsrB which will facilitate our understanding of EsrB coordinating virulence and basal metabolism gene expression. [ABSTRACT FROM AUTHOR]

Cervical cancer (CC) is the most common cancer in women. This study aims to explore the molecular mechanism of lactate secreted by CC cells modulating macrophage polarization in CC via histone lactylation. Normal cervical epithelium (NCE), low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), and cervical squamous cell carcinoma (CESC) were collected to assess H3K18la level and macrophage infiltration. Macrophages were incubated with SiHa cell-derived conditioned medium to detect M1 and M2 markers. NCE, HSIL, and CESC samples were used for ChIP-seq of H3K18la. Histone lactylation-dirven GPD2 was knocked down in macrophages. Compared to NCE, H3K18la level and M2 macrophage abundance were increased in LSIL, HSIL, and CESC. Lactate secreted by CC cells upregulated H3K18la and M2 markers but downregulated M1 markers in macrophages. ChIP-seq revealed that upregulated pathways in HSIL vs. NCE and CESC vs. HSIL were commonly enriched in lipid metabolism. Notably, lactate upregulated H3K18la-modified GPD2 expression in macrophages, and GPD2 knockdown reversed lactate induction to M2 macrophages. Collectively, lactate secreted by CC cells upregulates GPD2 via histone lactylation, thereby promoting M2 macrophage polarization in CC. This study provides new insights into the role of histone lactylation in macrophage polarization in the malignant transformation of CC. [ABSTRACT FROM AUTHOR]

Forkhead box (Fox) transcription factors (TFs) mediate multiple conserved cardiogenic processes in both mammals and Drosophila. Our prior work identified the roles of two Drosophila Fox genes, jumeau (jumu) and Checkpoint suppressor 1-like (CHES-1-like), in cardiac progenitor cell specification and division, and in the proper positioning of cardiac cell subtypes. Fox TF binding sites are also significantly enriched in the enhancers of genes expressed in the heart, suggesting that these genes may play a core regulatory role in one or more of these cardiogenic processes. We identified downstream targets of Jumu by comparing transcriptional expression profiles of flow cytometry-sorted mesodermal cells from wild-type embryos and embryos completely lacking the jumu gene and found that genes with functional annotation and ontological features suggesting roles in cell division were overrepresented among Jumu targets. Phenotypic analysis of a subset of these targets identified 21 jumu-regulated genes that mediate cardiac progenitor cell division, one of which, Retinal Homeobox (Rx), was characterized in more detail. Finally, the observation that many of these 21 genes and/or their orthologs exhibit genetic or physical interactions among themselves indicates that Jumu is a master regulator acting as a hub of a cardiac progenitor cell division-mediating network. [ABSTRACT FROM AUTHOR]

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]

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]

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]

Context: Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the central system in epigenomic exploration. Chromatin immunoprecipitation coupled with sequencing (ChIP-seq) is an important technology to identify the genome-wide location of DNA-binding proteins such as histones proteins, transcription factors, RNA polymerase, or any protein of interest. ChIP-seq has been used to study the binding sites and efficacy of drugs in cancer cell lines etc. Aims: In current research, breast cancer cell line data have been used to study the effect PADI2 (peptidyl arginine deiminase) gene in the progression of breast cancer. Further, this ChIP-seq data have also been used to study the binding site of Amanitin drug in breast cancer. Settings and Design: Breast cancer ChIP-seq data have been retrieved from the European Nucleotide Archive database with project Id PRJNA415426 short read archive. Four samples of FASTQ files were used and analyzed for the genome-wide analysis. Materials and Methods: Galaxy server (https://usegalaxy.org/) was used for complete ChIP-seq data analysis; different tools such as fast-quality control (QC), multi-QC, Bowtie2, model-based analysis of ChIP-sequencing, and ChIPseeker tools were used for motif enrichment and functional analysis. Motif analysis was done through the Multiple Expectation maximizations for Motif Elicitation database (https://meme-suite.org/meme/db/motifs). Results: Computational investigation demonstrates the binding sequences of the T47-D breast cancer cell line as TTTTGTATTTTTAGT, and this motif occurs 2123 times in the Homo Sapiens reference genome that is hg19. Conclusions: This research classifies the binding site and affinity of the T47-D human breast cancer cell line. Further, wet laboratory studies are required to verify the function of the predicted motifs and their importance in drug development or research in breast cancer.

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

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.

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.

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]

Context: Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the central system in epigenomic exploration. Chromatin immunoprecipitation coupled with sequencing (ChIP-seq) is an important technology to identify the genome-wide location of DNA-binding proteins such as histones proteins, transcription factors, RNA polymerase, or any protein of interest. ChIP-seq has been used to study the binding sites and efficacy of drugs in cancer cell lines etc. Aims: In current research, breast cancer cell line data have been used to study the effect PADI2 (peptidyl arginine deiminase) gene in the progression of breast cancer. Further, this ChIP-seq data have also been used to study the binding site of Amanitin drug in breast cancer. Settings and Design: Breast cancer ChIP-seq data have been retrieved from the European Nucleotide Archive database with project Id PRJNA415426 short read archive. Four samples of FASTQ files were used and analyzed for the genome-wide analysis. Materials and Methods: Galaxy server (https://usegalaxy.org/) was used for complete ChIP-seq data analysis; different tools such as fast-quality control (QC), multi-QC, Bowtie2, model-based analysis of ChIP-sequencing, and ChIPseeker tools were used for motif enrichment and functional analysis. Motif analysis was done through the Multiple Expectation maximizations for Motif Elicitation database (https://meme-suite.org/meme/db/motifs). Results: Computational investigation demonstrates the binding sequences of the T47-D breast cancer cell line as TTTTGTATTTTTAGT, and this motif occurs 2123 times in the Homo Sapiens reference genome that is hg19. Conclusions: This research classifies the binding site and affinity of the T47-D human breast cancer cell line. Further, wet laboratory studies are required to verify the function of the predicted motifs and their importance in drug development or research in breast cancer. [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]

Mitochondrial biogenesis requires the expression of genes encoded by both the nuclear and mitochondrial genomes. However, aside from a handful transcription factors regulating specific subsets of mitochondrial genes, the overall architecture of the transcriptional control of mitochondrial biogenesis remains to be elucidated. The mechanisms coordinating these two genomes are largely unknown. We performed a targeted RNAi screen in developing eyes with reduced mitochondrial DNA content, anticipating a synergistic disruption of tissue development due to impaired mitochondrial biogenesis and mitochondrial DNA (mtDNA) deficiency. Among 638 transcription factors annotated in the Drosophila genome, 77 were identified as potential regulators of mitochondrial biogenesis. Utilizing published ChIP-seq data of positive hits, we constructed a regulatory network revealing the logic of the transcription regulation of mitochondrial biogenesis. Multiple transcription factors in core layers had extensive connections, collectively governing the expression of nearly all mitochondrial genes, whereas factors sitting on the top layer may respond to cellular cues to modulate mitochondrial biogenesis through the underlying network. CG1603, a core component of the network, was found to be indispensable for the expression of most nuclear mitochondrial genes, including those required for mtDNA maintenance and gene expression, thus coordinating nuclear genome and mtDNA activities in mitochondrial biogenesis. Additional genetic analyses validated YL-1, a transcription factor upstream of CG1603 in the network, as a regulator controlling CG1603 expression and mitochondrial biogenesis.

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.

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.

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.

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.

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]

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]

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

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.

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

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]

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]

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 FROM AUTHOR]