High Chromatin Accessibility Research Articles

The single-molecule accessibility landscape of newly replicated mammalian chromatin.

We present replication-aware single-molecule accessibility mapping (RASAM), a method to nondestructively measure replication status and protein-DNA interactions on chromatin genome-wide. Using RASAM, we uncover a genome-wide state of single-molecule "hyperaccessibility" post-replication that resolves over several hours. Combining RASAM with cellular models for rapid protein degradation, we demonstrate that histone chaperone CAF-1 reduces nascent chromatin accessibility by filling single-molecular "gaps" and generating closely spaced dinucleosomes on replicated DNA. At cis-regulatory elements, we observe unique modes by which nascent chromatin hyperaccessibility resolves: at CCCTC-binding factor (CTCF)-binding sites, CTCF and nucleosomes compete, reducing CTCF occupancy and motif accessibility post-replication; at active transcription start sites, high chromatin accessibility is maintained, implying rapid re-establishment of nucleosome-free regions. Our study introduces a new paradigm for studying replicated chromatin fiber organization. More broadly, we uncover a unique organization of newly replicated chromatin that must be reset by active processes, providing a substrate for epigenetic reprogramming.

Single-cell multi-modal chromatin profiles revealing epigenetic regulations of cells in hepatocellular carcinoma.

Various epigenetic regulations systematically govern gene expression in cells involving various biological processes. Dysregulation of the epigenome leads to aberrant transcriptional programs and subsequently results in diseases, such as cancer. Therefore, comprehensive profiling epigenomics is essential for exploring the mechanisms underlying gene expression regulation during development and disease. In this study, we developed single-cell chromatin proteins and accessibility tagmentation (scCPA-Tag), a multi-modal single-cell epigenetic profile capturing technique based on barcoded Tn5 transposases and a droplet microfluidics platform. scCPA-Tag enables the simultaneous capture of DNA profiles of histone modification and chromatin accessibility in the same cell. By applying scCPA-Tag to K562 cells and a hepatocellular carcinoma (HCC) sample, we found that the silence of several chromatin-accessible genes can be attributed to lysine-27-trimethylation of the histone H3 tail (H3K27me3) modification. We characterized the epigenetic features of the tumour cells and different immune cell types in the HCC tumour tissue by scCPA-Tag. Besides, a tumour cell subtype (C2) with more aggressive features was identified and characterized by high chromatin accessibility and a lower abundance of H3K27me3 on tumour-promoting genes. Our multi-modal scCPA-Tag provides a comprehensive approach for exploring the epigenetic landscapes of heterogeneous cell types and revealing the mechanisms of gene expression regulation during developmental and pathological processes at the single-cell level. scCPA-Tag offers a highly efficient and high throughput technique to simultaneously profile histone modification and chromatin accessibility within a single cell. scCPA-Tag enables to uncover multiple epigenetic modification features of cellular compositions within tumor tissues. scCPA-Tag facilitates the exploration of the epigenetic landscapes of heterogeneous cell types and provides the mechanisms governing gene expression regulation.

Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish

BackgroundCis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq).ResultsUsing the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions.ConclusionsThus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

Abstract 1695: Genomic clustering tendency of transcription factors reflects phase-separated transcriptional condensates at cancer super-enhancers

