ATACseqQC: a Bioconductor package for post-alignment quality assessment of ATAC-seq data

A comparative analysis of planarian genomes reveals regulatory conservation in the face of rapid structural divergence

The planarian Schmidtea mediterranea is being studied as a model species for regeneration, but the assembly of planarian genomes remains challenging. Here, we report a high-quality haplotype-phased, chromosome-scale genome assembly of the sexual S2 strain of S. mediterranea and high-quality chromosome-scale assemblies of its three close relatives, S. polychroa, S. nova, and S. lugubris. Using hybrid gene annotations and optimized ATAC-seq and ChIP-seq protocols for regulatory element annotation, we provide valuable genome resources for the planarian research community and a first comparative perspective on planarian genome evolution. Our analyses reveal substantial divergence in protein-coding sequences and regulatory regions but considerable conservation within promoter and enhancer annotations. We also find frequent retrotransposon-associated chromosomal inversions and interchromosomal translocations within the genus Schmidtea and, remarkably, independent and nearly complete losses of ancestral metazoan synteny in Schmidtea and two other flatworm groups. Overall, our results suggest that platyhelminth genomes can evolve without syntenic constraints.

Arid1a-dependent canonical BAF complex suppresses inflammatory programs to drive efficient germinal center B cell responses

The mammalian Brg1/Brm-associated factor (BAF) complexes are major regulators of nucleosomal remodeling that are commonly mutated in several cancers, including germinal center (GC)-derived B cell lymphomas. However, the specific roles of different BAF complexes in GC B cell biology are not well understood. Here we show that the AT-rich interaction domain 1a (Arid1a) containing canonical BAF (cBAF) complex is required for maintenance of GCs and high-affinity antibody responses. While Arid1a-deficient B cells undergo initial activation, they fail to sustain the GC program. Arid1a establishes permissive chromatin landscapes for B cell activation and is concomitantly required to suppress inflammatory gene programs. The inflammatory signatures instigated by Arid1a deficiency promoted the recruitment of neutrophils and inflammatory monocytes. Dampening of inflammatory cues through interleukin-1β blockade or glucocorticoid receptor agonist partially rescued Arid1a-deficient GCs, highlighting a critical role for inflammation in impeding GCs. Our work reveals essential functions of Arid1a-dependent cBAF in promoting efficient GC responses.

Mechanism for controlled assembly of transcriptional condensates by Aire

Transcriptional condensates play a crucial role in gene expression and regulation, yet their assembly mechanisms remain poorly understood. Here, we report a multi-layered mechanism for condensate assembly by autoimmune regulator (Aire), an essential transcriptional regulator that orchestrates gene expression reprogramming for central T cell tolerance. Aire condensates assemble on enhancers, stimulating local transcriptional activities and connecting disparate inter-chromosomal loci. This functional condensate formation hinges upon the coordination between three Aire domains: polymerization domain caspase activation recruitment domain (CARD), histone-binding domain (first plant homeodomain (PHD1)), and C-terminal tail (CTT). Specifically, CTT binds coactivators CBP/p300, recruiting Aire to CBP/p300-rich enhancers and promoting CARD-mediated condensate assembly. Conversely, PHD1 binds to the ubiquitous histone mark H3K4me0, keeping Aire dispersed throughout the genome until Aire nucleates on enhancers. Our findings showed that the balance between PHD1-mediated suppression and CTT-mediated stimulation of Aire polymerization is crucial to form transcriptionally active condensates at target sites, providing new insights into controlled polymerization of transcriptional regulators.

Muscle growth differences in Lijiang pigs revealed by ATAC-seq multi-omics

As one of the largest tissues in the animal body, skeletal muscle plays a pivotal role in the production and quality of pork. Consequently, it is of paramount importance to investigate the growth and developmental processes of skeletal muscle. Lijiang pigs, which naturally have two subtypes, fast-growing and slow-growing, provide an ideal model for such studies by eliminating breed-related influences. In this study, we selected three fast-growing and three slow-growing 6-month-old Lijiang pigs as subjects. We utilized assay for transposase-accessible chromatin with sequencing (ATAC-seq) combined with genomics, RNA sequencing, and proteomics to screen for differentially expressed genes and transcription factors linked to increased longissimus dorsi muscle volume in Lijiang pigs. We identified 126 genes through ATAC-seq, including PPARA, TNRC6B, NEDD1, and FKBP5, that exhibited differential expression patterns during muscle growth. Additionally, we identified 59 transcription factors, including Foxh1, JunB, Mef2 family members (Mef2a/b/c/d), NeuroD1, and TEAD4. By examining open chromatin regions (OCRs) with significant genetic differentiation, genes such as SAV1, CACNA1H, PRKCG, and FGFR4 were found. Integrating ATAC-seq with transcriptomics and transcriptomics with proteomics, we identified differences in open chromatin regions, transcription, and protein levels of FKBP5 and SCARB2 genes in fast-growing and slow-growing Lijiang pigs. Utilizing multi-omics analysis with R packages, we jointed ATAC-seq, transcriptome, and proteome datasets, identifying enriched pathways related to glycogen metabolism and skeletal muscle cell differentiation. We pinpointed genes such as MYF6 and HABP2 that exhibit strong correlations across these diverse data types. This study provides a multi-faceted understanding of the molecular mechanisms that lead to differences in pig muscle fiber growth.

PHF6 cooperates with SWI/SNF complexes to facilitate transcriptional progression

Genes encoding subunits of SWI/SNF (BAF) chromatin remodeling complexes are mutated in nearly 25% of cancers. To gain insight into the mechanisms by which SWI/SNF mutations drive cancer, we contributed ten rhabdoid tumor (RT) cell lines mutant for SWI/SNF subunit SMARCB1 to a genome-scale CRISPR-Cas9 depletion screen performed across 896 cell lines. We identify PHF6 as specifically essential for RT cell survival and demonstrate that dependency on Phf6 extends to Smarcb1-deficient cancers in vivo. As mutations in either SWI/SNF or PHF6 can cause the neurodevelopmental disorder Coffin-Siris syndrome, our findings of a dependency suggest a previously unrecognized functional link. We demonstrate that PHF6 co-localizes with SWI/SNF complexes at promoters, where it is essential for maintenance of an active chromatin state. We show that in the absence of SMARCB1, PHF6 loss disrupts the recruitment and stability of residual SWI/SNF complex members, collectively resulting in the loss of active chromatin at promoters and stalling of RNA Polymerase II progression. Our work establishes a mechanistic basis for the shared syndromic features of SWI/SNF and PHF6 mutations in CSS and the basis for selective dependency on PHF6 in SMARCB1-mutant cancers.

