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

Unraveling the role of sperm retained histones in bull fertility and daughter fertility

During spermatogenesis, a substantial proportion of histones are substituted by protamine to condense the genome within the sperm head. Studies indicate that a minority of histones, typically ranging from 1 to 15 %, persist in mammalian sperm post-substitution. The persistence of histones in the zygote facilitates chromatin accessibility to transcription factors in regions crucial for early embryonic development. Nevertheless, the potential causal relationship between retained histones and fertility phenotypes remains uncertain. This study seeks to investigate this relationship. The results indicate that in mature bovine sperm, regions of DNA associated with fertility that bind to histones are primarily concentrated in promoters and transcription start sites, potentially impacting bull fertility and offspring fertility through the regulation of relevant genes. Furthermore, microRNAs and estradiol/ESR are suggested to be the main regulators of the canonical pathways identified, highlighting the need for additional research to investigate their potential utility as biomarkers.

Gene regulatory dynamics during the development of a paleopteran insect, the mayfly Cloeon dipterum

The evolution of insects has been marked by the appearance of key body plan innovations that promoted the outstanding ability of this lineage to adapt to new habitats, boosting the most successful radiation in animals. To understand the evolution of these new structures, it is essential to investigate which genes and gene regulatory networks participate during the embryonic development of insects. Great efforts have been made to fully understand gene expression and gene regulation during the development of holometabolous insects, in particular Drosophila melanogaster. Conversely, functional genomics resources and databases in other insect lineages are scarce. To provide a new platform to study gene regulation in insects, we generated ATAC-seq for the first time during the development of the mayfly Cloeon dipterum, which belongs to Paleoptera, the sister group to all other winged insects. With these comprehensive datasets along six developmental stages, we characterized pronounced changes in accessible chromatin between early and late embryogenesis. The application of ATAC-seq in mayflies provides a fundamental resource to understand the evolution of gene regulation in insects.

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

The tea cultivar 'Chungui' with jasmine-like aroma: From genome and epigenome to quality

'Chungui' is a newly promoted tea cultivar in China, renowned for producing oolong tea with a distinctive jasmine-like aroma. However, the genetic basis of this unique aroma remains unclear. In this study, the 'Chungui' genome, one of the most complete and well-annotated tea genomes, was assembled using PacBio HiFi reads and Hi-C sequencing. Through comparative analysis with typical jasmine flower volatiles, eight core compounds responsible for this aroma were identified. Further research revealed that the jasmine-like aroma in 'Chungui' is regulated by a coordinated mechanism involving a significant increase in chromatin accessibility and the demethylation of CHH and CHG in the promoter regions of key aroma-related genes during oolong tea processing. The study proposes that the formation of this unique aroma is driven by the synergistic effect of enhanced chromatin accessibility and reduced methylation, which together lead to the robust upregulation of genes involved in the biosynthesis of these core aroma components. These results provide a molecular foundation for understanding the unique jasmine-like aroma of 'Chungui' tea and sets the stage for future studies to explore the roles of these regulatory mechanisms in aroma formation.

Translation of genome-wide association study: from genomic signals to biological insights

Since the turn of the 21st century, genome-wide association study (GWAS) have successfully identified genetic signals associated with a myriad of common complex traits and diseases. As we transition from establishing robust genetic associations with diverse phenotypes, the central challenge is now focused on characterizing the underlying functional mechanisms driving these signals. Previous GWAS efforts have revealed multiple variants, each conferring relatively subtle susceptibility, collectively contributing to the pathogenesis of various common diseases. Such variants can further exhibit associations with multiple other traits and differ across ancestries, plus disentangling causal variants from non-causal due to linkage disequilibrium complexities can lead to challenges in drawing direct biological conclusions. Combined with cellular context considerations, such challenges can reduce the capacity to definitively elucidate the biological significance of GWAS signals, limiting the potential to define mechanistic insights. This review will detail current and anticipated approaches for functional interpretation of GWAS signals, both in terms of characterizing the underlying causal variants and the corresponding effector genes.

ARID1A-BAF coordinates ZIC2 genomic occupancy for epithelial-to-mesenchymal transition in cranial neural crest specification

The BAF chromatin remodeler regulates lineage commitment including cranial neural crest cell (CNCC) specification. Variants in BAF subunits cause Coffin-Siris syndrome (CSS), a congenital disorder characterized by coarse craniofacial features and intellectual disability. Approximately 50% of individuals with CSS harbor variants in one of the mutually exclusive BAF subunits, ARID1A/ARID1B. While Arid1a deletion in mouse neural crest causes severe craniofacial phenotypes, little is known about the role of ARID1A in CNCC specification. Using CSS-patient-derived ARID1A+/- induced pluripotent stem cells to model CNCC specification, we discovered that ARID1A-haploinsufficiency impairs epithelial-to-mesenchymal transition (EMT), a process necessary for CNCC delamination and migration from the neural tube. Furthermore, wild-type ARID1A-BAF regulates enhancers associated with EMT genes. ARID1A-BAF binding at these enhancers is impaired in heterozygotes while binding at promoters is unaffected. At the sequence level, these EMT enhancers contain binding motifs for ZIC2, and ZIC2 binding at these sites is ARID1A-dependent. When excluded from EMT enhancers, ZIC2 relocates to neuronal enhancers, triggering aberrant neuronal gene activation. In mice, deletion of Zic2 impairs NCC delamination, while ZIC2 overexpression in chick embryos at post-migratory neural crest stages elicits ectopic delamination from the neural tube. These findings reveal an essential ARID1A-ZIC2 axis essential for EMT and CNCC delamination.

TET activity safeguards pluripotency throughout embryonic dormancy

Dormancy is an essential biological process for the propagation of many life forms through generations and stressful conditions. Early embryos of many mammals are preservable for weeks to months within the uterus in a dormant state called diapause, which can be induced in vitro through mTOR inhibition. Cellular strategies that safeguard original cell identity within the silent genomic landscape of dormancy are not known. Here we show that the protection of cis-regulatory elements from silencing is key to maintaining pluripotency in the dormant state. We reveal a TET-transcription factor axis, in which TET-mediated DNA demethylation and recruitment of methylation-sensitive transcription factor TFE3 drive transcriptionally inert chromatin adaptations during dormancy transition. Perturbation of TET activity compromises pluripotency and survival of mouse embryos under dormancy, whereas its enhancement improves survival rates. Our results reveal an essential mechanism for propagating the cellular identity of dormant cells, with implications for regeneration and disease.

RNF20-mediated transcriptional pausing and VEGFA splicing orchestrate vessel growth

