Differential Chromatin Accessibility Research Articles

Abstract PR016: Dynamic chromatin landscapes of colorectal cancer development and cell-free DNA fragmentation

Abstract The stepwise progression of colorectal cancer (CRC) from healthy epithelium, to premalignant adenoma, to cancer is accompanied by genome-wide epigenetic reprogramming which may be reflected in circulating cell-free DNA (cfDNA). However, current analyses of these processes have been hampered by complex mixtures of cells in normal, adenoma and cancer tissue samples, as well as absence of suitable adenoma model systems. Cultured colorectal organoids, comprised of pure epithelial cells, could enable insights on the epigenetic changes leading to progression from healthy tissue to cancer. Such features could be instrumental in the early cancer detection using genome-wide analyses of cfDNA fragmentation. Using 43 patient-derived colorectal organoids derived from healthy, adenoma, and cancer tissues, as well as 10 freshly collected peripheral blood mononuclear cell (PBMC) control samples, we examined the chromatin landscape of tumor progression using transposase accessible chromatin analyses with next-generation sequencing. We defined a consensus set of nucleosome depleted regions for each disease state and identified loci with differential accessibility between PBMCs, colon healthy, colon adenoma, and colon cancer tissues. We analyzed these differential loci for fragmentation characteristics in the circulating cfDNA of 250 healthy individuals and 51 patients with metastatic CRC. Across all organoid samples, we identified >35 million nucleosome depleted regions that we coalesced into a consensus set of 62,708 regions. As 78% of these regions could be linked to genes, we assessed whether these signatures identify known pathways involved in colorectal cancer. Within these regions, we identified unique colon tumor-related signatures comprising genes of known and novel pathways, including those involved in WNT, hippo, and RAS signaling. While the majority of chromatin accessibility differences were shared between adenomas and cancers, we identified 895 chromatin changes in adenomas that were not present in either healthy or malignant tissues. In healthy individuals, we found that the cfDNA coverage was negatively correlated with accessibility of white blood cell accessibility loci (Spearman ρ= - 0.60) while in individuals with cancer, changes in chromatin accessibility were associated with independent circulating tumor DNA abundance by droplet digital PCR (R2=0.69, p<0.001). Chromatin accessibility profiles were correlated with genome-wide cell-free DNA fragmentation features and could be used to distinguish healthy individuals from those with colorectal cancer (AUC=0.97). These analyses are consistent with dynamic chromatin remodeling during colorectal tumor development and may provide an avenue for evaluating these through genome-wide analyses of cfDNA fragmentation. Citation Format: Nicholas A. Vulpescu, Zachariah H. Foda, Pieter H.A. Wisse, Christopher Cherry, Jaime E. Medina, Vilmos Adleff, Remond J.A. Fijneman, Robert B. Scharpf, Gerrit A. Meijer, Beatriz Carvalho, Victor E. Velculescu. Dynamic chromatin landscapes of colorectal cancer development and cell-free DNA fragmentation [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr PR016.

EPCO-49. SEX-SPECIFIC CHROMATIN REMODELING DRIVES TUMORIGENESIS IN GLIOBLASTOMA

Abstract Glioblastoma (GBM) is the most prevalent primary malignant brain cancer and has higher incidence and lower survival in males compared to females. The SWI/SNF chromatin remodeling complexes, including canonical BAF (cBAF), polybromo-associated BAF (PBAF), and non-canonical BAF (ncBAF), play a critical role in multiple cancers including gliomas. However, their role in gliomagenesis and progression remains largely unknown. We analyzed 18 human glioma specimens for BAF complex mRNA expression by CNS WHO grade and used these data to guide CRISPRi knockdown of BAF components in syngeneic murine male and female GBM astrocytes. BAF complex knockdown revealed marked sex-specific differences in proliferation. Specifically, Smarce1 or Smarcb1 knockdown (cBAF & PBAF) resulted in increased growth in male astrocytes but decreased growth in females. Brd7 (PBAF) knockdown enhanced growth in both sexes, while Brd9 (ncBAF) had minimal impact. We utilized extreme limiting dilution analysis (ELDA) to assess the clonogenic potential of GBM cells with depleted BAF components. Notably, cells with BRD7 knockdown displayed enhanced sphere-forming capacity in both male and female cells. Male cells with Smarce1or Smarcb1 knockdown showed increased clonogenic frequency, whereas females with Smarce1 knockdown had a marked reduction in sphere formation. Brd9 knockdown did not affect the sphere-forming capacity in male or female cells, in part due to compensatory plasticity of an alternate BAF complex. Epitranscriptomic analysis is ongoing and demonstrates sex-specific differences in chromatin accessibility and downstream gene expression programs based on BAF complex knockdown. Immunoprecipitation coupled with mass spectrometry was performed to elucidate the structure and assembly of the BAF complex in GBM cells. Understanding sex-specific mechanisms of BAF chromatin remodeling complex in glioma pathogenesis may yield critical insights for the development of precision therapeutics that have the potential to improve patient outcomes.

Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages.

Breast cancer (BC) is defined by distinct molecular subtypes with different cells of origin. The transcriptional networks that characterize the subtype-specific tumor-normal lineages are not established. In this work, we applied bulk, single-cell and single-nucleus multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 patients with BC to show characteristic links in gene expression and chromatin accessibility between BC subtypes and their putative cells of origin. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal BC and luminal mature cells and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like (SOX6 and KCNQ3) and luminal A/B (FAM155A and LRP1B) lineages. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like BC, suggesting an altered means of immune dysfunction. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single-cell level is a powerful tool for investigating cancer lineage and highlight transcriptional networks that define basal and luminal BC lineages.

Chromatin Accessibility Plays an Important Epigenetic Role on Antibody Expression From CMV Promoter and DNA Elements Flanking the CHO TI Host Landing-Pad.

Targeted integration (TI) Chinese hamster ovary (CHO) platforms are commonly used for protein expression. However, the impact of epigenetic modifications on protein expression in TI cell lines remains elusive since almost all the epigenetic studies focus on random integration (RI) of the gene of interest and only within the promoter region. To address the impact of epigenetic modifications on TI CHO cells, we utilized a standard mAb-1 to identify and characterize TI clones with the same transgene copy numbers but different levels of transgene transcription and titer. Surprisingly, while CMV promoters were not methylated and histone acetylation/methylation was present, these epigenetic markers did not trend with mRNA transcription and protein expression in our TI model. Instead, ATAC-seq data analysis revealed that differences in chromatin accessibility within the TI site could be a major factor impacting these observed differences. However, neither chromatin accessibility nor histone acetylation/methylation profiles in early cultures were predictive of high-expressing clones early during the CLD process. Finally, modulation of the histone profiles (H3K27ac and H3K4me3) at the CMV promoters within the TI integration site using dCas9 fusion proteins was not effective in further increasing mAb titers which could have been likely due to interference of the dCas9 fusion proteins with transcription from the CMV promoters. Overall, our data suggests increasing chromatin accessibility at the TI site is the most effective way to increase mRNA transcription and hence, productivity in TI cell lines.

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.

Single-cell transcriptome and chromatin accessibility mapping of upper lip and primary palate fusion.

Cleft lip and/or primary palate (CL/P) represent a prevalent congenital malformation, the aetiology of which is highly intricate. Although it is generally accepted that the condition arises from failed fusion between the upper lip and primary palate, the precise mechanism underlying this fusion process remains enigmatic. In this study, we utilized transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) to interrogate lambdoidal junction tissue derived from C57BL/6J mouse embryos at critical stages of embryogenesis (10.5, 11.5 and 12.5 embryonic days). We successfully identified distinct subgroups of mesenchymal and ectodermal cells involved in the fusion process and characterized their unique transcriptional profiles. Furthermore, we conducted cell differentiation trajectory analysis, revealing a dynamic repertoire of genes that are sequentially activated or repressed during pseudotime, facilitating the transition of relevant cell types. Additionally, we employed scATAC data to identify key genes associated with the fusion process and demonstrated differential chromatin accessibility across major cell types. Finally, we constructed a dynamic intercellular communication network and predicted upstream transcriptional regulators of critical genes involved in important signalling pathways. Our findings provide a valuable resource for future studies on upper lip and primary palate development, as well as congenital defects.

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

