ATAC-seq: A Method for Assaying Chromatin Accessibility Genome-Wide

Spatial chromatin accessibility sequencing resolves high-order spatial interactions of epigenomic markers

As the genome is organized into a three-dimensional structure in intracellular space, epigenomic information also has a complex spatial arrangement. However, most epigenetic studies describe locations of methylation marks, chromatin accessibility regions, and histone modifications in the horizontal dimension. Proper spatial epigenomic information has rarely been obtained. In this study, we designed spatial chromatin accessibility sequencing (SCA-seq) to resolve the genome conformation by capturing the epigenetic information in single-molecular resolution while simultaneously resolving the genome conformation. Using SCA-seq, we are able to examine the spatial interaction of chromatin accessibility (e.g. enhancer-promoter contacts), CpG island methylation, and spatial insulating functions of the CCCTC-binding factor. We demonstrate that SCA-seq paves the way to explore the mechanism of epigenetic interactions and extends our knowledge in 3D packaging of DNA in the nucleus.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis)

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Epigenetic reshaping through damage: promoting cell fate transition by BrdU and IdU incorporation

BACKGROUND

Thymidine analogs have long been recognized for their ability to randomly incorporate into DNA. However, the precise mechanisms through which thymidine analogs facilitate cell fate transition remains unclear.

RESULTS

Here, we discovered a strong correlation between the dosage dependence of thymidine analogs and their ability to overcome reprogramming barrier. The extraembryonic endoderm (XEN) state seems to be a cell's selective response to DNA damage repair (DDR), offering a shortcut to overcome reprogramming barriers. Meanwhile, we found that homologous recombination repair (HRR) pathway causes an overall epigenetic reshaping of cells and enabling them to overcome greater barriers. This response leads to the creation of a hypomethylated environment, which facilitates the transition of cell fate in various reprogramming systems. We term this mechanism as Epigenetic Reshaping through Damage (ERD).

CONCLUSION

Overall, our study finds that BrdU/IdU can activate the DNA damage repair pathway (HRR), leading to increased histone acetylation and genome-wide DNA demethylation, regulating somatic cell reprogramming. This offers valuable insights into mechanisms underlying cell fate transition.

Epstein-Barr virus induces germinal center light zone chromatin architecture and promotes survival through enhancer looping at the BCL2A1 locus

IMPORTANCE

Epstein-Barr virus has evolved with its human host leading to an intimate relationship where infection of antibody-producing B cells mimics the process by which these cells normally recognize foreign antigens and become activated. Virtually everyone in the world is infected by adulthood and controls this virus pushing it into life-long latency. However, immune-suppressed individuals are at high risk for EBV+ cancers. Here, we isolated B cells from tonsils and compare the underlying molecular genetic differences between these cells and those infected with EBV. We find similar regulatory mechanism for expression of an important cellular protein that enables B cells to survive in lymphoid tissue. These findings link an underlying relationship at the molecular level between EBV-infected B cells in vitro with normally activated B cells in vivo. Our studies also characterize the role of a key viral control mechanism for B cell survival involved in long-term infection.

Chromatin activity identifies differential gene regulation across human ancestries

BACKGROUND

Current evidence suggests that cis-regulatory elements controlling gene expression may be the predominant target of natural selection in humans and other species. Detecting selection acting on these elements is critical to understanding evolution but remains challenging because we do not know which mutations will affect gene regulation.

RESULTS

To address this, we devise an approach to search for lineage-specific selection on three critical steps in transcriptional regulation: chromatin activity, transcription factor binding, and chromosomal looping. Applying this approach to lymphoblastoid cells from 831 individuals of either European or African descent, we find strong signals of differential chromatin activity linked to gene expression differences between ancestries in numerous contexts, but no evidence of functional differences in chromosomal looping. Moreover, we show that enhancers rather than promoters display the strongest signs of selection associated with sites of differential transcription factor binding.

CONCLUSIONS

Overall, our study indicates that some cis-regulatory adaptation may be more easily detected at the level of chromatin than DNA sequence. This work provides a vast resource of genomic interaction data from diverse human populations and establishes a novel selection test that will benefit future study of regulatory evolution in humans and other species.

Activity-Dependent Transcriptional Program in NGN2+ Neurons Enriched for Genetic Risk for Brain-Related Disorders

BACKGROUND

Converging evidence from large-scale genetic and postmortem studies highlights the role of aberrant neurotransmission and genetic regulation in brain-related disorders. However, identifying neuronal activity-regulated transcriptional programs in the human brain and understanding how changes contribute to disease remain challenging.

METHODS

To better understand how the activity-dependent regulome contributes to risk for brain-related disorders, we profiled the transcriptomic and epigenomic changes following neuronal depolarization in human induced pluripotent stem cell-derived glutamatergic neurons (NGN2) from 6 patients with schizophrenia and 5 control participants.

RESULTS

Multiomic data integration associated global patterns of chromatin accessibility with gene expression and identified enhancer-promoter interactions in glutamatergic neurons. Within 1 hour of potassium chloride-induced depolarization, independent of diagnosis, glutamatergic neurons displayed substantial activity-dependent changes in the expression of genes regulating synaptic function. Depolarization-induced changes in the regulome revealed significant heritability enrichment for schizophrenia and Parkinson's disease, adding to mounting evidence that sequence variation within activation-dependent regulatory elements contributes to the genetic risk for brain-related disorders. Gene coexpression network analysis elucidated interactions among activity-dependent and disease-associated genes and pointed to a key driver (NAV3) that interacted with multiple genes involved in axon guidance.

CONCLUSIONS

Overall, we demonstrated that deciphering the activity-dependent regulome in glutamatergic neurons reveals novel targets for advanced diagnosis and therapy.

Polygenic risk for schizophrenia converges on alternative polyadenylation as molecular mechanism underlying synaptic impairment

Schizophrenia (SCZ) is a genetically heterogenous psychiatric disorder of highly polygenic nature. Correlative evidence from genetic studies indicate that the aggregated effects of distinct genetic risk factor combinations found in each patient converge onto common molecular mechanisms. To prove this on a functional level, we employed a reductionistic cellular model system for polygenic risk by differentiating induced pluripotent stem cells (iPSCs) from 104 individuals with high polygenic risk load and controls into cortical glutamatergic neurons (iNs). Multi-omics profiling identified widespread differences in alternative polyadenylation (APA) in the 3' untranslated region of many synaptic transcripts between iNs from SCZ patients and healthy donors. On the cellular level, 3'APA was associated with a reduction in synaptic density of iNs. Importantly, differential APA was largely conserved between postmortem human prefrontal cortex from SCZ patients and healthy donors, and strongly enriched for transcripts related to synapse biology. 3'APA was highly correlated with SCZ polygenic risk and affected genes were significantly enriched for SCZ associated common genetic variation. Integrative functional genomic analysis identified the RNA binding protein and SCZ GWAS risk gene PTBP2 as a critical trans-acting factor mediating 3'APA of synaptic genes in SCZ subjects. Functional characterization of PTBP2 in iNs confirmed its key role in 3'APA of synaptic transcripts and regulation of synapse density. Jointly, our findings show that the aggregated effects of polygenic risk converge on 3'APA as one common molecular mechanism that underlies synaptic impairments in SCZ.