Signal-responsive gene expression is essential for vascular development, yet the mechanisms integrating signaling inputs with transcriptional activities are largely unknown. Here we show that RNF20, the primary E3 ubiquitin ligase for histone H2B, plays a multifaceted role in sprouting angiogenesis. RNF20 mediates RNA polymerase (Pol II) promoter-proximal pausing at genes highly paused in endothelial cells, involved in VEGFA signaling, stress response, cell cycle control and mRNA splicing. It also orchestrates large-scale mRNA processing events that alter the bioavailability and function of critical pro-angiogenic factors, such as VEGFA. Mechanistically, RNF20 restricts ERG-dependent Pol II pause release at highly paused genes while binding to Notch1 to promote H2B monoubiquitination at Notch target genes and Notch-dependent gene expression. This balance is crucial, as loss of Rnf20 leads to uncontrolled tip cell specification. Our findings highlight the pivotal role of RNF20 in regulating VEGF-Notch signaling circuits during vessel growth, underscoring its potential for therapeutic modulation of angiogenesis.

Differential gene expression and chromatin accessibility in the medial prefrontal cortex associated with individual differences in rat behavioral models of opioid use disorder

Opioid use disorder (OUD) is a neuropsychological disease that has a devastating impact on public health. Substantial individual differences in vulnerability exist, the neurobiological substrates of which remain unclear. To address this question, we investigated genome-wide gene transcription (RNA-seq) and chromatin accessibility (ATAC-seq) in the medial prefrontal cortex (mPFC) of male and female rats exhibiting differential vulnerability in behavioral paradigms modeling different phases of OUD: Withdrawal-Induced Anhedonia (WIA), Demand, and Reinstatement. Ingenuity Pathway Analysis (IPA) of RNA-seq revealed greater changes in canonical pathways in Resilient (vs. Saline) rats in comparison to Vulnerable (vs. Saline) rats across 3 paradigms, suggesting brain adaptations that might contribute to resilience to OUD across its trajectory. Analyses of gene networks and upstream regulators implicated processes involved in oligodendrocyte maturation and myelination in WIA, neuroinflammation in Demand, and metabolism in Reinstatement. Motif analysis of ATAC-seq showed changes in chromatin accessibility to a small set of transcription factor (TF) binding sites as a function either of opioid exposure (i.e., morphine versus saline) generally or of individual vulnerability specifically. Some of these were shared across the 3 paradigms and others were unique to each. In conclusion, we have identified changes in biological pathways, TFs, and their binding motifs that vary with paradigm and OUD vulnerability. These findings point to the involvement of distinct transcriptional and epigenetic mechanisms in response to opioid exposure, vulnerability to OUD, and different stages of the disorder.

Agricultural insect pests as models for studying stress-induced evolutionary processes

Agricultural insect pests (AIPs) are widely successful in adapting to natural and anthropogenic stressors, repeatedly overcoming population bottlenecks and acquiring resistance to intensive management practices. Although they have been largely overlooked in evolutionary studies, AIPs are ideal systems for understanding rapid adaptation under novel environmental conditions. Researchers have identified several genomic mechanisms that likely contribute to adaptive stress responses, including positive selection on de novo mutations, polygenic selection on standing allelic variation and phenotypic plasticity (e.g., hormesis). However, new theory suggests that stress itself may induce epigenetic modifications, which may confer heritable physiological changes (i.e., stress-resistant phenotypes). In this perspective, we discuss how environmental stress from agricultural management generates the epigenetic and genetic modifications that are associated with rapid adaptation in AIPs. We summarise existing evidence for stress-induced evolutionary processes in the context of insecticide resistance. Ultimately, we propose that studying AIPs offers new opportunities and resources for advancing our knowledge of stress-induced evolution.

XIST dampens X chromosome activity in a SPEN-dependent manner during early human development

XIST (X-inactive specific transcript) long noncoding RNA (lncRNA) is responsible for X chromosome inactivation (XCI) in placental mammals, yet it accumulates on both X chromosomes in human female preimplantation embryos without triggering X chromosome silencing. The XACT (X-active coating transcript) lncRNA coaccumulates with XIST on active X chromosomes and may antagonize XIST function. Here, we used human embryonic stem cells in a naive state of pluripotency to assess the function of XIST and XACT in shaping the X chromosome chromatin and transcriptional landscapes during preimplantation development. We show that XIST triggers the deposition of polycomb-mediated repressive histone modifications and dampens the transcription of most X-linked genes in a SPEN-dependent manner, while XACT deficiency does not significantly affect XIST activity or X-linked gene expression. Our study demonstrates that XIST is functional before XCI, confirms the existence of a transient process of X chromosome dosage compensation and reveals that XCI and dampening rely on the same set of factors.

The stress response regulator HSF1 modulates natural killer cell anti-tumour immunity

Diverse cellular insults converge on activation of the heat shock factor 1 (HSF1), which regulates the proteotoxic stress response to maintain protein homoeostasis. HSF1 regulates numerous gene programmes beyond the proteotoxic stress response in a cell-type- and context-specific manner to promote malignancy. However, the role(s) of HSF1 in immune populations of the tumour microenvironment remain elusive. Here, we leverage an in vivo model of HSF1 activation and single-cell transcriptomic tumour profiling to show that augmented HSF1 activity in natural killer (NK) cells impairs cytotoxicity, cytokine production and subsequent anti-tumour immunity. Mechanistically, HSF1 directly binds and regulates the expression of key mediators of NK cell effector function. This work demonstrates that HSF1 regulates the immune response under the stress conditions of the tumour microenvironment. These findings have important implications for enhancing the efficacy of adoptive NK cell therapies and for designing combinatorial strategies including modulators of NK cell-mediated tumour killing.

Dissecting the Single-Cell Diversity and Heterogeneity Underlying Cervical Precancerous Lesions and Cancer Tissues

The underlying cellular diversity and heterogeneity from cervix precancerous lesions to cervical squamous cell carcinoma (CSCC) is investigated. Four single-cell datasets including normal tissues, normal adjacent tissues, precancerous lesions, and cervical tumors were integrated to perform disease stage analysis. Single-cell compositional data analysis (scCODA) was utilized to reveal the compositional changes of each cell type. Differentially expressed genes (DEGs) among cell types were annotated using BioCarta. An assay for transposase-accessible chromatin sequencing (ATAC-seq) analysis was performed to correlate epigenetic alterations with gene expression profiles. Lastly, a logistic regression model was used to assess the similarity between the original and new cohort data (HRA001742). After global annotation, seven distinct cell types were categorized. Eight consensus-upregulated DEGs were identified in B cells among different disease statuses, which could be utilized to predict the overall survival of CSCC patients. Inferred copy number variation (CNV) analysis of epithelial cells guided disease progression classification. Trajectory and ATAC-seq integration analysis identified 95 key transcription factors (TF) and one immunohistochemistry (IHC) testified key-node TF (YY1) involved in epithelial cells from CSCC initiation to progression. The consistency of epithelial cell subpopulation markers was revealed with single-cell sequencing, bulk sequencing, and RT-qPCR detection. KRT8 and KRT15, markers of Epi6, showed progressively higher expression with disease progression as revealed by IHC detection. The logistic regression model testified the robustness of the resemblance of clusters among the various datasets utilized in this study. Valuable insights into CSCC cellular diversity and heterogeneity provide a foundation for future targeted therapy.