Abstract Many transcription factors (TFs) have been shown to bind to super-enhancers, forming transcriptional condensates to activate transcription in various cellular systems including cancer. However, the genomic and epigenomic determinants of phase-separated transcriptional condensate formation remain poorly understood. Questions regarding which TFs tend to associate with transcriptional condensates and what factors influence their association are largely unanswered. Here we systematically analyzed 528 DNA sequence motifs across the human genome and 6,650 TF binding profiles across different cell types to identify the molecular features contributing to the formation of transcriptional condensates. We found that the genomic distributions of sequence motifs for different TFs exhibit distinct clustering tendencies. Notably, TF motifs with a high genomic clustering tendency are significantly associated with super-enhancers. TF binding profiles showing a high genomic clustering tendency are further enriched at cell-type-specific super-enhancers. TFs with a high binding clustering tendency also possess high liquid-liquid phase separation abilities. Compared to non-clustered TF binding, densely clustered TF binding sites are more enriched at cell-type-specific super-enhancers with higher chromatin accessibility, elevated chromatin interaction, and stronger association with cancer outcomes. Our results indicate that the clustered genomic binding patterns and the phase separation properties of TFs collectively contribute to the formation of transcriptional condensates at cancer super-enhancers. Citation Format: Zhenjia Wang, Shengyuan Wang, Chongzhi Zang. Genomic clustering tendency of transcription factors reflects phase-separated transcriptional condensates at cancer super-enhancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1695.

DNA methylation dynamics during yak adipocyte differentiation

In mammals, epigenetic modifications involving DNA methylation are necessary for the completion of the cell differentiation process. However, the global DNA methylation landscape and its dynamics during yak adipocyte differentiation remain unexplored. Here, we performed whole-genome bisulfite sequencing (WGBS) to asses DNA methylation in yak preadipocytes and adipocytes, combining these results with those of our previous studies on changes in chromatin accessibility and gene expression during yak adipogenesis. The results showed that CG methylation levels were lower in promoter than in exon and intron, and gradually decreasing from the distal regions to transcription start site (TSS). There was a significant negative correlation between CG methylation levels located in promoter and gene expression levels. The 46 genes shared by CG differentially methylated regions (DMRs) and differential chromatin accessibility were significantly enriched in Hedgehog and PI3K-Akt signaling pathways. ATAC-seq peaks with high chromatin accessibility located in both promoter (≤ 2 kb from TSS) and distal (> 2 kb from TSS) regions corresponded to low methylation levels. Taken together, these findings demonstrated that DNA methylation and its interactions with chromatin accessibility regulate gene expression during yak adipocyte differentiation, contributing to the understanding of mechanisms of various epigenetic modifications and their interactions in adipogenesis.

Chromatin accessibility analysis suggested vascular induction of the biliary epithelium via the Notch signaling pathway in the human liver

The biliary epithelial cells (cholangiocytes) in the liver originate from undifferentiated liver parenchymal cells (hepatoblasts) that are located adjacent to the portal vein. This differentiation process is driven by Notch signaling, which is recognized for generating salt-and-pepper (fine-grained) patterns, in contrast to one- or two-cell layer (spatially confined) patterning in cholangiocyte differentiation. It is unclear how Notch signaling acts and localizes only in cholangiocytes. A computer simulation study suggested that low production rates of the ligands or receptors of Notch signaling are crucial for the spatially confined patterning, although biochemical examination is lacking. Here, we analyzed a publicly available single-cell ATAC-sequencing dataset from human fetal liver samples. We showed high chromatin accessibility for the ligands only in vascular cells, while that for the receptor is limited to a small population of hepatoblasts. This finding strengthens the previously proposed idea that low production rates of the ligands or receptors of Notch signaling enable vascular induction of cholangiocytes.

Multi-omics integrative analysis revealed characteristic changes in blood cell immunity and amino acid metabolism in a silkworm model of hyperproteinemia

Hyperproteinemia is a serious metabolic disease of both humans and animals characterized by an abnormally high plasma protein concentration (HPPC). Although hyperproteinemia can cause an imbalance in blood cell homeostasis, the functional changes to blood cells remain unclear. Here, a HPPC silkworm model was used to assess changes to the chromatin accessibility and transcript levels of genes related to blood cell metabolism and immune function. The results showed that HPPC enhanced phagocytosis of blood cells, increased chromatin accessibility and transcript levels of genes involved in cell phagocytosis, proliferation, stress, and programmed death, while genes associated with aromatic amino acid metabolism, and antibacterial peptide synthesis were inhibited in blood cells. Further analysis of the chromatin accessibility of the promoter region found that the high chromatin accessibility of genes sensitive to HPPC, was related to histone modifications, including tri-methylation of lysine residue 4 of histone H3 and acetylation of lysine residue 27 of histone H3. Changes to the chromatin accessibility and transcript levels of genes related to immune function and amino acid metabolism in the blood cells of the HPPC silkworm model provided useful references for future studies of the mechanisms underlying epigenomic regulation mediated by hyperproteinemia.

