Di- and tri-methylation of histone H3K36 play distinct roles in DNA double-strand break repair.

Histone H3 Lys36 (H3K36) methylation and its associated modifiers are crucial for DNA double-strand break (DSB) repair, but the mechanism governing whether and how different H3K36 methylation forms impact repair pathways is unclear. Here, we unveil the distinct roles of H3K36 dimethylation (H3K36me2) and H3K36 trimethylation (H3K36me3) in DSB repair via non-homologous end joining (NHEJ) or homologous recombination (HR). Yeast cells lacking H3K36me2 or H3K36me3 exhibit reduced NHEJ or HR efficiency. yKu70 and Rfa1 bind H3K36me2- or H3K36me3-modified peptides and chromatin, respectively. Disrupting these interactions impairs yKu70 and Rfa1 recruitment to damaged H3K36me2- or H3K36me3-rich loci, increasing DNA damage sensitivity and decreasing repair efficiency. Conversely, H3K36me2-enriched intergenic regions and H3K36me3-enriched gene bodies independently recruit yKu70 or Rfa1 under DSB stress. Importantly, human KU70 and RPA1, the homologs of yKu70 and Rfa1, exclusively associate with H3K36me2 and H3K36me3 in a conserved manner. These findings provide valuable insights into how H3K36me2 and H3K36me3 regulate distinct DSB repair pathways, highlighting H3K36 methylation as a critical element in the choice of DSB repair pathway.

Role of Mod(mdg4)-67.2 Protein in Interactions between Su(Hw)-Dependent Complexes and Their Recruitment to Chromatin.

Su(Hw) belongs to the class of proteins that organize chromosome architecture, determine promoter activity, and participate in formation of the boundaries/insulators between the regulatory domains. This protein contains a cluster of 12 zinc fingers of the C2H2 type, some of which are responsible for binding to the consensus site. The Su(Hw) protein forms complex with the Mod(mdg4)-67.2 and the CP190 proteins, where the last one binds to all known Drosophila insulators. To further study functioning of the Su(Hw)-dependent complexes, we used the previously described su(Hw)E8 mutation with inactive seventh zinc finger, which produces mutant protein that cannot bind to the consensus site. The present work shows that the Su(Hw)E8 protein continues to directly interact with the CP190 and Mod(mdg4)-67.2 proteins. Through interaction with Mod(mdg4)-67.2, the Su(Hw)E8 protein can be recruited into the Su(Hw)-dependent complexes formed on chromatin and enhance their insulator activity. Our results demonstrate that the Su(Hw) dependent complexes without bound DNA can be recruited to the Su(Hw) binding sites through the specific protein-protein interactions that are stabilized by Mod(mdg4)-67.2.

Modification of the Hi-C Technology for Molecular Genetic Analysis of Formalin-Fixed Paraffin-Embedded Sections of Tumor Tissues.

Molecular genetic analysis of tumor tissues is the most important step towards understanding the mechanisms of cancer development; it is also necessary for the choice of targeted therapy. The Hi-C (high-throughput chromatin conformation capture) technology can be used to detect various types of genomic variants, including balanced chromosomal rearrangements, such as inversions and translocations. We propose a modification of the Hi-C method for the analysis of chromatin contacts in formalin-fixed paraffin-embedded (FFPE) sections of tumor tissues. The developed protocol allows to generate high-quality Hi-C data and detect all types of chromosomal rearrangements. We have analyzed various databases to compile a comprehensive list of translocations that hold clinical importance for the targeted therapy selection. The practical value of molecular genetic testing is its ability to influence the treatment strategies and to provide prognostic insights. Detecting specific chromosomal rearrangements can guide the choice of the targeted therapies, which is a critical aspect of personalized medicine in oncology.

Towards Development of the 4C-Based Method Detecting Interactions of Plasmid DNA with Host Genome.