Neutrophils Recruited by NKX2-1 Suppression via Activation of CXCLs/CXCR2 Axis Promote Lung Adenocarcinoma Progression

NK2 Homeobox 1 (NKX2-1) is a well-characterized pathological marker that delineates lung adenocarcinoma (LUAD) progression. The advancement of LUAD is influenced by the immune tumor microenvironment through paracrine signaling. However, the involvement of NKX2-1 in modeling the tumor immune microenvironment is still unclear. Here, the downregulation of NKX2-1 is observed in high-grade LUAD. Meanwhile, single-cell RNA sequencing and Visium in situ capturing profiling revealed the recruitment and infiltration of neutrophils in orthotopic syngeneic tumors exhibiting strong cell-cell communication through the activation of CXCLs/CXCR2 signaling. The depletion of NKX2-1 triggered the expression and secretion of CXCL1, CXCL2, CXCL3, and CXCL5 in LUAD cells. Chemokine secretion is analyzed by chemokine array and validated by qRT-PCR. ATAC-seq revealed the restrictive regulation of NKX2-1 on the promoters of CXCL1, CXCL2, and CXCL5 genes. This phenomenon led to increased tumor growth, and conversely, tumor growth decreased when inhibited by the CXCR2 antagonist SB225002. This study unveils how NKX2-1 modulates the infiltration of tumor-promoting neutrophils by inhibiting CXCLs/CXCR2-dependent mechanisms. Hence, targeting CXCR2 in NKX2-1-low tumors is a potential antitumor therapy that may improve LUAD patient outcomes.

Differential chromatin accessibility and Gene Expression Associated with Backfat Deposition in pigs

BACKGROUND

Backfat serves as a vital fat reservoir in pigs, and its excessive accumulation will adversely impact pig growth performance, farming efficiency, and pork quality. The aim of this research is to integrate assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing (RNA-seq) to explore the molecular mechanisms underlying porcine backfat deposition.

RESULTS

ATAC-seq analysis identified 568 genes originating from 698 regions exhibiting differential accessibility, which were significantly enriched in pathways pertinent to adipocyte differentiation and lipid metabolism. Besides, a total of 283 transcription factors (TFs) were identified by motif analysis. RNA-seq analysis revealed 978 differentially expressed genes (DEGs), which were enriched in pathways related to energy metabolism, cell cycle and signal transduction. The integration of ATAC-seq and RNA-seq data indicates that DEG expression levels are associated with chromatin accessibility. This comprehensive study highlights the involvement of critical pathways, including the Wnt signaling pathway, Jak-STAT signaling pathway, and fatty acid degradation, in the regulation of backfat deposition. Through rigorous analysis, we identified several candidate genes (LEP, CTBP2, EHHADH, OSMR, TCF7L2, BCL2, FGF1, UCP2, CCND1, TIMP1, and VDR) as potentially significant contributors to backfat deposition. Additionally, we constructed TF-TF and TF-target gene regulatory networks and identified a series of potential TFs related to backfat deposition (FOS, STAT3, SMAD3, and ESR1).

CONCLUSIONS

This study represents the first application of ATAC-seq and RNA-seq, affording a novel perspective into the mechanisms underlying backfat deposition and providing invaluable resources for the enhancement of pig breeding programs.

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 generate 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. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers 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.

CD70 identifies alloreactive T cells and represents a potential target for prevention and treatment of acute GVHD

ABSTRACT

Graft-versus-host disease (GVHD) remains a major challenge after allogeneic hematopoietic stem cell transplantation (allo-HSCT), and further understanding of its immunopathology is crucial for developing new treatments. CD70 interacts with CD27 and is upregulated transiently on T cells after recent T-cell receptor (TCR) engagement. Here, we investigated the functional and clinical significance of CD70 expression on T cells during the early posttransplantation period. CD70 was expressed on a subset of highly activated memory T cells within the first 2 weeks after transplant, which then gradually declined in most patients. CD70+ T cells exhibited an open chromatin landscape and a transcriptional profile indicative of intense Myelocytomatosis oncogene (MYC)-driven glycolysis and proliferation. CD4+ and CD8+CD70+ T-cell numbers increased by ninefold and fourfold, respectively, during acute GVHD (aGVHD) and displayed an oligoclonal TCR repertoire. These cells expressed CCR4 and CCR6 chemokine receptors and were markedly increased in aGVHD tissue samples. Furthermore, CD70+ T cells demonstrated alloreactive specificity in vitro, and proliferative and inflammatory cytokine responses were markedly attenuated by CD70 blockade. These findings identify CD70 as a marker of highly activated alloreactive T cells and reveal the potential therapeutic importance of inhibiting CD27-CD70 costimulation in both the prophylaxis and treatment of aGVHD.

Analysis of chromatin accessibility in peripheral blood mononuclear cells from patients with early-stage breast cancer

Objective: This study was aimed at exploring a specific open region of chromatin in the peripheral blood mononuclear cells (PBMCs) of patients with breast cancer and evaluating its feasibility as a biomarker for diagnosing and predicting breast cancer prognosis. Methods: We obtained PBMCs from breast cancer patients and healthy people for the assay for transposase-accessible chromatin (ATAC) sequencing (n = 3) and obtained the GSE27562 chip sequencing data for secondary analyses. Through bioinformatics analysis, we mined the pattern changes for chromatin accessibility in the PBMCs of breast cancer patients. Results: A total of 1,906 differentially accessible regions (DARs) and 1,632 differentially expressed genes (DEGs) were identified via ATAC sequencing. The upregulated DEGs in the disease group were mainly distributed in the cells, organelles, and cell-intima-related structures and were mainly responsible for biological functions such as cell nitrogen complex metabolism, macromolecular metabolism, and cell communication, in addition to functions such as nucleic acid binding, enzyme binding, hydrolase reaction, and transferase activity. Combined with microarray data analysis, the following set of nine DEGs showed intersection between the ATAC and microarray data: JUN, MSL2, CDC42, TRIB1, SERTAD3, RAB14, RHOB, RAB40B, and PRKDC. HOMER predicted and identified five transcription factors that could potentially bind to these peak sites, namely NFY, Sp 2, GFY, NRF, and ELK 1. Conclusion: Chromatin accessibility analysis of the PBMCs from patients with early-stage breast cancer underscores its potential as a significant avenue for biomarker discovery in breast cancer diagnostics and treatment. By screening the transcription factors and DEGs related to breast cancer, this study provides a comprehensive theoretical foundation that is expected to guide future clinical applications and therapeutic developments.