BackgroundBackfat 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.ResultsATAC-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).ConclusionsThis 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.

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

IntroductionOvarian 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. MethodsRNA 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. ResultsThe 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. ConclusionOur 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.

Abstract P393: Isolevuglandins Regulate Chromatin Accessibility in Hypertensive Human Monocytes

Background: Isolevuglandins (isoLGs) are products of lipid peroxidation that are formed in the presence of reactive oxygen species. IsoLGs are increased in circulating CD14 + monocytes in patients with hypertension and contribute to inflammatory activation. When adducted to histones, isoLGs prevent nucleosome assembly. Moreover, isoLGs induce disassembly of pre-formed nucleosomes. IsoLGs can be induced in cell culture with tert-butyl hydroperoxide (tBHP) and are efficiently scavenged by ethyl-2-hydroxybenzylamine (Et-2HOBA). We hypothesized that isoLGs regulate chromatin accessibility in hypertensive (HTN) monocytes. Methods: Human CD14 + monocytes were isolated from 3 normotensive (NT) and 5 HTN subjects. HTN monocytes were treated with vehicle. NT monocytes were treated with vehicle, tBHP, or co-treated with tBHP and Et-2HOBA. Cells were cultured for 24 hours and submitted for bulk ATAC sequencing to assess regions of open chromatin across the genome. Data were analyzed using Partek Flow software. Reads were aligned to the GRCh38 human genome. Differences in chromatin accessibility were determined by DESeq2 analysis and were considered significant at an adjusted P-value of 0.05 and fold change of 1.5. Reads were filtered to focus on accessible chromatin at transcription start sites (TSS). First, chromatin accessibility of HTN monocytes was compared to NT monocytes. Second, to determine chromatin regions with isoLG-regulated accessibility that are associated with hypertension, we compared control and tBHP+Et2HOBA treated cells to a combined HTN and tBHP treated group. Results: HTN monocytes exhibited increased chromatin accessibility at 131 TSS and decreased accessibility at 56 TSS compared to NT monocytes. Chromatin accessibility was increased at 205 TSS and decreased at 183 TSS in response to both hypertension and isoLGs. Gene set enrichment revealed a significant increase in chromatin accessibility in genes regulating chemokine signaling, RAS signaling, lysosomal activity, and focal adhesion. Conclusions: We determined that significant differences in chromatin accessibility are present in HTN monocytes compared to NT monocytes. IsoLGs regulate chromatin accessibility in HTN monocytes. Moreover, isoLGs regulate accessibility of the TSS of genes involved in inflammation and cell migration. These findings implicate isoLGs as novel post-translational histone modifications that regulate chromatin accessibility and contribute to inflammation in hypertension.

Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders.

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Profiling of chromatin accessibility and regulatory elements toward understanding the mechanisms underlying growth regulation in the Pacific oyster (Crassostrea gigas)