WRN loss accelerates abnormal adipocyte metabolism in Werner syndrome

BACKGROUND

Metabolic dysfunction is one of the main symptoms of Werner syndrome (WS); however, the underlying mechanisms remain unclear. Here, we report that loss of WRN accelerates adipogenesis at an early stage both in vitro (stem cells) and in vivo (zebrafish). Moreover, WRN depletion causes a transient upregulation of late-stage of adipocyte-specific genes at an early stage.

METHODS

In an in vivo study, we generated wrn-/- mutant zebrafish and performed histological stain and Oil Red O staining to assess the fat metabolism. In an in vitro study, we used RNA-seq and ATAC-seq to profile the transcriptional features and chromatin accessibility in WRN depleted adipocytes. Moreover, we performed ChIP-seq to further study the regulatory mechanisms of metabolic dysfunction in WS.

RESULTS

Our findings show that mechanistically WRN deficiency causes SMARCA5 upregulation. SMARCA5 is crucial in chromatin remodeling and gene regulation. Additionally, rescuing WRN could normalize SMARCA5 expression and adipocyte differentiation. Moreover, we find that nicotinamide riboside (NR) supplementation restores adipocyte metabolism in both stem cells and zebrafish models.

CONCLUSIONS

Our findings unravel a new mechanism for the influence of WRN in the early stage of adipogenesis and provide a possible treatment for metabolic dysfunction in WS. These data provide promising insights into potential therapeutics for ageing and ageing-related diseases.

Disease-associated astrocyte epigenetic memory promotes CNS pathology

Astrocytes play important roles in the central nervous system (CNS) physiology and pathology. Indeed, astrocyte subsets defined by specific transcriptional activation states contribute to the pathology of neurologic diseases, including multiple sclerosis (MS) and its pre-clinical model experimental autoimmune encephalomyelitis (EAE) 1-8 . However, little is known about the stability of these disease-associated astrocyte subsets, their regulation, and whether they integrate past stimulation events to respond to subsequent challenges. Here, we describe the identification of an epigenetically controlled memory astrocyte subset which exhibits exacerbated pro-inflammatory responses upon re-challenge. Specifically, using a combination of single-cell RNA sequencing (scRNA-seq), assay for transposase-accessible chromatin with sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing (ChIP-seq), focused interrogation of cells by nucleic acid detection and sequencing (FIND-seq), and cell-specific in vivo CRISPR/Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) used by the histone acetyltransferase p300 to control chromatin accessibility. ACLY + p300 + memory astrocytes are increased in acute and chronic EAE models; the genetic targeting of ACLY + p300 + astrocytes using CRISPR/Cas9 ameliorated EAE. We also detected responses consistent with a pro-inflammatory memory phenotype in human astrocytes in vitro ; scRNA-seq and immunohistochemistry studies detected increased ACLY + p300 + astrocytes in chronic MS lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, MS. These findings may guide novel therapeutic approaches for MS and other neurologic diseases.

CellSTAR: a comprehensive resource for single-cell transcriptomic annotation

Large-scale studies of single-cell sequencing and biological experiments have successfully revealed expression patterns that distinguish different cell types in tissues, emphasizing the importance of studying cellular heterogeneity and accurately annotating cell types. Analysis of gene expression profiles in these experiments provides two essential types of data for cell type annotation: annotated references and canonical markers. In this study, the first comprehensive database of single-cell transcriptomic annotation resource (CellSTAR) was thus developed. It is unique in (a) offering the comprehensive expertly annotated reference data for annotating hundreds of cell types for the first time and (b) enabling the collective consideration of reference data and marker genes by incorporating tens of thousands of markers. Given its unique features, CellSTAR is expected to attract broad research interests from the technological innovations in single-cell transcriptomics, the studies of cellular heterogeneity & dynamics, and so on. It is now publicly accessible without any login requirement at: https://idrblab.org/cellstar.

CLIC3 interacts with NAT10 to inhibit N4-acetylcytidine modification of p21 mRNA and promote bladder cancer progression

Chromatin accessibility plays important roles in revealing the regulatory networks of gene expression, while its application in bladder cancer is yet to be fully elucidated. Chloride intracellular channel 3 (CLIC3) protein has been reported to be associated with the progression of some tumors, whereas the specific mechanism of CLIC3 in tumor remains unclear. Here, we screened for key genes in bladder cancer through the identification of transcription factor binding site clustered region (TFCR) on the basis of chromatin accessibility and TF motif. CLIC3 was identified by joint profiling of chromatin accessibility data with TCGA database. Clinically, CLIC3 expression was significantly elevated in bladder cancer and was negatively correlated with patient survival. CLIC3 promoted the proliferation of bladder cancer cells by reducing p21 expression in vitro and in vivo. Mechanistically, CLIC3 interacted with NAT10 and inhibited the function of NAT10, resulting in the downregulation of ac4C modification and stability of p21 mRNA. Overall, these findings uncover an novel mechanism of mRNA ac4C modification and CLIC3 may act as a potential therapeutic target for bladder cancer.

Leveraging Tissue-Specific Enhancer-Target Gene Regulatory Networks Identifies Enhancer Somatic Mutations That Functionally Impact Lung Cancer

Enhancers are noncoding regulatory DNA regions that modulate the transcription of target genes, often over large distances along with the genomic sequence. Enhancer alterations have been associated with various pathological conditions, including cancer. However, the identification and characterization of somatic mutations in noncoding regulatory regions with a functional effect on tumorigenesis and prognosis remain a major challenge. Here, we present a strategy for detecting and characterizing enhancer mutations in a genome-wide analysis of patient cohorts, across three lung cancer subtypes. Lung tissue-specific enhancers were defined by integrating experimental data and public epigenomic profiles, and the genome-wide enhancer-target gene regulatory network of lung cells was constructed by integrating chromatin three-dimensional architecture data. Lung cancers possessed a similar mutation burden at tissue-specific enhancers and exons but with differences in their mutation signatures. Functionally relevant alterations were prioritized on the basis of the pathway-level integration of the effect of a mutation and the frequency of mutations on individual enhancers. The genes enriched for mutated enhancers converged on the regulation of key biological processes and pathways relevant to tumor biology. Recurrent mutations in individual enhancers also affected the expression of target genes, with potential relevance for patient prognosis. Together, these findings show that noncoding regulatory mutations have a potential relevance for cancer pathogenesis and can be exploited for patient classification.