Cohesin-mediated chromatin remodeling controls the differentiation and function of conventional dendritic cells

The cohesin protein complex extrudes chromatin loops, stopping at CTCF-bound sites, to organize chromosomes into topologically associated domains, yet the biological implications of this process are poorly understood. We show that cohesin is required for the post-mitotic differentiation and function of antigen-presenting dendritic cells (DCs), particularly for antigen cross-presentation and IL-12 secretion by type 1 conventional DCs (cDC1s) in vivo. The chromatin organization of DCs was shaped by cohesin and the DC-specifying transcription factor IRF8, which controlled chromatin looping and chromosome compartmentalization, respectively. Notably, optimal expression of IRF8 itself required CTCF/cohesin-binding sites demarcating the Irf8 gene. During DC activation, cohesin was required for the induction of a subset of genes with distal enhancers. Accordingly, the deletion of CTCF sites flanking the Il12b gene reduced IL-12 production by cDC1s. Our data reveal an essential role of cohesin-mediated chromatin regulation in cell differentiation and function in vivo, and its bi-directional crosstalk with lineage-specifying transcription factors.

Super-silencer perturbation by EZH2 and REST inhibition leads to large loss of chromatin interactions and reduction in cancer growth

Human silencers have been shown to regulate developmental gene expression. However, the functional importance of human silencers needs to be elucidated, such as whether they can form 'super-silencers' and whether they are linked to cancer progression. Here, we show two silencer components of the FGF18 gene can cooperate through compensatory chromatin interactions to form a super-silencer. Double knockout of two silencers exhibited synergistic upregulation of FGF18 expression and changes in cell identity. To perturb the super-silencers, we applied combinational treatment of an enhancer of zeste homolog 2 inhibitor GSK343, and a repressor element 1-silencing transcription factor inhibitor, X5050 ('GR'). Interestingly, GR led to severe loss of topologically associated domains and loops, which were associated with reduced CTCF and TOP2A mRNA levels. Moreover, GR synergistically upregulated super-silencer-controlled genes related to cell cycle, apoptosis and DNA damage, leading to anticancer effects in vivo. Overall, our data demonstrated a super-silencer example and showed that GR can disrupt super-silencers, potentially leading to cancer ablation.

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.

RNA Polymerase II coordinates histone deacetylation at active promoters

Nucleosomes at actively transcribed promoters have specific histone post-transcriptional modifications and histone variants. These features are thought to contribute to the formation and maintenance of a permissive chromatin environment. Recent reports have drawn conflicting conclusions about whether these histone modifications depend on transcription. We used triptolide to inhibit transcription initiation and degrade RNA Polymerase II and interrogated the effect on histone modifications. Transcription initiation was dispensable for de novo and steady-state histone acetylation at transcription start sites (TSSs) and enhancers. However, at steady state, blocking transcription initiation increased the levels of histone acetylation and H2AZ incorporation at active TSSs. These results demonstrate that deposition of specific histone modifications at TSSs is not dependent on transcription and that transcription limits the maintenance of these marks.

PRDM3/16 regulate chromatin accessibility required for NKX2-1 mediated alveolar epithelial differentiation and function

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. Combined deletion of the histone methyl transferases Prdm3 and Prdm16 in early lung endoderm causes perinatal lethality due to respiratory failure from loss of AT2 cells and the accumulation of partially differentiated AT1 cells. Combination of single-cell RNA-seq, bulk ATAC-seq, and CUT&RUN data demonstrate that PRDM3 and PRDM16 regulate chromatin accessibility at NKX2-1 transcriptional targets critical for perinatal AT2 cell differentiation and surfactant homeostasis. Lineage specific deletion of PRDM3/16 in AT2 cells leads 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.

Review and Evaluate the Bioinformatics Analysis Strategies of ATAC-seq and CUT&Tag Data

Efficient and reliable profiling methods are essential to study epigenetics. Tn5, one of the first identified prokaryotic transposases with high DNA-binding and tagmentation efficiency, is widely adopted in different genomic and epigenomic protocols for high-throughputly exploring the genome and epigenome. Based on Tn5, the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) and the Cleavage Under Targets and Tagmentation (CUT&Tag) were developed to measure chromatin accessibility and detect DNA-protein interactions. These methodologies can be applied to large amounts of biological samples with low-input levels, such as rare tissues, embryos, and sorted single cells. However, fast and proper processing of these epigenomic data has become a bottleneck because massive data production continues to increase quickly. Furthermore, inappropriate data analysis can generate biased or misleading conclusions. Therefore, it is essential to evaluate the performance of Tn5-based ATAC-seq and CUT&Tag data processing bioinformatics tools, many of which were developed mostly for analyzing chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here, we conducted a comprehensive benchmarking analysis to evaluate the performance of eight popular software for processing ATAC-seq and CUT&Tag data. We compared the sensitivity, specificity, and peak width distribution for both narrow-type and broad-type peak calling. We also tested the influence of the availability of control IgG input in CUT&Tag data analysis. Finally, we evaluated the differential analysis strategies commonly used for analyzing the CUT&Tag data. Our study provided comprehensive guidance for selecting bioinformatics tools and recommended analysis strategies, which were implemented into Docker/Singularity images for streamlined data analysis.

Predictive Prioritization of Enhancers Associated with Pancreas Disease Risk

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.

CYP1B1 promotes PARPi-resistance via histone H1.4 interaction and increased chromatin accessibility in ovarian cancer

INTRODUCTION

Ovarian cancer is the most lethal gynecological cancer and presents significant therapeutic challenges. The discovery of synthetic lethality between PARP inhibitors (PARPi) and homologous recombination deficiency marked a new era in treating BRCA1/2-mutated tumors. However, PARPi resistance remains a major clinical challenge.

METHODS

RNA sequencing was used to identify genes altered by PARPi treatment and LC-MS was used to detect proteins interacting with CYP1B1. Resistance mechanisms were explored through ATAC-seq and gene expression manipulation. Additional techniques, including micrococcal nuclease digestion assays, DAPI staining, and fluorescence microscopy, were used to assess changes in nuclear morphology and chromatin accessibility.

RESULTS

The gradual exposure of Olaparib has developed a PARPi-resistant cell line, A2780-OlaR, which exhibits significant upregulation of CYP1B1 at both RNA and protein levels. Down-regulating CYP1B1 expression or using specific inhibitors decreased the cellular response to Olaparib. Linker histone H1.4 was identified as associated with CYP1B1. ATAC-seq showed differential chromatin accessibility between A2780-OlaR and parental cells, indicating that the downregulation of H1.4 was associated with increased chromatin accessibility and higher cell viability after Olaparib treatment.