Changes in phenotype are precisely controlled by transcriptional regulators, which depend on chromatin remodeling. Although artificial genetic breeding has significantly improved the growth performance of the Pacific oyster (Crassostrea gigas), the identification of gene regulatory elements and transcription factors related to growth is obscure. In this study, we performed a genome-wide characterization of chromatin accessibility using two strains of C. gigas that exhibit significant differences in growth. The results showed that accessible chromatin regions were mainly distributed in the noncoding regions of the genome, with a greater proportion of distal intergenic regions than that of promoters and introns. A total of 2044 differentially accessible regions and 2574 differentially expressed genes were identified through comparative analysis. Genes annotated from regions with differential chromatin accessibility were enriched in cell adhesion molecules, JAK-STAT, insulin resistance, HIF-1α, and platelet activation signaling pathways. Correlation analysis of ATAC-seq and RNA-seq data showed that differential gene expression was positively correlated with differential chromatin accessibility. Furthermore, motif screening and identification of accessible chromatin regions revealed a set of transcription factors, including LHX3, TCF7L2, and LEF1, that may play important roles in the growth of C. gigas. Integrated analysis with public Hi-C datasets allowed the identification of a differential peak located 47 kb upstream of angiopoietin-2, which may act as a potential enhancer that affects the growth of C. gigas. In summary, our work contributes to a better understanding of growth-related regulation in C. gigas and provides valuable information for deciphering the economic traits underlying growth, which is conducive to the development of the breeding industry of C. gigas in the future.

