Accessible Chromatin Regions Research Articles

Binding of PtoRAP2.12 to demethylated and accessible chromatin regions in the PtoGntK promoter stimulates growth of poplar.

DNA methylation is an essential epigenetic modification for gene regulation in plant growth and development. However, the precise mechanisms of DNA methylation remain poorly understood, especially in woody plants. We employed whole-genome bisulfite sequencing (WGBS), assays for transposase-accessible chromatin using sequencing (ATAC-seq), and RNA-Seq to investigate epigenetic regulatory relationships in Populus tomentosa treated with DNA methylation inhibitor 5-azacitidine. Expression-quantitative trait methylation analysis (eQTM), epigenome-wide association study (EWAS), and joint linkage-linkage disequilibrium mapping were used to explore the epigenetic regulatory genes, and using CRISPR/Cas9 to identify the role of candidate genes. Plant developmental abnormalities occurred when DNA methylation levels were substantially reduced. DNA methylation regulated 112 expressed genes via chromatin accessibility, of which 61 genes were significantly influenced by DNA methylation variation at the population level. One DNA methylation-regulated gene, PtoGntK, was located in a major quantitative trait locus (QTL) for poplar growth. Overexpression and CRISPR/Cas9 of PtoGntK revealed it affected poplar height and stem diameter. The PtoRAP2.12 was found to bind to the demethylated accessible region in the PtoGntK promoter, thereby promoting growth in poplar. This study identified key genes with epigenetic regulation for plant growth and provides insights into epigenetic regulation mechanisms in woody plants.

Chromatin structure and 3D architecture define the differential functions of PU.1 regulatory elements in blood cell lineages

The precise spatiotemporal expression of the hematopoietic ETS transcription factor PU.1, a key determinant of hematopoietic cell fates, is tightly regulated at the chromatin level. However, how chromatin signatures are linked to this dynamic expression pattern across different blood cell lineages remains uncharacterized. Here, we performed an in-depth analysis of the relationships between gene expression, chromatin structure, 3D architecture, and trans-acting factors at PU.1 cis-regulatory elements (PCREs). By identifying phylogenetically conserved DNA elements within chromatin-accessible regions in primary human blood lineages, we discovered multiple novel candidate PCREs within the upstream region of the human PU.1 locus. A subset of these elements localizes within an 8-kb-wide cluster exhibiting enhancer features, including open chromatin, demethylated DNA, enriched enhancer histone marks, present enhancer RNAs, and PU.1 occupation, presumably mediating PU.1 autoregulation. Importantly, we revealed the presence of a common 35-kb-wide CTCF-flanked insulated neighborhood that contains the PCRE cluster (PCREC), forming a chromatin territory for lineage-specific and PCRE-mediated chromatin interactions. These include functional PCRE-promoter interactions in myeloid and B cells that are absent in erythroid and T cells. By correlating chromatin structure and 3D architecture with PU.1 expression in various lineages, we were able to attribute enhancer versus silencer functions to individual elements. Our findings provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in the chromatin regulation of PU.1 expression. This study lays crucial groundwork for additional experimental studies that validate and dissect the role of PCREs in epigenetic regulation of normal and malignant hematopoiesis.

On the Hunt for the Histone Code

Our genome is not made of naked DNA, but a fiber (chromatin) composed of DNA and proteins packaged into our chromosomes. The basic building block of chromatin is the nucleosome, which has two copies of each of the proteins called histones (H2A, H2B, H3, and H4) wrapped by 146 base pairs of DNA. Regions of our genetic material are found between the more open (euchromatin) and more compact (heterochromatin) regions of the genome that can be variably accessible to the underlying genes. Furthermore, post-translational modifications (PTMs) on histones, such as on H3, are critical for regulating chromatin accessibility and gene expression. While site specific antibodies were the tool of choice for histone PTM analysis in the early days (pre-2000s), enter Don Hunt changing the histone PTM field forever. Don’s clever thinking brought new innovative mass spectrometry-based approaches to the epigenetics field. His lab’s effort led to the discovery of many new histone modifications and methods to facilitate the detection and quantification of histone PTMs, which are still considered state of the art in the proteomics field today. Due to Don’s pioneering work in this area, many labs have been able to jump into the epigenetics field and “Hunt” down their own histone targets. A walkthrough of those early histone years in the Hunt Lab are described by three of us who were fortunate enough to be at the right place, at the right time.

