Nucleosome positioning and chromatin organization

Pan-cancer analysis of m1A writer gene RRP8: implications for immune infiltration and prognosis in human cancers

BACKGROUND

Ribosomal RNA Processing 8 (RRP8) is a gene associated with RNA modification and has been implicated in the development of several types of tumors in recent research. Nevertheless, the biological importance of RRP8 in pan-cancer has not yet been thoroughly and comprehensively investigated.

METHODS

In this study, we conducted an analysis of various public databases to investigate the biological functions of RRP8. Our analysis included examining its correlation with pan-cancer prognosis, heterogeneity, stemness, immune checkpoint genes, and immune cell infiltration. Furthermore, we utilized the GDSC and CTRP databases to assess the sensitivity of RRP8 to small molecule drugs.

RESULTS

Our findings indicate that RRP8 exhibits differential expression between tumor and normal samples, particularly impacting the prognosis of various cancers such as Adrenocortical carcinoma (ACC) and Kidney Chromophobe (KICH). The expression of RRP8 is intricately linked to tumor heterogeneity and stemness markers. Additionally, RRP8 shows a positive correlation with the presence of tumor-infiltrating cells, with TP53 being the predominant mutated gene in these malignancies.

CONCLUSION

Our findings suggest that RRP8 may serve as a potential prognostic marker and therapeutic target in a variety of cancer types.

Unraveling the mechanisms of PAMless DNA interrogation by SpRY-Cas9

CRISPR-Cas9 is a powerful tool for genome editing, but the strict requirement for an NGG protospacer-adjacent motif (PAM) sequence immediately next to the DNA target limits the number of editable genes. Recently developed Cas9 variants have been engineered with relaxed PAM requirements, including SpG-Cas9 (SpG) and the nearly PAM-less SpRY-Cas9 (SpRY). However, the molecular mechanisms of how SpRY recognizes all potential PAM sequences remains unclear. Here, we combine structural and biochemical approaches to determine how SpRY interrogates DNA and recognizes target sites. Divergent PAM sequences can be accommodated through conformational flexibility within the PAM-interacting region, which facilitates tight binding to off-target DNA sequences. Nuclease activation occurs ~1000-fold slower than for Streptococcus pyogenes Cas9, enabling us to directly visualize multiple on-pathway intermediate states. Experiments with SpG position it as an intermediate enzyme between Cas9 and SpRY. Our findings shed light on the molecular mechanisms of PAMless genome editing.

Decoding Chromatin Ubiquitylation: A Chemical Biology Perspective

Since Strahl and Allis proposed the "language of covalent histone modifications", a host of experimental studies have shed light on the different facets of chromatin regulation by epigenetic mechanisms. Initially proposed as a concept for controlling gene transcription, the regulation of deposition and removal of histone post-translational modifications (PTMs), such as acetylation, methylation, and phosphorylation, have been implicated in many chromatin regulation pathways. However, large PTMs such as ubiquitylation challenge research on many levels due to their chemical complexity. In recent years, chemical tools have been developed to generate chromatin in defined ubiquitylation states in vitro. Chemical biology approaches are now used to link specific histone ubiquitylation marks with downstream chromatin regulation events on the molecular level. Here, we want to highlight how chemical biology approaches have empowered the mechanistic study of chromatin ubiquitylation in the context of gene regulation and DNA repair with attention to future challenges.

Antitumor efficacy of a sequence-specific DNA-targeted γPNA-based c-Myc inhibitor

Targeting oncogenes at the genomic DNA level can open new avenues for precision medicine. Significant efforts are ongoing to target oncogenes using RNA-targeted and protein-targeted platforms, but no progress has been made to target genomic DNA for cancer therapy. Here, we introduce a gamma peptide nucleic acid (γPNA)-based genomic DNA-targeted platform to silence oncogenes in vivo. γPNAs efficiently invade the mixed sequences of genomic DNA with high affinity and specificity. As a proof of concept, we establish that γPNA can inhibit c-Myc transcription in multiple cell lines. We evaluate the in vivo efficacy and safety of genomic DNA targeting in three pre-clinical models. We also establish that anti-transcription γPNA in combination with histone deacetylase inhibitors and chemotherapeutic drugs results in robust antitumor activity in cell-line- and patient-derived xenografts. Overall, this strategy offers a unique therapeutic platform to target genomic DNA to inhibit oncogenes for cancer therapy.

Diffusion controls local versus dispersed inheritance of histones during replication and shapes epigenomic architecture