Loss of cold tolerance is conferred by absence of the WRKY34 promoter fragment during tomato evolution

Natural evolution has resulted in reduced cold tolerance in cultivated tomato (Solanum lycopersicum). Herein, we perform a combined analysis of ATAC-Seq and RNA-Seq in cold-sensitive cultivated tomato and cold-tolerant wild tomato (S. habrochaites). We identify that WRKY34 has the most significant association with differential chromatin accessibility and expression patterns under cold stress. We find that a 60 bp InDel in the WRKY34 promoter causes differences in its transcription and cold tolerance among 376 tomato accessions. This 60 bp fragment contains a GATA cis-regulatory element that binds to SWIBs and GATA29, which synergistically suppress WRKY34 expression under cold stress. Moreover, WRKY34 interferes with the CBF cold response pathway through regulating transcription and protein levels. Our findings emphasize the importance of polymorphisms in cis-regulatory regions and their effects on chromatin structure and gene expression during crop evolution.

ScaDA: A novel statistical method for differential analysis of single-cell chromatin accessibility sequencing data.

Single-cell ATAC-seq sequencing data (scATAC-seq) has been widely used to investigate chromatin accessibility on the single-cell level. One important application of scATAC-seq data analysis is differential chromatin accessibility (DA) analysis. However, the data characteristics of scATAC-seq such as excessive zeros and large variability of chromatin accessibility across cells impose a unique challenge for DA analysis. Existing statistical methods focus on detecting the mean difference of the chromatin accessible regions while overlooking the distribution difference. Motivated by real data exploration that distribution difference exists among cell types, we introduce a novel composite statistical test named "scaDA", which is based on zero-inflated negative binomial model (ZINB), for performing differential distribution analysis of chromatin accessibility by jointly testing the abundance, prevalence and dispersion simultaneously. Benefiting from both dispersion shrinkage and iterative refinement of mean and prevalence parameter estimates, scaDA demonstrates its superiority to both ZINB-based likelihood ratio tests and published methods by achieving the highest power and best FDR control in a comprehensive simulation study. In addition to demonstrating the highest power in three real sc-multiome data analyses, scaDA successfully identifies differentially accessible regions in microglia from sc-multiome data for an Alzheimer's disease (AD) study that are most enriched in GO terms related to neurogenesis and the clinical phenotype of AD, and AD-associated GWAS SNPs.

Century-old chromatin architecture revealed in formalin-fixed vertebrates

Gene expression is regulated by changes in chromatin architecture intrinsic to cellular differentiation and as an active response to environmental stimuli. Chromatin dynamics are a major driver of phenotypic diversity, regulation of development, and manifestation of disease. Remarkably, we know little about the evolutionary dynamics of chromatin reorganisation through time, data essential to characterise the impact of environmental stress during the ongoing biodiversity extinction crisis (20th–21st century). Linking the disparate fields of chromatin biology and museum science through their common use of the preservative formaldehyde (a constituent of formalin), we have generated historical chromatin profiles in museum specimens up to 117 years old. Historical chromatin profiles are reproducible, tissue-specific, sex-specific, and environmental condition-dependent in vertebrate specimens. Additionally, we show that over-fixation modulates differential chromatin accessibility to enable semi-quantitative estimates of relative gene expression in vertebrates and a yeast model. Our approach transforms formalin-fixed biological collections into an accurate, comprehensive, and global record of environmental impact on gene expression and phenotype.

Elucidating allergic reaction mechanisms in response to SARS-CoV-2 mRNA vaccination in adults.