ATACAmp: a tool for detecting ecDNA/HSRs from bulk and single-cell ATAC-seq data

BackgroundHigh oncogene expression in cancer cells is a major cause of rapid tumor progression and drug resistance. Recent cancer genome research has shown that oncogenes as well as regulatory elements can be amplified in the form of extrachromosomal DNA (ecDNA) or subsequently integrated into chromosomes as homogeneously staining regions (HSRs). These genome-level variants lead to the overexpression of the corresponding oncogenes, resulting in poor prognosis. Most existing detection methods identify ecDNA using whole genome sequencing (WGS) data. However, these techniques usually detect many false positive regions owing to chromosomal DNA interference.ResultsIn the present study, an algorithm called “ATACAmp” that can identify ecDNA/HSRs in tumor genomes using ATAC-seq data has been described. High chromatin accessibility, one of the characteristics of ecDNA, makes ATAC-seq naturally enriched in ecDNA and reduces chromosomal DNA interference. The algorithm was validated using ATAC-seq data from cell lines that have been experimentally determined to contain ecDNA regions. ATACAmp accurately identified the majority of validated ecDNA regions. AmpliconArchitect, the widely used ecDNA detecting tool, was used to detect ecDNA regions based on the WGS data of the same cell lines. Additionally, the Circle-finder software, another tool that utilizes ATAC-seq data, was assessed. The results showed that ATACAmp exhibited higher accuracy than AmpliconArchitect and Circle-finder. Moreover, ATACAmp supported the analysis of single-cell ATAC-seq data, which linked ecDNA to specific cells.ConclusionsATACAmp, written in Python, is freely available on GitHub under the MIT license: https://github.com/chsmiss/ATAC-amp. Using ATAC-seq data, ATACAmp offers a novel analytical approach that is distinct from the conventional use of WGS data. Thus, this method has the potential to reduce the cost and technical complexity associated ecDNA analysis.

Genome-wide epigenetic dynamics during postnatal skeletal muscle growth in Hu sheep

Hypertrophy and fiber transformation are two prominent features of postnatal skeletal muscle development. However, the role of epigenetic modifications is less understood. ATAC-seq, whole genome bisulfite sequencing, and RNA-seq were applied to investigate the epigenetic dynamics of muscle in Hu sheep at 3 days, 3 months, 6 months, and 12 months after birth. All 6865 differentially expressed genes were assigned into three distinct tendencies, highlighting the balanced protein synthesis, accumulated immune activities, and restrained cell division in postnatal development. We identified 3742 differentially accessible regions and 11799 differentially methylated regions that were associated with muscle-development-related pathways in certain stages, like D3-M6. Transcription factor network analysis, based on genomic loci with high chromatin accessibility and low methylation, showed that ARID5B, MYOG, and ENO1 were associated with muscle hypertrophy, while NR1D1, FADS1, ZFP36L2, and SLC25A1 were associated with muscle fiber transformation. Taken together, these results suggest that DNA methylation and chromatin accessibility contributed toward regulating the growth and fiber transformation of postnatal skeletal muscle in Hu sheep.

Abnormal genetic and epigenetic patterns of m6A regulators associated with tumor microenvironment in colorectal cancer.