Role of the circadian nuclear receptor REV-ERBα in dorsal raphe serotonin synthesis in mood regulation

Affective disorders are frequently associated with disrupted circadian rhythms. The existence of rhythmic secretion of central serotonin (5-hydroxytryptamine, 5-HT) pattern has been reported; however, the functional mechanism underlying the circadian control of 5-HTergic mood regulation remains largely unknown. Here, we investigate the role of the circadian nuclear receptor REV-ERBα in regulating tryptophan hydroxylase 2 (Tph2), the rate-limiting enzyme of 5-HT synthesis. We demonstrate that the REV-ERBα expressed in dorsal raphe (DR) 5-HTergic neurons functionally competes with PET-1-a nuclear activator crucial for 5-HTergic neuron development. In mice, genetic ablation of DR 5-HTergic REV-ERBα increases Tph2 expression, leading to elevated DR 5-HT levels and reduced depression-like behaviors at dusk. Further, pharmacological manipulation of the mice DR REV-ERBα activity increases DR 5-HT levels and affects despair-related behaviors. Our findings provide valuable insights into the molecular and cellular link between the circadian rhythm and the mood-controlling DR 5-HTergic systems.

Quantitative analysis of cis-regulatory elements in transcription with KAS-ATAC-seq

Cis-regulatory elements (CREs) are pivotal in orchestrating gene expression throughout diverse biological systems. Accurate identification and in-depth characterization of functional CREs are crucial for decoding gene regulation networks during cellular processes. In this study, we develop Kethoxal-Assisted Single-stranded DNA Assay for Transposase-Accessible Chromatin with Sequencing (KAS-ATAC-seq) to quantitatively analyze the transcriptional activity of CREs. A main advantage of KAS-ATAC-seq lies in its precise measurement of ssDNA levels within both proximal and distal ATAC-seq peaks, enabling the identification of transcriptional regulatory sequences. This feature is particularly adept at defining Single-Stranded Transcribing Enhancers (SSTEs). SSTEs are highly enriched with nascent RNAs and specific transcription factors (TFs) binding sites that define cellular identity. Moreover, KAS-ATAC-seq provides a detailed characterization and functional implications of various SSTE subtypes. Our analysis of CREs during mouse neural differentiation demonstrates that KAS-ATAC-seq can effectively identify immediate-early activated CREs in response to retinoic acid (RA) treatment. Our findings indicate that KAS-ATAC-seq provides more precise annotation of functional CREs in transcription. Future applications of KAS-ATAC-seq would help elucidate the intricate dynamics of gene regulation in diverse biological processes.

Epigenomics coverage data extraction and aggregation in R with tidyCoverage

SUMMARY

The tidyCoverage R package provides a framework for intuitive investigation of collections of genomic tracks over genomic features, relying on the principle of tidy data manipulation. It defines two data structures, CoverageExperiment and AggregatedCoverage classes, directly extending the SummarizedExperiment fundamental class, and introduces a principled approach to exploring genome-wide data. This infrastructure facilitates the extraction and manipulation of genomic coverage track data across individual or multiple sets of thousands of genomic loci. This allows the end user to rapidly visualize track coverage at individual genomic loci or aggregated coverage profiles over sets of genomic loci. tidyCoverage seamlessly combines with the existing Bioconductor ecosystem to accelerate the integration of genome-wide track data in epigenomic analysis workflows. tidyCoverage emerges as a valuable tool, contributing to the advancement of epigenomics research by promoting consistency, reproducibility, and accessibility in data analysis.

AVAILABILITY AND IMPLEMENTATION

