Identification of molecular subtypes based on chromatin regulator-related genes and experimental verification of the role of ASCL1 in conferring chemotherapy resistance to breast cancer.

Objective: The aim of this study was to identify the molecular subtypes of breast cancer based on chromatin regulator-related genes. Methods: The RNA sequencing data of The Cancer Genome Atlas-Breast Cancer cohort were obtained from the official website, while the single-cell data were downloaded from the Gene Expression Omnibus database (GSE176078). Validation was performed using the Molecular Taxonomy of Breast Cancer International Consortium dataset. Furthermore, the immune characteristics, tumor stemness, heterogeneity, and clinical characteristics of these molecular subtypes were analyzed. The correlation between chromatin regulators and chemotherapy resistance was examined in vitro using the quantitative real-time polymerase chain reaction (qRT-PCR) and Cell Counting Kit-8 (CCK8) assays. Results: This study identified three stable molecular subtypes with different prognostic and pathological features. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and protein-protein interaction analyses revealed that the differentially expressed genes were associated with disease processes, such as mitotic nuclear division, chromosome segregation, condensed chromosome, and specific chromosome region. The T stage and subtypes were correlated with the clinical features. Tumor heterogeneity (mutant-allele tumor heterogeneity, tumor mutational burden, purity, and homologous recombination deficiency) and tumor stemness (RNA expression-based stemness score, epigenetically regulated RNA expression-based stemness score, DNA methylation-based stemness score, and epigenetically regulated DNA methylation-based stemness score) significantly varied between the three subtypes. Furthermore, Western blotting, qRT-PCR, and CCK8 assays demonstrated that the expression of ASCL1 was positively correlated with chemotherapy resistance in breast cancer. Conclusion: This study identified the subtypes of breast cancer based on chromatin regulators and analyzed their clinical features, gene mutation status, immunophenotype, and drug sensitivity. The results of this study provide effective strategies for assessing clinical prognosis and developing personalized treatment strategies.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs.

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

The landscape of transcriptional profiles in human oocytes with different chromatin configurations.

With increasingly used assisted reproductive technology (ART), the acquisition of high-quality oocytes and early embryos has become the focus of much attention. Studies in mice have found that the transition of chromatin conformation from non-surrounded nucleolus (NSN) to surrounded nucleolus (SN) is essential for oocyte maturation and early embryo development, and similar chromatin transition also exists in human oocytes. In this study, we collected human NSN and SN oocytes and investigated their transcriptome. The analysis of differentially expressed genes showed that epigenetic functions, cyclin-dependent kinases and transposable elements may play important roles in chromatin transition during human oocyte maturation. Our findings provide new insights into the molecular mechanism of NSN-to-SN transition of human oocyte and obtained new clues for improvement of oocyte in vitro maturation technique.

Development of a New Model System to Study Long-Distance Interactions Supported by Architectural Proteins.

Chromatin architecture is critical for the temporal and tissue-specific activation of genes that determine eukaryotic development. The functional interaction between enhancers and promoters is controlled by insulators and tethering elements that support specific long-distance interactions. However, the mechanisms of the formation and maintenance of long-range interactions between genome regulatory elements remain poorly understood, primarily due to the lack of convenient model systems. Drosophila became the first model organism in which architectural proteins that determine the activity of insulators were described. In Drosophila, one of the best-studied DNA-binding architectural proteins, Su(Hw), forms a complex with Mod(mdg4)-67.2 and CP190 proteins. Using a combination of CRISPR/Cas9 genome editing and attP-dependent integration technologies, we created a model system in which the promoters and enhancers of two reporter genes are separated by 28 kb. In this case, enhancers effectively stimulate reporter gene promoters in cis and trans only in the presence of artificial Su(Hw) binding sites (SBS), in both constructs. The expression of the mutant Su(Hw) protein, which cannot interact with CP190, and the mutation inactivating Mod(mdg4)-67.2, lead to the complete loss or significant weakening of enhancer-promoter interactions, respectively. The results indicate that the new model system effectively identifies the role of individual subunits of architectural protein complexes in forming and maintaining specific long-distance interactions in the D. melanogaster model.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity.

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

PAX6 promotes neuroendocrine phenotypes of prostate cancer via enhancing MET/STAT5A-mediated chromatin accessibility.