Chromosome conformation capture techniques have revolutionized our understanding of chromatin architecture and dynamics at the genome-wide scale. In recent years, these methods have been applied to a diverse array of species, revealing fundamental principles of chromosomal organization. However, structural organization of the extrachromosomal entities, like viral genomes or plasmids, and their interactions with the host genome, remain relatively underexplored. In this work, we introduce an enhanced 4C-protocol tailored for probing plasmid DNA interactions. We design specific plasmid vector and optimize protocol to allow high detection rate of contacts between the plasmid and host DNA.

Age-dependent regulation of ELP1 exon 20 splicing in Familial Dysautonomia by RNA Polymerase II kinetics and chromatin structure.

Familial Dysautonomia (FD) is a rare disease caused by ELP1 exon 20 skipping. Here we clarify the role of RNA Polymerase II (RNAPII) and chromatin on this splicing event. A slow RNAPII mutant and chromatin-modifying chemicals that reduce the rate of RNAPII elongation induce exon skipping whereas chemicals that create a more relaxed chromatin exon inclusion. In the brain of a mouse transgenic for the human FD-ELP1 we observed on this gene an age-dependent decrease in the RNAPII density profile that was most pronounced on the alternative exon, a robust increase in the repressive marks H3K27me3 and H3K9me3 and a decrease of H3K27Ac, together with a progressive reduction in ELP1 exon 20 inclusion level. In HEK 293T cells, selective drug-induced demethylation of H3K27 increased RNAPII elongation on ELP1 and SMN2, promoted the inclusion of the corresponding alternative exons, and, by RNA-sequencing analysis, induced changes in several alternative splicing events. These data suggest a co-transcriptional model of splicing regulation in which age-dependent changes in H3K27me3/Ac modify the rate of RNAPII elongation and affect processing of ELP1 alternative exon 20.

The importance of DNA sequence for nucleosome positioning in transcriptional regulation.

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

The PCNA-Pol δ complex couples lagging strand DNA synthesis to parental histone transfer for epigenetic inheritance.

Inheritance of epigenetic information is critical for maintaining cell identity. The transfer of parental histone H3-H4 tetramers, the primary carrier of epigenetic modifications on histone proteins, represents a crucial yet poorly understood step in the inheritance of epigenetic information. Here, we show the lagging strand DNA polymerase, Pol δ, interacts directly with H3-H4 and that the interaction between Pol δ and the sliding clamp PCNA regulates parental histone transfer to lagging strands, most likely independent of their roles in DNA synthesis. When combined, mutations at Pol δ and Mcm2 that compromise parental histone transfer result in a greater reduction in nucleosome occupancy at nascent chromatin than mutations in either alone. Last, PCNA contributes to nucleosome positioning on nascent chromatin. On the basis of these results, we suggest that the PCNA-Pol δ complex couples lagging strand DNA synthesis to parental H3-H4 transfer, facilitating epigenetic inheritance.

The ARK2N-CK2 complex initiates transcription-coupled repair through enhancing the interaction of CSB with lesion-stalled RNAPII.

Transcription is extremely important for cellular processes but can be hindered by RNA polymerase II (RNAPII) pausing and stalling. Cockayne syndrome protein B (CSB) promotes the progression of paused RNAPII or initiates transcription-coupled nucleotide excision repair (TC-NER) to remove stalled RNAPII. However, the specific mechanism by which CSB initiates TC-NER upon damage remains unclear. In this study, we identified the indispensable role of the ARK2N-CK2 complex in the CSB-mediated initiation of TC-NER. The ARK2N-CK2 complex is recruited to damage sites through CSB and then phosphorylates CSB. Phosphorylation of CSB enhances its binding to stalled RNAPII, prolonging the association of CSB with chromatin and promoting CSA-mediated ubiquitination of stalled RNAPII. Consistent with this finding, Ark2n-/- mice exhibit a phenotype resembling Cockayne syndrome. These findings shed light on the pivotal role of the ARK2N-CK2 complex in governing the fate of RNAPII through CSB, bridging a critical gap necessary for initiating TC-NER.

Cell surface CD55 traffics to the nucleus leading to cisplatin resistance and stemness by inducing PRC2 and H3K27 trimethylation on chromatin in ovarian cancer.