CONCLUSION

Our findings reveal a novel role for CYP1B1 in driving PARPi resistance through distinct molecular mechanisms in A2780-OlaR. This study highlights the importance of chromatin accessibility in PARPi efficacy and suggests the CYP1B1/H1.4 axis as a promising therapeutic target for overcoming drug resistance in ovarian cancer, offering potentially therapeutic benefits.

SMARCC2 silencing suppresses oncogenic activation through modulation of chromatin accessibility in breast cancer

SWI/SNF chromatin remodeling complexes play a key role in gene transcription as epigenetic regulators and are typically considered to act as tumor suppressors in cancers. Compared to other cancer-related components of the SWI/SNF complex, research on SMARCC2, a component of the initial BAF core, has been relatively limited. This study aimed to elucidate the role of SMARCC2 in breast cancer by employing various in vitro and in vivo methods including cell proliferation assays, mammosphere formation, and xenograft models, complemented by RNA-seq, ATAC-seq, and ChIP analyses. The results showed that SMARCC2 silencing surprisingly led to the suppression of breast tumorigenesis, indicating a pro-tumorigenic function for SMARCC2 in breast cancer, which contrasts with the roles of other SWI/SNF subunits. In addition, SMARCC2 depletion reduces cancer stem cell features of breast cancer cells. Mechanistic study showed that SMARCC2 silencing downregulated the oncogenic Ras-PI3K signaling pathway, likely by directly regulating the chromatin accessibility of the enhancers of the key genes such as PIK3CB. Together, these results expand our understanding of the SWI/SNF complex's role in cancer development and identify SMARCC2 as a promising new target for breast cancer therapies.

Role of the SAF-A SAP domain in X inactivation, transcription, splicing, and cell proliferation

SAF-A is conserved throughout vertebrates and has emerged as an important factor regulating a multitude of nuclear functions, including lncRNA localization, gene expression, and splicing. SAF-A has several functional domains, including an N-terminal SAP domain that binds directly to DNA. Phosphorylation of SAP domain serines S14 and S26 are important for SAF-A localization and function during mitosis, however whether these serines are involved in interphase functions of SAF-A is not known. In this study we tested for the role of the SAP domain, and SAP domain serines S14 and S26 in X chromosome inactivation, protein dynamics, gene expression, splicing, and cell proliferation. Here we show that the SAP domain serines S14 and S26 are required to maintain XIST RNA localization and polycomb-dependent histone modifications on the inactive X chromosome in female cells. In addition, we present evidence that an Xi localization signal resides in the SAP domain. We found that that the SAP domain is not required to maintain gene expression and plays only a minor role in mRNA splicing. In contrast, the SAF-A SAP domain, in particular serines S14 and S26, are required for normal protein dynamics, and to maintain normal cell proliferation. We propose a model whereby dynamic phosphorylation of SAF-A serines S14 and S26 mediates rapid turnover of SAF-A interactions with DNA during interphase.

Systematic identification of post-transcriptional regulatory modules

In our cells, a limited number of RNA binding proteins (RBPs) are responsible for all aspects of RNA metabolism across the entire transcriptome. To accomplish this, RBPs form regulatory units that act on specific target regulons. However, the landscape of RBP combinatorial interactions remains poorly explored. Here, we perform a systematic annotation of RBP combinatorial interactions via multimodal data integration. We build a large-scale map of RBP protein neighborhoods by generating in vivo proximity-dependent biotinylation datasets of 50 human RBPs. In parallel, we use CRISPR interference with single-cell readout to capture transcriptomic changes upon RBP knockdowns. By combining these physical and functional interaction readouts, along with the atlas of RBP mRNA targets from eCLIP assays, we generate an integrated map of functional RBP interactions. We then use this map to match RBPs to their context-specific functions and validate the predicted functions biochemically for four RBPs. This study provides a detailed map of RBP interactions and deconvolves them into distinct regulatory modules with annotated functions and target regulons. This multimodal and integrative framework provides a principled approach for studying post-transcriptional regulatory processes and enriches our understanding of their underlying mechanisms.

Quantifying allele-specific CRISPR editing activity with CRISPECTOR2.0

Off-target effects present a significant impediment to the safe and efficient use of CRISPR-Cas genome editing. Since off-target activity is influenced by the genomic sequence, the presence of sequence variants leads to varying on- and off-target profiles among different alleles or individuals. However, a reliable tool that quantifies genome editing activity in an allelic context is not available. Here, we introduce CRISPECTOR2.0, an extended version of our previously published software tool CRISPECTOR, with an allele-specific editing activity quantification option. CRISPECTOR2.0 enables reference-free, allele-aware, precise quantification of on- and off-target activity, by using de novo sample-specific single nucleotide variant (SNV) detection and statistical-based allele-calling algorithms. We demonstrate CRISPECTOR2.0 efficacy in analyzing samples containing multiple alleles and quantifying allele-specific editing activity, using data from diverse cell types, including primary human cells, plants, and an original extensive human cell line database. We identified instances where an SNV induced changes in the protospacer adjacent motif sequence, resulting in allele-specific editing. Intriguingly, differential allelic editing was also observed in regions carrying distal SNVs, hinting at the involvement of additional epigenetic factors. Our findings highlight the importance of allele-specific editing measurement as a milestone in the adaptation of efficient, accurate, and safe personalized genome editing.

Central role of SUMOylation in the regulation of chromatin interactions and transcriptional outputs of the androgen receptor in prostate cancer cells

The androgen receptor (AR) is pivotal in prostate cancer (PCa) progression and represents a critical therapeutic target. AR-mediated gene regulation involves intricate interactions with nuclear proteins, with many mediating and undergoing post-translational modifications that present alternative therapeutic avenues. Through chromatin proteomics in PCa cells, we identified SUMO ligases together with nuclear receptor coregulators and pioneer transcription factors within the AR's protein network. Intriguingly, this network displayed a significant association with SUMO2/3. To elucidate the influence of SUMOylation on AR chromatin interactions and subsequent gene regulation, we inhibited SUMOylation using ML-792 (SUMOi). While androgens generally facilitated the co-occupancy of SUMO2/3 and AR on chromatin, SUMOi induced divergent effects dependent on the type of AR-binding site (ARB). SUMOi augmented AR's pioneer-like binding on inaccessible chromatin regions abundant in androgen response elements (AREs) and diminished its interaction with accessible chromatin regions sparse in AREs yet rich in pioneer transcription factor motifs. The SUMOi-impacted ARBs divergently influenced AR-regulated genes; those associated with AR-mediated activation played roles in negative regulation of cell proliferation, while those with AR-mediated repression were involved in pattern formation. In conclusion, our findings underscore the pervasive influence of SUMOylation in shaping AR's role in PCa cells, potentially unveiling new therapeutic strategies.