tidyCoverage is an R package freely available from Bioconductor ≥ 3.19 (https://www.bioconductor.org/packages/tidyCoverage) for R ≥ 4.4. The software is distributed under the MIT License and is accompanied by example files and data.

The CD8+ T cell tolerance checkpoint triggers a distinct differentiation state defined by protein translation defects

Peripheral CD8+ T cell tolerance is a checkpoint in both autoimmune disease and anti-cancer immunity. Despite its importance, the relationship between tolerance-induced states and other CD8+ T cell differentiation states remains unclear. Using flow cytometric phenotyping, single-cell RNA sequencing (scRNA-seq), and chromatin accessibility profiling, we demonstrated that in vivo peripheral tolerance to a self-antigen triggered a fundamentally distinct differentiation state separate from exhaustion, memory, and functional effector cells but analogous to cells defectively primed against tumors. Tolerant cells diverged early and progressively from effector cells, adopting a transcriptionally and epigenetically distinct state within 60 h of antigen encounter. Breaching tolerance required the synergistic actions of strong T cell receptor (TCR) signaling and inflammation, which cooperatively induced gene modules that enhanced protein translation. Weak TCR signaling during bystander infection failed to breach tolerance due to the uncoupling of effector gene expression from protein translation. Thus, tolerance engages a distinct differentiation trajectory enforced by protein translation defects.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Single-Cell Sequencing Technology in Ruminant Livestock: Challenges and Opportunities

Advancements in single-cell sequencing have transformed the genomics field by allowing researchers to delve into the intricate cellular heterogeneity within tissues at greater resolution. While single-cell omics are more widely applied in model organisms and humans, their use in livestock species is just beginning. Studies in cattle, sheep, and goats have already leveraged single-cell and single-nuclei RNA-seq as well as single-cell and single-nuclei ATAC-seq to delineate cellular diversity in tissues, track changes in cell populations and gene expression over developmental stages, and characterize immune cell populations important for disease resistance and resilience. Although challenges exist for the use of this technology in ruminant livestock, such as the precise annotation of unique cell populations and spatial resolution of cells within a tissue, there is vast potential to enhance our understanding of the cellular and molecular mechanisms underpinning traits essential for healthy and productive livestock. This review intends to highlight the insights gained from published single-cell omics studies in cattle, sheep, and goats, particularly those with publicly accessible data. Further, this manuscript will discuss the challenges and opportunities of this technology in ruminant livestock and how it may contribute to enhanced profitability and sustainability of animal agriculture in the future.

Chromatin remodeling protein BPTF regulates transcriptional stability in planarian stem cells

Trimethylation of histone H3 lysine 4 (H3K4me3) correlates strongly with gene expression in many different organisms, yet the question of whether it plays a causal role in transcriptional activity remains unresolved. Although H3K4me3 does not directly affect chromatin accessibility, it can indirectly affect genome accessibility by recruiting the ATP-dependent chromatin remodeling complex NuRF (Nucleosome Remodeling Factor). The largest subunit of NuRF, BPTF/NURF301, binds H3K4me3 specifically and recruits the NuRF complex to loci marked by this modification. Studies have shown that the strength and duration of BPTF binding likely also depends on additional chromatin features at these loci, such as lysine acetylation and variant histone proteins. However, the exact details of this recruitment mechanism vary between studies and have largely been tested in vitro. Here, we use stem cells isolated directly from live planarian animals to investigate the role of BPTF in regulating chromatin accessibility in vivo. We find that BPTF operates at gene promoters and is most effective at facilitating transcription at genes marked by Set1-dependent H3K4me3 peaks, which are significantly broader than those added by the lysine methyltransferase MLL1/2. Moreover, BPTF is essential for planarian stem cell biology and its loss of function phenotype mimics that of Set1 knockdown. Together, these data suggest that BPTF and H3K4me3 are important mediators of both transcription and in vivo stem cell function.

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk

Osteoarthritis (OA) poses a significant healthcare burden with limited treatment options. While genome-wide association studies (GWAS) have identified over 100 OA-associated loci, translating these findings into therapeutic targets remains challenging. Integrating expression quantitative trait loci (eQTL), 3D chromatin structure, and other genomic approaches with OA GWAS data offers a promising approach to elucidate disease mechanisms; however, comprehensive eQTL maps in OA-relevant tissues and conditions remain scarce. We mapped gene expression, chromatin accessibility, and 3D chromatin structure in primary human articular chondrocytes in both resting and OA-mimicking conditions. We identified thousands of differentially expressed genes, including those associated with differences in sex and age. RNA-seq in chondrocytes from 101 donors across two conditions uncovered 3782 unique eGenes, including 420 that exhibited strong and significant condition-specific effects. Colocalization with OA GWAS signals revealed 13 putative OA risk genes, 10 of which have not been previously identified. Chromatin accessibility and 3D chromatin structure provided insights into the mechanisms and conditional specificity of these variants. Our findings shed light on OA pathogenesis and highlight potential targets for therapeutic development.

Highlights

∘ Comprehensive analysis of sex- and age-related global gene expression in human chondrocytes revealed differences that correlate with osteoarthritis ∘ First response eQTLs in chondrocytes treated with an OA-related stimulus ∘ Deeply sequenced Hi-C in resting and activated chondrocytes helps connect OA risk variants to their putative causal genes ∘ Colocalization analysis reveals 13 (including 10 novel) putative OA risk genes.

YY1 binding is a gene-intrinsic barrier to Xist-mediated gene silencing

X chromosome inactivation (XCI) in mammals is mediated by Xist RNA which functions in cis to silence genes on a single X chromosome in XX female cells, thereby equalising levels of X-linked gene expression relative to XY males. XCI progresses over a period of several days, with some X-linked genes silencing faster than others. The chromosomal location of a gene is an important determinant of silencing rate, but uncharacterised gene-intrinsic features also mediate resistance or susceptibility to silencing. In this study, we examine mouse embryonic stem cell lines with an inducible Xist allele (iXist-ChrX mESCs) and integrate allele-specific data of gene silencing and decreasing inactive X (Xi) chromatin accessibility over time courses of Xist induction with cellular differentiation. Our analysis reveals that motifs bound by the transcription factor YY1 are associated with persistently accessible regulatory elements, including many promoters and enhancers of slow-silencing genes. We further show that YY1 is evicted relatively slowly from target sites on Xi, and that silencing of X-linked genes is increased upon YY1 degradation. Together our results suggest that YY1 acts as a barrier to Xist-mediated silencing until the late stages of the XCI process.

Clock-dependent chromatin accessibility rhythms regulate circadian transcription

Chromatin organization plays a crucial role in gene regulation by controlling the accessibility of DNA to transcription machinery. While significant progress has been made in understanding the regulatory role of clock proteins in circadian rhythms, how chromatin organization affects circadian rhythms remains poorly understood. Here, we employed ATAC-seq (Assay for Transposase-Accessible Chromatin with Sequencing) on FAC-sorted Drosophila clock neurons to assess genome-wide chromatin accessibility at dawn and dusk over the circadian cycle. We observed significant oscillations in chromatin accessibility at promoter and enhancer regions of hundreds of genes, with enhanced accessibility either at dusk or dawn, which correlated with their peak transcriptional activity. Notably, genes with enhanced accessibility at dusk were enriched with E-box motifs, while those more accessible at dawn were enriched with VRI/PDP1-box motifs, indicating that they are regulated by the core circadian feedback loops, PER/CLK and VRI/PDP1, respectively. Further, we observed a complete loss of chromatin accessibility rhythms in per01 null mutants, with chromatin consistently accessible at both dawn and dusk, underscoring the critical role of Period protein in driving chromatin compaction during the repression phase at dawn. Together, this study demonstrates the significant role of chromatin organization in circadian regulation, revealing how the interplay between clock proteins and chromatin structure orchestrates the precise timing of biological processes throughout the day. This work further implies that variations in chromatin accessibility might play a central role in the generation of diverse circadian gene expression patterns in clock neurons.

Multi-omics panoramic analysis of HBV integration, transcriptional regulation, and epigenetic modifications in PLC/PRF/5 cell line

The clearance or transcriptional silencing of integrated HBV DNA is crucial for achieving a functional cure in patients with chronic hepatitis B and reducing the risk of hepatocellular carcinoma development. The PLC/PRF/5 cell line is commonly used as an in vitro model for studying HBV integration. In this study, we employed a range of multi-omics techniques to gain a panoramic understanding of the characteristics of HBV integration in PLC/PRF/5 cells and to reveal the transcriptional regulatory mechanisms of integrated HBV DNA. Transcriptome long-read sequencing (ONT) was conducted to analyze and characterize the transcriptional activity of different HBV DNA integration sites in PLC/PRF/5 cells. Additionally, we collected data related to epigenetic regulation, including whole-genome bisulfite sequencing (WGBS), histone chromatin immunoprecipitation sequencing (ChIP-seq), and assays for transposase-accessible chromatin using sequencing (ATAC-seq), to explore the potential mechanisms involved in the transcriptional regulation of integrated HBV DNA. Long-read RNA sequencing analysis revealed significant transcriptional differences at various integration sites in the PLC/PRF/5 cell line, with higher HBV DNA transcription levels at integration sites on chr11, chr13, and the chr13/chr5 fusion chromosome t (13:5). Combining long-read DNA and RNA sequencing results, we found that transcription of integrated HBV DNA generally starts downstream of the SP1, SP2, or XP promoters. ATAC-seq data confirmed that chromatin accessibility has limited influence on the transcription of integrated HBV DNA in the PLC/PRF/5 cell line. Analysis of WGBS data showed that the methylation intensity of integrated HBV DNA was highly negatively correlated with its transcription level (r = -0.8929, p = 0.0123). After AzaD treatment, the transcription level of integrated HBV DNA significantly increased, especially for the integration chr17, which had the highest level of methylation. Through ChIP-seq data, we observed the association between histone modification of H3K4me3 and H3K9me3 with the transcription of integrated HBV DNA. Our findings suggest that the SP1, SP2 and XP in integrated HBV DNA, methylation level of surrounding host chromosome, and histone modifications affect the transcription of integrated HBV DNA in PLC/PRF/5 cells. This provides important clues for future studies on the expression and regulatory mechanisms of integrated HBV.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.

TRPS1 modulates chromatin accessibility to regulate estrogen receptor alpha (ER) binding and ER target gene expression in luminal breast cancer cells

Common genetic variants in the repressive GATA-family transcription factor (TF) TRPS1 locus are associated with breast cancer risk, and luminal breast cancer cell lines are particularly sensitive to TRPS1 knockout. We introduced an inducible degron tag into the native TRPS1 locus within a luminal breast cancer cell line to identify the direct targets of TRPS1 and determine how TRPS1 mechanistically regulates gene expression. We acutely deplete over 80 percent of TRPS1 from chromatin within 30 minutes of inducing degradation. We find that TRPS1 regulates transcription of hundreds of genes, including those related to estrogen signaling. TRPS1 directly regulates chromatin structure, which causes estrogen receptor alpha (ER) to redistribute in the genome. ER redistribution leads to both repression and activation of dozens of ER target genes. Downstream from these primary effects, TRPS1 depletion represses cell cycle-related gene sets and reduces cell doubling rate. Finally, we show that high TRPS1 activity, calculated using a gene expression signature defined by primary TRPS1-regulated genes, is associated with worse breast cancer patient prognosis. Taken together, these data suggest a model in which TRPS1 modulates the genomic distribution of ER, both activating and repressing transcription of genes related to cancer cell fitness.

Aging differentially alters the transcriptome and landscape of chromatin accessibility in the male and female mouse hippocampus

Aging-related memory impairment and pathological memory disorders such as Alzheimer's disease differ between males and females, and yet little is known about how aging-related changes in the transcriptome and chromatin environment differ between sexes in the hippocampus. To investigate this question, we compared the chromatin accessibility landscape and gene expression/alternative splicing pattern of young adult and aged mouse hippocampus in both males and females using ATAC-seq and RNA-seq. We detected significant aging-dependent changes in the expression of genes involved in immune response and synaptic function and aging-dependent changes in the alternative splicing of myelin sheath genes. We found significant sex-bias in the expression and alternative splicing of hundreds of genes, including aging-dependent female-biased expression of myelin sheath genes and aging-dependent male-biased expression of genes involved in synaptic function. Aging was associated with increased chromatin accessibility in both male and female hippocampus, especially in repetitive elements, and with an increase in LINE-1 transcription. We detected significant sex-bias in chromatin accessibility in both autosomes and the X chromosome, with male-biased accessibility enriched at promoters and CpG-rich regions. Sex differences in gene expression and chromatin accessibility were amplified with aging, findings that may shed light on sex differences in aging-related and pathological memory loss.

Deterministic reprogramming of neutrophils within tumors

Neutrophils are increasingly recognized as key players in the tumor immune response and are associated with poor clinical outcomes. Despite recent advances characterizing the diversity of neutrophil states in cancer, common trajectories and mechanisms governing the ontogeny and relationship between these neutrophil states remain undefined. Here, we demonstrate that immature and mature neutrophils that enter tumors undergo irreversible epigenetic, transcriptional, and proteomic modifications to converge into a distinct, terminally differentiated dcTRAIL-R1+ state. Reprogrammed dcTRAIL-R1+ neutrophils predominantly localize to a glycolytic and hypoxic niche at the tumor core and exert pro-angiogenic function that favors tumor growth. We found similar trajectories in neutrophils across multiple tumor types and in humans, suggesting that targeting this program may provide a means of enhancing certain cancer immunotherapies.

Supergene evolution via gain of auto-regulation

The development of complex phenotypes requires the coordinated action of many genes across space and time, yet many species have evolved the ability to develop multiple discrete, alternate phenotypes1-3. Such polymorphisms are often controlled by supergenes, sets of tightly-linked mutations in one or more loci that function together to produce a complex phenotype4. Although theories of supergene evolution are well-established, the mutations that cause functional differences between supergene alleles remain essentially unknown. doublesex is the master regulator of insect sexual differentiation but functions as a supergene in multiple Papilio swallowtail butterflies, where divergent dsx alleles control development of discrete non-mimetic or mimetic female wing color patterns5-7. Here we demonstrate that the functional elements of the mimetic allele in Papilio alphenor are six new cis-regulatory elements (CREs) spread across 150 kb that are bound by DSX itself. Our findings provide experimental support to classic supergene theory and suggest that the evolution of auto-regulation may provide a simple route to supergene origination and to the co-option of pleiotropic genes into new developmental roles.

Repression of latent NF-κB enhancers by PDX1 regulates β cell functional heterogeneity

Interactions between lineage-determining and activity-dependent transcription factors determine single-cell identity and function within multicellular tissues through incompletely known mechanisms. By assembling a single-cell atlas of chromatin state within human islets, we identified β cell subtypes governed by either high or low activity of the lineage-determining factor pancreatic duodenal homeobox-1 (PDX1). β cells with reduced PDX1 activity displayed increased chromatin accessibility at latent nuclear factor κB (NF-κB) enhancers. Pdx1 hypomorphic mice exhibited de-repression of NF-κB and impaired glucose tolerance at night. Three-dimensional analyses in tandem with chromatin immunoprecipitation (ChIP) sequencing revealed that PDX1 silences NF-κB at circadian and inflammatory enhancers through long-range chromatin contacts involving SIN3A. Conversely, Bmal1 ablation in β cells disrupted genome-wide PDX1 and NF-κB DNA binding. Finally, antagonizing the interleukin (IL)-1β receptor, an NF-κB target, improved insulin secretion in Pdx1 hypomorphic islets. Our studies reveal functional subtypes of single β cells defined by a gradient in PDX1 activity and identify NF-κB as a target for insulinotropic therapy.

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.

Epigenomic mapping reveals distinct B cell acute lymphoblastic leukemia chromatin architectures and regulators

B cell lineage acute lymphoblastic leukemia (B-ALL) is composed of diverse molecular subtypes, and while transcriptional and DNA methylation profiling has been extensively examined, the chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq in primary B-ALL cells from 156 patients spanning ten molecular subtypes and present this dataset as a resource. Differential chromatin accessibility and transcription factor (TF) footprint profiling were employed and identified B-ALL cell of origin, TF-target gene interactions enriched in B-ALL, and key TFs associated with accessible chromatin sites preferentially active in B-ALL. We further identified over 20% of accessible chromatin sites exhibiting strong subtype enrichment and candidate TFs that maintain subtype-specific chromatin architectures. Over 9,000 genetic variants were uncovered, contributing to variability in chromatin accessibility among patient samples. Our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants that promote unique gene-regulatory networks.

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.

PAK1-dependent mechanotransduction enables myofibroblast nuclear adaptation and chromatin organization during fibrosis

Myofibroblasts are responsible for scarring during fibrosis. The scar propagates mechanical signals inducing a radical transformation in myofibroblast cell state and increasing profibrotic phenotype. Here, we show mechanical stress from progressive scarring induces nuclear softening and de-repression of heterochromatin. The parallel loss of H3K9Me3 enables a permissive state for distinct chromatin accessibility and profibrotic gene regulation. Integrating chromatin accessibility profiles with RNA expression provides insight into the transcription network underlying the switch in profibrotic myofibroblast states, emphasizing mechanoadaptive regulation of PAK1 as key drivers. Through genetic manipulation in liver and lung fibrosis, loss of PAK1-dependent signaling impairs the mechanoadaptive response in vitro and dramatically improves fibrosis in vivo. Moreover, we provide human validation for mechanisms underpinning PAK1-mediated mechanotransduction in liver and lung fibrosis. Collectively, these observations provide insight into the nuclear mechanics driving the profibrotic chromatin landscape in fibrosis, highlighting actomyosin-dependent mechanisms as potential therapeutic targets in fibrosis.

The MYC axis in advanced prostate cancer is impacted through concurrent targeting of ERβ and AR using a novel ERβ-selective ligand alongside Enzalutamide

We have dissected the role of Estrogen receptor beta (ERβ) in prostate cancer (PCa) with a novel ERβ ligand, OSU-ERb-12. Drug screens revealed additive interactions between OSU-ERB-12 and either epigenetic inhibitors or the androgen receptor antagonist, Enzalutamide (Enza). Clonogenic and cell biolody studies supported the potent additive effects of OSU-ERB-12 (100nM) and Enza (1μM). The cooperative behavior was in PCa cell lines treated with either OSU-ERB-12 plus Enza or combinations involving 17β-estradiol (E2). OSU-ERb-12 plus Enza uniquely impacted the transcriptiome, accessible chromatin, and the AR, MYC and H3K27ac cistromes. This included skewed transcriptional responses including suppression of the androgen and MYC transcriptomes, and repressed MYC protein. OSU-ERb-12 plus Enza uniquely impacted chromatin accessibility at approximately 3000 nucleosome-free sites, enriched at enhancers, enriched for basic Helix-Loop-Helix motifs. CUT&RUN experiments revealed combination treatment targeting of MYC, AR, and H3K27ac again shaping enhancer accessibility. Specifically, it repressed MYC interactions at enhancer regions enriched for bHLH motifs, and overlapped with publicly-available bHLH cistromes. Finally, cistrome-transcriptome analyses identified ~200 genes that distinguished advanced PCa tumors in the SU2C cohort with high androgen and low neuroendocrine scores.

cisDynet: An integrated platform for modeling gene-regulatory dynamics and networks

Chromatin accessibility sequencing has been widely used for uncovering genetic regulatory mechanisms and inferring gene regulatory networks. However, effectively integrating large-scale chromatin accessibility datasets has posed a significant challenge. This is due to the lack of a comprehensive end-to-end solution, as many existing tools primarily emphasize data preprocessing and overlook downstream analyses. To bridge this gap, we have introduced cisDynet, a holistic solution that combines streamlined data preprocessing using Snakemake and R functions with advanced downstream analysis capabilities. cisDynet excels in conventional data analyses, encompassing peak statistics, peak annotation, differential analysis, motif enrichment analysis, and more. Additionally, it allows to perform sophisticated data exploration, such as tissue-specific peak identification, time course data modeling, integration of RNA-seq data to establish peak-to-gene associations, constructing regulatory networks, and conducting enrichment analysis of genome-wide association study (GWAS) variants. As a proof of concept, we applied cisDynet to reanalyze comprehensive ATAC-seq datasets across various tissues from the Encyclopedia of DNA Elements (ENCODE) project. The analysis successfully delineated tissue-specific open chromatin regions (OCRs), established connections between OCRs and target genes, and effectively linked these discoveries with 1861 GWAS variants. Furthermore, cisDynet was instrumental in dissecting the time course open chromatin data of mouse embryonic development, revealing the dynamic behavior of OCRs over developmental stages and identifying key transcription factors governing differentiation trajectories. In summary, cisDynet offers researchers a user-friendly solution that minimizes the need for extensive coding, ensures the reproducibility of results, and greatly simplifies the exploration of epigenomic data.

Chromatin accessibility profiling of targeted cell populations with laser capture microdissection coupled to ATAC-seq

Spatially resolved omics technologies reveal context-dependent cellular regulatory networks in tissues of interest. Beyond transcriptome analysis, information on epigenetic traits and chromatin accessibility can provide further insights on gene regulation in health and disease. Nevertheless, compared to the enormous advancements in spatial transcriptomics technologies, the field of spatial epigenomics is much younger and still underexplored. In this study, we report laser capture microdissection coupled to ATAC-seq (LCM-ATAC-seq) applied to fresh frozen samples for the spatial characterization of chromatin accessibility. We first demonstrate the efficient use of LCM coupled to in situ tagmentation and evaluate its performance as a function of cell number, microdissected areas, and tissue type. Further, we demonstrate its use for the targeted chromatin accessibility analysis of discrete contiguous or scattered cell populations in tissues via single-nuclei capture based on immunostaining for specific cellular markers.

Inflammatory stress-mediated chromatin changes underlie dysfunction in endothelial cells

Inflammatory stresses underlie endothelial dysfunction and contribute to the development of chronic cardiovascular disorders such as atherosclerosis and vascular fibrosis. The initial transcriptional response of endothelial cells to pro-inflammatory cytokines such as TNF-alpha is well established. However, very few studies uncover the effects of inflammatory stresses on chromatin architecture. We used integrative analysis of ATAC-seq and RNA-seq data to investigate chromatin alterations in human endothelial cells in response to TNF-alpha and febrile-range heat stress exposure. Multi-omics data analysis suggests a correlation between the transcription of stress-related genes and endothelial dysfunction drivers with chromatin regions exhibiting differential accessibility. Moreover, microscopy identified the dynamics in the nuclear organization, specifically, the changes in a subset of heterochromatic nucleoli-associated chromatin domains, the centromeres. Upon inflammatory stress exposure, the centromeres decreased association with nucleoli in a p38-dependent manner and increased the number of transcripts from pericentromeric regions. Overall, we provide two lines of evidence that suggest chromatin alterations in vascular endothelial cells during inflammatory stresses.

A comparative atlas of single-cell chromatin accessibility in the human brain

Recent advances in single-cell transcriptomics have illuminated the diverse neuronal and glial cell types within the human brain. However, the regulatory programs governing cell identity and function remain unclear. Using a single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), we explored open chromatin landscapes across 1.1 million cells in 42 brain regions from three adults. Integrating this data unveiled 107 distinct cell types and their specific utilization of 544,735 candidate cis-regulatory DNA elements (cCREs) in the human genome. Nearly a third of the cCREs demonstrated conservation and chromatin accessibility in the mouse brain cells. We reveal strong links between specific brain cell types and neuropsychiatric disorders including schizophrenia, bipolar disorder, Alzheimer's disease (AD), and major depression, and have developed deep learning models to predict the regulatory roles of noncoding risk variants in these disorders.