BACKGROUND: Neuroendocrine prostate cancer (NEPC) is a lethal subset of prostate cancer which is characterized by neuroendocrine differentiation and loss of androgen receptor (AR) signaling. Growing evidence reveals that cell lineage plasticity is crucial in the failure of NEPC therapies. Although studies suggest the involvement of the neural transcription factor PAX6 in drug resistance, its specific role in NEPC remains unclear. METHODS: The expression of PAX6 in NEPC was identified via bioinformatics and immunohistochemistry. CCK8 assay, colony formation assay, tumorsphere formation assay and apoptosis assay were used to illustrate the key role of PAX6 in the progression of in vitro. ChIP and Dual-luciferase reporter assays were conducted to confirm the binding sequences of AR in the promoter region of PAX6, as well as the binding sequences of PAX6 in the promoter regions of STAT5A and MET. For in vivo validation, the xenograft model representing NEPC subtype underwent pathological analysis to verify the significant role of PAX6 in disease progression. Complementary diagnoses were established through public clinical datasets and transcriptome sequencing of specific cell lines. ATAC-seq was used to detect the chromatin accessibility of specific cell lines. RESULTS: PAX6 expression was significantly elevated in NEPC and negatively regulated by AR signaling. Activation of PAX6 in non-NEPC cells led to NE trans-differentiation, while knock-down of PAX6 in NEPC cells inhibited the development and progression of NEPC. Importantly, loss of AR resulted in an enhanced expression of PAX6, which reprogramed the lineage plasticity of prostate cancer cells to develop NE phenotypes through the MET/STAT5A signaling pathway. Through ATAC-seq, we found that a high expression level of PAX6 elicited enhanced chromatin accessibility, mainly through attenuation of H4K20me3, which typically causes chromatin silence in cancer cells. CONCLUSION: This study reveals a novel neural transcription factor PAX6 could drive NEPC progression and suggest that it might serve as a potential therapeutic target for the management of NEPC.

The miR-144/Hmgn2 regulatory axis orchestrates chromatin organization during erythropoiesis.

Differentiation of stem and progenitor cells is a highly regulated process that involves the coordinated action of multiple layers of regulation. Here we show how the post-transcriptional regulatory layer instructs the level of chromatin regulation via miR-144 and its targets to orchestrate chromatin condensation during erythropoiesis. The loss of miR-144 leads to impaired chromatin condensation during erythrocyte maturation. Among the several targets of miR-144 that influence chromatin organization, the miR-144-dependent regulation of Hmgn2 is conserved from fish to humans. Our genetic probing of the miR-144/Hmgn2 regulatory axis establish that intact miR-144 target sites in the Hmgn2 3'UTR are necessary for the proper maturation of erythrocytes in both zebrafish and human iPSC-derived erythroid cells while loss of Hmgn2 rescues in part the miR-144 null phenotype. Altogether, our results uncover miR-144 and its target Hmgn2 as the backbone of the genetic regulatory circuit that controls the terminal differentiation of erythrocytes in vertebrates.

The single-cell chromatin landscape in gonadal cell lineage specification.

Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions.

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Annotation of nuclear lncRNAs based on chromatin interactions.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Disease related changes in ATAC-seq of iPSC-derived motor neuron lines from ALS patients and controls.

Amyotrophic Lateral Sclerosis (ALS), like many other neurodegenerative diseases, is highly heritable, but with only a small fraction of cases explained by monogenic disease alleles. To better understand sporadic ALS, we report epigenomic profiles, as measured by ATAC-seq, of motor neuron cultures derived from a diverse group of 380 ALS patients and 80 healthy controls. We find that chromatin accessibility is heavily influenced by sex, the iPSC cell type of origin, ancestry, and the inherent variance arising from sequencing. Once these covariates are corrected for, we are able to identify ALS-specific signals in the data. Additionally, we find that the ATAC-seq data is able to predict ALS disease progression rates with similar accuracy to methods based on biomarkers and clinical status. These results suggest that iPSC-derived motor neurons recapitulate important disease-relevant epigenomic changes.

Drosophila TET acts with PRC1 to activate gene expression independently of its catalytic activity.

Enzymes of the ten-eleven translocation (TET) family play a key role in the regulation of gene expression by oxidizing 5-methylcytosine (5mC), a prominent epigenetic mark in many species. Yet, TET proteins also have less characterized noncanonical modes of action, notably in Drosophila, whose genome is devoid of 5mC. Here, we show that Drosophila TET activates the expression of genes required for larval central nervous system (CNS) development mainly in a catalytic-independent manner. Genome-wide profiling shows that TET is recruited to enhancer and promoter regions bound by Polycomb group complex (PcG) proteins. We found that TET interacts and colocalizes on chromatin preferentially with Polycomb repressor complex 1 (PRC1) rather than PRC2. Furthermore, PRC1 but not PRC2 is required for the activation of TET target genes. Last, our results suggest that TET and PRC1 binding to activated genes is interdependent. These data highlight the importance of TET noncatalytic function and the role of PRC1 for gene activation in the Drosophila larval CNS.

Mono-methylated histones control PARP-1 in chromatin and transcription.

PARP-1 is central to transcriptional regulation under both normal and stress conditions, with the governing mechanisms yet to be fully understood. Our biochemical and ChIP-seq-based analyses showed that PARP-1 binds specifically to active histone marks, particularly H4K20me1. We found that H4K20me1 plays a critical role in facilitating PARP-1 binding and the regulation of PARP-1-dependent loci during both development and heat shock stress. Here, we report that the sole H4K20 mono-methylase, pr-set7, and parp-1 Drosophila mutants undergo developmental arrest. RNA-seq analysis showed an absolute correlation between PR-SET7- and PARP-1-dependent loci expression, confirming co-regulation during developmental phases. PARP-1 and PR-SET7 are both essential for activating hsp70 and other heat shock genes during heat stress, with a notable increase of H4K20me1 at their gene body. Mutating pr-set7 disrupts monomethylation of H4K20 along heat shock loci and abolish PARP-1 binding there. These data strongly suggest that H4 monomethylation is a key triggering point in PARP-1 dependent processes in chromatin.