BACKGROUND: Platinum resistance is the primary cause of poor survival in ovarian cancer (OC) patients. Targeted therapies and biomarkers of chemoresistance are critical for the treatment of OC patients. Our previous studies identified cell surface CD55, a member of the complement regulatory proteins, drives chemoresistance and maintenance of cancer stem cells (CSCs). CSCs are implicated in tumor recurrence and metastasis in multiple cancers. METHODS: Protein localization assays including immunofluorescence and subcellular fractionation were used to identify CD55 at the cell surface and nucleus of cancer cells. Protein half-life determinations were used to compare cell surface and nuclear CD55 stability. CD55 deletion mutants were generated and introduced into cancer cells to identify the nuclear trafficking code, cisplatin sensitivity, and stem cell frequency that were assayed using in vitro and in vivo models. Detection of CD55 binding proteins was analyzed by immunoprecipitation followed by mass spectrometry. Target pathways activated by CD55 were identified by RNA sequencing. RESULTS: CD55 localizes to the nucleus of a subset of OC specimens, ascites from chemoresistant patients, and enriched in chemoresistant OC cells. We determined that nuclear CD55 is glycosylated and derived from the cell surface pool of CD55. Nuclear localization is driven by a trafficking code containing the serine/threonine (S/T) domain of CD55. Nuclear CD55 is necessary for cisplatin resistance, stemness, and cell proliferation in OC cells. CD55 S/T domain is necessary for nuclear entry and inducing chemoresistance to cisplatin in both in vitro and in vivo models. Deletion of the CD55 S/T domain is sufficient to sensitize chemoresistant OC cells to cisplatin. In the nucleus, CD55 binds and attenuates the epigenetic regulator and tumor suppressor ZMYND8 with a parallel increase in H3K27 trimethylation and members of the Polycomb Repressive Complex 2. CONCLUSIONS: For the first time, we show CD55 localizes to the nucleus in OC and promotes CSC and chemoresistance. Our studies identify a therapeutic mechanism for treating platinum resistant ovarian cancer by blocking CD55 nuclear entry.

Chromatin profiling identifies chondrocyte-specific Sox9 enhancers important for skeletal development.

The transcription factor SRY-related HMG box 9 (Sox9) is essential for chondrogenesis. Mutations in and around SOX9 cause campomelic dysplasia (CD) characterized by skeletal malformations. Although the function of Sox9 in this context is well studied, the mechanisms that regulate Sox9 expression in chondrocytes remain to be elucidated. Here, we have used genome-wide profiling to identify 2 Sox9 enhancers located in a proximal breakpoint cluster responsible for CD. Enhancer activity of E308 (located 308 kb 5' upstream) and E160 (located 160 kb 5' upstream) correlated with Sox9 expression levels, and both enhancers showed a synergistic effect in vitro. While single deletions in mice had no apparent effect, simultaneous deletion of both E308 and E160 caused a dwarf phenotype, concomitant with a reduction of Sox9 expression in chondrocytes. Moreover, bone morphogenetic protein 2-dependent chondrocyte differentiation of limb bud mesenchymal cells was severely attenuated in E308/E160 deletion mice. Finally, we found that an open chromatin region upstream of the Sox9 gene was reorganized in the E308/E160 deletion mice to partially compensate for the loss of E308 and E160. In conclusion, our findings reveal a mechanism of Sox9 gene regulation in chondrocytes that might aid in our understanding of the pathophysiology of skeletal disorders.

Genome wide nucleosome landscape shapes 3D chromatin organization.