Chromatin Accessibility Landscape of Human Triple-negative Breast Cancer Cell Lines Reveals Variation by Patient Donor Ancestry

African American (AA) women have an excessive risk of developing triple-negative breast cancer (TNBC). We employed Assay for Transposase-Accessible Chromatin using sequencing to characterize differences in chromatin accessibility between nine commonly used TNBC cell lines derived from patients of European and African ancestry. Principal component and chromosome mapping analyses of accessibility peaks with the most variance revealed separation of chromatin profiles by patient group. Motif enrichment and footprinting analyses of disparate open chromatin regions revealed differences in transcription factor activity, identifying 79 with ancestry-associated binding patterns (FDR < 0.01). AA TNBC cell lines exhibited increased accessibility for 62 transcription factors associated with epithelial-to-mesenchymal transition, cancer stemness/chemotherapeutic resistance, proliferation, and aberrant p53 regulation, as well as KAISO, which has been previously linked to aggressive tumor characteristics in AA patients with cancer. Differential Assay for Transposase-Accessible Chromatin signal analysis identified 1,596 genes located within promoters of differentially open chromatin regions in AA-derived TNBC, identifying DNA methyltransferase 1 as the top upregulated gene associated with African ancestry. Pathway analyses with these genes revealed enrichment in several pathways, including hypoxia. Culturing cells under hypoxia showed ancestry-specific stress responses that led to the identification of a core set of AA-associated transcription factors, which included members of the Kruppel-like factor and Sp subfamilies, as well as KAISO, and identified ZDHHC1, a gene previously implicated in immunity and STING activation, as the top upregulated AA-specific gene under hypoxia. Together, these data reveal a differential chromatin landscape in TNBC associated with donor ancestry. The open chromatin structure of AA TNBC may contribute to a more lethal disease.