INO80 regulates chromatin accessibility to facilitate suppression of sex-linked gene expression during mouse spermatogenesis.

The INO80 protein is the main catalytic subunit of the INO80-chromatin remodeling complex, which is critical for DNA repair and transcription regulation in murine spermatocytes. In this study, we explored the role of INO80 in silencing genes on meiotic sex chromosomes in male mice. INO80 immunolocalization at the XY body in pachytene spermatocytes suggested a role for INO80 in the meiotic sex body. Subsequent deletion of Ino80 resulted in high expression of sex-linked genes. Furthermore, the active form of RNA polymerase II at the sex chromosomes of Ino80-null pachytene spermatocytes indicates incomplete inactivation of sex-linked genes. A reduction in the recruitment of initiators of meiotic sex chromosome inhibition (MSCI) argues for INO80-facilitated recruitment of DNA repair factors required for silencing sex-linked genes. This role of INO80 is independent of a common INO80 target, H2A.Z. Instead, in the absence of INO80, a reduction in chromatin accessibility at DNA repair sites occurs on the sex chromosomes. These data suggest a role for INO80 in DNA repair factor localization, thereby facilitating the silencing of sex-linked genes during the onset of pachynema.

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules—e.g., enhancers and silencers—that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

Regulatory genome annotation of 33 insect species.

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

The HSV-1 encoded CCCTC-binding factor, CTRL2, impacts the nature of viral chromatin during HSV-1 lytic infection.

HSV-1 genomes are rapidly heterochromatinized following entry by host cells to limit viral gene expression. Efficient HSV-1 genome replication requires mechanisms that de-repress chromatin associated with the viral genome. CCCTC-binding factors, or CTCF insulators play both silencing and activating roles in cellular transcriptional regulation. Importantly, the HSV-1 genome encodes several CTCF insulators that flank IE genes, implying that individual HSV-1 encoded CTCF insulators regulate IE transcription during all stages of the HSV-1 life cycle. We previously reported that the HSV-1 encoded CTCF insulator located downstream of the LAT (CTRL2) controlled IE gene silencing during latency. To further characterize the role of this insulator during the lytic infection we leveraged a ΔCTRL2 recombinant virus to show that there was a genome replication defect that stemmed from decreased IE gene expression in fibroblasts and epithelial cells at early times following initiation of infection. Further experiments indicated that the defect in gene expression resulted from chromatin inaccessibility in the absence of the insulator. To elucidate how chromatin accessibility was altered in the absence of the CTRL2 insulator, we showed that enrichment of Alpha-thalassemia/mental retardation, X-linked chromatin remodeler (ATRX), and the histone variant H3.3, both of which are known for their roles in maintaining repressive histone markers on the HSV-1 viral genome were increased on IE regions of HSV-1. Finally, both H3K27me3 and H3K9me3 repressive histone marks remained enriched by 4 hours post infection in the absence of the CTRL2 insulator, confirming that the CTRL2 insulator is required for de-repression of IE genes of viral genomes. To our knowledge these are the first data that show that a specific CTCF insulator in the HSV-1 genome (CTRL2) regulates chromatin accessibility during the lytic infection.

Evolution of plant cell-type-specific cis-regulatory elements.

Cis-regulatory elements (CREs) are critical in regulating gene expression, and yet understanding of CRE evolution remains challenging. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa, integrating data from 103,911 nuclei representing 126 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species (Zea mays, Sorghum bicolor, Panicum miliaceum, and Urochloa fusca). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1-derived epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially silencer-like CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of plant cell-type-specific CRE evolution.

Chromatin Accessibility and Gene Expression Vary Between a New and Evolved Autopolyploid of Arabidopsis arenosa.