N6-methyladenosine (m6A) has a critical role in the development and progression of cancer. However, the genetic and epigenetic patterns, as well as tumor microenvironment (TME) infiltration characteristics of m6A regulators in colorectal cancer (CRC) remain largely unknown. Molecular patterns of m6A modifications of 24 m6A regulators in CRC samples were evaluated using data from The Cancer Genome Atlas (TCGA). Mutations, copy number variations (CNVs), DNA methylation, and chromatin accessibility were examined to investigate the underlying mechanisms of the aberrant expression of m6A regulators. Correlations between m6A-related genes and TME cell-infiltrating characteristics were evaluated using Tumor Immune Estimation Resource (TIMER). The m6A regulators were frequently dysregulated in CRC, with two downregulated and 16 upregulated. All the m6A regulators had mutations (frequency ranging from 0.9% to 7%), with active mutations tending to occur in RBM15 and inactive mutations in ZC3H13. Only five m6A regulators had CNV frequency greater than 1%: YTHDC2 (2.4%), YTHDF1 (7.0%), YTHDF3 (1.9%), VIRMA (1.7%), and ZC3H13 (3.0%). The copy numbers of these five genes were positively correlated with their expression levels. The m6A regulators frequently showed imbalanced methylation in CRC, with hypomethylation of YTHDF2, IGF2BP3, FTO, and hypermethylation of HNRNPC, METTL3, and WTAP. Most m6A regulators had high chromatin accessibility, which was positively correlated with their gene expression. IGF2BP1 was identified as an independent prognostic factor for overall survival. Moreover, the expression of most m6A regulators was positively correlated with the infiltration of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. Aberrant expression of m6A regulators is associated with mutation, CNV, and chromatin accessibility, owing to both genetic and epigenetic modifications. The TME infiltration characterization of m6A regulators could guide the development of more effective immunotherapy strategies in CRC.

Application of single cell multiomics points to changes in chromatin accessibility near calcitonin receptor like receptor and a possible role for adrenomedullin in the post-shock lung.

Acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) is a commonly occurring sequelae of traumatic injury resulting from indirect insults like hypovolemic shock and/or extrapulmonary sepsis. The high lethality rate associated with these pathologies outlines the importance of clarifying the "priming" effects seen in the post-shock lung microenvironment, which are understood to bring about a dysregulated or overt immune response when triggered by a secondary systemic infectious/septic challenge culminating in ALI. In this pilot project, we test the hypothesis that application of a single cell multiomics approach can elucidate novel phenotype specific pathways potentially contributing to shock-induced ALI/ARDS. Hypovolemic shock was induced in C57BL/6 (wild-type), PD-1, PD-L1, or VISTA gene deficient male mice, 8-12 weeks old. Wild-type sham surgeries function as negative controls. A total of 24-h post-shock rodents were sacrificed, their lungs harvested and sectioned, with pools prepared from 2 mice per background, and flash frozen on liquid nitrogen. N = 2 biological replicates (representing 4 mice total) were achieved for all treatment groups across genetic backgrounds. Samples were received by the Boas Center for Genomics and Human Genetics, where single cell multiomics libraries were prepared for RNA/ATAC sequencing. The analysis pipeline Cell Ranger ARC was implemented to attain feature linkage assessments across genes of interest. Sham (pre-shock) results suggest high chromatin accessibility around calcitonin receptor like receptor (CALCRL) across cellular phenotypes with 17 and 18 feature links, exhibiting positive correlation with gene expression between biological replicates. Similarity between both sample chromatin profiles/linkage arcs is evident. Post-shock wild-type accessibility is starkly reduced across replicates where the number of feature links drops to 1 and 3, again presenting similar replicate profiles. Samples from shocked gene deficient backgrounds displayed high accessibility and similar profiles to the pre-shock lung microenvironment. High pre-shock availability of DNA segments and their positive correlation with CALCRL gene expression suggests an apparent regulatory capacity on transcription. Post-shock gene deficient chromatin profiles presented similar results to that of pre-shock wild-type samples, suggesting an influence on CALCRL accessibility. Key changes illustrated in the pre-ALI context of shock may allow for additional resolution of "priming" and "cellular pre-activation/pre-disposition" processes within the lung microenvironment.

Transcriptomic and Chromatin Landscape Analysis Reveals That Involvement of Pituitary Level Transcription Factors Modulate Incubation Behaviors of Magang Geese.