SIGNIFICANCE

We identify an ancestry-associated open chromatin landscape and related transcription factors that may contribute to aggressive TNBC in AA women. Furthermore, this study advocates for the inclusion of diversely sourced cell lines in experimental in vitro studies to advance health equity at all levels of scientific research.

Ribosome biogenesis disruption mediated chromatin structure changes revealed by SRAtac, a customizable end to end analysis pipeline for ATAC-seq

The nucleolus is a large nuclear body that serves as the primary site for ribosome biogenesis. Recent studies have suggested that it also plays an important role in organizing chromatin architecture. However, to establish a causal relationship between nucleolar ribosome assembly and chromatin architecture, genetic tools are required to disrupt nucleolar ribosome biogenesis. In this study, we used ATAC-seq to investigate changes in chromatin accessibility upon specific depletion of two ribosome biogenesis components, RPOA-2 and GRWD-1, in the model organism Caenorhabditis elegans. To facilitate the analysis of ATAC-seq data, we introduced two tools: SRAlign, an extensible NGS data processing workflow, and SRAtac, a customizable end-to-end ATAC-seq analysis pipeline. Our results revealed highly comparable changes in chromatin accessibility following both RPOA-2 and GRWD-1 perturbations. However, we observed a weak correlation between changes in chromatin accessibility and gene expression. While our findings corroborate the idea of a feedback mechanism between ribosomal RNA synthesis, nucleolar ribosome large subunit biogenesis, and chromatin structure during the L1 stage of C. elegans development, they also prompt questions regarding the functional impact of these alterations on gene expression.