During the COVID-19 pandemic, novel nanoparticle-based mRNA vaccines were developed. A small number of individuals developed allergic reactions to these vaccines although the mechanisms remain undefined. To understand COVID-19 vaccine-mediated allergic reactions, we enrolled 19 participants who developed allergic events within 2 h of vaccination and 13 controls, nonreactors. Using standard hemolysis assays, we demonstrated that sera from allergic participants induced stronger complement activation compared to nonallergic subjects following exvivo vaccine exposure. Vaccine-mediated complement activation correlated with anti-polyethelyne glycol (PEG) IgG (but not IgM) levels while anti-PEG IgE was undetectable in all subjects. Depletion of total IgG suppressed complement activation in select individuals. To investigate the effects of vaccine excipients on basophil function, we employed a validated indirect basophil activation test that stratified the allergic populations into high and low responders. Complement C3a and C5a receptor blockade in this system suppressed basophil response, providing strong evidence for complement involvement in vaccine-mediated basophil activation. Single-cell multiome analysis revealed differential expression of genes encoding the cytokine response and Toll-like receptor (TLR) pathways within the monocyte compartment. Differential chromatin accessibility for IL-13 and IL-1B genes was found in allergic and nonallergic participants, suggesting that invivo, epigenetic modulation of mononuclear phagocyte immunophenotypes determines their subsequent functional responsiveness, contributing to the overall physiologic manifestation of vaccine reactions. These findings provide insights into the mechanisms underlying allergic reactions to COVID-19 mRNA vaccines, which may be used for future vaccine strategies in individuals with prior history of allergies or reactions and reduce vaccine hesitancy.

Transcription Factor Activity Regulating Macrophage Heterogeneity during Skin Wound Healing.

Monocytes and macrophages (Mos/Mϕs) play diverse roles in wound healing by adopting a spectrum of functional phenotypes; however, the regulation of such heterogeneity remains poorly defined. We enhanced our previously published Bayesian inference TF activity model, incorporating both single-cell RNA sequencing and single-cell ATAC sequencing data to infer transcription factor (TF) activity in Mos/Mϕs during skin wound healing. We found that wound Mos/Mϕs clustered into early-stage Mos/Mϕs, late-stage Mϕs, and APCs, and that each cluster showed differential chromatin accessibility and differential predicted TF activity that did not always correlate with mRNA or protein expression. Network analysis revealed two highly connected large communities involving a total of 19 TFs, highlighting TF cooperation in regulating wound Mos/Mϕs. This analysis also revealed a small community populated by NR4A1 and NFKB1, supporting a proinflammatory link between these TFs. Importantly, we validated a proinflammatory role for NR4A1 activity during wound healing, showing that Nr4a1 knockout mice exhibit decreased inflammatory gene expression in early-stage wound Mos/Mϕs, along with delayed wound re-epithelialization and impaired granulation tissue formation. In summary, our study provides insight into TF activity that regulates Mo/Mϕ heterogeneity during wound healing and provides a rational basis for targeting Mo/Mϕ TF networks to alter phenotypes and improve healing.

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

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

ScBlood: A comprehensive single-cell accessible chromatin database of blood cells