The hierarchical chromatin organization begins with formation of nucleosomes, which fold into chromatin domains punctuated by boundaries and ultimately chromosomes. In a hierarchal organization, lower levels shape higher levels. However, the dependence of higher-order 3D chromatin organization on the nucleosome-level organization has not been studied in cells. We investigated the relationship between nucleosome-level organization and higher-order chromatin organization by perturbing nucleosomes across the genome by deleting Imitation SWItch (ISWI) and Chromodomain Helicase DNA-binding (CHD1) chromatin remodeling factors in budding yeast. We find that changes in nucleosome-level properties are accompanied by changes in 3D chromatin organization. Short-range chromatin contacts up to a few kilo-base pairs decrease, chromatin domains weaken, and boundary strength decreases. Boundary strength scales with accessibility and moderately with width of nucleosome-depleted region. Change in nucleosome positioning seems to alter the stiffness of chromatin, which can affect formation of chromatin contacts. Our results suggest a biomechanical "bottom-up" mechanism by which nucleosome distribution across genome shapes 3D chromatin organization.

FOXA2 rewires AP-1 for transcriptional reprogramming and lineage plasticity in prostate cancer.

FOXA family proteins act as pioneer factors by remodeling compact chromatin structures. FOXA1 is crucial for the chromatin binding of the androgen receptor (AR) in both normal prostate epithelial cells and the luminal subtype of prostate cancer (PCa). Recent studies have highlighted the emergence of FOXA2 as an adaptive response to AR signaling inhibition treatments. However, the role of the FOXA1 to FOXA2 transition in regulating cancer lineage plasticity remains unclear. Our study demonstrates that FOXA2 binds to distinct classes of developmental enhancers in multiple AR-independent PCa subtypes, with its binding depending on LSD1. Moreover, we reveal that FOXA2 collaborates with JUN at chromatin and promotes transcriptional reprogramming of AP-1 in lineage-plastic cancer cells, thereby facilitating cell state transitions to multiple lineages. Overall, our findings underscore the pivotal role of FOXA2 as a pan-plasticity driver that rewires AP-1 to induce the differential transcriptional reprogramming necessary for cancer cell lineage plasticity.

Ki67 deficiency impedes chromatin accessibility and BCR gene rearrangement.

The proliferation marker Ki67 has been attributed critical functions in maintaining mitotic chromosome morphology and heterochromatin organization during the cell cycle, indicating a potential role in developmental processes requiring rigid cell-cycle control. Here, we discovered that despite normal fecundity and organogenesis, germline deficiency in Ki67 resulted in substantial defects specifically in peripheral B and T lymphocytes. This was not due to impaired cell proliferation but rather to early lymphopoiesis at specific stages where antigen-receptor gene rearrangements occurred. We identified that Ki67 was required for normal global chromatin accessibility involving regulatory regions of genes critical for checkpoint stages in B cell lymphopoiesis. In line with this, mRNA expression of Rag1 was diminished and gene rearrangement was less efficient in the absence of Ki67. Transgenes encoding productively rearranged immunoglobulin heavy and light chains complemented Ki67 deficiency, completely rescuing early B cell development. Collectively, these results identify a unique contribution from Ki67 to somatic antigen-receptor gene rearrangement during lymphopoiesis.

Single-cell 3D genome structure reveals distinct human pluripotent states.

BACKGROUND: Pluripotent states of embryonic stem cells (ESCs) with distinct transcriptional profiles affect ESC differentiative capacity and therapeutic potential. Although single-cell RNA sequencing has revealed additional subpopulations and specific features of naive and primed human pluripotent stem cells (hPSCs), the underlying mechanisms that regulate their specific transcription and that control their pluripotent states remain elusive. RESULTS: By single-cell analysis of high-resolution, three-dimensional (3D) genomic structure, we herein demonstrate that remodeling of genomic structure is highly associated with the pluripotent states of human ESCs (hESCs). The naive pluripotent state is featured with specialized 3D genomic structures and clear chromatin compartmentalization that is distinct from the primed state. The naive pluripotent state is achieved by remodeling the active euchromatin compartment and reducing chromatin interactions at the nuclear center. This unique genomic organization is linked to enhanced chromatin accessibility on enhancers and elevated expression levels of naive pluripotent genes localized to this region. In contradistinction, the primed state exhibits intermingled genomic organization. Moreover, active euchromatin and primed pluripotent genes are distributed at the nuclear periphery, while repressive heterochromatin is densely concentrated at the nuclear center, reducing chromatin accessibility and the transcription of naive genes. CONCLUSIONS: Our data provide insights into the chromatin structure of ESCs in their naive and primed states, and we identify specific patterns of modifications in transcription and chromatin structure that might explain the genes that are differentially expressed between naive and primed hESCs. Thus, the inversion or relocation of heterochromatin to euchromatin via compartmentalization is related to the regulation of chromatin accessibility, thereby defining pluripotent states and cellular identity.

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.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci.