SIGNIFICANCE

Mapping enhancer-target gene regulatory interactions and analyzing enhancer mutations at the level of their target genes and pathways reveal convergence of recurrent enhancer mutations on biological processes involved in tumorigenesis and prognosis.

Transcriptional and epigenetic dysregulation impairs generation of proliferative neural stem and progenitor cells during brain aging

The decline in stem cell function during aging may affect the regenerative capacity of mammalian organisms; however, the gene regulatory mechanism underlying this decline remains unclear. Here we show that the aging of neural stem and progenitor cells (NSPCs) in the male mouse brain is characterized by a decrease in the generation efficacy of proliferative NSPCs rather than the changes in lineage specificity of NSPCs. We reveal that the downregulation of age-dependent genes in NSPCs drives cell aging by decreasing the population of actively proliferating NSPCs while increasing the expression of quiescence markers. We found that epigenetic deregulation of the MLL complex at promoters leads to transcriptional inactivation of age-dependent genes, highlighting the importance of the dynamic interaction between histone modifiers and gene regulatory elements in regulating transcriptional program of aging cells. Our study sheds light on the key intrinsic mechanisms driving stem cell aging through epigenetic regulators and identifies potential rejuvenation targets that could restore the function of aging stem cells.

Epigenetic Regulation of Ameloblast Differentiation by HMGN Proteins

Dental enamel formation is coordinated by ameloblast differentiation, production of enamel matrix proteins, and crystal growth. The factors regulating ameloblast differentiation are not fully understood. Here we show that the high mobility group N (HMGN) nucleosomal binding proteins modulate the rate of ameloblast differentiation and enamel formation. We found that HMGN1 and HMGN2 proteins are downregulated during mouse ameloblast differentiation. Genetically altered mice lacking HMGN1 and HMGN2 proteins show faster ameloblast differentiation and a higher rate of enamel deposition in mice molars and incisors. In vitro differentiation of induced pluripotent stem cells to dental epithelium cells showed that HMGN proteins modulate the expression and chromatin accessibility of ameloblast-specific genes and affect the binding of transcription factors epiprofin and PITX2 to ameloblast-specific genes. Our results suggest that HMGN proteins regulate ameloblast differentiation and enamel mineralization by modulating lineage-specific chromatin accessibility and transcription factor binding to ameloblast regulatory sites.

Single-Nucleus ATAC-seq for Mapping Chromatin Accessibility in Individual Cells of Murine Hearts

During the last decade a wide range of single-cell and single-nucleus next-generation sequencing techniques have been developed, which revolutionized detection of rare cell populations, enabling creation of comprehensive cell atlases of complex organs and tissues. State-of-the-art methods do not only allow classical transcriptomics of individual cells but also comprise a number of epigenetic approaches, including assessment of chromatin accessibility by single-nucleus Assay for Transposase Accessible Chromatin ATAC-seq (snATAC-seq). The snATAC-seq assay detects "open chromatin," a term for low nucleosome occupancy of genomic regions, which is a prerequisite for effective transcription factor binding. Information about open chromatin at the single-nucleus level helps to recognize epigenetic changes, sometimes before transcription of respective genes occurs. snATAC-seq detects cellular heterogeneity in otherwise still transcriptionally and/or morphologically homogeneous cell populations. Chromatin accessibility assays may be used to detect epigenetic changes in cardiac lineages during heart development, chromatin landscape changes during aging, and epigenetic alterations in heart diseases. Here, we provide an optimized protocol for snATAC-seq of murine hearts. We describe isolation of single nuclei from snap-frozen hearts, provide hints for preparation of libraries suitable for snATAC-seq next-generation sequencing (NGS) using the Chromium 10× platform, and give general recommendations for downstream analysis using conventional bioinformatic pipelines and packages. The protocol should serve as a beginner's guide to generate high-quality snATAC-seq datasets and to perform chromatin accessibility analysis of individual heart-derived cell nuclei.

Chromatin and DNA Dynamics in Mouse Models of Liver Cancers

In recent years, important efforts have been made to understand how the expression of a specific gene repertoire correlates with chromatin accessibility, histone mark deposition, as well as with chromatin looping establishing connectivity with regulatory regions. The emergence of new techniques for genome-wide analyses and their progressive optimization to work on low amounts of material allows the scientific community to obtain an integrated view of transcriptional landscapes in physiology and disease. Here, we describe our own experience aiming at correlating the TCF-4/β-catenin cistrome during liver tumorigenesis with chromatin remodeling, histone mark modifications, and long-distance DNA looping.

Transcriptional regulation of FACT involves Coordination of chromatin accessibility and CTCF binding

Histone chaperone FACT (facilitates chromatin transcription) is well known to promote chromatin recovery during transcription. However, the mechanism how FACT regulates genome-wide chromatin accessibility and transcription factor binding has not been fully elucidated. Through loss-of-function studies, we show here that FACT component Ssrp1 is required for DNA replication and DNA damage repair and is also essential for progression of cell phase transition and cell proliferation in mouse embryonic fibroblast cells. On the molecular level, absence of the Ssrp1 leads to increased chromatin accessibility, enhanced CTCF binding, and a remarkable change in dynamic range of gene expression. Our study thus unequivocally uncovers a unique mechanism by which FACT complex regulates transcription by coordinating genome-wide chromatin accessibility and CTCF binding.

Chromatin remodeling analysis reveals the RdDM pathway responds to low-phosphorus stress in maize

Chromatin in eukaryotes folds into a complex three-dimensional (3D) structure that is essential for controlling gene expression and cellular function and is dynamically regulated in biological processes. Studies on plant phosphorus signaling have concentrated on single genes and gene interactions. It is critical to expand the existing signaling pathway in terms of its 3D structure. In this study, low-Pi treatment led to greater chromatin volume. Furthermore, low-Pi stress increased the insulation score and the number of TAD-like domains, but the effects on the A/B compartment were not obvious. The methylation levels of target sites (hereafter as RdDM levels) peaked at specific TAD-like boundaries, whereas RdDM peak levels at conserved TAD-like boundaries shifted and decreased sharply. The distribution pattern of RdDM sites originating from the Helitron transposons matched that of genome-wide RdDM sites near TAD-like boundaries. RdDM pathway genes were upregulated in the middle or early stages and downregulated in the later stages under low-Pi conditions. The RdDM pathway mutant ddm1a showed increased tolerance to low-Pi stress, with shortened and thickened roots contributing to higher Pi uptake from the shallow soil layer. ChIP-seq results revealed that ZmDDM1A could bind to Pi- and root development-related genes. Strong associations were found between interacting genes in significantly different chromatin-interaction regions and root traits. These findings not only expand the mechanisms by which plants respond to low-Pi stress through the RdDM pathway but also offer a crucial framework for the analysis of biological issues using 3D genomics.

