ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation

ATF4:p52 complex activates oncogenic enhancers in multiple myeloma via p300/CBP recruitment to regulate BACH1

Multiple myeloma (MM) is a B-cell malignancy accounting for 20 % of all blood-associated cancers. MM patients with a poorer prognosis and high-risk stratification were previously observed to be causally linked to the constitutive activation of non-canonical NF-κB (ncNF-κB) pathway. Consistent with this, the ncNF-κB p52 transcription factor was earlier found to regulate the enhancer landscape of MM to potentiate oncogenic transcription. However, the mechanism by which aberrant p52 expression is involved in coordinating enhancer activity has not been well explored. In this study, we analysed H3K27ac ChIP-seq and ATAC-seq data from MM cell lines and patient samples to screen for putative transcription factors that cooperate with p52 to regulate enhancers activated in MM. We report that ATF4 interacts with p52 and together, this complex mediates the activity of a subset of MM-associated enhancers through the recruitment of histone acetyltransferases (HATs), p300 and CBP (CREB-binding protein). We also identified a ATF4:p52 regulated target gene BACH1 under the regulation of a proximal super-enhancer, which was found to drive oncogenesis in MM by promoting cell cycle progression and proliferation. Together, our findings provide further mechanistic insights into how aberrant enhancer activation observed in MM tumours could lead to disease progression.

Repeat-element RNAs integrate a neuronal growth circuit

Neuronal growth and regeneration are regulated by local translation of mRNAs in axons. We examined RNA polyadenylation changes upon sensory neuron injury and found upregulation of a subset of polyadenylated B2-SINE repeat elements, hereby termed GI-SINEs (growth-inducing B2-SINEs). GI-SINEs are induced from ATF3 and other AP-1 promoter-associated extragenic loci in injured sensory neurons but are not upregulated in lesioned retinal ganglion neurons. Exogenous GI-SINE expression elicited axonal growth in injured sensory, retinal, and corticospinal tract neurons. GI-SINEs interact with ribosomal proteins and nucleolin, an axon-growth-regulating RNA-binding protein, to regulate translation in neuronal cytoplasm. Finally, antisense oligos against GI-SINEs perturb sensory neuron outgrowth and nucleolin-ribosome interactions. Thus, a specific subfamily of transposable elements is integral to a physiological circuit linking AP-1 transcription with localized RNA translation.

Temporally discordant chromatin accessibility and DNA demethylation define short- and long-term enhancer regulation during cell fate specification

Chromatin and DNA modifications mediate the transcriptional activity of lineage-specifying enhancers, but recent work challenges the dogma that joint chromatin accessibility and DNA demethylation are prerequisites for transcription. To understand this paradox, we established a highly resolved timeline of their dynamics during neural progenitor cell differentiation. We discovered that, while complete demethylation appears delayed relative to shorter-lived chromatin changes for thousands of enhancers, DNA demethylation actually initiates with 5-hydroxymethylation before appreciable accessibility and transcription factor occupancy is observed. The extended timeline of DNA demethylation creates temporal discordance appearing as heterogeneity in enhancer regulatory states. Few regions ever gain methylation, and resulting enhancer hypomethylation persists long after chromatin activities have dissipated. We demonstrate that the temporal methylation status of CpGs (mC/hmC/C) predicts past, present, and future chromatin accessibility using machine learning models. Thus, chromatin and DNA methylation collaborate on different timescales to shape short- and long-term enhancer regulation during cell fate specification.

ProOvErlap: Assessing feature proximity/overlap and testing statistical significance from genomic intervals

Feature overlap is a critical concept in bioinformatics and occurs when two genomic intervals, usually represented as BED files, are located in the same genomic regions. Instead, feature proximity refers to the spatial proximity of genomic elements. For example, promoters typically overlap or are close to the genes they regulate. Overlap and proximity are also important in epigenetic studies. Here, the overlap of regions enriched for specific epigenetic modifications or accessible chromatin can elucidate complex molecular phenotypes. Consequently, the ability to analyze and interpret feature overlap and proximity is essential for understanding the biological processes that contribute to a given phenotype. To address this need, we present a computational method capable of analyzing data represented in the BED format. This method aims to quantitatively assess the degree of proximity or overlap between genomic features and to determine the statistical significance of these events in the context of a non-parametric randomization test. The aim is to understand whether the observed state differs from what would be expected by chance. The method is designed to be easy to use, requiring only a single command line to run, allowing straightforward overlap and proximity analysis. It also provides clear visualizations and publication-quality figures. In conclusion, this study highlights the importance of feature overlap and proximity in epigenetic studies and presents a method to improve the systematic assessment and interpretation of these features. A new resource for identifying biologically significant interactions between genomic features in both healthy and disease states.

Chromatin loops are an ancestral hallmark of the animal regulatory genome

In bilaterian animals, gene regulation is shaped by a combination of linear and spatial regulatory information. Regulatory elements along the genome are integrated into gene regulatory landscapes through chromatin compartmentalization1,2, insulation of neighbouring genomic regions3,4 and chromatin looping that brings together distal cis-regulatory sequences5. However, the evolution of these regulatory features is unknown because the three-dimensional genome architecture of most animal lineages remains unexplored6,7. To trace the evolutionary origins of animal genome regulation, here we characterized the physical organization of the genome in non-bilaterian animals (sponges, ctenophores, placozoans and cnidarians)8,9 and their closest unicellular relatives (ichthyosporeans, filastereans and choanoflagellates)10 by combining high-resolution chromosome conformation capture11,12 with epigenomic marks and gene expression data. Our comparative analysis showed that chromatin looping is a conserved feature of genome architecture in ctenophores, placozoans and cnidarians. These sequence-determined distal contacts involve both promoter-enhancer and promoter-promoter interactions. By contrast, chromatin loops are absent in the unicellular relatives of animals. Our findings indicate that spatial genome regulation emerged early in animal evolution. This evolutionary innovation introduced regulatory complexity, ultimately facilitating the diversification of animal developmental programmes and cell type repertoires.

HDAC7 promotes cardiomyocyte proliferation by suppressing myocyte enhancer factor 2

Postnatal mammalian cardiomyocytes (CMs) rapidly lose proliferative capacity and exit the cell cycle to undergo further differentiation and maturation. Cell cycle activation has been a major strategy to stimulate postnatal CM proliferation, albeit achieving modest effects. One impediment is that postnatal CMs may need to undergo dedifferentiation before proliferation, if not simultaneously. Here, we report that overexpression of Hdac7 in neonatal mouse CMs results in significant CM dedifferentiation and proliferation. Mechanistically, we show that histone deacetylase 7 (HDAC7)-mediated CM proliferation is contingent on dedifferentiation, which is accomplished by suppressing myocyte enhance factor 2 (MEF2). Hdac7 overexpression in CM shifts the chromatin state from binding with MEF2, which favors the transcriptional program toward differentiation, to binding with AP-1, which favors the transcriptional program toward proliferation. Furthermore, we found that HDAC7 interacts with minichromosome maintenance complex components to initiate cell cycle progression. Our findings reveal that HDAC7 promotes CM proliferation by its dual action on CM dedifferentiation and proliferation, uncovering a potential new strategy for heart regeneration/repair.

Histone chaperones coupled to DNA replication and transcription control divergent chromatin elements to maintain cell fate