EmptyDropsMultiome discriminates real cells from background in single-cell multiomics assays.

Multiomic droplet-based technologies allow different molecular modalities, such as chromatin accessibility and gene expression (scATAC-seq and scRNA-seq), to be probed in the same nucleus. We develop EmptyDropsMultiome, an approach that distinguishes true nuclei-containing droplets from background. Using simulations, we show that EmptyDropsMultiome has higher statistical power and accuracy than existing approaches, including CellRanger-arc and EmptyDrops. On real datasets, we observe that CellRanger-arc misses more than half of the nuclei identified by EmptyDropsMultiome and, moreover, is biased against certain cell types, some of which have a retrieval rate lower than 20%.

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment.

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

VRK1 Regulates Sensitivity to Oxidative Stress by Altering Histone Epigenetic Modifications and the Nuclear Phosphoproteome in Tumor Cells.

The chromatin organization and its dynamic remodeling determine its accessibility and sensitivity to DNA damage oxidative stress, the main source of endogenous DNA damage. We studied the role of the VRK1 chromatin kinase in the response to oxidative stress. which alters the nuclear pattern of histone epigenetic modifications and phosphoproteome pathways. The early effect of oxidative stress on chromatin was studied by determining the levels of 8-oxoG lesions and the alteration of the epigenetic modification of histones. Oxidative stress caused an accumulation of 8-oxoG DNA lesions that were increased by VRK1 depletion, causing a significant accumulation of DNA strand breaks detected by labeling free 3'-DNA ends. In addition, oxidative stress altered the pattern of chromatin epigenetic marks and the nuclear phosphoproteome pathways that were impaired by VRK1 depletion. Oxidative stress induced the acetylation of H4K16ac and H3K9 and the loss of H3K4me3. The depletion of VRK1 altered all these modifications induced by oxidative stress and resulted in losses of H4K16ac and H3K9ac and increases in the H3K9me3 and H3K4me3 levels. All these changes were induced by the oxidative stress in the epigenetic pattern of histones and impaired by VRK1 depletion, indicating that VRK1 plays a major role in the functional reorganization of chromatin in the response to oxidative stress. The analysis of the nuclear phosphoproteome in response to oxidative stress detected an enrichment of the phosphorylated proteins associated with the chromosome organization and chromatin remodeling pathways, which were significantly decreased by VRK1 depletion. VRK1 depletion alters the histone epigenetic pattern and nuclear phosphoproteome pathways in response to oxidative stress. The enzymes performing post-translational epigenetic modifications are potential targets in synthetic lethality strategies for cancer therapies.

Chromatin context-dependent regulation and epigenetic manipulation of prime editing.

We set out to exhaustively characterize the impact of the cis-chromatin environment on prime editing, a precise genome engineering tool. Using a highly sensitive method for mapping the genomic locations of randomly integrated reporters, we discover massive position effects, exemplified by editing efficiencies ranging from ∼0% to 94% for an identical target site and edit. Position effects on prime editing efficiency are well predicted by chromatin marks, e.g., positively by H3K79me2 and negatively by H3K9me3. Next, we developed a multiplex perturbational framework to assess the interaction of trans-acting factors with the cis-chromatin environment on editing outcomes. Applying this framework to DNA repair factors, we identify HLTF as a context-dependent repressor of prime editing. Finally, several lines of evidence suggest that active transcriptional elongation enhances prime editing. Consistent with this, we show we can robustly decrease or increase the efficiency of prime editing by preceding it with CRISPR-mediated silencing or activation, respectively.

Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components.

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

The variation and evolution of complete human centromeres.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

scGHOST: identifying single-cell 3D genome subcompartments.

Single-cell Hi-C (scHi-C) technologies allow for probing of genome-wide cell-to-cell variability in three-dimensional (3D) genome organization from individual cells. Computational methods have been developed to reveal single-cell 3D genome features based on scHi-C, including A/B compartments, topologically associating domains and chromatin loops. However, no method exists for annotating single-cell subcompartments, which is important for understanding chromosome spatial localization in single cells. Here we present scGHOST, a single-cell subcompartment annotation method using graph embedding with constrained random walk sampling. Applications of scGHOST to scHi-C data and contact maps derived from single-cell 3D genome imaging demonstrate reliable identification of single-cell subcompartments, offering insights into cell-to-cell variability of nuclear subcompartments. Using scHi-C data from complex tissues, scGHOST identifies cell-type-specific or allele-specific subcompartments linked to gene transcription across various cell types and developmental stages, suggesting functional implications of single-cell subcompartments. scGHOST is an effective method for annotating single-cell 3D genome subcompartments in a broad range of biological contexts.