Visualization of the landscape of the read alignment shape of ATAC-seq data using Hellinger distance metric

Assay for Transposase-Accessible Chromatin using high-throughput sequencing (ATAC-seq) is the popular technique using next-generation sequencing to measure chromatin accessibility and identify open chromatin regions. While read alignment shape information of next-generation sequencing data with intensity information has been used in various bioinformatics methods, few studies have focused on pure shape information alone. In this study, we investigated what types of ATAC-seq read alignment shapes are observed for the promoter region and whether the pure shape information was related or unrelated to other gene features. We introduced a novel concept and pipeline for handling the pure shape information of NGS data as probability distributions and quantifying their dissimilarities by information theory. Based on this concept, we demonstrate that the pure shape information of ATAC-seq data is correlated with chromatin openness and some gene characteristics. On the other hand, it is suggested that the pure information of ATAC-seq read alignment shape is unlikely to contain additional information to explain differences in RNA expression. Our study suggests that viewing the read alignment shape of NGS data as probability distributions enables us to capture the characteristics of the genome-wide landscape of such data in a non-parametric manner.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Revisiting chromatin packaging in mouse sperm

Mammalian sperm show an unusual and heavily compacted genomic packaging state. In addition to its role in organizing the compact and hydrodynamic sperm head, it has been proposed that sperm chromatin architecture helps to program gene expression in the early embryo. Scores of genome-wide surveys in sperm have reported patterns of chromatin accessibility, nucleosome localization, histone modification, and chromosome folding. Here, we revisit these studies in light of recent reports that sperm obtained from the mouse epididymis are contaminated with low levels of cell-free chromatin. In the absence of proper sperm lysis, we readily recapitulate multiple prominent genome-wide surveys of sperm chromatin, suggesting that these profiles primarily reflect contaminating cell-free chromatin. Removal of cell-free DNA, and appropriate lysis conditions, are together required to reveal a sperm chromatin state distinct from most previous reports. Using ATAC-seq to explore relatively accessible genomic loci, we identify a landscape of open loci associated with early development and transcriptional control. Histone modification and chromosome folding profiles also strongly support the hypothesis that prior studies suffer from contamination, but technical challenges associated with reliably preserving the architecture of the compacted sperm head prevent us from confidently assaying true localization patterns for these epigenetic marks. Together, our studies show that our knowledge of chromosome packaging in mammalian sperm remains largely incomplete, and motivate future efforts to more accurately characterize genome organization in mature sperm.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads

Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.

Selection on the promoter regions plays an important role in complex traits during duck domestication

BACKGROUND

Identifying the key factors that underlie complex traits during domestication is a great challenge for evolutionary and biological studies. In addition to the protein-coding region differences caused by variants, a large number of variants are located in the noncoding regions containing multiple types of regulatory elements. However, the roles of accumulated variants in gene regulatory elements during duck domestication and economic trait improvement are poorly understood.

RESULTS

We constructed a genomics, transcriptomics, and epigenomics map of the duck genome and assessed the evolutionary forces that have been in play across the whole genome during domestication. In total, 304 (42.94%) gene promoters have been specifically selected in Pekin duck among all selected genes. Joint multi-omics analysis reveals that 218 genes (72.01%) with selected promoters are located in open and active chromatin, and 267 genes (87.83%) with selected promoters were highly and differentially expressed in domestic trait-related tissues. One important candidate gene ELOVL3, with a strong signature of differentiation on the core promoter region, is known to regulate fatty acid elongation. Functional experiments showed that the nearly fixed variants in the top selected ELOVL3 promoter in Pekin duck decreased binding ability with HLF and increased gene expression, with the overexpression of ELOVL3 able to increase lipid deposition and unsaturated fatty acid enrichment.

CONCLUSIONS

This study presents genome resequencing, RNA-Seq, Hi-C, and ATAC-Seq data of mallard and Pekin duck, showing that selection of the gene promoter region plays an important role in gene expression and phenotypic changes during domestication and highlights that the variants of the ELOVL3 promoter may have multiple effects on fat and long-chain fatty acid content in ducks.

Super-enhancers include classical enhancers and facilitators to fully activate gene expression

Super-enhancers are compound regulatory elements that control expression of key cell identity genes. They recruit high levels of tissue-specific transcription factors and co-activators such as the Mediator complex and contact target gene promoters with high frequency. Most super-enhancers contain multiple constituent regulatory elements, but it is unclear whether these elements have distinct roles in activating target gene expression. Here, by rebuilding the endogenous multipartite α-globin super-enhancer, we show that it contains bioinformatically equivalent but functionally distinct element types: classical enhancers and facilitator elements. Facilitators have no intrinsic enhancer activity, yet in their absence, classical enhancers are unable to fully upregulate their target genes. Without facilitators, classical enhancers exhibit reduced Mediator recruitment, enhancer RNA transcription, and enhancer-promoter interactions. Facilitators are interchangeable but display functional hierarchy based on their position within a multipartite enhancer. Facilitators thus play an important role in potentiating the activity of classical enhancers and ensuring robust activation of target genes.

The role of enhancers in psoriasis and atopic dermatitis

Regulatory elements, particularly enhancers, play a crucial role in disease susceptibility and progression. Enhancers are DNA sequences that activate gene expression and can be affected by epigenetic modifications, interactions with transcription factors (TFs) or changes to the enhancer DNA sequence itself. Altered enhancer activity impacts gene expression and contributes to disease. In this review, we define enhancers and the experimental techniques used to identify and characterize them. We also discuss recent studies that examine how enhancers contribute to atopic dermatitis (AD) and psoriasis. Articles in the PubMed database were identified (from 1 January 2010 to 28 February 2023) that were relevant to enhancer variants, enhancer-associated TFs and enhancer histone modifications in psoriasis or AD. Most enhancers associated with these conditions regulate genes affecting epidermal homeostasis or immune function. These discoveries present potential therapeutic targets to complement existing treatment options for AD and psoriasis.

PRDM3/16 Regulate Chromatin Accessibility Required for NKX2-1 Mediated Alveolar Epithelial Differentiation and Function