The manipulation of DNA replication and transcription can be harnessed to control cell fate. Central to the regulation of these DNA-templated processes are histone chaperones, which in turn are emerging as cell fate regulators. Histone chaperones are a group of proteins with diverse functions that are primarily involved in escorting histones to assemble nucleosomes and maintain the chromatin landscape. Whether distinct histone chaperone pathways control cell fate and whether they function using related mechanisms remain unclear. To address this, we performed a screen to assess the requirement of diverse histone chaperones in the self-renewal of hematopoietic stem and progenitor cells. Remarkably, all candidates were required to maintain cell fate to differing extents, with no clear correlation with their specific histone partners or DNA-templated process. Among all the histone chaperones, the loss of the transcription-coupled histone chaperone SPT6 most strongly promoted differentiation, even more than the major replication-coupled chromatin assembly factor complex CAF-1. To directly compare how DNA replication- and transcription-coupled histone chaperones maintain stem cell self-renewal, we generated an isogenic dual-inducible system to perturb each pathway individually. We found that SPT6 and CAF-1 perturbations required cell division to induce differentiation but had distinct effects on cell cycle progression, chromatin accessibility, and lineage choice. CAF-1 depletion led to S-phase accumulation, increased heterochromatic accessibility (particularly at H3K27me3 sites), and aberrant multilineage gene expression. In contrast, SPT6 loss triggered cell cycle arrest, altered accessibility at promoter elements, and drove lineage-specific differentiation, which is in part influenced by AP-1 transcription factors. Thus, CAF-1 and SPT6 histone chaperones maintain cell fate through distinct mechanisms, highlighting how different chromatin assembly pathways can be leveraged to alter cell fate.

IL-10 targets IRF transcription factors to suppress IFN and inflammatory response genes by epigenetic mechanisms

Interleukin-10 (IL-10) is pivotal in suppressing innate immune activation, in large part by suppressing induction of inflammatory genes. Despite decades of research, the molecular mechanisms underlying this inhibition have not been resolved. Here we utilized an integrated epigenomic analysis to investigate IL-10-mediated suppression of LPS and TNF responses in primary human monocytes. Instead of inhibiting core TLR4-activated pathways such as NF-κB, MAPK-AP-1 and TBK1-IRF3 signaling, IL-10 targeted IRF transcription factor activity and DNA binding, particularly IRF5 and an IRF1-mediated amplification loop. This resulted in suppression of inflammatory NF-κB target genes and near-complete suppression of interferon-stimulated genes. Mechanisms of gene inhibition included downregulation of chromatin accessibility, de novo enhancer formation and IRF1-associated H3K27ac activating histone marks. These results provide a mechanism by which IL-10 suppresses inflammatory NF-κB target genes, highlight the role of IRF1 in inflammatory gene expression and describe the suppression of IFN responses by epigenetic mechanisms.

Super RNA Pol II domains enhance minor ZGA through 3D interaction to ensure the integrity of major transcriptional waves in late-ZGA mammals

Zygotic genome activation (ZGA) occurs at distinct stages across mammals, with mice initiating ZGA at the 2-cell stage and bovines and humans activating the process in the 4- to 8-cell stages. RNA polymerase II (RNA Pol II) gradually initiates ZGA in mice, but regulation in late-ZGA species remains unclear. Here, RNA Pol II profiling in bovine embryos identified strong intergenic clusters that boost minor ZGA gene expression via chromatin interactions and are named super RNA Pol II domains (SPDs). CRISPRi perturbation of SPDs in bovine embryos decreases the expression of minor ZGA genes, whereas the knockdown of these genes disrupts major ZGA and embryogenesis. Rapid enhancement of minor ZGA genes also occurs in human embryos. Alternatively, mouse and porcine oocytes precociously express these minor ZGA genes without SPDs. Thus, SPDs appear to be an adaptation in bovine embryos, promoting minor ZGA gene expression to comparable levels as early-ZGA species, illuminating species-specific regulation of ZGA timing.

Uncovering the transcriptional regulatory network underlying selenium tolerance in maize seedlings

Selenium (Se) plays a dual role in plant growth, functioning as both an essential micronutrient and a potential toxin. Understanding the regulatory mechanisms of Se tolerance is crucial for enhancing crop resilience and biofortification. In this study, we integrated transcriptomics (RNA-seq), chromatin accessibility (ATAC-seq), and genome-wide association studies (GWAS) to elucidate the regulatory networks governing Se responses in maize seedlings. Low Se concentrations (≤ 0.05 mM) enhanced plant growth and biomass accumulation, whereas high Se concentrations (≥ 0.1 mM) induced toxicity and suppressed growth. Different treatment groups exhibited dose-dependent transcriptional reprogramming, with significant upregulation of genes involved in glutathione biosynthesis, Se metabolism, and jasmonic acid (JA) signaling. Concurrent chromatin accessibility remodeling in promoter regions orchestrated the transcriptional responses of these key genes. Characterization of the ZmGSTs gene family revealed subfamily-specific expression patterns and regulatory mechanisms under Se stress. Integration of high-confidence transcriptional regulatory networks with GWAS data led to the identification of a key metabolic gene (ZmGSR2) in the selenium metabolism pathway and three important transcription factors (ZmWRKY48, ZmbZIP123, and ZmKNOX6) that specifically activated distinct ZmGSTs genes. This study provides novel insights into the genetic and epigenetic mechanisms underlying selenium tolerance and identifies potential targets for improving crop selenium adaptability.

Variant-specific disruption to notch signalling in PAX6 microphthalmia and aniridia patient-derived hiPSC optic cup-like organoids

The homeobox-containing transcription factor PAX6 is a key regulator of eye development. Pathogenic heterozygous PAX6 variants lead to variable ocular phenotypes, most commonly haploinsufficiency-induced aniridia. Missense variants are typically associated with milder ocular conditions, although variants in the DNA-binding paired domain which alter target binding lead to severe ocular phenotypes including bilateral microphthalmia, similar to SOX2-anophthalmia syndrome. However, the variant-specific pathway disruption resulting in phenotypic heterogeneity is not well understood. To investigate pathogenic mechanisms of PAX6 variants, transcriptomic and chromatin accessibility analysis was performed on hiPSC derived 3D optic cup-like organoids generated from patients with variants (i) PAX6N124K displaying combined microphthalmia, aniridia and optic nerve coloboma, and (ii) PAX6R261X exhibiting typical aniridia. Total RNA sequencing analysis revealed downregulation of SOX2 in missense PAX6N124K cups compared to both wildtype and PAX6R261X haploinsufficient aniridia controls, along with Notch signalling components and markers of proliferation and differentiation. Transcription factor binding motifs of Notch-related genes were also found to be differentially bound in PAX6N124K cups through ATACseq footprinting analysis. Our analysis of PAX6-related oculopathies using in vitro models reveals disruption to DNA binding perturbs SOX2 and Notch signalling, contributing to severe ocular phenotypes in patients with missense changes in the paired domain. This work reveals a previously unestablished role for PAX6 in SOX2 and Notch signalling regulation during early oculogenesis, as well as illuminating disease mechanisms underlying variant-specific ocular phenotypes and genotype-phenotype correlations. These novel insights can influence clinical care, and provide valuable data on potential therapeutic targets, which can guide future translational research.

Orchestration of pluripotent stem cell genome reactivation during mitotic exit

Cell identity maintenance faces many challenges during mitosis, as most DNA-binding proteins are evicted from DNA and transcription is virtually abolished. How cells maintain their identity through division and faithfully re-initiate gene expression during mitotic exit is unclear. Here, we develop a novel reporter system enabling cell cycle synchronization-free separation of pluripotent stem cells in temporal bins of <30 min during mitotic exit. This allows us to quantify genome-wide reactivation of transcription, sequential changes in chromatin accessibility and transcription factor footprints, and re-binding of the pluripotency transcription factors OCT4, SOX2, and NANOG (OSN). We find that transcriptional activity progressively ramps up after mitosis and that OSN rapidly reoccupy the genome during the anaphase-telophase transition. We also demonstrate transcription factor-specific, dynamic relocation patterns and a hierarchical reorganization of the OSN binding landscape governed by OCT4 and SOX2. Our study sheds light on the dynamic orchestration of transcriptional reactivation after mitosis.