The incubation behavior of geese seriously affects their egg production performance. Studies on incubation behavior have identified functional genes, but the regulatory architecture relationship between functional genes and chromatin accessibility remains poorly understood. Here, we present an integrated analysis of open chromatin profiles and transcriptome to identify the cis-regulatory element and their potential transcription factors involved in regulating incubation behavior in goose pituitary. Assay for transposase-accessible chromatin sequencing (ATAC-seq) revealed that open chromatin regions increased in the pituitary during the transition from incubation behavior to laying. We identified 920 significant differential accessible regions (DARs) in the pituitary. Compared to the laying stage, most DARs had higher chromatin accessibility in the brooding stage. Motif analysis of open DARs showed that the most significant transcription factor (TF) occupied sites predominantly enriched in motifs binding to the RFX family (RFX5, RFX2, and RFX1). While the majority of TF motifs enriched under sites of the nuclear receptor (NR) family (ARE, GRE, and PGR) in closed DARs at the incubation behavior stage. Footprint analysis indicated that the transcription factor RFX family exhibited higher binding on chromatin at the brooding stage. To further elucidate the effect of changes in chromatin accessibility on gene expression levels, a comparison of the transcriptome revealed 279 differentially expressed genes (DEGs). The transcriptome changes were associated with processes of steroid biosynthesis. By integrating ATAC-seq and RNA-seq, few DARs directly affect incubation behavior by regulating the transcription levels of genes. Five DAR-related DEGs were found to be closely related to maintaining the incubation behavior in geese. Footprinting analysis revealed a set of transcription factors (RFX1, RFX2, RFX3, RFX5, BHLHA15, SIX1, and DUX) which displayed the highest activity at the brooding stage. SREBF2 was predicted to be the unique differentially expressed transcription factor whose mRNA level was down-regulated and enriched in hyper-accessible regions of PRL in the broody stage. In the present study, we comprehensively profiled the transcriptome and chromatin accessibility in the pituitary related to incubation behavior. Our findings provided insight into the identification and analysis of regulatory elements in goose incubation behavior. The epigenetic alterations profiled here can help decipher the epigenetic mechanisms that contribute to the regulation of incubation behavior in birds.

Epigenetic and chromosomal features drive transposon insertion in Drosophila melanogaster.

Transposons are mobile genetic elements prevalent in the genomes of most species. The distribution of transposons within a genome reflects the actions of two opposing processes: initial insertion site selection, and selective pressure from the host. By analyzing whole-genome sequencing data from transposon-activated Drosophila melanogaster, we identified 43 316 de novo and 237 germline insertions from four long-terminal-repeat (LTR) transposons, one LINE transposon (I-element), and one DNA transposon (P-element). We found that all transposon types favored insertion into promoters de novo, but otherwise displayed distinct insertion patterns. De novo and germline P-element insertions preferred replication origins, often landing in a narrow region around transcription start sites and in regions of high chromatin accessibility. De novo LTR transposon insertions preferred regions with high H3K36me3, promoters and exons of active genes; within genes, LTR insertion frequency correlated with gene expression. De novo I-element insertion density increased with distance from the centromere. Germline I-element and LTR transposon insertions were depleted in promoters and exons, suggesting strong selective pressure to remove transposons from functional elements. Transposon movement is associated with genome evolution and disease; therefore, our results can improve our understanding of genome and disease biology.

Upregulated YTHDF1 associates with tumor immune microenvironment in head and neck squamous cell carcinomas.