Clock-dependent chromatin accessibility rhythms regulate circadian transcription

Chromatin organization plays a crucial role in gene regulation by controlling the accessibility of DNA to transcription machinery. While significant progress has been made in understanding the regulatory role of clock proteins in circadian rhythms, how chromatin organization affects circadian rhythms remains poorly understood. Here, we employed ATAC-seq (Assay for Transposase-Accessible Chromatin with Sequencing) on FAC-sorted Drosophila clock neurons to assess genome-wide chromatin accessibility over the circadian cycle. We observed significant circadian oscillations in chromatin accessibility at promoter and enhancer regions of hundreds of genes, with enhanced accessibility either at dusk or dawn, which correlated with their peak transcriptional activity. Notably, genes with enhanced accessibility at dusk were enriched with E-box motifs, while those more accessible at dawn were enriched with VRI/PDP1-box motifs, indicating that they are regulated by the core circadian feedback loops, PER/CLK and VRI/PDP1, respectively. Further, we observed a complete loss of chromatin accessibility rhythms in per01 null mutants, with chromatin consistently accessible throughout the circadian cycle, underscoring the critical role of Period protein in driving chromatin compaction during the repression phase. Together, this study demonstrates the significant role of chromatin organization in circadian regulation, revealing how the interplay between clock proteins and chromatin structure orchestrates the precise timing of biological processes throughout the day. This work further implies that variations in chromatin accessibility might play a central role in the generation of diverse circadian gene expression patterns in clock neurons.

Best practices for single-cell analysis across modalities

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

The HSV-1 ICP22 protein selectively impairs histone repositioning upon Pol II transcription downstream of genes

Herpes simplex virus 1 (HSV-1) infection and stress responses disrupt transcription termination by RNA Polymerase II (Pol II). In HSV-1 infection, but not upon salt or heat stress, this is accompanied by a dramatic increase in chromatin accessibility downstream of genes. Here, we show that the HSV-1 immediate-early protein ICP22 is both necessary and sufficient to induce downstream open chromatin regions (dOCRs) when transcription termination is disrupted by the viral ICP27 protein. This is accompanied by a marked ICP22-dependent loss of histones downstream of affected genes consistent with impaired histone repositioning in the wake of Pol II. Efficient knock-down of the ICP22-interacting histone chaperone FACT is not sufficient to induce dOCRs in ΔICP22 infection but increases dOCR induction in wild-type HSV-1 infection. Interestingly, this is accompanied by a marked increase in chromatin accessibility within gene bodies. We propose a model in which allosteric changes in Pol II composition downstream of genes and ICP22-mediated interference with FACT activity explain the differential impairment of histone repositioning downstream of genes in the wake of Pol II in HSV-1 infection.