The genetic architecture of cell type-specific cis regulation in maize

Gene expression and complex phenotypes are determined by the activity of cis-regulatory elements. However, an understanding of how extant genetic variants affect cis regulation remains limited. Here, we investigated the consequences of cis-regulatory diversity using single-cell genomics of more than 0.7 million nuclei across 172 Zea mays (maize) inbreds. Our analyses pinpointed cis-regulatory elements distinct to domesticated maize and revealed how historical transposon activity has shaped the cis-regulatory landscape. Leveraging population genetics principles, we fine-mapped about 22,000 chromatin accessibility-associated genetic variants with widespread cell type-specific effects. Variants in TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR-binding sites were the most prevalent determinants of chromatin accessibility. Finally, integrating chromatin accessibility-associated variants, organismal trait variation, and population differentiation revealed how local adaptation has rewired regulatory networks in unique cellular contexts to alter maize flowering.

A Sox2 enhancer cluster regulates region-specific neural fates from mouse embryonic stem cells

Sex-determining region Y box 2 (Sox2) is a critical transcription factor for embryogenesis and neural stem and progenitor cell (NSPC) maintenance. While distal enhancers control Sox2 in embryonic stem cells (ESCs), enhancers closer to the gene are implicated in Sox2 transcriptional regulation in neural development. We hypothesize that a downstream enhancer cluster, termed Sox2 regulatory regions 2-18 (SRR2-18), regulates Sox2 transcription in neural stem cells and we investigate this in NSPCs derived from mouse ESCs. Using functional genomics and CRISPR-Cas9-mediated deletion analyses, we investigate the role of SRR2-18 in Sox2 regulation during neural differentiation. Transcriptome analyses demonstrate that the loss of even 1 copy of SRR2-18 disrupts the region-specific identity of NSPCs, reducing the expression of genes associated with more anterior regions of the embryonic nervous system. Homozygous deletion of this Sox2 neural enhancer cluster causes reduced SOX2 protein, less frequent interaction with transcriptional machinery, and leads to perturbed chromatin accessibility genome-wide further affecting the expression of neurodevelopmental and anterior-posterior regionalization genes. Furthermore, homozygous NSPC deletants exhibit self-renewal defects and impaired differentiation into cell types found in the brain. Altogether, our data define a cis-regulatory enhancer cluster controlling Sox2 transcription in NSPCs and highlight the sensitivity of neural differentiation processes to decreased Sox2 transcription, which causes differentiation into posterior neural fates, specifically the caudal neural tube. This study highlights the importance of precise Sox2 regulation by SRR2-18 in neural differentiation.

Hypoxia-induced phase separation of ZHX2 alters chromatin looping to drive cancer metastasis

Hypoxia and dysregulated phase separation can both activate oncogenic transcriptomic profiles. However, whether hypoxia regulates transcription-associated phase separation remains unknown. Here, we find that zinc fingers and homeoboxes 2 (ZHX2) undergoes phase separation in response to hypoxia, promoting their occupancy on chromatin and activating a cluster of oncogene transcription that is enriched by metastatic genes distinct from the targets of hypoxia-inducible factor (HIF) and pathologically relevant to breast cancer. Hypoxia induces ZHX2 phase separation via a proline-rich intrinsically disordered region (IDR), enhancing phosphorylation of ZHX2 at S625 and S628 that incorporates CCCTC-binding factor (CTCF) in condensates to alter chromatin looping, consequently driving metastatic gene transcription and cancer metastasis. Our findings provide significant insight into oncogene activation and suggest a phase-separation-based therapeutic strategy for cancer.

Inferring differential protein binding from time-series chromatin accessibility data

Motivation

Due to internal and external factors, the epigenomic landscape is constantly changing in ways that are linked to changes in gene expression. Chromatin accessibility data, such as MNase-seq, provide valuable insights into this landscape and have been used to compute chromatin occupancy profiles. Multiple datasets generated over time or under different conditions can thus be used to study dynamic changes in chromatin occupancy across the genome.

Results

Our existing model, RoboCOP, computes a genome-wide chromatin occupancy profile for nucleosomes and hundreds of transcription factors. Here, we present a new method called DynaCOP that takes multiple chromatin occupancy profiles and uses them to generate a series of nucleosome-guided difference profiles. These profiles identify differentially binding transcription factors and reveal changes in nucleosome occupancy and positioning. We apply DynaCOP to chromatin occupancy profiles derived from deeply sequenced time-series MNase-seq data to study differential chromatin occupancy in the yeast genome under cadmium stress. We find strong correlations between the observed chromatin changes and changes in transcription.

Availability and implementation

https://github.com/HarteminkLab/RoboCOP.

Characterization of non-coding variants associated with transcription-factor binding through ATAC-seq-defined footprint QTLs in liver

Non-coding variants discovered by genome-wide association studies (GWASs) are enriched in regulatory elements harboring transcription factor (TF) binding motifs, strongly suggesting a connection between disease association and the disruption of cis-regulatory sequences. Occupancy of a TF inside a region of open chromatin can be detected in ATAC-seq where bound TFs block the transposase Tn5, leaving a pattern of relatively depleted Tn5 insertions known as a "footprint." Here, we sought to identify variants associated with TF binding, or "footprint quantitative trait loci" (fpQTLs), in ATAC-seq data generated from 170 human liver samples. We used computational tools to scan the ATAC-seq reads to quantify TF binding likelihood as "footprint scores" at variants derived from whole-genome sequencing generated in the same samples. We tested for association between genotype and footprint score and observed 809 fpQTLs associated with footprint-inferred TF binding (FDR < 5%). Given that Tn5 insertion sites are measured with base-pair resolution, we show that fpQTLs can aid GWAS and QTL fine-mapping by precisely pinpointing TF activity within broad trait-associated loci where the underlying causal variant is unknown. Liver fpQTLs were strongly enriched across ChIP-seq peaks, liver expression QTLs (eQTLs), and liver-related GWAS loci, and their inferred effect on TF binding was concordant with their effect on underlying sequence motifs in 78% of cases. We conclude that fpQTLs can reveal causal GWAS variants, define the role of TF binding-site disruption in complex traits, and provide functional insights into non-coding variants, ultimately informing novel treatments for common diseases.

Binding domain mutations provide insight into CTCF's relationship with chromatin and its contribution to gene regulation

Here we used a series of CTCF mutations to explore CTCF's relationship with chromatin and its contribution to gene regulation. CTCF's impact depends on the genomic context of bound sites and the unique binding properties of WT and mutant CTCF proteins. Specifically, CTCF's signal strength is linked to changes in accessibility, and the ability to block cohesin is linked to its binding stability. Multivariate modeling reveals that both CTCF and accessibility contribute independently to cohesin binding and insulation, but CTCF signal strength has a stronger effect. CTCF and chromatin have a bidirectional relationship such that at CTCF sites, accessibility is reduced in a cohesin-dependent, mutant-specific fashion. In addition, each mutant alters TF binding and accessibility in an indirect manner, changes which impart the most influence on rewiring transcriptional networks and the cell's ability to differentiate. Collectively, the mutant perturbations provide a rich resource for determining CTCF's site-specific effects.

LMX1B missense-perturbation of regulatory element footprints disrupts serotonergic forebrain axon arborization