Polyploids arise from whole-genome duplication (WGD) events, which have played important roles in genome evolution across eukaryotes. WGD can increase genome complexity, yield phenotypic novelty, and influence adaptation. Neo-polyploids have been reported to often show seemingly stochastic epigenetic and transcriptional changes, but this leaves open the question whether these changes persist in evolved polyploids. A powerful approach to address this is to compare diploids, neo-polyploids, and evolved polyploids of the same species. Arabidopsis arenosa is a species that allows us to do this-natural diploid and autotetraploid populations exist, while neo-tetraploids can be artificially generated. Here, we use ATAC-seq to assay local chromatin accessibility, and RNA-seq to study gene expression on matched leaf and petal samples from diploid, neo-tetraploid and evolved tetraploid A. arenosa. We found over 8,000 differentially accessible chromatin regions across all samples. These are largely tissue specific and show distinct trends across cytotypes, with roughly 70% arising upon WGD. Interestingly, only a small proportion is associated with expression changes in nearby genes. However, accessibility variation across cytotypes associates strongly with the number of nearby transposable elements. Relatively few genes were differentially expressed upon genome duplication, and ∼60% of these reverted to near-diploid levels in the evolved tetraploid, suggesting that most initial perturbations do not last. Our results provide new insights into how epigenomic and transcriptional mechanisms jointly respond to genome duplication and subsequent evolution of autopolyploids, and importantly, show that one cannot be directly predicted from the other.

Untargeted CUT&Tag and BG4 CUT&Tag are both enriched at G-quadruplexes and accessible chromatin.

G-quadruplex DNA structures (G4s) form within single-stranded DNA in nucleosome-free chromatin. As G4s modulate gene expression and genomic stability, genome-wide mapping of G4s has generated strong research interest. Recently, the Cleavage Under Targets and Tagmentation (CUT&Tag) method was performed with the G4-specific BG4 antibody to target Tn5 transposase to G4s. While this method generated a novel high-resolution map of G4s, we unexpectedly observed a strong correlation between the genome-wide signal distribution of BG4 CUT&Tag and accessible chromatin. To examine whether untargeted Tn5 cutting at accessible chromatin contributes to BG4 CUT&Tag signal, we examined the genome-wide distribution of signal from untargeted (i.e. negative control) CUT&Tag datasets. We observed that untargeted CUT&Tag signal distribution was highly similar to both that of accessible chromatin and of BG4 CUT&Tag. We also observed that BG4 CUT&Tag signal increased at mapped G4s, but this increase was accompanied by a concomitant increase in untargeted CUT&Tag at the same loci. Consequently, enrichment of BG4 CUT&Tag over untargeted CUT&Tag was not increased at mapped G4s. These results imply that either the vast majority of accessible chromatin regions contain mappable G4s or that the presence of G4s within accessible chromatin cannot reliably be determined using BG4 CUT&Tag alone.

Large-scale multi-omic analysis identifies noncoding somatic driver mutations and nominates ZFP36L2 as a driver gene for pancreatic ductal adenocarcinoma.

Identification of somatic driver mutations in the noncoding genome remains challenging. To comprehensively characterize noncoding driver mutations for pancreatic ductal adenocarcinoma (PDAC), we first created genome-scale maps of accessible chromatin regions (ACRs) and histone modification marks (HMMs) in pancreatic cell lines and purified pancreatic acinar and duct cells. Integration with whole-genome mutation calls from 506 PDACs revealed 314 ACRs/HMMs significantly enriched with 3,614 noncoding somatic mutations (NCSMs). Functional assessment using massively parallel reporter assays (MPRA) identified 178 NCSMs impacting reporter activity (19.45% of those tested). Focused luciferase validation confirmed negative effects on gene regulatory activity for NCSMs near CDKN2A and ZFP36L2 . For the latter, CRISPR interference (CRISPRi) further identified ZFP36L2 as a target gene (16.0 - 24.0% reduced expression, P = 0.023-0.0047) with disrupted KLF9 binding likely mediating the effect. Our integrative approach provides a catalog of potentially functional noncoding driver mutations and nominates ZFP36L2 as a PDAC driver gene.