Differential chromatin accessibility accompanies and mediates transcriptional control of diverse cell fates and their differentiation during embryogenesis. While the critical role of NKX2-1 and its transcriptional targets in lung morphogenesis and pulmonary epithelial cell differentiation is increasingly known, mechanisms by which chromatin accessibility alters the epigenetic landscape and how NKX2-1 interacts with other co-activators required for alveolar epithelial cell differentiation and function are not well understood. Here, we demonstrate that the paired domain zinc finger transcriptional regulators PRDM3 and PRDM16 regulate chromatin accessibility to mediate cell differentiation decisions during lung morphogenesis. Combined deletion of Prdm3 and Prdm16 in early lung endoderm caused perinatal lethality due to respiratory failure from loss of AT2 cell function. Prdm3/16 deletion led to the accumulation of partially differentiated AT1 cells and loss of AT2 cells. Combination of single cell RNA-seq, bulk ATAC-seq, and CUT&RUN demonstrated that PRDM3 and PRDM16 enhanced chromatin accessibility at NKX2-1 transcriptional targets in peripheral epithelial cells, all three factors binding together at a multitude of cell-type specific cis-active DNA elements. Network analysis demonstrated that PRDM3/16 regulated genes critical for perinatal AT2 cell differentiation, surfactant homeostasis, and innate host defense. Lineage specific deletion of PRDM3/16 in AT2 cells led to lineage infidelity, with PRDM3/16 null cells acquiring partial AT1 fate. Together, these data demonstrate that NKX2-1-dependent regulation of alveolar epithelial cell differentiation is mediated by epigenomic modulation via PRDM3/16.

HIV-1 active and latent infections induce disparate chromatin reorganization and transcriptional regulation of mRNAs and lncRNAs in SupT1 cells

IMPORTANCE

HIV-1 infection of T-lymphocytes depends on co-opting cellular transcriptional and translational machineries for viral replication. This requires significant changes in the cellular microenvironment. We have characterized and compared the changes in cellular chromatin structures as well as gene expression landscapes in T cells that are either actively or latently infected with HIV-1. Our results reveal that chromatin accessibility and expression of both protein-coding mRNAs and non-coding lncRNAs are uniquely regulated in HIV-1-infected T cells, depending on whether the virus is actively transcribing or remains in a transcriptionally silent, latent state. HIV-1 latent infection elicits more robust changes in the cellular chromatin organization than active viral infection. Our analysis also identifies the effects of such epigenomic changes on the cellular gene expression and subsequent biological pathways. This study comprehensively characterizes the cellular epigenomic and transcriptomic states that support active and latent HIV-1 infection in an in vitro model of SupT1 cells.

Dynamics of chromatin accessibility and genome wide control of desiccation tolerance in the resurrection plant Haberlea rhodopensis

BACKGROUND

Drought is one of the main consequences of global climate change and this problem is expected to intensify in the future. Resurrection plants evolved the ability to withstand the negative impact of long periods of almost complete desiccation and to recover at rewatering. In this respect, many physiological, transcriptomic, proteomic and genomic investigations have been performed in recent years, however, few epigenetic control studies have been performed on these valuable desiccation-tolerant plants so far.

RESULTS

In the present study, for the first time for resurrection plants we provide evidences about the differential chromatin accessibility of Haberlea rhodopensis during desiccation stress by ATAC-seq (Assay for Transposase Accessible Chromatin with high-throughput sequencing). Based on gene similarity between species, we used the available genome of the closely related resurrection plant Dorcoceras hygrometricum to identify approximately nine hundred transposase hypersensitive sites (THSs) in H. rhodopensis. The majority of them corresponds to proximal and distal regulatory elements of different genes involved in photosynthesis, carbon metabolism, synthesis of secondary metabolites, cell signalling and transcriptional regulation, cell growth, cell wall, stomata conditioning, chaperons, oxidative stress, autophagy and others. Various types of binding motifs recognized by several families of transcription factors have been enriched from the THSs found in different stages of drought. Further, we used the previously published RNA-seq data from H. rhodopensis to evaluate the expression of transcription factors putatively interacting with the enriched motifs, and the potential correlation between the identified THS and the expression of their corresponding genes.

CONCLUSIONS

These results provide a blueprint for investigating the epigenetic regulation of desiccation tolerance in resurrection plant H. rhodopensis and comparative genomics between resurrection and non-resurrection species with available genome information.

Chromatin state transitions in the Drosophila intestinal lineage identify principles of cell-type specification

Tissue homeostasis relies on rewiring of stem cell transcriptional programs into those of differentiated cells. Here, we investigate changes in chromatin occurring in a bipotent adult stem cells. Combining mapping of chromatin-associated factors with statistical modeling, we identify genome-wide transitions during differentiation in the adult Drosophila intestinal stem cell (ISC) lineage. Active, stem-cell-enriched genes transition to a repressive heterochromatin protein-1-enriched state more prominently in enteroendocrine cells (EEs) than in enterocytes (ECs), in which the histone H1-enriched Black state is preeminent. In contrast, terminal differentiation genes associated with metabolic functions follow a common path from a repressive, primed, histone H1-enriched Black state in ISCs to active chromatin states in EE and EC cells. Furthermore, we find that lineage priming has an important function in adult ISCs, and we identify histone H1 as a mediator of this process. These data define underlying principles of chromatin changes during adult multipotent stem cell differentiation.

Aged intestinal stem cells propagate cell-intrinsic sources of inflammaging in mice

Low-grade chronic inflammation is a hallmark of ageing, associated with impaired tissue function and disease development. However, how cell-intrinsic and -extrinsic factors collectively establish this phenotype, termed inflammaging, remains poorly understood. We addressed this question in the mouse intestinal epithelium, using mouse organoid cultures to dissect stem cell-intrinsic and -extrinsic sources of inflammaging. At the single-cell level, we found that inflammaging is established differently along the crypt-villus axis, with aged intestinal stem cells (ISCs) strongly upregulating major histocompatibility complex class II (MHC-II) genes. Importantly, the inflammaging phenotype was stably propagated by aged ISCs in organoid cultures and associated with increased chromatin accessibility at inflammation-associated loci in vivo and ex vivo, indicating cell-intrinsic inflammatory memory. Mechanistically, we show that the expression of inflammatory genes is dependent on STAT1 signaling. Together, our data identify that intestinal inflammaging in mice is promoted by a cell-intrinsic mechanism, stably propagated by ISCs, and associated with a disbalance in immune homeostasis.

Maternal antibiotic exposure enhances ILC2 activation in neonates via downregulation of IFN1 signaling

Microbiota have an important function in shaping and priming neonatal immunity, although the cellular and molecular mechanisms underlying these effects remain obscure. Here we report that prenatal antibiotic exposure causes significant elevation of group 2 innate lymphoid cells (ILC2s) in neonatal lungs, in both cell numbers and functionality. Downregulation of type 1 interferon signaling in ILC2s due to diminished production of microbiota-derived butyrate represents the underlying mechanism. Mice lacking butyrate receptor GPR41 (Gpr41-/-) or type 1 interferon receptor IFNAR1 (Ifnar1-/-) recapitulate the phenotype of neonatal ILC2s upon maternal antibiotic exposure. Furthermore, prenatal antibiotic exposure induces epigenetic changes in ILC2s and has a long-lasting deteriorative effect on allergic airway inflammation in adult offspring. Prenatal supplementation of butyrate ameliorates airway inflammation in adult mice born to antibiotic-exposed dams. These observations demonstrate an essential role for the microbiota in the control of type 2 innate immunity at the neonatal stage, which suggests a therapeutic window for treating asthma in early life.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