Pathogenic coding mutations are prevalent in human neuronal transcription factors (TFs) but how they disrupt development is poorly understood. Lmx1b is a master transcriptional regulator of postmitotic Pet1 neurons that give rise to mature serotonin (5-HT) neurons; over two hundred pathogenic heterozygous mutations have been discovered in human LMX1B, yet their impact on brain development has not been investigated. Here, we developed mouse models with different LMX1B DNA-binding missense mutations. Missense heterozygosity broadly altered Pet1 neuron transcriptomes, but expression changes converged on axon and synapse genes. Missense heterozygosity effected highly specific deficits in the postnatal maturation of forebrain serotonin axon arbors, primarily in the hippocampus and motor cortex, which was associated with spatial memory defects. Digital genomic footprinting (DGF) revealed that missense heterozygosity caused complete loss of Lmx1b motif protection and chromatin accessibility at sites enriched for a distal active enhancer/active promoter histone signature and homeodomain binding motifs; at other bound Lmx1b motifs, varying levels of losses, gains, or no change in motif binding and accessibility were found. The spectrum of footprint changes was strongly associated with synapse and axon genes. Further, Lmx1b missense heterozygosity caused wide disruption of Lmx1b-dependent GRNs comprising diverse TFs expressed in Pet1 neurons. These findings reveal an unanticipated continuum of Lmx1b missense-forced perturbations on Pet1 neuron regulatory element TF binding and accessibility. Our work illustrates DGF's utility for gaining unique insight into how expressed TF missense mutations interfere with developing neuronal GRNs.

Menin-MLL1 complex cooperates with NF-Y to promote HCC survival

Identification of new therapeutic targets in hepatocellular carcinoma (HCC) remains critical. Chromatin regulating complexes are frequently mutated or aberrantly expressed in HCC, suggesting dysregulation of chromatin environments is a key feature driving liver cancer. To investigate whether the altered chromatin state in HCC cells could be targeted, we designed and utilized an epigenome-focused CRISPR library that targets genes involved in chromatin regulation. This focused approach allowed us to test multiple HCC cell lines in both 2D and 3D growth conditions, which revealed striking differences in the essentiality of genes involved in ubiquitination and multiple chromatin regulators vital for HCC cell survival in 2D but whose loss promoted growth in 3D. We found the core subunits of the menin-MLL1 complex among the strongest essential genes for HCC survival in all screens and thoroughly characterized the mechanism through which the menin-MLL1 complex promotes HCC cell growth. Inhibition of the menin-MLL1 interaction led to global changes in occupancy of the complex with concomitant decreases in H3K4me3 and expression of genes involved in PI3K/AKT/mTOR signaling pathway. Menin inhibition affected chromatin accessibility in HCC cells, revealing that increased chromatin accessibility at sites not bound by menin-MLL1 was associated with the recruitment of the pioneer transcription factor complex NF-Y. A CRISPR/Cas9 screen of chromatin regulators in the presence of menin inhibitor SNDX-5613 revealed a significantly increased cell death when combined with NFYB knockout. Together these data show that menin-MLL1 is necessary for HCC cell survival and cooperates with NF-Y to regulate oncogenic gene transcription.

Amplified dosage of the NKX2-1 lineage transcription factor controls its oncogenic role in lung adenocarcinoma

Amplification-mediated oncogene overexpression is a critical and widespread driver event in cancer, yet our understanding of how amplification and dosage mediate oncogene regulation is limited. Here, we find that the most significant focal amplification event in lung adenocarcinoma (LUAD) targets a lineage "super-enhancer" near the NKX2-1 lineage transcription factor. The NKX2-1 super-enhancer is targeted by focal and co-amplification with NKX2-1 and controls NKX2-1 expression and regulation. We find that NKX2-1 directly controls enhancer accessibility to drive a lineage-addicted state in LUAD. We precisely map the effects of NKX2-1 dosage modulation upon both overexpression and knockdown and identify both linear and non-linear regulation by NKX2-1 dosage. We find that NKX2-1 is a widespread dependency in LUAD cell lines and that NKX2-1 confers persistence to EGFR inhibitors. Our data suggest a defining role for dosage in the oncogenic regulation of amplified NKX2-1 and that amplified NKX2-1 lineage addiction defines LUAD tumors.

Time-course transcriptome and chromatin accessibility analyses reveal the dynamic transcriptional regulation shaping spikelet hull size

The grains of rice (Oryza sativa) are enclosed by a spikelet hull comprising the lemma and palea. Development of the spikelet hull determines the storage capacity of the grain, thus affecting grain yield and quality. Although multiple signaling pathways controlling grain size have been identified, the transcriptional regulatory mechanisms underlying grain development remain limited. Here, we used RNA-seq and ATAC-seq to characterize the transcription and chromatin accessibility dynamics during the development of spikelet hulls. A time-course analysis showed that more than half of the genes were sequentially expressed during hull development and that the accessibility of most open chromatin regions (OCRs) changed moderately, although some regions positively or negatively affected the expression of their closest genes. We revealed a crucial role of GROWTH-REGULATING FACTORs in shaping grain size by influencing multiple metabolic and signaling pathways, and a coordinated transcriptional regulation in response to auxin and cytokinin signaling. We also demonstrated the function of SCL6-IIb, a member of the GRAS family transcription factors, in regulating grain size, with SCL6-IIb expression being activated by SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 18 (OsSPL18). When we edited the DNA sequences within OCRs upstream of the start codon of BRASSINAZOLE-RESISTANT 1 (BZR1) and SCL6-IIb, we generated multiple mutant lines with longer grains. These findings offer a comprehensive overview of the cis-regulatory landscape involved in forming grain capacity and a valuable resource for exploring the regulatory network behind grain development.

Evolutionary Process Underlying Receptor Gene Expansion and Cellular Divergence of Olfactory Sensory Neurons in Honeybees

Olfaction is crucial for animals' survival and adaptation. Unlike the strict singular expression of odorant receptor (OR) genes in vertebrate olfactory sensory neurons (OSNs), insects exhibit complex OR gene expression patterns. In honeybees (Apis mellifera), a significant expansion of OR genes implies a selection preference for the olfactory demands of social insects. However, the mechanisms underlying receptor expression specificity and their contribution to OSN divergence remain unclear. In this study, we used single-nucleus multiomics profiling to investigate the transcriptional regulation of OR genes and the cellular identity of OSNs in A. mellifera. We identified three distinct OR expression patterns, singular OR expression, co-expression of multiple OR genes with a single active promoter, and co-expression of multiple OR genes with multiple active promoters. Notably, ∼50% of OSNs co-expressed multiple OR genes, driven by polycistronic transcription of tandemly duplicated OR genes via a single active promoter. In these OSNs, their identity was determined by the first transcribed receptor. The divergent activation of the promoter for duplicated OR genes ensures the coordinated increased divergence of OSN population. By integrating multiomics data with genomic architecture, we illustrate how fundamental genetic mechanisms drive OR gene expansion and influence flanking regulatory elements, ultimately contributing to the cellular divergence of OSNs. Our findings highlight the interplay between gene duplication and regulatory evolution in shaping OSN diversity, providing new insights into the evolution and adaptation of olfaction in social insects. This study also sheds light on how genetic innovations contribute to the evolution of complex traits.

Resilience in Resistance: The Role of Cell Wall Integrity in Multidrug-Resistant Candida

The Candida species cell wall plays a pivotal role as a structural and functional barrier against external aggressors and as an intermediary in host-pathogen interactions. Candida species exhibit unique adaptations in their cell wall composition, with varying proportions of chitin, mannans, and β-glucans influenced by the environmental conditions and the morphological states. These components not only maintain cellular viability under osmotic, thermal, and chemical stress, but also serve as the key targets for novel antifungal strategies. MAPK signaling pathways, like the cell wall integrity pathway and the high-osmolarity glycerol pathway, play a crucial role in responding to cell wall stressors. Due to the rise of antifungal resistance and its clinical challenges, there is a need to identify new antifungal targets. This review discusses the recent advances in understanding the mechanisms underlying cell wall integrity, their impact on antifungal resistance and virulence, and their potential as therapeutic targets of C. albicans, N. glabratus, and C. auris.