The dynamics of inheritance of histones and their associated modifications across cell divisions can have major consequences on maintenance of the cellular epigenomic state. Recent experiments contradict the long-held notion that histone inheritance during replication is always local, suggesting that active and repressed regions of the genome exhibit fundamentally different histone dynamics independent of transcription-coupled turnover. Here we develop a stochastic model of histone dynamics at the replication fork and demonstrate that differential diffusivity of histones in active versus repressed chromatin is sufficient to quantitatively explain these recent experiments. Further, we use the model to predict patterns in histone mark similarity between pairs of genomic loci that should be developed as a result of diffusion, but cannot originate from either PRC2 mediated mark spreading or transcriptional processes. Interestingly, using a combination of CHIP-seq, replication timing and Hi-C datasets we demonstrate that all the computationally predicted patterns are consistently observed for both active and repressive histone marks in two different cell lines. While direct evidence for histone diffusion remains controversial, our results suggest that dislodged histones in euchromatin and facultative heterochromatin may exhibit some level of diffusion within "Diffusion-Accessible-Domains" (DADs), leading to redistribution of epigenetic marks within and across chromosomes. Preservation of the epigenomic state across cell divisions therefore might be achieved not by passing on strict positional information of histone marks, but by maintaining the marks in somewhat larger DADs of the genome.

Multiple Roles of dXNP and dADD1-Drosophila Orthologs of ATRX Chromatin Remodeler

The Drosophila melanogaster dADD1 and dXNP proteins are orthologues of the ADD and SNF2 domains of the vertebrate ATRX (Alpha-Thalassemia with mental Retardation X-related) protein. ATRX plays a role in general molecular processes, such as regulating chromatin status and gene expression, while dADD1 and dXNP have similar functions in the Drosophila genome. Both ATRX and dADD1/dXNP interact with various protein partners and participate in various regulatory complexes. Disruption of ATRX expression in humans leads to the development of α-thalassemia and cancer, especially glioma. However, the mechanisms that allow ATRX to regulate various cellular processes are poorly understood. Studying the functioning of dADD1/dXNP in the Drosophila model may contribute to understanding the mechanisms underlying the multifunctional action of ATRX and its connection with various cellular processes. This review provides a brief overview of the currently available information in mammals and Drosophila regarding the roles of ATRX, dXNP, and dADD1. It discusses possible mechanisms of action of complexes involving these proteins.

Base pair compositional variability influences DNA structural stability and tunes hydration thermodynamics and dynamics

DNA deformability and differential hydration are crucial determinants of biological processes ranging from genetic material packaging to gene expression; their associative details, however, remain inadequately understood. Herein, we report investigations of the dynamic and thermodynamic responses of the local hydration of a variety of base pair sequences. Leveraging in silico sampling and our in-house analyses, we first report the local conformational propensity of sequences that are either predisposed toward the canonical A- or B-conformations or are restrained to potential transitory pathways. It is observed that the transition from the unrestrained A-form to the B-form leads to lengthwise structural deformation. The insertion of intermittent -(CG)- base pairs in otherwise homogeneous -(AT)- sequences bears dynamical consequences for the vicinal hydration layer. Calculation of the excess (pair) entropy suggests substantially higher values of hydration water surrounding A conformations over the B- conformations. Applying the Rosenfeld approximation, we project that the diffusivity of water molecules proximal to canonical B conformation is least for the minor groove of the canonical B-conformation. We determine that structure, composition, and conformation specific groove dimension together influence the local hydration characteristics and, therefore, are expected to be important determinants of biological processes.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining

Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.

How Do Nucleosome Dynamics Regulate Protein Search on DNA?

Nearly three-fourths of all eukaryotic DNA is occupied by nucleosomes, protein-DNA complexes comprising octameric histone core proteins and ∼150 base pairs of DNA. In addition to acting as a DNA compaction vehicle, the dynamics of nucleosomes regulate the DNA site accessibility for the nonhistone proteins, thereby controlling regulatory processes involved in determining the cell identity and cell fate. Here, we propose an analytical framework to analyze the role of nucleosome dynamics on the target search process of transcription factors through a simple discrete-state stochastic description of the search process. By considering the experimentally determined kinetic rates associated with protein and nucleosome dynamics as the only inputs, we estimate the target search time of a protein via first-passage probability calculations separately during nucleosome breathing and sliding dynamics. Although both the nucleosome dynamics permit transient access to the DNA sites that are otherwise occluded by the histone proteins, our result suggests substantial differences between the protein search mechanism on a nucleosome performing breathing and sliding dynamics. Furthermore, we identify the molecular factors that influence the search efficiency and demonstrate how these factors together portray a highly dynamic landscape of gene regulation. Our analytical results are validated using extensive Monte Carlo simulations.

Disruption of polyhomeotic polymerization decreases nucleosome occupancy and alters genome accessibility