Investigating the cis-regulatory basis of C3 and C4 photosynthesis in grasses at single-cell resolution

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-dependent malic enzyme), Panicum miliaceum (NAD-dependent malic enzyme), Urochloa fusca (phosphoenolpyruvate carboxykinase), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of cis-regulatory elements evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

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.

Transcription factor networks disproportionately enrich for heritability of blood cell phenotypes.

Most phenotype-associated genetic variants map to non-coding regulatory regions of the human genome. Moreover, variants associated with blood cell phenotypes are enriched in regulatory regions active during hematopoiesis. To systematically explore the nature of these regions, we developed a highly efficient strategy, Perturb-multiome, that makes it possible to simultaneously profile both chromatin accessibility and gene expression in single cells with CRISPR-mediated perturbation of a range of master transcription factors (TFs). This approach allowed us to examine the connection between TFs, accessible regions, and gene expression across the genome throughout hematopoietic differentiation. We discovered that variants within the TF-sensitive accessible chromatin regions, while representing less than 0.3% of the genome, show a ~100-fold enrichment in heritability across certain blood cell phenotypes; this enrichment is strikingly higher than for other accessible chromatin regions. Our approach facilitates large-scale mechanistic understanding of phenotype-associated genetic variants by connecting key cis-regulatory elements and their target genes within gene regulatory networks.

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.

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.

A STAR for the ages: a 30-year historical perspective of the role of transcription factors in the regulation of steroidogenic acute regulatory gene expression.

The steroidogenic acute regulatory (STAR) protein is an essential cholesterol transporter that shuttles cholesterol from the outer to the inner mitochondrial membrane in the major steroidogenic endocrine organs. It is a key player in the acute regulation of steroid hormone biosynthesis in response to tropic hormone stimulation. Its discovery 30 years ago sparked immediate interest in understanding how STAR action is controlled. Since increased STAR gene expression is a classic feature of the acute regulation of steroidogenesis, a special emphasis was placed on defining the transcriptional regulatory mechanisms that underlie its rapid induction in response to tropic hormone stimulation. These mechanisms include the effects of enhancers, the regulation of chromatin accessibility, the impact of epigenetic factors, and the role of transcription factors. Over the past three decades, understanding the transcription factors that regulate STAR gene expression has been the subject of more than 170 independent scientific publications, making it one of, and if not the best, studied genes in the steroidogenic pathway. This intense research effort has led to the identification of dozens of transcription factors and their related binding sites in STAR 5' flanking (promoter) sequences across multiple species. STAR gene transcription appears to be complex in that a large number of transcription factors have been proposed to interact with either isolated or overlapping regulatory sequences that are tightly clustered over a relatively short promoter region upstream of the STAR transcription start site. Many of these transcription factors appear to work in unique combinatorial codes and are impacted by diverse hormonal and intracellular signaling pathways. This review provides a retrospective overview of the transcription factors proposed to regulate both basal and acute (hormonal) STAR gene expression, and how insights in this area have evolved over the past 30 years.

ΔNp63 bookmarks and creates an accessible epigenetic environment for TGFβ-induced cancer cell stemness and invasiveness