BACKGROUND: Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space. RESULTS: To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4. CONCLUSIONS: Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

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

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.

Drought response revealed by chromatin organization variation and transcriptional regulation in cotton.

BACKGROUND: Cotton is a major world cash crop and an important source of natural fiber, oil, and protein. Drought stress is becoming a restrictive factor affecting cotton production. To facilitate the development of drought-tolerant cotton varieties, it is necessary to study the molecular mechanism of drought stress response by exploring key drought-resistant genes and related regulatory factors. RESULTS: In this study, two cotton varieties, ZY007 (drought-sensitive) and ZY168 (drought-tolerant), showing obvious phenotypic differences under drought stress, were selected. A total of 25,898 drought-induced genes were identified, exhibiting significant enrichment in pathways related to plant stress responses. Under drought induction, At subgenome expression bias was observed at the whole-genome level, which may be due to stronger inhibition of Dt subgenome expression. A gene co-expression module that was significantly associated with drought resistance was identified. About 90% of topologically associating domain (TAD) boundaries were stable, and 6613 TAD variation events were identified between the two varieties under drought. We identified 92 genes in ZY007 and 98 in ZY168 related to chromatin 3D structural variation and induced by drought stress. These genes are closely linked to the cotton response to drought stress through canonical hormone-responsive pathways, modulation of kinase and phosphatase activities, facilitation of calcium ion transport, and other related molecular mechanisms. CONCLUSIONS: These results lay a foundation for elucidating the molecular mechanism of the cotton drought response and provide important regulatory locus and gene resources for the future molecular breeding of drought-resistant cotton varieties.

Epigenetically active chromatin in neonatal iWAT reveals GABPα as a potential regulator of beige adipogenesis.

Background: Thermogenic beige adipocytes, which dissipate energy as heat, are found in neonates and adults. Recent studies show that neonatal beige adipocytes are highly plastic and contribute to >50% of beige adipocytes in adults. Neonatal beige adipocytes are distinct from recruited beige adipocytes in that they develop independently of temperature and sympathetic innervation through poorly defined mechanisms. Methods: We characterized the neonatal beige adipocytes in the inguinal white adipose tissue (iWAT) of C57BL6 postnatal day 3 and 20 mice (P3 and P20) by imaging, genome-wide RNA-seq analysis, ChIP-seq analysis, qRT-PCR validation, and biochemical assays. Results: We found an increase in acetylated histone 3 lysine 27 (H3K27ac) on the promoter and enhancer regions of beige-specific gene UCP1 in iWAT of P20 mice. Furthermore, H3K27ac ChIP-seq analysis in the iWAT of P3 and P20 mice revealed strong H3K27ac signals at beige adipocyte-associated genes in the iWAT of P20 mice. The integration of H3K27ac ChIP-seq and RNA-seq analysis in the iWAT of P20 mice reveal epigenetically active signatures of beige adipocytes, including oxidative phosphorylation and mitochondrial metabolism. We identify the enrichment of GA-binding protein alpha (GABPα) binding regions in the epigenetically active chromatin regions of the P20 iWAT, particularly on beige genes, and demonstrate that GABPα is required for beige adipocyte differentiation. Moreover, transcriptomic analysis and glucose oxidation assays revealed increased glycolytic activity in the neonatal iWAT from P20. Conclusions: Our findings demonstrate that epigenetic mechanisms regulate the development of peri-weaning beige adipocytes via GABPα. Further studies to better understand the upstream mechanisms that regulate epigenetic activation of GABPα and characterization of the metabolic identity of neonatal beige adipocytes will help us harness their therapeutic potential in metabolic diseases.