Chromatin attains its three-dimensional (3D) conformation by establishing contacts between different noncontiguous regions. Sterile Alpha Motif (SAM)-mediated polymerization of the polyhomeotic (PH) protein regulates subnuclear clustering of Polycomb Repressive Complex 1 (PRC1) and chromatin topology. The mutations that perturb the ability of the PH to polymerize, disrupt long-range chromatin contacts, alter Hox gene expression, and lead to developmental defects. To understand the underlying mechanism, we combined the experiments and theory to investigate the effect of this SAM domain mutation on nucleosome occupancy and accessibility on a genome wide scale. Our data show that disruption of PH polymerization because of SAM domain mutation decreases nucleosome occupancy and alters accessibility. Polymer simulations investigating the interplay between distant chromatin contacts and nucleosome occupancy, both of which are regulated by PH polymerization, suggest that nucleosome density increases when contacts between different regions of chromatin are established. Taken together, it appears that SAM domain-mediated PH polymerization biomechanically regulates the organization of chromatin at multiple scales from nucleosomes to chromosomes and we suggest that higher order organization can have a top-down causation effect on nucleosome occupancy.

Epigenetic regulation in metabolic diseases: mechanisms and advances in clinical study

Epigenetics regulates gene expression and has been confirmed to play a critical role in a variety of metabolic diseases, such as diabetes, obesity, non-alcoholic fatty liver disease (NAFLD), osteoporosis, gout, hyperthyroidism, hypothyroidism and others. The term 'epigenetics' was firstly proposed in 1942 and with the development of technologies, the exploration of epigenetics has made great progresses. There are four main epigenetic mechanisms, including DNA methylation, histone modification, chromatin remodelling, and noncoding RNA (ncRNA), which exert different effects on metabolic diseases. Genetic and non-genetic factors, including ageing, diet, and exercise, interact with epigenetics and jointly affect the formation of a phenotype. Understanding epigenetics could be applied to diagnosing and treating metabolic diseases in the clinic, including epigenetic biomarkers, epigenetic drugs, and epigenetic editing. In this review, we introduce the brief history of epigenetics as well as the milestone events since the proposal of the term 'epigenetics'. Moreover, we summarise the research methods of epigenetics and introduce four main general mechanisms of epigenetic modulation. Furthermore, we summarise epigenetic mechanisms in metabolic diseases and introduce the interaction between epigenetics and genetic or non-genetic factors. Finally, we introduce the clinical trials and applications of epigenetics in metabolic diseases.

Current understanding of the cGAS-STING signaling pathway: Structure, regulatory mechanisms, and related diseases

The innate immune system protects the host from external pathogens and internal damage in various ways. The cGAS-STING signaling pathway, comprised of cyclic GMP-AMP synthase (cGAS), stimulator of interferon genes (STING), and downstream signaling adaptors, plays an essential role in protective immune defense against microbial DNA and internal damaged-associated DNA and is responsible for various immune-related diseases. After binding with DNA, cytosolic cGAS undergoes conformational change and DNA-linked liquid-liquid phase separation to produce 2'3'-cGAMP for the activation of endoplasmic reticulum (ER)-localized STING. However, further studies revealed that cGAS is predominantly expressed in the nucleus and strictly tethered to chromatin to prevent binding with nuclear DNA, and functions differently from cytosolic-localized cGAS. Detailed delineation of this pathway, including its structure, signaling, and regulatory mechanisms, is of great significance to fully understand the diversity of cGAS-STING activation and signaling and will be of benefit for the treatment of inflammatory diseases and cancer. Here, we review recent progress on the above-mentioned perspectives of the cGAS-STING signaling pathway and discuss new avenues for further study.

Understanding the Target Search by Multiple Transcription Factors on Nucleosomal DNA

The association of multiple Transcription Factors (TFs) in the cis-regulatory region is imperative for developmental changes in eukaryotes. The underlying process is exceedingly complex, and it is not at all clear what orchestrates the overall search process by multiple TFs. In this study, by developing a theoretical model based on a discrete-state stochastic approach, we investigated the target search mechanism of multiple TFs on nucleosomal DNA. Experimental kinetic rate constants of different TFs are taken as input to estimate the Mean-First-Passage time to recognize the binding motifs by two TFs on a dynamic nucleosome model. The theory systematically analyzes when the TFs search their binding motifs hierarchically and when simultaneously by proceeding via the formation of a protein-protein complex. Our results, validated by extensive Monte Carlo simulations, elucidate the molecular basis of the complex target search phenomenon of multiple TFs on nucleosomal DNA.

Mapping Nucleosome Location Using FS-Seq

The organization of nucleosomes in eukaryotic chromatin is thought to play a critical role in the regulation of the biological function of the chromatin. Because of this potential role in regulation, a number of techniques have been developed, which combine chromatin fragmentation around nucleosomes with next-generation sequencing to map the location of nucleosomes in chromatin. In this section, a procedure using a kit from New England Biolabs (NEB NEXT Ultra II FS DNA library prep Kit) to fragment chromatin in preparation for next-generation sequencing is described and compared to other available procedures for mapping nucleosome location.

Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains

Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.