Genotype inference from aggregated chromatin accessibility data reveals genetic regulatory mechanisms

BACKGROUND

Understanding the genetic causes underlying variability in chromatin accessibility can shed light on the molecular mechanisms through which genetic variants may affect complex traits. Thousands of ATAC-seq samples have been collected that hold information about chromatin accessibility across diverse cell types and contexts, but most of these are not paired with genetic information and come from distinct projects and laboratories.

RESULTS

We report here joint genotyping, chromatin accessibility peak calling, and discovery of quantitative trait loci which influence chromatin accessibility (caQTLs), demonstrating the capability of performing caQTL analysis on a large scale in a diverse sample set without pre-existing genotype information. Using 10,293 profiling samples representing 1454 unique donor individuals across 653 studies from public databases, we catalog 24,159 caQTLs in total. After joint discovery analysis, we cluster samples based on accessible chromatin profiles to identify context-specific caQTLs. We find that caQTLs are strongly enriched for annotations of gene regulatory elements across diverse cell types and tissues and are often linked with genetic variation associated with changes in expression (eQTLs), indicating that caQTLs can mediate genetic effects on gene expression. We demonstrate sharing of causal variants for chromatin accessibility across human traits, enabling a more complete picture of the genetic mechanisms underlying complex human phenotypes.

CONCLUSIONS

Our work provides a proof of principle for caQTL calling from previously ungenotyped samples and represents one of the largest, most diverse caQTL resources currently available, informing mechanisms of genetic regulation of gene expression and contribution to disease.

Genome-wide chromatin accessibility and selective signals of meat rabbits reveal key Cis-regulatory elements and variants during postnatal development of skeletal muscles in rabbits

BACKGROUND

The development of skeletal muscles is intricately modulated by multiple genetic factors and significantly impacts the economic value of meat rabbits. However, our knowledge of epigenetics in rabbit skeletal muscles remains largely unknown.

RESULTS

In this study, we collected leg skeletal muscles of rabbits and performed assays for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) to detect open chromatin across three developmental stages: birth (D1), weaning (D35), and adulthood (D75). A total of 126,959 accessible chromatin regions (ACRs) were identified across samples, and a broad increase and decrease in chromatin accessibility were found from D1 to D35 and D35 to D75, respectively. Integrative analysis of chromatin accessibility and transcriptome data revealed ACRs that were nearly closed at D1 but highly accessible at D35 and D75 were significantly enriched in skeletal muscle development. Cis-regulation analysis further revealed that genes dominated by enhancers mainly play roles in the neuron development of rabbit skeletal muscles. Moreover, the detection of selection signals of meat rabbits and the footprinting analysis of transcription factor at open chromatin revealed that both base transversion (Chr13:12144967 A-> G) and the dynamics of chromatin accessibility at the PRDM1 binding site might regulate ZSWIM5 during the development of skeletal muscles in rabbits.

CONCLUSIONS

Our study provided a category of potential cis-regulatory elements for understanding the development of skeletal muscles at the tissue level and might facilitate potential insights into growth regulation in rabbits.

The role of transposon activity in shaping cis-regulatory element evolution after whole-genome duplication

Whole-genome duplications (WGDs) and transposable element (TE) activity can act synergistically in genome evolution. WGDs can increase TE activity directly through cellular stress or indirectly by relaxing selection against TE insertions in functionally redundant, duplicated regions. Because TEs can function as, or evolve into, TE-derived cis-regulatory elements (TE-CREs), bursts of TE activity following WGD are therefore likely to impact evolution of gene regulation. Yet, the role of TEs in genome regulatory evolution after WGDs is not well understood. Here we used Atlantic salmon as a model system to explore how TE activity after the salmonid WGD ∼100 MYA shaped CRE evolution. We identified 55,080 putative TE-CREs using chromatin accessibility data from the liver and brain. Retroelements were both the dominant source of TE-CREs and had higher regulatory activity in MPRA experiments compared with DNA elements. A minority of TE subfamilies (16%) accounted for 46% of TE-CREs, but these "CRE superspreaders" were mostly active prior to the WGD. Analysis of individual TE insertions, however, revealed enrichment of TE-CREs originating from WGD-associated TE activity, particularly for the DTT (Tc1-Mariner) DNA elements. Furthermore, coexpression analyses supported the presence of TE-driven gene regulatory network evolution, including DTT elements active at the time of WGD. In conclusion, our study supports a scenario in which TE activity has been important in genome regulatory evolution, either through relaxed selective constraints or through strong selection to recalibrate optimal gene expression phenotypes, during a transient period following genome doubling.

Multiomic QTL mapping reveals phenotypic complexity of GWAS loci and prioritizes putative causal variants

Most GWAS loci are presumed to affect gene regulation; however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we map eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental-like tissues. Through colocalization, we annotate 10.4% (n = 540) of GWAS loci in 15 traits by QTL phenotype, temporal specificity, and complexity. We show that integration of chromatin QTLs results in a 2.3-fold higher annotation rate of GWAS loci because they capture distal GWAS loci missed by eQTLs, and that 5.4% (n = 13) of GWAS colocalizing eQTLs are early developmental specific. Finally, we utilize the iPSCORE multiomic QTLs to prioritize putative causal variants overlapping transcription factor motifs to elucidate the potential genetic underpinnings of 296 GWAS-QTL colocalizations.

Multiomic single-cell profiling identifies critical regulators of postnatal brain

Human brain development spans from embryogenesis to adulthood, with dynamic gene expression controlled by cell-type-specific cis-regulatory element activity and three-dimensional genome organization. To advance our understanding of postnatal brain development, we simultaneously profiled gene expression and chromatin accessibility in 101,924 single nuclei from four brain regions across ten donors, covering five key postnatal stages from infancy to late adulthood. Using this dataset and chromosome conformation capture data, we constructed enhancer-based gene regulatory networks to identify cell-type-specific regulators of brain development and interpret genome-wide association study loci for ten main brain disorders. Our analysis connected 2,318 cell-specific loci to 1,149 unique genes, representing 41% of loci linked to the investigated traits, and highlighted 55 genes influencing several disease phenotypes. Pseudotime analysis revealed distinct stages of postnatal oligodendrogenesis and their regulatory programs. These findings provide a comprehensive dataset of cell-type-specific gene regulation at critical timepoints in postnatal brain development.

IFNγ activates an immune-like regulatory network in the cardiac vascular endothelium

The regulatory mechanisms underlying the response to pro-inflammatory cytokines in cardiac diseases are poorly understood. Here, we use iPSC-derived cardiovascular progenitor cells (CVPCs) to model the response to interferon gamma (IFNγ) in human cardiac tissue. We generate RNA-seq and ATAC-seq for four CVPCs that were treated with IFNγ and compare them with paired untreated controls. Transcriptional differences after treatment show that IFNγ initiates an innate immune cell-like response, shifts the CVPC transcriptome toward coronary artery and aorta profiles, and stimulates expression of endothelial cell-specific genes. Analysis of the accessible chromatin shows that IFNγ is a potent chromatin remodeler and establishes an IRF-STAT immune-cell like regulatory network. Finally, we show that 11 GWAS risk variants for 8 common cardiac diseases overlap IFNγ-upregulated ATAC-seq peaks. Our findings reveal insights into IFNγ-induced activation of an immune-like regulatory network in human cardiac tissue and the potential role that regulatory elements in this pathway play in common cardiac diseases.

Epigenomic landscape of the human dorsal root ganglion: sex differences and transcriptional regulation of nociceptive genes

ABSTRACT