The nature of the tumor immune microenvironment (TME) is essential for the head and neck squamous cell carcinomas (HNSCC) initiation, prognosis, and response to immunotherapy. However, its gene regulatory network remains to be elucidated. To identify N6-methyladenosine (m6A) regulators that are involved in regulating the HNSCC TME, a computational screen was applied to The Cancer Genome Atlas (TCGA) HNSCC patient samples. The effects of mutation, copy number variation (CNV), and transcriptional regulation on YTHDF1 expression were analyzed. We analyzed the TME infiltration, cancer-immunity cycle activities, and YTHDF1-related Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Among the 24 m6A regulators, 3 factors (YTHDF1, ELAVL1, and METTL3) were highly correlated with TME infiltration. As the top candidate, YTHDF1 was up-regulated and amplified in HNSCC. YTHDF1 promoter gains active histone marks and high chromatin accessibility, which might be transcriptionally activated by SOX2 and TP63. Moreover, YTHDF1 expression significantly associates with tumor malignant phenotype in HNSCC, which has a positive correlation with CD4+ T cells and a negative correlation with CD8+ T cells infiltration. Specifically, YTHDF1 expression is negatively associated with the cancer-immunity cycle and immune checkpoint inhibitors. In terms of the underlying biological mechanisms, YTHDF1 may interact with YTHDF2/3 to regulate several vital immune-related pathways. We identify YTHDF1 associated with TME and elucidate an underlying mechanism of immune escape in HNSCC, which might be used as a predictive marker in guiding immunotherapy.

Integrating ATAC-seq and RNA-seq Reveals the Dynamics of Chromatin Accessibility and Gene Expression in Apple Response to Drought.

Drought resistance in plants is influenced by multiple signaling pathways that involve various transcription factors, many target genes, and multiple types of epigenetic modifications. Studies on epigenetic modifications of drought focus on DNA methylation and histone modifications, with fewer on chromatin remodeling. Changes in chromatin accessibility can play an important role in abiotic stress in plants by affecting RNA polymerase binding and various regulatory factors. However, the changes in chromatin accessibility during drought in apples are not well understood. In this study, the landscape of chromatin accessibility associated with the gene expression of apple (GL3) under drought conditions was analyzed by Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq. Differential analysis between drought treatment and control identified 23,466 peaks of upregulated chromatin accessibility and 2447 peaks of downregulated accessibility. The drought-induced chromatin accessibility changed genes were mainly enriched in metabolism, stimulus, and binding pathways. By combining results from differential analysis of RNA-seq and ATAC-seq, we identified 240 genes with higher chromatin accessibility and increased gene expression under drought conditions that may play important functions in the drought response process. Among them, a total of nine transcription factor genes were identified, including ATHB7, HAT5, and WRKY26. These transcription factor genes are differentially expressed with different chromatin accessibility motif binding loci that may participate in apple response to drought by regulating downstream genes. Our study provides a reference for chromatin accessibility under drought stress in apples and the results will facilitate subsequent studies on chromatin remodelers and transcription factors.

Runx2 regulates chromatin accessibility to direct the osteoblast program at neonatal stages.

SUMMARYThe transcriptional regulator Runx2 (runt-related transcription factor 2) has essential but distinct roles in osteoblasts and chondrocytes in skeletal development. However, Runx2-mediated regulatory mechanisms underlying the distinctive programming of osteoblasts and chondrocytes are not well understood. Here, we perform an integrative analysis to investigate Runx2-DNA binding and chromatin accessibility ex vivo using neonatal osteoblasts and chondrocytes. We find that Runx2 engages with cell-type-distinct chromatin-accessible regions, potentially interacting with different combinations of transcriptional regulators, forming cell-type-specific hotspots, and potentiating chromatin accessibility. Genetic analysis and direct cellular reprogramming studies suggest that Runx2 is essential for establishment of chromatin accessibility in osteoblasts. Functional enhancer studies identify an Sp7 distal enhancer driven by Runx2-dependent binding and osteoblast-specific chromatin accessibility, contributing to normal osteoblast differentiation. Our findings provide a framework for understanding the regulatory landscape encompassing Runx2-mediated and cell-type-distinct enhancer networks that underlie the specification of osteoblasts.

Super-enhancer hijacking LINC01977 promotes malignancy of early-stage lung adenocarcinoma addicted to the canonical TGF-β/SMAD3 pathway