Differences in nanoscale organization of regulatory active and inactive human chromatin

Methodological advances in conformation capture techniques have fundamentally changed our understanding of chromatin architecture. However, the nanoscale organization of chromatin and its cell-to-cell variance are less studied. Analyzing genome-wide data from 733 human cell and tissue samples, we identified 2 prototypical regions that exhibit high or absent hypersensitivity to deoxyribonuclease I, respectively. These regulatory active or inactive regions were examined in the lymphoblast cell line K562 by using high-throughput super-resolution microscopy. In both regions, we systematically measured the physical distance of 2 fluorescence in situ hybridization spots spaced by only 5 kb of DNA. Unexpectedly, the resulting distance distributions range from very compact to almost elongated configurations of more than 200-nm length for both the active and inactive regions. Monte Carlo simulations of a coarse-grained model of these chromatin regions based on published data of nucleosome occupancy in K562 cells were performed to understand the underlying mechanisms. There was no parameter set for the simulation model that can explain the microscopically measured distance distributions. Obviously, the chromatin state given by the strength of internucleosomal interaction, nucleosome occupancy, or amount of histone H1 differs from cell to cell, which results in the observed broad distance distributions. This large variability was not expected, especially in inactive regions. The results for the mechanisms for different distance distributions on this scale are important for understanding the contacts that mediate gene regulation. Microscopic measurements show that the inactive region investigated here is expected to be embedded in a more compact chromatin environment. The simulation results of this region require an increase in the strength of internucleosomal interactions. It may be speculated that the higher density of chromatin is caused by the increased internucleosomal interaction strength.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

NucPosDB: a database of nucleosome positioning in vivo and nucleosomics of cell-free DNA

Nucleosome positioning is involved in many gene regulatory processes happening in the cell, and it may change as cells differentiate or respond to the changing microenvironment in a healthy or diseased organism. One important implication of nucleosome positioning in clinical epigenetics is its use in the “nucleosomics” analysis of cell-free DNA (cfDNA) for the purpose of patient diagnostics in liquid biopsies. The rationale for this is that the apoptotic nucleases that digest chromatin of the dying cells mostly cut DNA between nucleosomes. Thus, the short pieces of DNA in body fluids reflect the positions of nucleosomes in the cells of origin. Here, we report a systematic nucleosomics database — NucPosDB — curating published nucleosome positioning datasets in vivo as well as datasets of sequenced cell-free DNA (cfDNA) that reflect nucleosome positioning in situ in the cells of origin. Users can select subsets of the database by a number of criteria and then obtain raw or processed data. NucPosDB also reports the originally determined regions with stable nucleosome occupancy across several individuals with a given condition. An additional section provides a catalogue of computational tools for the analysis of nucleosome positioning or cfDNA experiments and theoretical algorithms for the prediction of nucleosome positioning preferences from DNA sequence. We provide an overview of the field, describe the structure of the database in this context, and demonstrate data variability using examples of different medical conditions. NucPosDB is useful both for the analysis of fundamental gene regulation processes and the training of computational models for patient diagnostics based on cfDNA. The database currently curates ~ 400 publications on nucleosome positioning in cell lines and in situ as well as cfDNA from > 10,000 patients and healthy volunteers. For open-access cfDNA datasets as well as key MNase-seq datasets in human cells, NucPosDB allows downloading processed mapped data in addition to the regions with stable nucleosome occupancy. NucPosDB is available at https://generegulation.org/nucposdb/.

Nucleosome-induced homology recognition in chromatin

One of the least understood properties of chromatin is the ability of its similar regions to recognize each other through weak interactions. Theories based on electrostatic interactions between helical macromolecules suggest that the ability to recognize sequence homology is an innate property of the non-ideal helical structure of DNA. However, this theory does not account for the nucleosomal packing of DNA. Can homologous DNA sequences recognize each other while wrapped up in the nucleosomes? Can structural homology arise at the level of nucleosome arrays? Here, we present a theoretical model for the recognition potential well between chromatin fibres sliding against each other. This well is different from the one predicted for bare DNA; the minima in energy do not correspond to literal juxtaposition, but are shifted by approximately half the nucleosome repeat length. The presence of this potential well suggests that nucleosome positioning may induce mutual sequence recognition between chromatin fibres and facilitate the formation of chromatin nanodomains. This has implications for nucleosome arrays enclosed between CTCF-cohesin boundaries, which may form stiffer stem-like structures instead of flexible entropically favourable loops. We also consider switches between chromatin states, e.g. through acetylation/deacetylation of histones, and discuss nucleosome-induced recognition as a precursory stage of genetic recombination.