Cell states are influenced by the regulation of gene expression orchestrated by transcription factors capable of binding to accessible DNA regions. To uncover if sex differences exist in chromatin accessibility in the human dorsal root ganglion (hDRG), where nociceptive neurons innervating the body are found, we performed bulk and spatial assays for transposase-accessible chromatin technology followed by sequencing (ATAC-seq) from organ donors without a history of chronic pain. Using bulk ATAC-seq, we detected abundant sex differences in the hDRG. In women, differentially accessible regions (DARs) mapped mostly to the X chromosome, whereas in men, they mapped to autosomal genes. Hormone-responsive transcription factor binding motifs such as EGR1/3 were abundant within DARs in women, while JUN, FOS, and other activating protein 1 factor motifs were enriched in men, suggesting a higher activation state of cells compared with women. These observations were consistent with spatial ATAC-seq data. Furthermore, we validated that EGR1 expression is biased to female hDRG using RNAscope. In neurons, spatial ATAC-seq revealed higher chromatin accessibility in GABAergic, glutamatergic, and interferon-related genes in women and in Ca2+-signaling-related genes in men. Strikingly, XIST, responsible for inactivating 1 X chromosome by compacting it and maintaining at the periphery of the nucleus, was found to be highly dispersed in female neuronal nuclei. This is likely related to the higher chromatin accessibility in X in female hDRG neurons observed using both ATAC-seq approaches. We have documented baseline epigenomic sex differences in the hDRG which provide important descriptive information to test future hypotheses.

The chromatin accessibility landscape during early maize seed development

Cis-regulatory elements (CREs) are enriched in accessible chromatin regions (ACRs) of eukaryotes. Despite extensive research on genome-wide ACRs in various plant tissues, the global impact of these changes on developmental processes in maize seeds remains poorly understood. In this study, we employed the assay for transposase-accessible chromatin sequencing (ATAC-seq) to reveal the chromatin accessibility profile throughout the genome during the early stages of maize seed development. We identified a total of 37 952 to 59 887 high-quality ACRs in maize seeds at 0 to 8 days after pollination (DAP). Furthermore, we examined the correlation between the identified ACRs and gene expression. We observed a positive correlation between the open degree of promoter-ACRs and the expression of most genes. Moreover, we identified binding footprints of numerous transcription factors (TFs) within chromatin accessibility regions and revealed key TF families involved in different stages. Through the footprints of accessible chromatin regions, we predicted transcription factor regulatory networks during early maize embryo development. Additionally, we discovered that DNA sequence diversity was notably reduced at ACRs, yet trait-associated SNPs were more likely to be located within ACRs. We edited the ACR containing the trait-associated SNP of NKD1. Both NKD1pro-1 and NKD1pro-2 showed phenotypes corresponding to the trait-associated SNP. Our results suggest that alterations in chromatin accessibility play a crucial role in maize seed development and highlight the potential contribution of open chromatin regions to advancements in maize breeding.

Dysregulated intestinal nutrient absorption in obesity is associated with epigenomic alterations in epithelia

Obesity is an epidemic with myriad health effects, but little is understood regarding individual obese phenotypes and how they may respond to therapy. Epigenetic changes associated with obesity have been detected in blood, liver, pancreas, and adipose tissues. Previous work using human organoids found that dietary glucose hyperabsorption is a steadfast trait in cultures derived from some obese subjects, but detailed transcriptional or epigenomic features of the intestinal epithelia associated with this persistent phenotype are unknown. This study evaluated differentially expressed genes and relative chromatin accessibility in intestinal organoids established from donors classified as non-obese, obese, or obese hyperabsorptive by body mass index and glucose transport assays. Transcriptomic analysis indicated that obese hyperabsorptive subject organoids have significantly upregulated dietary nutrient absorption transcripts and downregulated type I interferon targets. Chromatin accessibility and transcription factor footprinting predicted that enhanced HNF4G binding may promote the obese hyperabsorption phenotype. Quantitative RT-PCR assessment in organoids representing a larger subject cohort suggested that intestinal epithelial expression of CUBN, GIP, SLC5A11, and SLC2A5 were highly correlated with hyperabsorption. Thus, the obese hyperabsorption phenotype was characterized by transcriptional changes that support increased nutrient uptake by intestinal epithelia, potentially driven by differentially accessible chromatin. Recognizing unique intestinal phenotypes in obesity provides a new perspective in considering therapeutic targets and options to manage the disease.

IL-23 tunes inflammatory functions of human mucosal-associated invariant T cells

IL-23 signaling plays a key role in the pathogenesis of chronic inflammatory and infectious diseases, yet the cellular targets and signaling pathways affected by this cytokine remain poorly understood. We show that IL-23 receptors are expressed on the large majority of human mucosal-associated invariant T (MAIT), but not of conventional T cells. Protein and transcriptional profiling at the population and single cell level demonstrates that stimulation with IL-23 or the structurally related cytokine IL-12 drives distinct functional profiles, revealing a high level of plasticity of MAIT cells. IL-23, in particular, affects key molecules and pathways related to autoimmunity and cytotoxic functions. Integrated analysis of transcriptomes and chromatin accessibility, supported by CRISPR-Cas9 mediated deletion, shows that AP-1 transcription factors constitute a key regulatory node of the IL-23 pathway in MAIT cells. In conclusion, our findings indicate that MAIT cells are key mediators of IL-23 functions in immunity to infections and chronic inflammatory diseases.

Binding domain mutations provide insight into CTCF's relationship with chromatin and its contribution to gene regulation

Here we used a series of CTCF mutations to explore CTCF's relationship with chromatin and its contribution to gene regulation. CTCF's impact depends on the genomic context of bound sites and the unique binding properties of WT and mutant CTCF proteins. Specifically, CTCF's signal strength is linked to changes in accessibility, and the ability to block cohesin is linked to its binding stability. Multivariate modelling reveals that both CTCF and accessibility contribute independently to cohesin binding and insulation, however CTCF signal strength has a stronger effect. CTCF and chromatin have a bidirectional relationship such that at CTCF sites, accessibility is reduced in a cohesin-dependent, mutant specific fashion. In addition, each mutant alters TF binding and accessibility in an indirect manner, changes which impart the most influence on rewiring transcriptional networks and the cell's ability to differentiate. Collectively, the mutant perturbations provide a rich resource for determining CTCF's site-specific effects.

Dynamic Chromatin Accessibility and Transcriptional Regulation in the Eyes of Red Tilapia (Oreochromis sp.) in Response to Wintering Stress

During wintering, red tilapia may develop variable black spots on their bodies, significantly reducing their market value. Understanding the mechanisms driving this phenomenon is essential for molecular improvements in body color. In this study, we investigated chromatin accessibility landscapes in the eyes of red tilapia with two distinct phenotypes (normal pure red and black spot) under wintering stress using ATAC-seq and RNA-seq analyses. We observed that approximately 32.7% of chromatin accessibility peaks were located in promoter regions, followed by intergenic regions (32.4%) and intronic regions (26.7%). One thousand two hundred twenty-nine differentially accessible regions (DARs) and 1448 differentially expressed genes (DEGs) were identified between the RS and DS groups. Notably, DEGs associated with melanin synthesis, including tyrp1, tyr, tyrp1b, pmela, slc24a5, and mlph, were significantly upregulated in the DS group, which aligns with the observed 1.85-fold increase in melanin content, compared to the RS group. 92 DEGs were associated with significant changes in chromatin accessibility between groups (R2 = 0.8059; p < 0.0001), indicating potential regulatory relationships. Interestingly, 23.92% of the DARs were located on the chromosome 3. Specifically, a 2.5-fold difference in average peak height on LG3: 11,215,273-11,217,225 were observed between DS and RS tilapia. In the region, transcription factors including HSF1 and HSF2 were identified as key regulators of chromatin structure and gene expression under wintering stress. Our findings reveal that dynamic chromatin accessibility in the eyes of red tilapia facilitates adaptation to wintering stress by regulating visual signaling, melanin production, and downstream pigmentation.