Differential gene regulation by a synthetic vitamin D receptor ligand and active vitamin D in human cells

Vitamin D (VD) exerts a wide variety of biological functions including calcemic activity. VD nutritional status is closely associated with the onset and development of chronic diseases. To develop a VD analog with the desired VD activity but without calcemic activity, we screened synthetic VDR antagonists. We identified 1α,25-dihydroxyvitamin D3-26-23-lactams (DLAM)-2a-d (DLAM-2s) as nuclear vitamin D receptor (VDR) ligands in a competitive VDR binding assay for 1α,25(OH)2 vitamin D3 (1α,25(OH)2D3), and DLAM-2s showed an antagonistic effect on 1α,25(OH)2 D3-induced cell differentiation in HL60 cells. In a luciferase reporter assay in which human VDR was exogenously expressed in cultured COS-1 cells, DLAM-2s acted as transcriptional antagonists. Consistently, DLAM-2s had an antagonistic effect on the 1α,25(OH)2D3-induced expression of a known VD target gene [Cytochrome P450 24A1 (CYP24A1)], and VDR bound DLAM-2s was recruited to an endogenous VD response element in chromatin in human keratinocytes (HaCaT cells) endogenously expressing VDR. In an ATAC-seq assay, the effects of 1α,25(OH)2 D3 and DLAM-2b on chromatin reorganization were undetectable in HaCaT cells, while the effect of an androgen receptor (AR) antagonist (bicalutamide) was confirmed in prostate cancer cells (LNCaP) expressing endogenous AR. However, whole genome analysis using RNA-seq and ATAC (Assay for Transposase Accessible Chromatin)-seq revealed differential gene expression profiles regulated by DLAM-2b versus 1α,25(OH)2D3. The upregulated and downregulated genes only partially overlapped between cells treated with 1α,25(OH)2D3 and those treated with DLAM-2b. Thus, the present findings illustrate a novel VDR ligand with gene regulatory activity differing from that of 1α,25(OH)2D3.

Marine Invertebrates One Cell at A Time: Insights from Single-Cell Analysis

Over the past decade, single-cell RNA-sequencing (scRNA-seq) has made it possible to study the cellular diversity of a broad range of organisms. Technological advances in single-cell isolation and sequencing have expanded rapidly, allowing the transcriptomic profile of individual cells to be captured. As a result, there has been an explosion of cell type atlases created for many different marine invertebrate species from across the tree of life. Our focus in this review is to synthesize current literature on marine invertebrate scRNA-seq. Specifically, we provide perspectives on key insights from scRNA-seq studies, including descriptive studies of cell type composition, how cells respond in dynamic processes such as development and regeneration, and the evolution of new cell types. Despite these tremendous advances, there also lie several challenges ahead. We discuss the important considerations that are essential when making comparisons between experiments, or between datasets from different species. Finally, we address the future of single-cell analyses in marine invertebrates, including combining scRNA-seq data with other 'omics methods to get a fuller understanding of cellular complexities. The full diversity of cell types across marine invertebrates remains unknown and understanding this diversity and evolution will provide rich areas for future study.

Loss of cohesin regulator PDS5A reveals repressive role of Polycomb loops

Polycomb Repressive Complexes 1 and 2 (PRC1, PRC2) are conserved epigenetic regulators that promote transcriptional gene silencing. PRC1 and PRC2 converge on shared targets, catalyzing repressive histone modifications. Additionally, a subset of PRC1/PRC2 targets engage in long-range interactions whose functions in gene silencing are poorly understood. Using a CRISPR screen in mouse embryonic stem cells, we found that the cohesin regulator PDS5A links transcriptional silencing by Polycomb and 3D genome organization. PDS5A deletion impairs cohesin unloading and results in derepression of a subset of endogenous PRC1/PRC2 target genes. Importantly, derepression is not linked to loss of Polycomb chromatin domains. Instead, PDS5A removal causes aberrant cohesin activity leading to ectopic insulation sites, which disrupt the formation of ultra-long Polycomb loops. We show that these loops are important for robust silencing at a subset of PRC1/PRC2 target genes and that maintenance of cohesin-dependent genome architecture is critical for Polycomb regulation.

Regulation of transcription patterns, poly-ADP-ribose, and RNA-DNA hybrids by the ATM protein kinase

The ATM protein kinase is a master regulator of the DNA damage response and also an important sensor of oxidative stress. Analysis of gene expression in Ataxia-telangiectasia patient brain tissue shows that large-scale transcriptional changes occur in patient cerebellum that correlate with expression level and GC content of transcribed genes. In human neuron-like cells in culture we map locations of poly-ADP-ribose and RNA-DNA hybrid accumulation genome-wide with ATM inhibition and find that these marks also coincide with high transcription levels, active transcription histone marks, and high GC content. Antioxidant treatment reverses the accumulation of R-loops in transcribed regions, consistent with the central role of ROS in promoting these lesions. Based on these results we postulate that transcription-associated lesions accumulate in ATM-deficient cells and that the single-strand breaks and PARylation at these sites ultimately generate changes in transcription that compromise cerebellum function and lead to neurodegeneration over time in A-T patients.

Interactions between BRD4S, LOXL2, and MED1 drive cell cycle transcription in triple-negative breast cancer

Triple-negative breast cancer (TNBC) often develops resistance to single-agent treatment, which can be circumvented using targeted combinatorial approaches. Here, we demonstrate that the simultaneous inhibition of LOXL2 and BRD4 synergistically limits TNBC proliferation in vitro and in vivo. Mechanistically, LOXL2 interacts in the nucleus with the short isoform of BRD4 (BRD4S), MED1, and the cell cycle transcriptional regulator B-MyB. These interactions sustain the formation of BRD4 and MED1 nuclear transcriptional foci and control cell cycle progression at the gene expression level. The pharmacological co-inhibition of LOXL2 and BRD4 reduces BRD4 nuclear foci, BRD4-MED1 colocalization, and the transcription of cell cycle genes, thus suppressing TNBC cell proliferation. Targeting the interaction between BRD4S and LOXL2 could be a starting point for the development of new anticancer strategies for the treatment of TNBC.

Inhibiting stromal Class I HDACs curbs pancreatic cancer progression