A primer on single-cell RNA-seq analysis using dendritic cells as a case study

Recent advances in single-cell (sc) transcriptomics have revolutionized our understanding of dendritic cells (DCs), pivotal players of the immune system. ScRNA-sequencing (scRNA-seq) has unraveled a previously unrecognized complexity and heterogeneity of DC subsets, shedding light on their ontogeny and specialized roles. However, navigating the rapid technological progress and computational methods can be daunting for researchers unfamiliar with the field. This review aims to provide immunologists with a comprehensive introduction to sc transcriptomic analysis, offering insights into recent developments in DC biology. Addressing common analytical queries, we guide readers through popular tools and methodologies, supplemented with references to benchmarks and tutorials for in-depth understanding. By examining findings from pioneering studies, we illustrate how computational techniques have expanded our knowledge of DC biology. Through this synthesis, we aim to equip researchers with the necessary tools and knowledge to navigate and leverage scRNA-seq for unraveling the intricacies of DC biology and advancing immunological research.

Flexible use of conserved motif vocabularies constrains genome access in cell type evolution

Cell types evolve into a hierarchy with related types grouped into families. How cell type diversification is constrained by the stable separation between families over vast evolutionary times remains unknown. Here, integrating single-nucleus multiomic sequencing and deep learning, we show that hundreds of sequence features (motifs) divide into distinct sets associated with accessible genomes of specific cell type families. This division is conserved across highly divergent, early-branching animals including flatworms and cnidarians. While specific interactions between motifs delineate cell type relationships within families, surprisingly, these interactions are not conserved between species. Consistently, while deep learning models trained on one species can predict accessibility of other species' sequences, their predictions frequently rely on distinct, but synonymous, motif combinations. We propose that long-term stability of cell type families is maintained through genome access specified by conserved motif sets, or 'vocabularies', whereas cell types diversify through flexible use of motifs within each set.

Genotype inference from aggregated chromatin accessibility data reveals genetic regulatory mechanisms

Background

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

Results

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

Conclusions

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

Inhibition of novel human-HPV hybrid ecDNA enhancers reduces oncogene expression and tumor growth in oropharyngeal cancer

Extrachromosomal circular DNA (ecDNA) have been found in most types of human cancers, and ecDNA incorporating viral genomes has recently been described, specifically in human papillomavirus (HPV)-mediated oropharyngeal cancer (OPC). However, the molecular mechanisms of human-viral hybrid ecDNA (hybrid ecDNA) for carcinogenesis remains elusive. We characterized the epigenetic status of hybrid ecDNA using HPVOPC cell lines and patient-derived tumor xenografts, identifying HPV oncogenes E6/E7 in hybrid ecDNA were flanked by novel somatic DNA enhancers and HPV L1 enhancers, with strong cis-interaction. Targeting of these enhancers by clustered regularly interspaced short palindromic repeats interference or hybrid ecDNA by bromodomain and extra-terminal inhibitor reduced E6/E7 expression, and significantly inhibited in vitro and/or in vivo growth only in ecDNA(+) models. HPV DNA in hybrid ecDNA structures are associated with novel somatic and HPV enhancers in hybrid ecDNA that drive HPV ongogene expression and carcinogenesis, and can be targeted with ecDNA disrupting therapeutics.

Answering open questions in biology using spatial genomics and structured methods

Genomics methods have uncovered patterns in a range of biological systems, but obscure important aspects of cell behavior: the shapes, relative locations, movement, and interactions of cells in space. Spatial technologies that collect genomic or epigenomic data while preserving spatial information have begun to overcome these limitations. These new data promise a deeper understanding of the factors that affect cellular behavior, and in particular the ability to directly test existing theories about cell state and variation in the context of morphology, location, motility, and signaling that could not be tested before. Rapid advancements in resolution, ease-of-use, and scale of spatial genomics technologies to address these questions also require an updated toolkit of statistical methods with which to interrogate these data. We present a framework to respond to this new avenue of research: four open biological questions that can now be answered using spatial genomics data paired with methods for analysis. We outline spatial data modalities for each open question that may yield specific insights, discuss how conflicting theories may be tested by comparing the data to conceptual models of biological behavior, and highlight statistical and machine learning-based tools that may prove particularly helpful to recover biological understanding.

Deletion of Trps1 regulatory elements recapitulates postnatal hip joint abnormalities and growth retardation of Trichorhinophalangeal syndrome in mice

Trichorhinophalangeal syndrome (TRPS) is a genetic disorder caused by point mutations or deletions in the gene-encoding transcription factor TRPS1. TRPS patients display a range of skeletal dysplasias, including reduced jaw size, short stature, and a cone-shaped digit epiphysis. Certain TRPS patients experience early onset coxarthrosis that leads to a devastating drop in their daily activities. The etiologies of congenital skeletal abnormalities of TRPS were revealed through the analysis of Trps1 mutant mouse strains. However, early postnatal lethality in Trps1 knockout mice has hampered the study of postnatal TRPS pathology. Here, through epigenomic analysis we identified two previously uncharacterized candidate gene regulatory regions in the first intron of Trps1. We deleted these regions, either individually or simultaneously, and examined their effects on skeletal morphogenesis. Animals that were deleted individually for either region displayed only modest phenotypes. In contrast, the Trps1Δint/Δint mouse strain with simultaneous deletion of both genomic regions exhibit postnatal growth retardation. This strain displayed delayed secondary ossification center formation in the long bones and misshaped hip joint development that resulted in acetabular dysplasia. Reducing one allele of the Trps1 gene in Trps1Δint mice resulted in medial patellar dislocation that has been observed in some patients with TRPS. Our novel Trps1 hypomorphic strain recapitulates many postnatal pathologies observed in human TRPS patients, thus positioning this strain as a useful animal model to study postnatal TRPS pathogenesis. Our observations also suggest that Trps1 gene expression is regulated through several regulatory elements, thus guaranteeing robust expression maintenance in skeletal cells.

Detecting Small Cell Transformation in Patients with Advanced EGFR Mutant Lung Adenocarcinoma through Epigenomic cfDNA Profiling

PURPOSE

Histologic transformation to small cell lung cancer (SCLC) is a mechanism of treatment resistance in patients with advanced oncogene-driven lung adenocarcinoma (LUAD) that currently requires histologic review for diagnosis. Herein, we sought to develop an epigenomic cell-free DNA (cfDNA)-based approach to noninvasively detect small cell transformation in patients with EGFR mutant (EGFRm) LUAD.