Phosphorylation of endothelial histone H3.3 serine 31 by PKN1 links flow-induced signaling to proatherogenic gene expression

Atherosclerotic lesions develop preferentially in arterial regions exposed to disturbed blood flow, where endothelial cells acquire an inflammatory phenotype. How disturbed flow induces endothelial cell inflammation is incompletely understood. Here we show that histone H3.3 phosphorylation at serine 31 (H3.3S31) regulates disturbed-flow-induced endothelial inflammation by allowing rapid induction of FOS and FOSB, required for inflammatory gene expression. We identified protein kinase N1 (PKN1) as the kinase responsible for disturbed-flow-induced H3.3S31 phosphorylation. Disturbed flow activates PKN1 in an integrin α5β1-dependent manner and induces its translocation into the nucleus, and PKN1 is also involved in the phosphorylation of the AP-1 transcription factor JUN. Mice with endothelium-specific PKN1 loss or endothelial expression of S31 phosphorylation-deficient H.3.3 mutants show reduced endothelial inflammation and disturbed-flow-induced vascular remodeling in vitro and in vivo. Together, we identified a pathway whereby disturbed flow through PKN1-mediated histone phosphorylation and FOS/FOSB induction promotes inflammatory gene expression and vascular inflammation.

Multiscale footprints reveal the organization of cis-regulatory elements

Cis-regulatory elements (CREs) control gene expression and are dynamic in their structure and function, reflecting changes in the composition of diverse effector proteins over time1. However, methods for measuring the organization of effector proteins at CREs across the genome are limited, hampering efforts to connect CRE structure to their function in cell fate and disease. Here we developed PRINT, a computational method that identifies footprints of DNA-protein interactions from bulk and single-cell chromatin accessibility data across multiple scales of protein size. Using these multiscale footprints, we created the seq2PRINT framework, which uses deep learning to allow precise inference of transcription factor and nucleosome binding and interprets regulatory logic at CREs. Applying seq2PRINT to single-cell chromatin accessibility data from human bone marrow, we observe sequential establishment and widening of CREs centred on pioneer factors across haematopoiesis. We further discover age-associated alterations in the structure of CREs in murine haematopoietic stem cells, including widespread reduction of nucleosome footprints and gain of de novo identified Ets composite motifs. Collectively, we establish a method for obtaining rich insights into DNA-binding protein dynamics from chromatin accessibility data, and reveal the architecture of regulatory elements across differentiation and ageing.

Oncofetal reprogramming drives phenotypic plasticity in WNT-dependent colorectal cancer

Targeting cancer stem cells (CSCs) is crucial for effective cancer treatment, yet resistance mechanisms to LGR5+ CSC depletion in WNT-driven colorectal cancer (CRC) remain elusive. In the present study, we revealed that mutant intestinal stem cells (SCs) depart from their canonical identity, traversing a dynamic phenotypic spectrum. This enhanced plasticity is initiated by oncofetal (OnF) reprogramming, driven by YAP and AP-1, with subsequent AP-1 hyperactivation promoting lineage infidelity. The retinoid X receptor serves as a gatekeeper of OnF reprogramming and its deregulation after adenomatous polyposis coli (APC) loss of function establishes an OnF 'memory' sustained by YAP and AP-1. Notably, the clinical significance of OnF and LGR5+ states in isolation is constrained by their functional redundancy. Although the canonical LGR5+ state is sensitive to the FOLFIRI regimen, an active OnF program correlates with resistance, supporting its role in driving drug-tolerant states. Targeting this program in combination with the current standard of care is pivotal for achieving effective and durable CRC treatment.

Chromatin accessibility and differentially expressed genes profiling in large yellow croaker (Larimichthys crocea) head kidney cells following iridovirus infection

Introduction

The large yellow croaker iridovirus (LYCIV) poses a significant threat to the aquaculture industry of Larimichthys crocea. Understanding the host defense response to LYCIV infection is crucial for developing effective strategies to mitigate its impact.

Methods

In this study, an epigenetic approach was employed to investigate dynamic changes in chromatin accessibility using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq). Additionally, RNA sequencing (RNA-seq) was used to analyze the expression pattern of immune response genes upon LYCIV infection.

Results

Substantial alterations in chromatin accessibility were observed, particularly in the regulatory regions of key immune-related genes. Significant changes in the expression of AP-1 transcription factors, including the Batf gene, were noted. CUT&Tag results revealed that AP-1 was significantly enriched in the open chromatin regions of cytokine genes, with Batf potentially regulating the cytokine genes LIF and CLCF1.

Discussion

These findings suggest that AP-1 may play a crucial role in the defense response against viral infection by modulating inflammatory cytokines and contributing to cellular inflammatory responses. This study provides a comprehensive analysis of the epigenomic landscape and gene expression regulation during iridovirus infection in L. crocea, offering valuable insights for breeding programs aimed at combating iridovirus infections.

Refining the cis-regulatory grammar learned by sequence-to-activity models by increasing model resolution

Chromatin accessibility can be measured genome-wide with ATAC-seq, enabling the discovery of regulatory regions that control gene expression and determine cell type. Deep genomic sequence-to-function (S2F) models link underlying genomic sequences to the measured chromatin state and identify motifs that regulate chromatin accessibility. Previously, we developed AI-TAC, a S2F model that predicts chromatin accessibility across 81 immune cell types and identifies sequence patterns that control their differential ATAC-seq signals. While AI-TAC provided valuable insights into the regulatory patterns that govern immune cell differentiation, later research established that ATAC-seq profiles (the distribution of Tn5 cuts) contain additional information about the exact location and strength of TF binding. To make use of this additional information, we developed bpAI-TAC, a multi-task neural network which models ATAC-seq at base-pair resolution across 90 immune cell types. We show that adding ATAC-profile information consistently improves predictions of differential chromatin accessibility. We also demonstrate that simultaneous learning of related cell types through multi-task modeling leads to better predictions than single task models. We then present a systematic framework for comparing how differences in model performance can be attributed to differences in what the model has learned. To understand what additional information bpAI-TAC gleans from ATAC-profiles, we use sequence attributions and identify motifs that have different effect sizes when trained on profiles. We conclude that modeling ATAC-seq at base-pair resolution enables the model to learn a more sensitive representation of the regulatory syntax that drives differences between immunocytes, and therefore will improve predictions of variant effects.

Genetic coupling of enhancer activity and connectivity in gene expression control

Gene enhancers often form long-range contacts with promoters, but it remains unclear if the activity of enhancers and their chromosomal contacts are mediated by the same DNA sequences and recruited factors. Here, we study the effects of expression quantitative trait loci (eQTLs) on enhancer activity and promoter contacts in primary monocytes isolated from 34 male individuals. Using eQTL-Capture Hi-C and a Bayesian approach considering both intra- and inter-individual variation, we initially detect 19 eQTLs associated with enhancer-eGene promoter contacts, most of which also associate with enhancer accessibility and activity. Capitalising on these shared effects, we devise a multi-modality Bayesian strategy, identifying 629 "trimodal QTLs" jointly associated with enhancer accessibility, eGene promoter contact, and gene expression. Causal mediation analysis and CRISPR interference reveal causal relationships between these three modalities. Many detected QTLs overlap disease susceptibility loci and influence the predicted binding of myeloid transcription factors, including SPI1, GABPB and STAT3. Additionally, a variant associated with PCK2 promoter contact directly disrupts a CTCF binding motif and impacts promoter insulation from downstream enhancers. Jointly, our findings suggest an inherent genetic coupling of enhancer activity and connectivity in gene expression control relevant to human disease and highlight the regulatory role of genetically determined chromatin boundaries.

Chromatin accessibility provides a window into the genetic etiology of human brain disease