The time-resolved genomic impact of Wnt/β-catenin signaling

Wnt signaling orchestrates gene expression via its effector, β-catenin. However, it is unknown whether β-catenin binds its target genomic regions simultaneously and how this impacts chromatin dynamics to modulate cell behavior. Using a combination of time-resolved CUT&RUN against β-catenin, ATAC-seq, and perturbation assays in different cell types, we show that Wnt/β-catenin physical targets are tissue-specific, β-catenin "moves" on different loci over time, and its association to DNA accompanies changing chromatin accessibility landscapes that determine cell behavior. In particular, Wnt/β-catenin progressively shapes the chromatin of human embryonic stem cells (hESCs) as they undergo mesodermal differentiation, a behavior that we define as "plastic." In HEK293T cells, on the other hand, Wnt/β-catenin drives a transient chromatin opening, followed by re-establishment of the pre-stimulation state, a response that we define as "elastic." Future experiments shall assess whether other cell communication mechanisms, in addition to Wnt signaling, are ruled by time, cellular idiosyncrasies, and chromatin constraints. A record of this paper's transparent peer review process is included in the supplemental information.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Signaling mechanisms in renal compensatory hypertrophy revealed by multi-omics

Loss of a kidney results in compensatory growth of the remaining kidney, a phenomenon of considerable clinical importance. However, the mechanisms involved are largely unknown. Here, we use a multi-omic approach in a unilateral nephrectomy model in male mice to identify signaling processes associated with renal compensatory hypertrophy, demonstrating that the lipid-activated transcription factor peroxisome proliferator-activated receptor alpha (PPARα) is an important determinant of proximal tubule cell size and is a likely mediator of compensatory proximal tubule hypertrophy.

Epigenetic programming defines haematopoietic stem cell fate restriction

Haematopoietic stem cells (HSCs) are multipotent, but individual HSCs can show restricted lineage output in vivo. Currently, the molecular mechanisms and physiological role of HSC fate restriction remain unknown. Here we show that lymphoid fate is epigenetically but not transcriptionally primed in HSCs. In multi-lineage HSCs that produce lymphocytes, lymphoid-specific upstream regulatory elements (LymUREs) but not promoters are preferentially accessible compared with platelet-biased HSCs that do not produce lymphoid cell types, providing transcriptionally silent lymphoid lineage priming. Runx3 is preferentially expressed in multi-lineage HSCs, and reinstating Runx3 expression increases LymURE accessibility and lymphoid-primed multipotent progenitor 4 (MPP4) output in old, platelet-biased HSCs. In contrast, platelet-biased HSCs show elevated levels of epigenetic platelet-lineage priming and give rise to MPP2 progenitors with molecular platelet bias. These MPP2 progenitors generate platelets with faster kinetics and through a more direct cellular pathway compared with MPP2s derived from multi-lineage HSCs. Epigenetic programming therefore predicts both fate restriction and differentiation kinetics in HSCs.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Optimized protocol for assay for transposase-accessible chromatin by sequencing (ATAC-seq) from Drosophila melanogaster brain tissue

Transposase-accessible chromatin by sequencing (ATAC-seq) has emerged as an advantageous technique to assess chromatin accessibility owing to the robustness of "tagmentation" process and a relatively faster library preparation. A comprehensive ATAC-seq protocol from Drosophila brain tissue is currently unavailable. Here, we have provided a detailed protocol of ATAC-seq assay from Drosophila brain tissue. Starting from dissection and transposition to amplification of libraries has been elaborated. Furthermore, a robust ATAC-seq analysis pipeline has been presented. The protocol can be easily adapted for other soft tissues.

Protocol for Bulk-ATAC sequencing in head and neck squamous cell carcinoma

The transposase-accessible chromatin using sequencing (ATAC-seq) offers a simplified approach to detect chromatin changes in cancer cells after genetic intervention and drug treatment. Here, we present an optimized ATAC-seq protocol to elucidate chromatin accessibility changes at the epigenetic level in head and neck squamous cell carcinoma cells. We describe steps for cell lysate preparation, transposition, and tagmentation, followed by library amplification and purification. We then detail next-generation sequencing and data analysis. For complete details on the use and execution of this protocol, please refer to Buenrostro et al.,¹ Chen et al.,².

Loss of H3K9 trimethylation alters chromosome compaction and transcription factor retention during mitosis

Recent studies have shown that repressive chromatin machinery, including DNA methyltransferases and polycomb repressor complexes, binds to chromosomes throughout mitosis and their depletion results in increased chromosome size. In the present study, we show that enzymes that catalyze H3K9 methylation, such as Suv39h1, Suv39h2, G9a and Glp, are also retained on mitotic chromosomes. Surprisingly, however, mutants lacking histone 3 lysine 9 trimethylation (H3K9me3) have unusually small and compact mitotic chromosomes associated with increased histone H3 phospho Ser10 (H3S10ph) and H3K27me3 levels. Chromosome size and centromere compaction in these mutants were rescued by providing exogenous first protein lysine methyltransferase Suv39h1 or inhibiting Ezh2 activity. Quantitative proteomic comparisons of native mitotic chromosomes isolated from wild-type versus Suv39h1/Suv39h2 double-null mouse embryonic stem cells revealed that H3K9me3 was essential for the efficient retention of bookmarking factors such as Esrrb. These results highlight an unexpected role for repressive heterochromatin domains in preserving transcription factor binding through mitosis and underscore the importance of H3K9me3 for sustaining chromosome architecture and epigenetic memory during cell division.

African American Prostate Cancer Displays Quantitatively Distinct Vitamin D Receptor Cistrome-transcriptome Relationships Regulated by BAZ1A