Oncogenic lesions in pancreatic ductal adenocarcinoma (PDAC) hijack the epigenetic machinery in stromal components to establish a desmoplastic and therapeutic resistant tumor microenvironment (TME). Here we identify Class I histone deacetylases (HDACs) as key epigenetic factors facilitating the induction of pro-desmoplastic and pro-tumorigenic transcriptional programs in pancreatic stromal fibroblasts. Mechanistically, HDAC-mediated changes in chromatin architecture enable the activation of pro-desmoplastic programs directed by serum response factor (SRF) and forkhead box M1 (FOXM1). HDACs also coordinate fibroblast pro-inflammatory programs inducing leukemia inhibitory factor (LIF) expression, supporting paracrine pro-tumorigenic crosstalk. HDAC depletion in cancer-associated fibroblasts (CAFs) and treatment with the HDAC inhibitor entinostat (Ent) in PDAC mouse models reduce stromal activation and curb tumor progression. Notably, HDAC inhibition (HDACi) enriches a lipogenic fibroblast subpopulation, a potential precursor for myofibroblasts in the PDAC stroma. Overall, our study reveals the stromal targeting potential of HDACi, highlighting the utility of this epigenetic modulating approach in PDAC therapeutics.

BATF represses BIM to sustain tolerant T cells in the periphery

T cells that encounter self-antigens after exiting the thymus avert autoimmunity through peripheral tolerance. Pathways for this include an unresponsive state known as anergy, clonal deletion, and T regulatory (Treg) cell induction. The transcription factor cues and kinetics that guide distinct peripheral tolerance outcomes remain unclear. Here, we found that anergic T cells are epigenetically primed for regulation by the non-classical AP-1 family member BATF. Tolerized BATF-deficient CD4+ T cells were resistant to anergy induction and instead underwent clonal deletion due to proapoptotic BIM (Bcl2l11) upregulation. During prolonged antigen exposure, BIM derepression resulted in fewer PD-1+ conventional T cells as well as loss of peripherally induced FOXP3+ Treg cells. Simultaneous Batf and Bcl2l11 knockdown meanwhile restored anergic T cell survival and Treg cell maintenance. The data identify the AP-1 nuclear factor BATF as a dominant driver of sustained T cell anergy and illustrate a mechanism for divergent peripheral tolerance fates.

Quantification of absolute transcription factor binding affinities in the native chromatin context using BANC-seq

Transcription factor binding across the genome is regulated by DNA sequence and chromatin features. However, it is not yet possible to quantify the impact of chromatin context on transcription factor binding affinities. Here, we report a method called binding affinities to native chromatin by sequencing (BANC-seq) to determine absolute apparent binding affinities of transcription factors to native DNA across the genome. In BANC-seq, a concentration range of a tagged transcription factor is added to isolated nuclei. Concentration-dependent binding is then measured per sample to quantify apparent binding affinities across the genome. BANC-seq adds a quantitative dimension to transcription factor biology, which enables stratification of genomic targets based on transcription factor concentration and prediction of transcription factor binding sites under non-physiological conditions, such as disease-associated overexpression of (onco)genes. Notably, whereas consensus DNA binding motifs for transcription factors are important to establish high-affinity binding sites, these motifs are not always strictly required to generate nanomolar-affinity interactions in the genome.

OCT4 regulates WNT/β-catenin signaling and prevents mesoendoderm differentiation by repressing EOMES in porcine pluripotent stem cells

The regulatory network between signaling pathways and transcription factors (TFs) is crucial for the maintenance of pluripotent stem cells. However, little is known about how the key TF OCT4 coordinates signaling pathways to regulate self-renewal and lineage differentiation of porcine pluripotent stem cells (pPSCs). Here, we explored the function of OCT4 in pPSCs by transcriptome and chromatin accessibility analysis. The TFs motif enrichment analysis revealed that, following OCT4 knockdown, the regions of increased chromatin accessibility were enriched with EOMES, GATA6, and FOXA1, indicating that pPSCs differentiated toward the mesoendoderm (ME) lineage. Besides, pPSCs rapidly differentiated into ME when the WNT/β-catenin inhibitor XAV939 was removed. However, the ME differentiation of pPSCs caused by OCT4 knockdown did not rely on the activation of WNT/β-catenin signaling because the target gene of WNT/β-catenin signaling, AXIN2 was not upregulated after OCT4 knockdown, despite significant upregulation of WLS and some WNT ligands. Importantly, OCT4 is directly bound to the promoter and enhancers of EOMES and repressed its transcription. Overexpression of EOMES was sufficient to induce ME differentiation in the presence of XAV939. These results demonstrate that OCT4 can regulate WNT/β-catenin signaling and prevent ME differentiation of pPSCs by repressing EOMES transcription.

ROCCO: a robust method for detection of open chromatin via convex optimization

MOTIVATION

Analysis of open chromatin regions across multiple samples from two or more distinct conditions can determine altered gene regulatory patterns associated with biological phenotypes and complex traits. The ATAC-seq assay allows for tractable genome-wide open chromatin profiling of large numbers of samples. Stable, broadly applicable genomic annotations of open chromatin regions are not available. Thus, most studies first identify open regions using peak calling methods for each sample independently. These are then heuristically combined to obtain a consensus peak set. Reconciling sample-specific peak results post hoc from larger cohorts is particularly challenging, and informative spatial features specific to open chromatin signals are not leveraged effectively.

RESULTS

We propose a novel method, ROCCO, that determines consensus open chromatin regions across multiple samples simultaneously. ROCCO employs robust summary statistics and solves a constrained optimization problem formulated to account for both enrichment and spatial dependence of open chromatin signal data. We show this formulation admits attractive theoretical and conceptual properties as well as superior empirical performance compared to current methodology.

AVAILABILITY AND IMPLEMENTATION

Source code, documentation, and usage demos for ROCCO are available on GitHub at: https://github.com/nolan-h-hamilton/ROCCO. ROCCO can also be installed as a stand-alone binary utility using pip/PyPI.

An injury-responsive mmp14b enhancer is required for heart regeneration

Mammals have limited capacity for heart regeneration, whereas zebrafish have extraordinary regeneration abilities. During zebrafish heart regeneration, endothelial cells promote cardiomyocyte cell cycle reentry and myocardial repair, but the mechanisms responsible for promoting an injury microenvironment conducive to regeneration remain incompletely defined. Here, we identify the matrix metalloproteinase Mmp14b as an essential regulator of heart regeneration. We identify a TEAD-dependent mmp14b endothelial enhancer induced by heart injury in zebrafish and mice, and we show that the enhancer is required for regeneration, supporting a role for Hippo signaling upstream of mmp14b. Last, we show that MMP-14 function in mice is important for the accumulation of Agrin, an essential regulator of neonatal mouse heart regeneration. These findings reveal mechanisms for extracellular matrix remodeling that promote heart regeneration.

Dose rate dependent reduction in chromatin accessibility at transcriptional start sites long time after exposure to gamma radiation