Neuropsychiatric and neurodegenerative diseases have a significant genetic component. Risk variants often affect the noncoding genome, altering cis-regulatory elements (CREs) and chromatin structure, ultimately impacting gene expression. Chromatin accessibility profiling methods, especially assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), have been used to pinpoint disease-associated SNPs and link them to affected genes and cell types in the brain. The integration of single-cell technologies with genome-wide association studies (GWAS) and transcriptomic data has further advanced our understanding of cell-specific chromatin dynamics. This review discusses recent findings regarding the role played by chromatin accessibility in brain disease, highlighting the need for high-quality data and rigorous computational tools. Future directions include spatial chromatin studies and CRISPR-based functional validation to bridge genetic discovery and clinical applications, paving the way for targeted gene-regulatory therapies.

KAS-ATAC reveals the genome-wide single-stranded accessible chromatin landscape of the human genome

Gene regulation in most eukaryotes involves two fundamental processes: alterations in genome packaging by nucleosomes, with active cis-regulatory elements (CREs) generally characterized by open-chromatin configuration, and transcriptional activation. Mapping these physical properties and biochemical activities, through profiling chromatin accessibility and active transcription, is a key tool for understanding the logic and mechanisms of transcription and its regulation. However, the relationship between these two states has not been accessible to simultaneous measurement. To this end, we developed KAS-ATAC, a combination of the kethoxal-assisted ssDNA sequencing (KAS-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) methods for mapping single-stranded DNA (and thus active transcription) and chromatin accessibility, respectively, enabling the genome-wide identification of DNA fragments that are simultaneously accessible and contain ssDNA. We use KAS-ATAC to evaluate levels of active transcription over different CRE classes, to estimate absolute levels of transcribed accessible DNA over CREs, to map nucleosomal configurations associated with RNA polymerase activities, and to assess transcription factor association with transcribed DNA through transcription factor binding site (TFBS) footprinting. We observe lower levels of transcription over distal enhancers compared with promoters and distinct nucleosomal configurations around transcription initiation sites associated with active transcription. We find that most TFs associate equally with transcribed and nontranscribed DNA, but a few factors specifically do not exhibit footprints over ssDNA-containing fragments. We anticipate KAS-ATAC to continue to derive useful insights into chromatin organization and transcriptional regulation in other contexts in the future.

Common and specific gene regulatory programs in zebrafish caudal fin regeneration at single-cell resolution

Following amputation, zebrafish regenerate their injured caudal fin through lineage-restricted reprogramming. Although previous studies have charted various genetic and epigenetic dimensions of this process, the intricate gene regulatory programs shared by, or unique to, different regenerating cell types remain underinvestigated. Here, we mapped the regulatory landscape of fin regeneration by applying paired snRNA-seq and snATAC-seq on uninjured and regenerating fins. This map delineates the regulatory dynamics of predominant cell populations at multiple stages of regeneration. We observe a marked increase in the accessibility of chromatin regions associated with regenerative and developmental processes at 1 dpa, followed by a gradual closure across major cell types at later stages. This pattern is distinct from that of transcriptomic dynamics, which is characterized by several waves of gene upregulation and downregulation. We identified and in vivo validated cell-type-specific and position-specific regeneration-responsive enhancers and constructed regulatory networks by cell type and stage. Our single-cell resolution transcriptomic and chromatin accessibility map across regenerative stages provides new insights into regeneration regulatory mechanisms and serves as a valuable resource for the community.

Uncovering the chromatin-mediated transcriptional regulatory network governing cold stress responses in fish immune cells

Temperature fluctuations challenge ectothermic species, particularly tropical fish dependent on external temperatures for physiological regulation. However, the molecular mechanisms through which low-temperature stress impacts immune responses in these species, especially in relation to chromatin accessibility and epigenetic regulation, remain poorly understood. In this study, we investigate chromatin and transcriptional changes in the head kidney and thymus tissues of Nile tilapia (Oreochromis niloticus), a tropical fish of significant economic importance, under cold stress. By analyzing cis-regulatory elements in open chromatin regions and their associated transcription factors (TFs), we construct a comprehensive transcriptional regulatory network (TRN) governing immune responses, including DNA damage-induced apoptosis. Our analysis identifies 119 TFs within the TRN, with Stat1 emerging as a central hub exhibiting distinct binding dynamics under cold stress, as revealed by footprint analysis. Overexpression of Stat1 in immune cells leads to apoptosis and increases the expression of apoptosis-related genes, many of which contain Stat1 binding sites in their regulatory regions, emphasizing its critical role in immune cell survival during cold stress. These results provide insights into the transcriptional and epigenetic regulation of immune responses to cold stress in tilapia and highlight Stat1 as a promising target for enhancing cold tolerance in tropical fish species.

Schistosoma mansoni antigen induced innate immune memory features mitochondrial biogenesis and can be inhibited by ovarian produced hormones

We have previously identified that S. mansoni infection induces a unique form of myeloid training that protects male but not female mice from high fat diet induced disease. Here we demonstrate that ovarian derived hormones account for this sex specific difference. Ovariectomy of females prior to infection permits metabolic reprogramming of the myeloid lineage, with BMDM exhibiting carbon source flexibility for cellular respiration, and mice protected from systemic metabolic disease. The innate training phenotype of infection can be replicated by in vivo injection of SEA, and by exposure of bone marrow to SEA in culture prior to macrophage differentiation (Day 0). This protective phenotype is linked to increased chromatin accessibility of lipid and mitochondrial pathways in BMDM including Nrf1 and Tfam, as well as mitochondrial biogenesis. This work provides evidence that S. mansoni antigens induce a unique form of innate training inhibited by ovarian-derived hormones in females.

Application of computational algorithms for single-cell RNA-seq and ATAC-seq in neurodegenerative diseases

Recent advancements in single-cell technologies, including single-cell RNA sequencing (scRNA-seq) and Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), have greatly improved our insight into the epigenomic landscapes across various biological contexts and diseases. This paper reviews key computational tools and machine learning approaches that integrate scRNA-seq and scATAC-seq data to facilitate the alignment of transcriptomic data with chromatin accessibility profiles. Applying these integrated single-cell technologies in neurodegenerative diseases, such as Alzheimer's disease and Parkinson's disease, reveals how changes in chromatin accessibility and gene expression can illuminate pathogenic mechanisms and identify potential therapeutic targets. Despite facing challenges like data sparsity and computational demands, ongoing enhancements in scATAC-seq and scRNA-seq technologies, along with better analytical methods, continue to expand their applications. These advancements promise to revolutionize our approach to medical research and clinical diagnostics, offering a comprehensive view of cellular function and disease pathology.

ChromBPNet: bias factorized, base-resolution deep learning models of chromatin accessibility reveal cis-regulatory sequence syntax, transcription factor footprints and regulatory variants

Despite extensive mapping of cis-regulatory elements (cREs) across cellular contexts with chromatin accessibility assays, the sequence syntax and genetic variants that regulate transcription factor (TF) binding and chromatin accessibility at context-specific cREs remain elusive. We introduce ChromBPNet, a deep learning DNA sequence model of base-resolution accessibility profiles that detects, learns and deconvolves assay-specific enzyme biases from regulatory sequence determinants of accessibility, enabling robust discovery of compact TF motif lexicons, cooperative motif syntax and precision footprints across assays and sequencing depths. Extensive benchmarks show that ChromBPNet, despite its lightweight design, is competitive with much larger contemporary models at predicting variant effects on chromatin accessibility, pioneer TF binding and reporter activity across assays, cell contexts and ancestry, while providing interpretation of disrupted regulatory syntax. ChromBPNet also helps prioritize and interpret regulatory variants that influence complex traits and rare diseases, thereby providing a powerful lens to decode regulatory DNA and genetic variation.