African American (AA) prostate cancer associates with vitamin D3 deficiency, but vitamin D receptor (VDR) genomic actions have not been investigated in this context. We undertook VDR proteogenomic analyses in European American (EA) and AA prostate cell lines and four clinical cohorts. Rapid immunoprecipitation mass spectrometry of endogenous protein (RIME) analyses revealed that nonmalignant AA RC43N prostate cells displayed the greatest dynamic protein content in the VDR complex. Likewise, in AA cells, Assay for Transposase-Accessible Chromatin using sequencing established greater 1α,25(OH)2D3-regulated chromatin accessibility, chromatin immunoprecipitation sequencing revealed significant enhancer-enriched VDR cistrome, and RNA sequencing identified the largest 1α,25(OH)2D3-dependent transcriptome. These VDR functions were significantly corrupted in the isogenic AA RC43T prostate cancer cells, and significantly distinct from EA cell models. We identified reduced expression of the chromatin remodeler, BAZ1A, in three AA prostate cancer cohorts as well as RC43T compared with RC43N. Restored BAZ1A expression significantly increased 1α,25(OH)2D3-regulated VDR-dependent gene expression in RC43T, but not HPr1AR or LNCaP cells. The clinical impact of VDR cistrome-transcriptome relationships were tested in three different clinical prostate cancer cohorts. Strikingly, only in AA patients with prostate cancer, the genes bound by VDR and/or associated with 1α,25(OH)2D3-dependent open chromatin (i) predicted progression from high-grade prostatic intraepithelial neoplasia to prostate cancer; (ii) responded to vitamin D3 supplementation in prostate cancer tumors; (iii) differentially responded to 25(OH)D3 serum levels. Finally, partial correlation analyses established that BAZ1A and components of the VDR complex identified by RIME significantly strengthened the correlation between VDR and target genes in AA prostate cancer only. Therefore, VDR transcriptional control is most potent in AA prostate cells and distorted through a BAZ1A-dependent control of VDR function.

Significance

Our study identified that genomic ancestry drives the VDR complex composition, genomic distribution, and transcriptional function, and is disrupted by BAZ1A and illustrates a novel driver for AA prostate cancer.

High level of novelty under the hood of convergent evolution

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

Genomic glucocorticoid action in embryonic mouse neural stem cells

Prenatal exposure to synthetic glucocorticoids (sGCs) reprograms brain development and predisposes the developing fetus towards potential adverse neurodevelopmental outcomes. Using a mouse model of sGC administration, previous studies show that these changes are accompanied by sexually dimorphic alterations in the transcriptome of neural stem and progenitor cells (NSPCs) derived from the embryonic telencephalon. Because cell type-specific gene expression profiles tightly regulate cell fate decisions and are controlled by a flexible landscape of chromatin domains upon which transcription factors and enhancer elements act, we multiplexed data from four genome-wide assays: RNA-seq, ATAC-seq (assay for transposase accessible chromatin followed by genome wide sequencing), dual cross-linking ChIP-seq (chromatin immunoprecipitation followed by genome wide sequencing), and microarray gene expression to identify novel relationships between gene regulation, chromatin structure, and genomic glucocorticoid receptor (GR) action in NSPCs. These data reveal that GR binds preferentially to predetermined regions of accessible chromatin to influence gene programming and cell fate decisions. In addition, we identify SOX2 as a transcription factor that impacts the genomic response of select GR target genes to sGCs (i.e., dexamethasone) in NSPCs.

Chicken chromatin accessibility atlas accelerates epigenetic annotation of birds and gene fine-mapping associated with growth traits

The development of epigenetic maps, such as the ENCODE project in humans, provides resources for gene regulation studies and a reference for research of disease-related regulatory elements. However, epigenetic information, such as a bird-specific chromatin accessibility atlas, is currently lacking for the thousands of bird species currently described. The major genomic difference between birds and mammals is their shorter introns and intergenic distances, which seriously hinders the use of humans and mice as a reference for studying the function of important regulatory regions in birds. In this study, using chicken as a model bird species, we systematically compiled a chicken chromatin accessibility atlas using 53 Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) samples across 11 tissues. An average of 50 796 open chromatin regions were identified per sample, cumulatively accounting for 20.36% of the chicken genome. Tissue specificity was largely reflected by differences in intergenic and intronic peaks, with specific functional regulation achieved by two mechanisms: recruitment of several sequence-specific transcription factors and direct regulation of adjacent functional genes. By integrating data from genome-wide association studies, our results suggest that chicken body weight is driven by different regulatory variants active in growth-relevant tissues. We propose CAB39L (active in the duodenum), RCBTB1 (muscle and liver), and novel long non-coding RNA ENSGALG00000053256 (bone) as candidate genes regulating chicken body weight. Overall, this study demonstrates the value of epigenetic data in fine-mapping functional variants and provides a compendium of resources for further research on the epigenetics and evolution of birds and mammals.

maxATAC: Genome-scale transcription-factor binding prediction from ATAC-seq with deep neural networks

Transcription factors read the genome, fundamentally connecting DNA sequence to gene expression across diverse cell types. Determining how, where, and when TFs bind chromatin will advance our understanding of gene regulatory networks and cellular behavior. The 2017 ENCODE-DREAM in vivo Transcription-Factor Binding Site (TFBS) Prediction Challenge highlighted the value of chromatin accessibility data to TFBS prediction, establishing state-of-the-art methods for TFBS prediction from DNase-seq. However, the more recent Assay-for-Transposase-Accessible-Chromatin (ATAC)-seq has surpassed DNase-seq as the most widely-used chromatin accessibility profiling method. Furthermore, ATAC-seq is the only such technique available at single-cell resolution from standard commercial platforms. While ATAC-seq datasets grow exponentially, suboptimal motif scanning is unfortunately the most common method for TFBS prediction from ATAC-seq. To enable community access to state-of-the-art TFBS prediction from ATAC-seq, we (1) curated an extensive benchmark dataset (127 TFs) for ATAC-seq model training and (2) built "maxATAC", a suite of user-friendly, deep neural network models for genome-wide TFBS prediction from ATAC-seq in any cell type. With models available for 127 human TFs, maxATAC is the largest collection of high-performance TFBS prediction models for ATAC-seq. maxATAC performance extends to primary cells and single-cell ATAC-seq, enabling improved TFBS prediction in vivo. We demonstrate maxATAC's capabilities by identifying TFBS associated with allele-dependent chromatin accessibility at atopic dermatitis genetic risk loci.

Nuclear Biophysical Changes during Human Melanoma Plasticity

Tumor plasticity is an emerging property of tumor cells which allows them to change their phenotype in dependence on the environment. The epithelial-mesenchymal transition plays a crucial role in helping cells acquire a more aggressive phenotype when they are in the mesenchymal state. Herein, we investigated the biophysical changes occurring during phenotypic switching in human melanoma cells, considering the blebbiness of the nuclei, their stiffness, and the involvement of polycombs with lamins. We show that the formation of cellular heterogeneity involves many crucial nuclear changes including the interaction between different types of polycombs with lamins and chromosome accessibility. Altogether, our results shed new light on the molecular mechanisms involved in the formation of a heterogeneous cell population during phenotypic switching. In particular, our results show that phenotypic switching in melanoma involves chromatin remodeling changing the transcriptional activity of cells and consequently their phenotype.