Ionizing radiation (IR) impact cellular and molecular processes that require chromatin remodelling relevant for cellular integrity. However, the cellular implications of ionizing radiation (IR) delivered per time unit (dose rate) are still debated. This study investigates whether the dose rate is relevant for inflicting changes to the epigenome, represented by chromatin accessibility, or whether it is the total dose that is decisive. CBA/CaOlaHsd mice were whole-body exposed to either chronic low dose rate (2.5 mGy/h for 54 d) or the higher dose rates (10 mGy/h for 14 d and 100 mGy/h for 30 h) of gamma radiation (60Co, total dose: 3 Gy). Chromatin accessibility was analysed in liver tissue samples using Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-Seq), both one day after and over three months post-radiation (>100 d). The results show that the dose rate contributes to radiation-induced epigenomic changes in the liver at both sampling timepoints. Interestingly, chronic low dose rate exposure to a high total dose (3 Gy) did not inflict long-term changes to the epigenome. In contrast to the acute high dose rate given to the same total dose, reduced accessibility at transcriptional start sites (TSS) was identified in genes relevant for the DNA damage response and transcriptional activity. Our findings link dose rate to essential biological mechanisms that could be relevant for understanding long-term changes after ionizing radiation exposure. However, future studies are needed to comprehend the biological consequence of these findings.

Cisplatin exposure alters tRNA-derived small RNAs but does not affect epimutations in C. elegans

BACKGROUND

The individual lifestyle and environment of an organism can influence its phenotype and potentially the phenotype of its offspring. The different genetic and non-genetic components of the inheritance system and their mutual interactions are key mechanisms to generate inherited phenotypic changes. Epigenetic changes can be transmitted between generations independently from changes in DNA sequence. In Caenorhabditis elegans, epigenetic differences, i.e. epimutations, mediated by small non-coding RNAs, particularly 22G-RNAs, as well as chromatin have been identified, and their average persistence is three to five generations. In addition, previous research showed that some epimutations had a longer duration and concerned genes that were enriched for multiple components of xenobiotic response pathways. These results raise the possibility that environmental stresses might change the rate at which epimutations occur, with potential significance for adaptation.

RESULTS

In this work, we explore this question by propagating C. elegans lines either in control conditions or in moderate or high doses of cisplatin, which introduces genotoxic stress by damaging DNA. Our results show that cisplatin has a limited effect on global small non-coding RNA epimutations and epimutations in gene expression levels. However, cisplatin exposure leads to increased fluctuations in the levels of small non-coding RNAs derived from tRNA cleavage. We show that changes in tRNA-derived small RNAs may be associated with gene expression changes.

CONCLUSIONS

Our work shows that epimutations are not substantially altered by cisplatin exposure but identifies transient changes in tRNA-derived small RNAs as a potential source of variation induced by genotoxic stress.

TGF-β broadly modifies rather than specifically suppresses reactivated memory CD8 T cells in a dose-dependent manner

Transforming growth factor β (TGF-β) directly acts on naive, effector, and memory T cells to control cell fate decisions, which was shown using genetic abrogation of TGF-β signaling. TGF-β availability is altered by infections and cancer; however, the dose-dependent effects of TGF-β on memory CD8 T cell (Tmem) reactivation are still poorly defined. We examined how activation and TGF-β signals interact to shape the functional outcome of Tmem reactivation. We found that TGF-β could suppress cytotoxicity in a manner that was inversely proportional to the strength of the activating TCR or proinflammatory signals. In contrast, even high doses of TGF-β had a comparatively modest effect on IFN-γ expression in the context of weak and strong reactivation signals. Since CD8 Tmem may not always receive TGF-β signals concurrently with reactivation, we also explored whether the temporal order of reactivation versus TGF-β signals is of importance. We found that exposure to TGF-β before or after an activation event were both sufficient to reduce cytotoxic effector function. Concurrent ATAC-seq and RNA-seq analysis revealed that TGF-β altered ~10% of the regulatory elements induced by reactivation and also elicited transcriptional changes indicative of broadly modulated functional properties. We confirmed some changes on the protein level and found that TGF-β-induced expression of CCR8 was inversely proportional to the strength of the reactivating TCR signal. Together, our data suggest that TGF-β is not simply suppressing CD8 Tmem but modifies functional and chemotactic properties in context of their reactivation signals and in a dose-dependent manner.

Integrative approaches based on genomic techniques in the functional studies on enhancers

With the development of sequencing technology and the dramatic drop in sequencing cost, the functions of noncoding genes are being characterized in a wide variety of fields (e.g. biomedicine). Enhancers are noncoding DNA elements with vital transcription regulation functions. Tens of thousands of enhancers have been identified in the human genome; however, the location, function, target genes and regulatory mechanisms of most enhancers have not been elucidated thus far. As high-throughput sequencing techniques have leapt forwards, omics approaches have been extensively employed in enhancer research. Multidimensional genomic data integration enables the full exploration of the data and provides novel perspectives for screening, identification and characterization of the function and regulatory mechanisms of unknown enhancers. However, multidimensional genomic data are still difficult to integrate genome wide due to complex varieties, massive amounts, high rarity, etc. To facilitate the appropriate methods for studying enhancers with high efficacy, we delineate the principles, data processing modes and progress of various omics approaches to study enhancers and summarize the applications of traditional machine learning and deep learning in multi-omics integration in the enhancer field. In addition, the challenges encountered during the integration of multiple omics data are addressed. Overall, this review provides a comprehensive foundation for enhancer analysis.

CEMIG: prediction of the cis-regulatory motif using the de Bruijn graph from ATAC-seq

Sequence motif discovery algorithms enhance the identification of novel deoxyribonucleic acid sequences with pivotal biological significance, especially transcription factor (TF)-binding motifs. The advent of assay for transposase-accessible chromatin using sequencing (ATAC-seq) has broadened the toolkit for motif characterization. Nonetheless, prevailing computational approaches have focused on delineating TF-binding footprints, with motif discovery receiving less attention. Herein, we present Cis rEgulatory Motif Influence using de Bruijn Graph (CEMIG), an algorithm leveraging de Bruijn and Hamming distance graph paradigms to predict and map motif sites. Assessment on 129 ATAC-seq datasets from the Cistrome Data Browser demonstrates CEMIG's exceptional performance, surpassing three established methodologies on four evaluative metrics. CEMIG accurately identifies both cell-type-specific and common TF motifs within GM12878 and K562 cell lines, demonstrating its comparative genomic capabilities in the identification of evolutionary conservation and cell-type specificity. In-depth transcriptional and functional genomic studies have validated the functional relevance of CEMIG-identified motifs across various cell types. CEMIG is available at https://github.com/OSU-BMBL/CEMIG, developed in C++ to ensure cross-platform compatibility with Linux, macOS and Windows operating systems.