EXPERIMENTAL DESIGN

To characterize the epigenomic landscape of transformed (t)SCLC relative to LUAD and de novo SCLC, we performed chromatin immunoprecipitation sequencing (ChIP-seq) to profile the histone modifications H3K27ac, H3K4me3, and H3K27me3; methylated DNA immunoprecipitation sequencing (MeDIP-seq); assay for transposase-accessible chromatin sequencing; and RNA sequencing on 26 lung cancer patient-derived xenograft (PDX) tumors. We then generated and analyzed H3K27ac ChIP-seq, MeDIP-seq, and whole genome sequencing cfDNA data from 1 mL aliquots of plasma from patients with EGFRm LUAD with or without tSCLC.

RESULTS

Analysis of 126 epigenomic libraries from the lung cancer PDXs revealed widespread epigenomic reprogramming between LUAD and tSCLC, with a large number of differential H3K27ac (n = 24,424), DNA methylation (n = 3,298), and chromatin accessibility (n = 16,352) sites between the two histologies. Tumor-informed analysis of each of these three epigenomic features in cfDNA resulted in accurate noninvasive discrimination between patients with EGFRm LUAD versus tSCLC [area under the receiver operating characteristic curve (AUROC) = 0.82-0.87]. A multianalyte cfDNA-based classifier integrating these three epigenomic features discriminated between EGFRm LUAD versus tSCLC with an AUROC of 0.94.

CONCLUSIONS

These data demonstrate the feasibility of detecting small cell transformation in patients with EGFRm LUAD through epigenomic cfDNA profiling of 1 mL of patient plasma.

Androgen receptor monomers and dimers regulate opposing biological processes in prostate cancer cells

Most prostate cancers express the androgen receptor (AR), and tumor growth and progression are facilitated by exceptionally low levels of systemic or intratumorally produced androgens. Thus, absolute inhibition of the androgen signaling axis remains the goal of current therapeutic approaches to treat prostate cancer (PCa). Paradoxically, high dose androgens also exhibit considerable efficacy as a treatment modality in patients with late-stage metastatic PCa. Here we show that low levels of androgens, functioning through an AR monomer, facilitate a non-genomic activation of the mTOR signaling pathway to drive proliferation. Conversely, high dose androgens facilitate the formation of AR dimers/oligomers to suppress c-MYC expression, inhibit proliferation and drive a transcriptional program associated with a differentiated phenotype. These findings highlight the inherent liabilities in current approaches used to inhibit AR action in PCa and are instructive as to strategies that can be used to develop new therapeutics for this disease and other androgenopathies.

Dynamic Foxp3-chromatin interaction controls tunable Treg cell function

Nuclear factor Foxp3 determines regulatory T (Treg) cell fate and function via mechanisms that remain unclear. Here, we investigate the nature of Foxp3-mediated gene regulation in suppressing autoimmunity and antitumor immune response. Contrasting with previous models, we find that Foxp3-chromatin binding is regulated by Treg activation states, tumor microenvironment, and antigen and cytokine stimulations. Proteomics studies uncover dynamic proteins within Foxp3 proximity upon TCR or IL-2 receptor signaling in vitro, reflecting intricate interactions among Foxp3, signal transducers, and chromatin. Pharmacological inhibition and genetic knockdown experiments indicate that NFAT and AP-1 protein Batf are required for enhanced Foxp3-chromatin binding in activated Treg cells and tumor-infiltrating Treg cells to modulate target gene expression. Furthermore, mutations at the Foxp3 DNA-binding domain destabilize the Foxp3-chromatin association. These representative settings delineate context-dependent Foxp3-chromatin interaction, suggesting that Foxp3 associates with chromatin by hijacking DNA-binding proteins resulting from Treg activation or differentiation, which is stabilized by direct Foxp3-DNA binding, to dynamically regulate Treg cell function according to immunological contexts.

Stem cells tightly regulate dead cell clearance to maintain tissue fitness

Billions of cells are eliminated daily from our bodies1-4. Although macrophages and dendritic cells are dedicated to migrating and engulfing dying cells and debris, many epithelial and mesenchymal tissue cells can digest nearby apoptotic corpses1-4. How these non-motile, non-professional phagocytes sense and eliminate dying cells while maintaining their normal tissue functions is unclear. Here we explore the mechanisms that underlie their multifunctionality by exploiting the cyclical bouts of tissue regeneration and degeneration during hair cycling. We show that hair follicle stem cells transiently unleash phagocytosis at the correct time and place through local molecular triggers that depend on both lipids released by neighbouring apoptotic corpses and retinoids released by healthy counterparts. We trace the heart of this dual ligand requirement to RARγ-RXRα, whose activation enables tight regulation of apoptotic cell clearance genes and provides an effective, tunable mechanism to offset phagocytic duties against the primary stem cell function of preserving tissue integrity during homeostasis. Finally, we provide functional evidence that hair follicle stem cell-mediated phagocytosis is not simply redundant with professional phagocytes but rather has clear benefits to tissue fitness. Our findings have broad implications for other non-motile tissue stem or progenitor cells that encounter cell death in an immune-privileged niche.

A multi-omics approach to reveal critical mechanisms of activator protein 1 (AP-1)

The Activator Protein 1 (AP-1) transcription factor complex plays a pivotal role in the regulation of cancer-related genes, influencing cancer cell proliferation, invasion, migration, angiogenesis, and apoptosis. Composed of multiple subunits, AP-1 has diverse roles across different cancer types and environmental contexts, but its specific mechanisms remain unclear. The advent of multi-omics approaches has shed light on a more comprehensive understanding of AP-1's role and mechanism in gene regulation. This review collates recent genome-wide data on AP-1 and provides an overview of its expression, structure, function, and interaction across different diseases. An examination of these findings can illuminate the intricate nature of AP-1 regulation and its significant involvement in the progression of different diseases. Moreover, we discuss the potential use of AP-1 as a target for individual therapy and explore the various challenges associated with such an approach. Ultimately, this review provides valuable insights into the biology of AP-1 and its potential as a therapeutic target for cancer and disease treatments.

MYBL2 Drives Prostate Cancer Plasticity: Inhibiting Its Transcriptional Target CDK2 for RB1-Deficient Neuroendocrine Prostate Cancer