BackgroundLung adenocarcinoma (LUAD) is the leading cause of death worldwide. However, the roles of long noncoding RNAs (lncRNAs) hijacked by super-enhancers (SEs), vital regulatory elements of the epigenome, remain elusive in the progression of LUAD metastasis.MethodsSE-associated lncRNA microarrays were used to identify the dysregulated lncRNAs in LUAD. ChIP-seq, Hi-C data analysis, and luciferase reporter assays were utilized to confirm the hijacking of LINC01977 by SE. The functions and mechanisms of LINC01977 in LUAD were explored by a series of in vitro and in vivo assays.ResultsWe found that LINC01977, a cancer-testis lncRNA, was hijacked by SE, which promoted proliferation and invasion both in vitro and in vivo. LINC01977 interacted with SMAD3 to induce its nuclear transport, which facilitated the interaction between SMAD3 and CBP/P300, thereby regulating the downstream target gene ZEB1. Additionally, SMAD3 up-regulated LINC09177 transcription by simultaneously binding the promoter and SE, which was induced by the infiltration of M2-like tumor-associated macrophages (TAM2), subsequently activating the TGF-β/SMAD3 pathway. Moreover, LINC01977 expression was positively correlated with TAM2 infiltration and SMAD3 expression, especially in early-stage LUAD. Higher chromatin accessibility in the SE region of LINC01977 was observed with high expression of TGF-β. Early-stage LUAD patients with high LIN01977 expression had a shorter disease-free survival.ConclusionsTAM2 infiltration induced a rich TGF-β microenvironment, activating SMAD3 to bind the promoter and the SE of LINC01977, which up-regulated LINC01977 expression. LINC01977 also promoted malignancy via the canonical TGF-β/SMAD3 pathway. LINC01977 hijacked by SE could be a valuable therapeutic target, especially for the treatment of early-stage LUAD.

Xist exerts gene-specific silencing during XCI maintenance and impacts lineage-specific cell differentiation and proliferation during hematopoiesis

X chromosome inactivation (XCI) is a dosage compensation phenomenon that occurs in females. Initiation of XCI depends on Xist RNA, which triggers silencing of one of the two X chromosomes, except for XCI escape genes that continue to be biallelically expressed. In the soma XCI is stably maintained with continuous Xist expression. How Xist impacts XCI maintenance remains an open question. Here we conditionally delete Xist in hematopoietic system of mice and report differentiation and cell cycle defects in female hematopoietic stem and progenitor cells (HSPCs). By utilizing female HSPCs and mouse embryonic fibroblasts, we find that X-linked genes show variable tolerance to Xist loss. Specifically, XCI escape genes exhibit preferential transcriptional upregulation, which associates with low H3K27me3 occupancy and high chromatin accessibility that accommodates preexisting binding of transcription factors such as Yin Yang 1 (YY1) at the basal state. We conclude that Xist is necessary for gene-specific silencing during XCI maintenance and impacts lineage-specific cell differentiation and proliferation during hematopoiesis.

Relationship between ATOH1 and tumor microenvironment in colon adenocarcinoma patients with different microsatellite instability status

BackgroundColon adenocarcinoma (COAD) is one of the major varieties of malignant tumors threatening human health today. Immune checkpoint inhibitors (ICIs) have recently begun to emerge as an effective option for the treatment of COAD patients, but not all patients can benefit from ICI treatment. Previous studies have suggested that ICIs boast significant clinical effects on patients with microsatellite instability-high (MSI-H), while conversely patients with microsatellite-stable/microsatellite instability-low (MSS/MSI-L) have shown limited response.MethodsWe used ATAC-seq, RNA-seq, and mutation data from The Cancer Genome Atlas Colon adenocarcinoma (TCGA-COAD) cohort to perform multi-omics differential analysis on COAD samples with different MSI statuses, then further screened genes by additionally combining these results with survival analysis. We analyzed the effects of the screened genes on the tumor microenvironment and immunogenicity of COAD patients, and subsequently determined their influence on the efficacy of ICIs in COAD patients using a series of predictive indexes.ResultsTwelve genes were screened in the TCGA-COAD cohort, and after the combined survival analysis, we identified ATOH1 as having significant effects. ATOH1 is characterized by high chromatin accessibility, high expression, and high mutation in COAD patients in the MSI-H group. COAD patients with high ATOH1 expression are associated with a better prognosis, unique immune microenvironment, and higher efficacy in ICI treatment. Enrichment analysis showed that COAD patients with high ATOH1 expression displayed significant upregulation in their humoral immunity and other related pathways.ConclusionsWe speculate that ATOH1 may influence the efficacy of ICIs therapy in patients with COAD by affecting the immune microenvironment and immunogenicity of the tumor.