The advent of single cell transposase-accessible chromatin sequencing (scATAC-seq) technology enables us to explore the genomic characteristics and chromatin accessibility of blood cells at the single-cell level. To fully make sense of the roles and regulatory complexities of blood cells, it is critical to collect and analyze these rapidly accumulating scATAC-seq datasets at a system level. Here, we present scBlood (https://bio.liclab.net/scBlood/), a comprehensive single-cell accessible chromatin database of blood cells. The current version of scBlood catalogs 770,907 blood cells and 452,247 non-blood cells from ∼400 high-quality scATAC-seq samples covering 30 tissues and 21 disease types. All data hosted on scBlood have undergone preprocessing from raw fastq files and multiple standards of quality control. Furthermore, we conducted comprehensive downstream analyses, including multi-sample integration analysis, cell clustering and annotation, differential chromatin accessibility analysis, functional enrichment analysis, co-accessibility analysis, gene activity score calculation, and transcription factor (TF) enrichment analysis. In summary, scBlood provides a user-friendly interface for searching, browsing, analyzing, visualizing, and downloading scATAC-seq data of interest. This platform facilitates insights into the functions and regulatory mechanisms of blood cells, as well as their involvement in blood-related diseases.

Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars.

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

Abstract PO4-13-11: Epigenetically defined sub-clonal heterogeneity drives therapy resistance in triple-negative breast cancer

Abstract Background Triple-negative breast cancer (TNBC) is a highly aggressive and heterogeneous disease with an average survival of less than 50%. Although immunotherapy is showing promising results in early TNBC, many patients do not respond, and chemotherapy remains the main treatment. High recurrence rates in TNBC are in part due to its inherent molecular heterogeneity and sub-clonal diversity, whereby cells present in minority sub-populations escape therapeutic pressure. Recurrent genomic alterations are not selected for upon therapy resistance in neoadjuvant chemotherapy (NAC) resistant TNBC, however sub-clonal transcriptionally defined cells are “primed” to drive resistance and are pre-existing, suggesting phenotypic heterogeneity is governed by inheritable activation of epigenetically defined regulatory regions. Here, we sought to characterise the sub-clonal and spatial diversity of chemotherapy resistant TNBC to identify the epigenetic hallmarks of lethal sub-populations of cells that may induce cell state transition, therapy resistance and ultimately metastatic progression. Methods Sixteen primary untreated triple-negative breast tumours (eight chemotherapy sensitive (ChemoS) and eight chemotherapy resistant (ChemoR, who died of their disease within three years)) were subjected to cell-matched high-depth multiome single nuclei RNA- and ATAC-sequencing. Driver DNA mutations and copy number alterations (CNA’s) were identified using a clinically validated targeted capture panel. Four tumours were subject to spatial RNA and ATAC profiling. Fourteen TNBC PDX models generated from patients with residual disease post NAC were subjected to single nuclei sequencing and used for validation. Meta-clustering of epithelial cells from the RNA profiles was used to identify sub-populations of cells harbouring distinct epigenetic and transcriptomic features enriched in ChemoR patients. Findings were validated in the BrighTNess TNBC NAC clinical trial and in chemotherapy treated TNBC patients from the SCAN-B study. Results Unsupervised meta-clustering (MC) of epithelial cells in G1 phase of the cell cycle identified nine meta-clusters across the dataset. A distinct immune-like MC was contributed by epithelial G1 cells, suggesting a cross-talk of epithelial cells with spatially proximal immune cells. Two MCs were dominated by cells from ChemoR patients and were defined by pathways related to hypoxia, EMT, and extra-cellular matrix signalling. Using scATAC profiles, we identified that the sub-populations of cells in the ChemoR dominant MCs were underpinned by both promoter and distal differential chromatin accessibility, that mapped to H3K27ac enhancer sites, suggesting a role of distal enhancer chromatin modification in regulating MC transcription. None of the meta-clusters were dominated by tumours with specific genetic mutations or CNA’s. Spatial ATAC and RNA profiling of a subset of tumours identified spatial epigenetically defined regions that showed distinct chromatin profiles, suggesting that the epigenetic regulation of genes that define the MCs are spatially distinct. ChemoR clusters were independently reproduced in a cohort of fourteen residual disease TNBC PDXs. Genes pertaining to the immune-like MC were associated with a good prognosis in independent cohorts of contemporary treated TNBC (SCAN-B) and associated with response to NAC in the BrighTNess trial. Marker genes from ChemoR dominant MCs were associated with poor prognosis and significant levels of residual disease following NAC. Conclusions Our multi-modal integrative analysis reveals unprecedented insight into the role of epigenetically defined spatially distinct, sub-clonal populations that pre-exist in treatment naïve TNBC and predict therapy response. This suggests pre-existing populations of cells have already acquired the epigenetic footprints and hence transcriptomic features to allow them to evade therapy. Citation Format: Fatemeh Ahmadimoghari, Yu Zhang, Ioanna Mavrommati, Shefali Thakur, Chris Starling, Ioannis Roxanis, Syed Haider, Rachael Natrajan. Epigenetically defined sub-clonal heterogeneity drives therapy resistance in triple-negative breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO4-13-11.