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.

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.

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.

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.

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.

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.

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.

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.

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.

The probability of chromatin to be at the nuclear lamina has no systematic effect on its transcription level in fruit flies.

BACKGROUND: Multiple studies have demonstrated a negative correlation between gene expression and positioning of genes at the nuclear envelope (NE) lined by nuclear lamina, but the exact relationship remains unclear, especially in light of the highly stochastic, transient nature of the gene association with the NE. RESULTS: In this paper, we ask whether there is a causal, systematic, genome-wide relationship between the expression levels of the groups of genes in topologically associating domains (TADs) of Drosophila nuclei and the probabilities of TADs to be found at the NE. To investigate the nature of this possible relationship, we combine a coarse-grained dynamic model of the entire Drosophila nucleus with genome-wide gene expression data; we analyze the TAD averaged transcription levels of genes against the probabilities of individual TADs to be in contact with the NE in the control and lamins-depleted nuclei. Our findings demonstrate that, within the statistical error margin, the stochastic positioning of Drosophila melanogaster TADs at the NE does not, by itself, systematically affect the mean level of gene expression in these TADs, while the expected negative correlation is confirmed. The correlation is weak and disappears completely for TADs not containing lamina-associated domains (LADs) or TADs containing LADs, considered separately. Verifiable hypotheses regarding the underlying mechanism for the presence of the correlation without causality are discussed. These include the possibility that the epigenetic marks and affinity to the NE of a TAD are determined by various non-mutually exclusive mechanisms and remain relatively stable during interphase. CONCLUSIONS: At the level of TADs, the probability of chromatin being in contact with the nuclear envelope has no systematic, causal effect on the transcription level in Drosophila. The conclusion is reached by combining model-derived time-evolution of TAD locations within the nucleus with their experimental gene expression levels.

Neuronal differentiation requires BRAT1 complex to remove REST from chromatin.

Repressor element-1 silencing transcription factor (REST) is required for the formation of mature neurons. REST dysregulation underlies a key mechanism of neurodegeneration associated with neurological disorders. However, the mechanisms leading to alterations of REST-mediated silencing of key neurogenesis genes are not known. Here, we show that BRCA1 Associated ATM Activator 1 (BRAT1), a gene linked to neurodegenerative diseases, is required for the activation of REST-responsive genes during neuronal differentiation. We find that INTS11 and INTS9 subunits of Integrator complex interact with BRAT1 as a distinct trimeric complex to activate critical neuronal genes during differentiation. BRAT1 depletion results in persistence of REST residence on critical neuronal genes disrupting the differentiation of NT2 cells into astrocytes and neuronal cells. We identified BRAT1 and INTS11 co-occupying the promoter region of these genes and pinpoint a role for BRAT1 in recruiting INTS11 to their promoters. Disease-causing mutations in BRAT1 diminish its association with INTS11/INTS9, linking the manifestation of disease phenotypes with a defect in transcriptional activation of key neuronal genes by BRAT1/INTS11/INTS9 complex. Finally, loss of Brat1 in mouse embryonic stem cells leads to a defect in neuronal differentiation assay. Importantly, while reconstitution with wild-type BRAT1 restores neuronal differentiation, the addition of a BRAT1 mutant is unable to associate with INTS11/INTS9 and fails to rescue the neuronal phenotype. Taken together, our study highlights the importance of BRAT1 association with INTS11 and INTS9 in the development of the nervous system.

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

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.

Effect of loops on the mean-square displacement of Rouse-model chromatin.

Chromatin polymer dynamics are commonly described using the classical Rouse model. The subsequent discovery, however, of intermediate-scale chromatin organization known as topologically associating domains (TADs) in experimental Hi-C contact maps for chromosomes across the tree of life, together with the success of loop extrusion factor (LEF) model in explaining TAD formation, motivates efforts to understand the effect of loops and loop extrusion on chromatin dynamics. This paper seeks to fulfill this need by combining LEF-model simulations with extended Rouse-model polymer simulations to investigate the dynamics of chromatin with loops and dynamic loop extrusion. We show that loops significantly suppress the averaged mean-square displacement (MSD) of a gene locus, consistent with recent experiments that track fluorescently labeled chromatin loci. We also find that loops reduce the MSD's stretching exponent from the classical Rouse-model value of 1/2 to a loop-density-dependent value in the 0.45-0.40 range. Remarkably, stretching exponent values in this range have also been observed in recent experiments [Weber et al., Phys. Rev. Lett. 104, 238102 (2010)0031-900710.1103/PhysRevLett.104.238102; Bailey et al., Mol. Biol. Cell 34, ar78 (2023)1059-152410.1091/mbc.E23-04-0119]. We also show that the dynamics of loop extrusion itself negligibly affects chromatin mobility. By studying static "rosette" loop configurations, we also demonstrate that chromatin MSDs and stretching exponents depend on the location of the locus in question relative to the position of the loops and on the local friction environment.

G-quadruplex folding in Xist RNA antagonizes PRC2 activity for stepwise regulation of X chromosome inactivation.

How Polycomb repressive complex 2 (PRC2) is regulated by RNA remains an unsolved problem. Although PRC2 binds G-tracts with the potential to form RNA G-quadruplexes (rG4s), whether rG4s fold extensively in vivo and whether PRC2 binds folded or unfolded rG4 are unknown. Using the X-inactivation model in mouse embryonic stem cells, here we identify multiple folded rG4s in Xist RNA and demonstrate that PRC2 preferentially binds folded rG4s. High-affinity rG4 binding inhibits PRC2's histone methyltransferase activity, and stabilizing rG4 in vivo antagonizes H3 at lysine 27 (H3K27me3) enrichment on the inactive X chromosome. Surprisingly, mutagenizing the rG4 does not affect PRC2 recruitment but promotes its release and catalytic activation on chromatin. H3K27me3 marks are misplaced, however, and gene silencing is compromised. Xist-PRC2 complexes become entrapped in the S1 chromosome compartment, precluding the required translocation into the S2 compartment. Thus, Xist rG4 folding controls PRC2 activity, H3K27me3 enrichment, and the stepwise regulation of chromosome-wide gene silencing.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome.

BACKGROUND: In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools. RESULTS: In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays. CONCLUSIONS: Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.