Abstract 6213: Chromatin accessibility landscape of pediatric high-risk B-cell acute lymphoblastic leukemia

Abstract B-cell acute lymphoblastic leukemia (B-ALL) is the most common pediatric malignancy. Although overall cure rate has exceeded 90%, the prognosis of high-risk B-ALL remains poor. Genome sequencing studies have revealed key genetic aberrations in B-ALL. But driver remains unclear in 8% of B-ALL for leukemogenesis, and 40% of high-risk B-ALL for drug resistance. Epigenetic mechanisms contributing to leukemia have recently been recognized, including histone modification and methylation. To date, a comprehensive landscape of chromatin accessibility in B-ALL is still lacking. Here, we performed Assays for Transposase Accessible Chromatin using sequencing (ATAC-seq) on a total of 61 high-risk B-ALLs treated at Shanghai Children’s Medical Center, including major B-ALL subtypes. We generated 144 high quality chromatin accessibility profiles from 79 tumors, including 18 diagnosis-relapse paired tumors as well as 11 diagnosis and 32 relapsed tumors. We observed significantly higher chromatin accessibility in B-ALL as compared to normal B-cells. Functional chromatin state annotation showed a median of 27.1% open chromatin regions (OCRs) in B-ALL were constituted from quiescent regions that were absent for known histone modifications, indicating unveiled role of quiescent regions in B-ALL. We further investigated the allelic imbalanced chromatin accessibility in 32 B-ALL with matched ATAC-seq and whole genome sequencing data. Unexpectedly, we found a median of 10.8% of OCRs showed imbalanced accessibility between two heterozygous alleles. Moreover, the adjacent imbalanced OCRs tend to be in the same topologically associating domain, suggesting the transcription regulation from chromatin accessibility was constrained within 3-D genome architecture. We next looked into the OCRs among B-ALL subtypes. Consistent with the regulatory role of chromatin accessibility on gene transcription, we found subtype specific pattern of OCRs in B-ALL, with open regions from enhancer and bivalent chromatin states showed more power in subtype discrimination. A total of 15516 out of 765788 (2.03%) peaks were identified as subtype specific OCRs. Motif analysis associated these regions to 212 transcription factors (TFs). Consistently, genes regulated by these TFs also exhibited subtype specific expression. Finally, we analyzed differences in OCRs between diagnosis (D) and relapse (R) in each B-ALL subtype. Only 2.44% D-R differential peaks were shared between any two subtypes, suggesting subtype specific role of chromatin accessibility in relapse. Potential target genes were further identified with quantitative trait loci analysis. With data from DepMap project, we verified that the expression of the target genes identified was associated with drug response. These data unveiled potential role of chromatin accessibility in high-risk B-ALL and the treatment response of this malignancy. Citation Format: Han Wang, Huiying Sun, Bilin Liang, Fang Zhang, Fan Yang, Bowen Cui, Lixia Ding, Ronghua Wang, Yanjing Tang, Jianan Rao, Wenting Hu, Shuang Zhao, Wenyan Wu, Benshang Li, Jingyan Tang, Shuhong Shen, Yu Liu. Chromatin accessibility landscape of pediatric high-risk B-cell acute lymphoblastic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6213.