Phenotypic plasticity is a recognized mechanism driving therapeutic resistance in patients with prostate cancer. Although underlying molecular causations driving phenotypic plasticity have been identified, therapeutic success is yet to be achieved. To identify putative master regulator transcription factors (MR-TF) driving phenotypic plasticity in prostate cancer, this work utilized a multiomic approach using genetically engineered mouse models of prostate cancer combined with patient data to identify MYB proto-oncogene like 2 (MYBL2) as a significantly enriched transcription factor in prostate cancer exhibiting phenotypic plasticity. Genetic inhibition of Mybl2 using independent murine prostate cancer cell lines representing phenotypic plasticity demonstrated Mybl2 loss significantly decreased in vivo growth as well as cell fitness and repressed gene expression signatures involved in pluripotency and stemness. Because MYBL2 is currently not druggable, a MYBL2 gene signature was employed to identify cyclin-dependent kinase-2 (CDK2) as a potential therapeutic target. CDK2 inhibition phenocopied genetic loss of Mybl2 and significantly decreased in vivo tumor growth associated with enrichment of DNA damage. Together, this work demonstrates MYBL2 as an important MR-TF driving phenotypic plasticity in prostate cancer. Furthermore, high MYBL2 activity identifies prostate cancer that would be responsive to CDK2 inhibition.

SIGNIFICANCE

Prostate cancers that escape therapy targeting the androgen receptor signaling pathways via phenotypic plasticity are currently untreatable. Our study identifies MYBL2 as a MR-TF in phenotypic plastic prostate cancer and implicates CDK2 inhibition as a novel therapeutic target for this most lethal subtype of prostate cancer.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Fibroblasts as an in vitro model of circadian genetic and genomic studies

Bipolar disorder (BD) is a heritable disorder characterized by shifts in mood that manifest in manic or depressive episodes. Clinical studies have identified abnormalities of the circadian system in BD patients as a hallmark of underlying pathophysiology. Fibroblasts are a well-established in vitro model for measuring circadian patterns. We set out to examine the underlying genetic architecture of circadian rhythm in fibroblasts, with the goal to assess its contribution to the polygenic nature of BD disease risk. We collected, from primary cell lines of 6 healthy individuals, temporal genomic features over a 48 h period from transcriptomic data (RNA-seq) and open chromatin data (ATAC-seq). The RNA-seq data showed that only a limited number of genes, primarily the known core clock genes such as ARNTL, CRY1, PER3, NR1D2 and TEF display circadian patterns of expression consistently across cell cultures. The ATAC-seq data identified that distinct transcription factor families, like those with the basic helix-loop-helix motif, were associated with regions that were increasing in accessibility over time. Whereas known glucocorticoid receptor target motifs were identified in those regions that were decreasing in accessibility. Further evaluation of these regions using stratified linkage disequilibrium score regression analysis failed to identify a significant presence of them in the known genetic architecture of BD, and other psychiatric disorders or neurobehavioral traits in which the circadian rhythm is affected. In this study, we characterize the biological pathways that are activated in this in vitro circadian model, evaluating the relevance of these processes in the context of the genetic architecture of BD and other disorders, highlighting its limitations and future applications for circadian genomic studies.

Integrated multi-omics analysis identifies features that predict human pluripotent stem cell-derived progenitor differentiation to cardiomyocytes

Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) are advancing cardiovascular development and disease modeling, drug testing, and regenerative therapies. However, hPSC-CM production is hindered by significant variability in the differentiation process. Establishment of early quality markers to monitor lineage progression and predict terminal differentiation outcomes would address this robustness and reproducibility roadblock in hPSC-CM production. An integrated transcriptomic and epigenomic analysis assesses how attributes of the cardiac progenitor cell (CPC) affect CM differentiation outcome. Resulting analysis identifies predictive markers of CPCs that give rise to high purity CM batches, including TTN, TRIM55, DGKI, MEF2C, MAB21L2, MYL7, LDB3, SLC7A11, and CALD1. Predictive models developed from these genes provide high accuracy in determining terminal CM purities at the CPC stage. Further, insights into mechanisms of batch failure and dominant non-CM cell types generated in failed batches are elucidated. Namely EMT, MAPK, and WNT signaling emerge as significant drivers of batch divergence, giving rise to off-target populations of fibroblasts/mural cells, skeletal myocytes, epicardial cells, and a non-CPC SLC7A11+ subpopulation. This study demonstrates how integrated multi-omic analysis of progenitor cells can identify quality attributes of that progenitor and predict differentiation outcomes, thereby improving differentiation protocols and increasing process robustness.

Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic β-cell identity and function. Elimination of Lkb1 from the β-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of β cell-selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from β-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors (TFs) important for β-cell identity, such as FOXA, MAFA or RFX6, and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates β-cell identity and function.

An autoimmune transcriptional circuit drives FOXP3+ regulatory T cell dysfunction

Autoimmune diseases, among the most common disorders of young adults, are mediated by genetic and environmental factors. Although CD4+FOXP3+ regulatory T cells (Tregs) play a central role in preventing autoimmunity, the molecular mechanism underlying their dysfunction is unknown. Here, we performed comprehensive transcriptomic and epigenomic profiling of Tregs in the autoimmune disease multiple sclerosis (MS) to identify critical transcriptional programs regulating human autoimmunity. We found that up-regulation of a primate-specific short isoform of PR domain zinc finger protein 1 (PRDM1-S) induces expression of serum and glucocorticoid-regulated kinase 1 (SGK1) independent from the evolutionarily conserved long PRDM1, which led to destabilization of forkhead box P3 (FOXP3) and Treg dysfunction. This aberrant PRDM1-S/SGK1 axis is shared among other autoimmune diseases. Furthermore, the chromatin landscape profiling in Tregs from individuals with MS revealed enriched activating protein-1 (AP-1)/interferon regulatory factor (IRF) transcription factor binding as candidate upstream regulators of PRDM1-S expression and Treg dysfunction. Our study uncovers a mechanistic model where the evolutionary emergence of PRDM1-S and epigenetic priming of AP-1/IRF may be key drivers of dysfunctional Tregs in autoimmune diseases.

Clonal evolution of the 3D chromatin landscape in patients with relapsed pediatric B-cell acute lymphoblastic leukemia

Relapsed pediatric B-cell acute lymphoblastic leukemia (B-ALL) remains one of the leading causes of cancer mortality in children. We performed Hi-C, ATAC-seq, and RNA-seq on 12 matched diagnosis/relapse pediatric leukemia specimens to uncover dynamic structural variants (SVs) and 3D chromatin rewiring that may contribute to relapse. While translocations are assumed to occur early in leukemogenesis and be maintained throughout progression, we discovered novel, dynamic translocations and confirmed several fusion transcripts, suggesting functional and therapeutic relevance. Genome-wide chromatin remodeling was observed at all organizational levels: A/B compartments, TAD interactivity, and chromatin loops, including some loci shared by 25% of patients. Shared changes were found to drive the expression of genes/pathways previously implicated in resistance as well as novel therapeutic candidates, two of which (ATXN1 and MN1) we functionally validated. Overall, these results demonstrate chromatin reorganization under the selective pressure of therapy and offer the potential for discovery of novel therapeutic interventions.