Backgroundp63 is a transcription factor with intrinsic pioneer factor activity and pleiotropic functions. Transforming growth factor β (TGFβ) signaling via activation and cooperative action of canonical, SMAD, and non-canonical, MAP-kinase (MAPK) pathways, elicits both anti- and pro-tumorigenic properties, including cell stemness and invasiveness. TGFβ activates the ΔNp63 transcriptional program in cancer cells; however, the link between TGFβ and p63 in unmasking the epigenetic landscape during tumor progression allowing chromatin accessibility and gene transcription, is not yet reported.MethodsSmall molecule inhibitors, including protein kinase inhibitors and RNA-silencing, provided loss of function analyses. Sphere formation assays in cancer cells, chromatin immunoprecipitation and mRNA expression assays were utilized in order to gain mechanistic evidence. Mass spectrometry analysis coupled to co-immunoprecipitation assays revealed novel p63 interactors and their involvement in p63-dependent transcription.ResultsThe sphere-forming capacity of breast cancer cells was enhanced upon TGFβ stimulation and significantly decreased upon ΔNp63 depletion. Activation of TGFβ signaling via p38 MAPK signaling induced ΔNp63 phosphorylation at Ser 66/68 resulting in stabilized ΔNp63 protein with enhanced DNA binding properties. TGFβ stimulation altered the ratio of H3K27ac and H3K27me3 histone modification marks, pointing towards higher H3K27ac and increased p300 acetyltransferase recruitment to chromatin. By silencing the expression of ΔNp63, the TGFβ effect on chromatin remodeling was abrogated. Inhibition of H3K27me3, revealed the important role of TGFβ as the upstream signal for guiding ΔNp63 to the TGFβ/SMAD gene loci, as well as the indispensable role of ΔNp63 in recruiting histone modifying enzymes, such as p300, to these genomic regions, regulating chromatin accessibility and gene transcription. Mechanistically, TGFβ through SMAD activation induced dissociation of ΔNp63 from NURD or NCOR/SMRT histone deacetylation complexes, while promoted the assembly of ΔNp63-p300 complexes, affecting the levels of histone acetylation and the outcome of ΔNp63-dependent transcription.ConclusionsΔNp63, phosphorylated and recruited by TGFβ to the TGFβ/SMAD/ΔNp63 gene loci, promotes chromatin accessibility and transcription of target genes related to stemness and cell invasion.

Functional mapping of epigenetic regulators uncovers coordinated tumor suppression by the HBO1 and MLL1 complexes.

Epigenetic dysregulation is widespread in cancer. However, the specific epigenetic regulators and the processes they control to drive cancer phenotypes are poorly understood. Here, we employed a novel, scalable and high-throughput in vivo method to perform iterative functional screens of over 250 epigenetic regulatory genes within autochthonous oncogenic KRAS-driven lung tumors. We identified multiple novel epigenetic tumor suppressor and tumor dependency genes. We show that a specific HBO1 complex and the MLL1 complex are among the most impactful tumor suppressive epigenetic regulators in lung. The histone modifications generated by the HBO1 complex are frequently absent or reduced in human lung adenocarcinomas. The HBO1 and MLL1 complexes regulate chromatin accessibility of shared genomic regions, lineage fidelity and the expression of canonical tumor suppressor genes. The HBO1 and MLL1 complexes are epistatic during lung tumorigenesis, and their functional correlation is conserved in human cancer cell lines. Together, these results demonstrate the value of quantitative methods to generate a phenotypic roadmap of epigenetic regulatory genes in tumorigenesis in vivo .

Dissociation of Human and Mouse Tumor Tissue Samples for Single-cell RNA Sequencing.

Human tumor samples hold a plethora of information about their microenvironment and immune repertoire. Effective dissociation of human tissue samples into viable cell suspensions is a required input for the single-cell RNA sequencing (scRNAseq) pipeline. Unlike bulk RNA sequencing approaches, scRNAseq enables us to infer the transcriptional heterogeneity in tumor specimens at the single-cell level. Incorporating this approach in recent years has led to many discoveries, such as identifying immune and tumor cellular states and programs associated with clinical responses to immunotherapies and other types of treatments. Moreover, single-cell technologies applied to dissociated tissues can be used to identify accessible chromatin regions T and B cell receptor repertoire, and the expression of proteins, using DNA barcoded antibodies (CITEseq). The viability and quality of the dissociated sample are critical variables when using these technologies, as these can dramatically affect the cross-contamination of single cells with ambient RNA, the quality of the data, and interpretation. Moreover, long dissociation protocols can lead to the elimination of sensitive cell populations and the upregulation of a stress response gene signature. To overcome these limitations, we devised a rapid universal dissociation protocol, which has been validated on multiple types of human and murine tumors. The process begins with mechanical and enzymatic dissociation, followed by filtration, red blood lysis, and live dead enrichment, suitable for samples with a low input of cells (e.g., needle core biopsies). This protocol ensures a clean and viable single-cell suspension paramount to the successful generation of Gel Bead-In Emulsions (GEMs), barcoding, and sequencing.