Nucleosome Positioning Research Articles

Investigating Nucleosome Stability via FRET

DNA's susceptibility to oxidative damage leads to the formation of the mutagenic Spiroiminodihydantoin (Sp) lesion. The Sp lesion contains a chiral center that disrupts the DNA double helix. When the Sp lesion is not repaired, it may lead to permanent mutations, apoptosis, cellular aging, or cancer. Previous work in our lab suggests that the Sp lesion can cause changes in nucleosome stability and DNA positioning. A nucleosome is the fundamental level of DNA packaging into chromosomes and forms when 146 base pair DNA duplex wraps around a histone octamer. This structure regulates access to DNA by repair proteins, transcription factors, and polymerases. With the use of Florescence Resonance Energy Transfer (FRET), we can examine how the Sp lesion affects nucleosome stability. We synthesized 146-mer duplexes with Cy3 and Cy5 fluorophores with Sp lesions at three potential locations. Then, we assembled nucleosome core particles with fluorescent donor-acceptor pairs on the DNA. FRET is a function of the proximity of the donor and acceptor fluorophores. Thus, when nucleosomes are gradually titrated with NaCl, DNA ends begin peeling and leads to a loss in acceptor fluorescence. With FRET, we are quantitatively measuring the correlation between nucleosome stability and location of the Sp lesion. We hope to advance our understanding of how the Sp lesion influences DNA compaction into nucleosomes.

Positive Torsional Stress on DNA Enhances Unwrapping of Nucleosomal DNA

Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. Some chromatin Remodelers are thought to change nucleosome positioning by adding torsional stress on DNA. To understand how such torsional stress affects the stability of the nucleosome, we examined the unwraping of the outer superhelical turn of the nucleosomal DNA under positive or negative torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates asymmetric unwrapping of the DNA while negative stress that undertwists DNA does not. Under positive stress, DNA was straightened and rigid, which facilitated a large-scale asymmetric unwrapping of the DNA without a large extension of two DNA. In the structure, one end of the DNA was disconnected from H3 and H2A/H2B while the other end of the DNA stably remained wrapped. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Our results suggest that the stability of individual nucleosomes in the chromatin can be changed by torsional stress, thereby DNA accessibility of regulatory proteins can be produced to maintain gene expression without large conformational changes even in heterochromatin.

Role of nucleosome positioning in 3D chromatin organization and loop formation

We present a physics-based polymer model that can investigate 3D organization of chromatin accounting for DNA elasticity, DNA-bending due to nucleosomes, and 1D organization of nucleosomes along DNA. We find that the packing density of chromatin oscillates between densities corresponding to highly folded and extended configurations as we change the nucleosome organization (length of linker DNA). We compute the looping probability of chromatin and show that the presence of nucleosomes increases the looping probability of the chain compared to that of a bare DNA. We also show that looping probability has a large variability depending on the nature of nucleosome organization and density of linker histones.

Construction of the gene regulatory network identifies MYC as a transcriptional regulator of SWI/SNF complex

Precise positioning of nucleosomes at the gene regulatory elements mediated by the SWI/SNF family of remodelling complex is important for the transcriptional regulation of genes. A wide set of genes are either positively or negatively regulated by SWI/SNF. In higher eukaryotes, around thirty genes were found to code for SWI/SNF subunits. The construction of a gene regulatory network of SWI/SNF subunits identifies MYC as a common regulator for many of the SWI/SNF subunit genes. A meta-analysis study was conducted to investigate the MYC dependent regulation of SWI/SNF remodelling complex. Subunit information and the promoter sequences of the subunit genes were used to find the canonical E-box motif and its variants. Detailed analysis of mouse and human ChIP-Seq at the SWI/SNF subunit loci indicates the presence of MYC binding peaks overlapping with E-boxes. The co-expression correlation and the differential expression analysis of wt vs. MYC perturbed MEFs indicate the MYC dependent regulation of some of the SWI/SNF subunits. The extension of the analysis was done on MYC proficient and MYC deficient embryonic fibroblast cell lines, TGR1 and HO15, and in one of the MYC amplified cancer types, Medulloblastoma. A transcriptional regulatory feedback loop between MYC and SWI/SNF could be a major factor contributing to the aggressiveness of MYC dependent cancers.

Analyzing RNA-DNA Triplex Formation in Chromatin.

A significant fraction of non-coding RNAs (ncRNAs) is associated with chromatin, shown to regulate gene expression and to organize nuclear architecture. Mechanisms of direct and indirect RNA-chromatin interactions have been described, including the sequence-specific formation of triple helix structures. Triplexes are formed by the sequence-specific binding of RNA to the bases located in the major groove of DNA. We recently showed that triplexes do exist in the context of cellular chromatin and that these structures are stabilized by the histone H3 tail of adjacent nucleosomes. The in vitro characterization of the specificity and binding affinity of triplex sequences next to nucleosomes are essential parameters to identify potential sites of RNA-chromatin interaction in vivo. Here we provide a detailed protocol to determine the influence of nucleosome positioning on triple helix formation. This assay allows the comparative quantification of triplex formation and specificity for triplex targeting sequences relative to the spatial nucleosome position.

Epigenetic biomarkers and preterm birth.

Preterm birth (PTB) is a major public health challenge, and novel, sensitive approaches to predict PTB are still evolving. Epigenomic markers are being explored as biomarkers of PTB because of their molecular stability compared to gene expression. This approach is also relatively new compared to gene-based diagnostics, which relies on mutations or single nucleotide polymorphisms. The fundamental principle of epigenome diagnostics is that epigenetic reprogramming in the target tissue (e.g. placental tissue) might be captured by more accessible surrogate tissue (e.g. blood) using biochemical epigenome assays on circulating DNA that incorporate methylation, histone modifications, nucleosome positioning, and/or chromatin accessibility. Epigenomic-based biomarkers may hold great potential for early identification of the majority of PTBs that are not associated with genetic variants or mutations. In this review, we discuss recent advances made in the development of epigenome assays focusing on its potential exploration for association and prediction of PTB. We also summarize population-level cohort studies conducted in the USA and globally that provide opportunities for genetic and epigenetic marker development for PTB. In addition, we summarize publicly available epigenome resources and published PTB studies. We particularly focus on ongoing genome-wide DNA methylation and epigenome-wide association studies. Finally, we review the limitations of current research, the importance of establishing a comprehensive biobank, and possible directions for future studies in identifying effective epigenome biomarkers to enhance health outcomes for pregnant women at risk of PTB and their infants.

Contact Sites of rDNA Clusters with FANK1 Gene Correspond to Repressed Chromatin

rDNA genes play an important role in epigenetic regulation and in differentiation of eukaryotic cells. Using the 4C (circular chromosome conformation capture) approach and model HEK293T cells, we analyzed the rDNA-contacting gene FANK1, using anchor located inside rDNA genes. At the 5' end of the FANK1 gene we detected frequent contacts with rDNA clusters. The contact sites coincide with the border where chromatin state changes and nucleosome positioning. The adjacent genes DHX32, BCCIP and UROS are located in the active chromatin and are transcribed, but do not contact with rDNA genes, while FANK1 gene is silenced, and is located in repressed chromatin. Heat shock treatment dramatically changes the pattern of rDNA contacts in the region and induces about 4-fold increase in activation of the FANK1 gene. We conclude that rDNA contacts may be involved in repression of the FANK1 gene.

Power to the protein: enhancing and combining activities using the Spy toolbox.

Proteins span an extraordinary range of shapes, sizes and functionalities. Therefore generic approaches are needed to overcome this diversity and stream-line protein analysis or application. Here we review SpyTag technology, now used in hundreds of publications or patents, and its potential for detecting and controlling protein behaviour. SpyTag forms a spontaneous and irreversible isopeptide bond upon binding its protein partner SpyCatcher, where both parts are genetically-encoded. New variants of this pair allow reaction at a rate approaching the diffusion limit, while reversible versions allow purification of SpyTagged proteins or tuned dynamic interaction inside cells. Anchoring of SpyTag-linked proteins has been established to diverse nanoparticles or surfaces, including gold, graphene and the air/water interface. SpyTag/SpyCatcher is mechanically stable, so is widely used for investigating protein folding and force sensitivity. A toolbox of scaffolds allows SpyTag-fusions to be assembled into defined multimers, from dimers to 180-mers, or unlimited 1D, 2D or 3D networks. Icosahedral multimers are being evaluated for vaccination against malaria, HIV and cancer. For enzymes, Spy technology has increased resilience, promoted substrate channelling, and assembled hydrogels for continuous flow biocatalysis. Combinatorial increase in functionality has been achieved through modular derivatisation of antibodies, light-emitting diodes or viral vectors. In living cells, SpyTag allowed imaging of protein trafficking, retargeting of CAR-T cell killing, investigation of heart contraction, and control of nucleosome position. The simple genetic encoding and rapid irreversible reaction provide diverse opportunities to enhance protein function. We describe limitations as well as future directions.

Tissue and Cell-Free DNA-Based Epigenomic Approaches for Cancer Detection.

Over 9 million people die of cancer each year worldwide, reflecting the unmet need for effective biomarkers for both cancer diagnosis and prognosis. Cancer diagnosis is complex because the majority of malignant tumors present with long periods of latency and lack of clinical presentation at early stages. During carcinogenesis, premalignant cells experience changes in their epigenetic landscapes, such as differential DNA methylation, histone modifications, nucleosome positioning, and higher orders of chromatin changes that confer growth advantage and contribute to determining the biologic phenotype of human cancers. Recent progress in microarray platforms and next-generation sequencing approaches has allowed the characterization of abnormal epigenetic patterns genome wide in a large number of cancer cases. The sizable amount of processed data also comes with challenges regarding data management and assessment for effective biomarker exploration to be further applied in prospective clinical trials. Epigenetics-based single or panel tests of genes are being explored for clinical management to fulfill unmet needs in oncology. The advance of these tests to the clinical routine will depend on rigorous, extensive, and independent validation in well-annotated cohort of patients and commercial development of clinical routine-friendly and adequate procedures. In this review we discuss the analytic validation of tissue and cell-free DNA-based epigenomic approaches for early cancer detection, diagnosis, and treatment monitoring and the clinical utility of candidate epigenetic alterations applied to colorectal, glioblastoma, breast, prostate, bladder, and lung cancer management.

Low-Phosphate Chromatin Dynamics Predict a Cell Wall Remodeling Network in Rice Shoots.

Phosphorus (P) is an essential plant macronutrient vital to fundamental metabolic processes. Plant-available P is low in most soils, making it a frequent limiter of growth. Declining P reserves for fertilizer production exacerbates this agricultural challenge. Plants modulate complex responses to fluctuating P levels via global transcriptional regulatory networks. Although chromatin structure plays a substantial role in controlling gene expression, the chromatin dynamics involved in regulating P homeostasis have not been determined. Here we define distinct chromatin states across the rice (Oryza sativa) genome by integrating multiple chromatin marks, including the H2A.Z histone variant, H3K4me3 modification, and nucleosome positioning. In response to P starvation, 40% of all protein-coding genes exhibit a transition from one chromatin state to another at their transcription start site. Several of these transitions are enriched in subsets of genes differentially expressed under P deficiency. The most prominent subset supports the presence of a coordinated signaling network that targets cell wall structure and is regulated in part via a decrease of H3K4me3 at transcription start sites. The P starvation-induced chromatin dynamics and correlated genes identified here will aid in enhancing P use efficiency in crop plants, benefitting global agriculture.

Q-Nuc: a bioinformatics pipeline for the quantitative analysis of nucleosomal profiles.

Nucleosomal profiling is an effective method to determine the positioning and occupancy of nucleosomes, which is essential to understand their roles in genomic processes. However, the positional randomness across the genome and its relationship with nucleosome occupancy remains poorly understood. Here we present a computational method that segments the profile into nucleosomal domains and quantifies their randomness and relative occupancy level. Applying this method to published data, we find on average ~ 3-fold differences in the degree of positional randomness between regions typically considered "well-ordered", as well as an unexpected predominance of only two types of domains of positional randomness in yeast cells. Further, we find that occupancy levels between domains actually differ maximally by ~ 2-3-fold in both cells, which has not been described before. We also developed a procedure by which one can estimate the sequencing depth that is required to identify nucleosomal positions even when regional positional randomness is high. Overall, we have developed a pipeline to quantitatively characterize domain-level features of nucleosome randomness and occupancy genome-wide, enabling the identification of otherwise unknown features in nucleosomal organization.

Genome-wide profiling of nucleosome position and chromatin accessibility in single cells using scMNase-seq.

Nucleosome organization is important for chromatin compaction and accessibility. Profiling nucleosome positioning genome-wide in single cells provides critical information to understand the cell-to-cell heterogeneity of chromatin states within a cell population. This protocol describes single-cell micrococcal nuclease sequencing (scMNase-seq), a method for detecting genome-wide nucleosome positioning and chromatin accessibility simultaneously from a small number of cells or single cells. To generate scMNase-seq libraries, single cells are isolated by FACS sorting, lysed and digested by MNase. DNA is purified, end-repaired and ligated to Y-shaped adaptors. Following PCR amplification with indexing primers, the subnucleosome-sized (fragments with a length of ≤80 bp) and mononucleosome-sized (fragments with a length between 140 and 180 bp) DNA fragments are recovered and sequenced on Illumina HiSeq platforms. On average, 0.5-1 million unique mapped reads are obtained for each single cell. The mononucleosome-sized DNA fragments precisely define genome-wide nucleosome positions in single cells, while the subnucleosome-sized DNA fragments provide information on chromatin accessibility. Library preparation of scMNase-seq takes only 2 d, requires only standard molecular biology techniques and does not require sophisticated laboratory equipment. Processing of high-throughput sequencing data requires basic bioinformatics skills and uses publicly available bioinformatics software.

A ligation-based single-stranded library preparation method to analyze cell-free DNA and synthetic oligos

BackgroundCell-free DNA (cfDNA), present in circulating blood plasma, contains information about prenatal health, organ transplant reception, and cancer presence and progression. Originally developed for the genomic analysis of highly degraded ancient DNA, single-stranded DNA (ssDNA) library preparation methods are gaining popularity in the field of cfDNA analysis due to their efficiency and ability to convert short, fragmented DNA into sequencing libraries without altering DNA ends. However, current ssDNA methods are costly and time-consuming.ResultsHere we present an efficient ligation-based single-stranded library preparation method that is engineered to produce complex libraries in under 2.5 h from as little as 1 nanogram of input DNA without alteration to the native ends of template molecules. Our method, called Single Reaction Single-stranded LibrarY or SRSLY, ligates uniquely designed Next-Generation Sequencing (NGS) adapters in a one-step combined phosphorylation/ligation reaction that foregoes end-polishing. Using synthetic DNA oligos and cfDNA, we demonstrate the efficiency and utility of this approach and compare with existing double-stranded and single-stranded approaches for library generation. Finally, we demonstrate that cfDNA NGS data generated from SRSLY can be used to analyze DNA fragmentation patterns to deduce nucleosome positioning and transcription factor binding.ConclusionsSRSLY is a versatile tool for converting short and fragmented DNA molecules, like cfDNA fragments, into sequencing libraries while retaining native lengths and ends.

Genomic methods in profiling DNA accessibility and factor localization

Recent advancements in next-generation sequencing technologies and accompanying reductions in cost have led to an explosion of techniques to examine DNA accessibility and protein localization on chromatin genome-wide. Generally, accessible regions of chromatin are permissive for factor binding and are therefore hotspots for regulation of gene expression; conversely, genomic regions that are highly occupied by histone proteins are not permissive for factor binding and are less likely to be active regulatory regions. Identifying regions of differential accessibility can be useful to uncover putative gene regulatory regions, such as enhancers, promoters, and insulators. In addition, DNA-binding proteins, such as transcription factors that preferentially bind certain DNA sequences and histone proteins that form the core of the nucleosome, play essential roles in all DNA-templated processes. Determining the genomic localization of chromatin-bound proteins is therefore essential in determining functional roles, sequence motifs important for factor binding, and regulatory networks controlling gene expression. In this review, we discuss techniques for determining DNA accessibility and nucleosome positioning (DNase-seq, FAIRE-seq, MNase-seq, and ATAC-seq) and techniques for detecting and functionally characterizing chromatin-bound proteins (ChIP-seq, DamID, and CUT&RUN). These methods have been optimized to varying degrees of resolution, specificity, and ease of use. Here, we outline some advantages and disadvantages of these techniques, their general protocols, and a brief discussion of their development. Together, these complimentary approaches have provided an unparalleled view of chromatin architecture and functional gene regulation.

TFIIIC Binding to Alu Elements Controls Gene Expression via Chromatin Looping and Histone Acetylation.

SummaryHow repetitive elements, epigenetic modifications, and architectural proteins interact ensuring proper genome expression remains poorly understood. Here, we report regulatory mechanisms unveiling a central role of Alu elements (AEs) and RNA polymerase III transcription factor C (TFIIIC) in structurally and functionally modulating the genome via chromatin looping and histone acetylation. Upon serum deprivation, a subset of AEs pre-marked by the activity-dependent neuroprotector homeobox Protein (ADNP) and located near cell-cycle genes recruits TFIIIC, which alters their chromatin accessibility by direct acetylation of histone H3 lysine-18 (H3K18). This facilitates the contacts of AEs with distant CTCF sites near promoter of other cell-cycle genes, which also become hyperacetylated at H3K18. These changes ensure basal transcription of cell-cycle genes and are critical for their re-activation upon serum re-exposure. Our study reveals how direct manipulation of the epigenetic state of AEs by a general transcription factor regulates 3D genome folding and expression.

Rap1-mediated nucleosome displacement can regulate gene expression in senescent cells without impacting the pace of senescence.

Cell senescence is accompanied, and in part mediated, by changes in chromatin, including histone losses, but underlying mechanisms are not well understood. We reported previously that during yeast cell senescence driven by telomere shortening, the telomeric protein Rap1 plays a major role in reprogramming gene expression by relocalizing hundreds of new target genes (called NRTS, for new Rap1 targets at senescence) to the promoters. This leads to two types of histone loss: Rap1 lowers histone level globally by repressing histone gene expression, and it also causes local nucleosome displacement at the promoters of upregulated NRTS. Here, we present evidence of direct binding between Rap1 and histone H3/H4 heterotetramers, and map amino acids involved in the interaction within the Rap1 SANT domain to amino acids 392–394 (SHY). Introduction of a point mutation within the native RAP1 locus that converts these residues to alanines (RAP1SHY), and thus disrupts Rap1‐H3/H4 interaction, does not interfere with Rap1 relocalization to NRTS at senescence, but prevents full nucleosome displacement and gene upregulation, indicating direct Rap1‐H3/H4 contacts are involved in nucleosome displacement. Consistent with this, the histone H3/H4 chaperone Asf1 is similarly unnecessary for Rap1 localization to NRTS but is required for full Rap1‐mediated nucleosome displacement and gene activation. Remarkably, RAP1SHY does not affect the pace of senescence‐related cell cycle arrest, indicating that some changes in gene expression at senescence are not coupled to this arrest.

Geometrical Heterogeneity Dominates Thermal Fluctuations in Facilitating Chromatin Contacts.

Within a living cell, the myriad of proteins that bind DNA introduce heterogeneously spaced kinks into an otherwise semiflexible DNA double helix. To investigate the effects of heterogeneous nucleosome binding on chromatin organization, we extend the wormlike chain model to include statistically spaced, rigid kinks. On timescales where nucleosome positions are fixed, we find that the probability of chromatin loop formation can vary by up to six orders of magnitude between two sets of nucleosome positions drawn from the same distribution. On longer timescales, we show that continuous rerandomization due to nucleosome turnover results in chromatin tracing out an effective WLC with a dramatically smaller Kuhn length than bare DNA. Together, these observations demonstrate that nucleosome spacing acts as the primary source of the structural heterogeneity that dominates local and global chromatin organization.

A distinct class of eukaryotic MT-A70 methyltransferases maintain symmetric DNA N6-adenine methylation at the ApT dinucleotides as an epigenetic mark associated with transcription.

Rediscovered as a potential eukaryotic epigenetic mark, DNA N6-adenine methylation (6mA) varies across species in abundance and its relationships with transcription. Here we characterize AMT1—representing a distinct MT-A70 family methyltransferase—in the ciliate Tetrahymena thermophila. AMT1 loss-of-function leads to severe defects in growth and development. Single Molecule, Real-Time (SMRT) sequencing reveals that AMT1 is required for the bulk of 6mA and all symmetric methylation at the ApT dinucleotides. The detection of hemi-methylated ApT sites suggests a semi-conservative mechanism for maintaining symmetric methylation. AMT1 affects expression of many genes; in particular, RAB46, encoding a Rab family GTPase involved in contractile vacuole function, is likely a direct target. The distribution of 6mA resembles H3K4 methylation and H2A.Z, two conserved epigenetic marks associated with RNA polymerase II transcription. Furthermore, strong 6mA and nucleosome positioning in wild-type cells is attenuated in ΔAMT1 cells. Our results support that AMT1-catalyzed 6mA is an integral part of the transcription-associated epigenetic landscape. AMT1 homologues are generally found in protists and basal fungi featuring ApT hyper-methylation associated with transcription, which are missing in animals, plants, and true fungi. This dichotomy of 6mA functions and the underlying molecular mechanisms may have implications in eukaryotic diversification.

Theory of Active Chromatin Remodeling.

Nucleosome positioning controls the accessible regions of chromatin and plays essential roles in DNA-templated processes. ATP driven remodeling enzymes are known to be crucial for its establishment invivo, but their nonequilibrium nature has hindered the development of a unified theoretical framework for nucleosome positioning. Using a perturbation theory, we show that the effect of these enzymes can be well approximated by effective equilibrium models with rescaled temperatures and interactions. Numerical simulations support the accuracy of the theory in predicting both kinetic and steady-state quantities, including the effective temperature and the radial distribution function, in biologically relevant regimes. The energy landscape view emerging from our study provides an intuitive understanding for the impact of remodeling enzymes in either reinforcing or overwriting intrinsic signals for nucleosome positioning, and may help improve the accuracy of computational models for its prediction insilico.

The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription.

Transcription start sites (TSS) in eukaryotes are characterized by a nucleosome-depleted region (NDR), which appears to be flanked upstream and downstream by strongly positioned nucleosomes incorporating the histone variant H2A.Z. H2A.Z associates with both active and repressed TSS and is important for priming genes for rapid transcriptional activation. However, the determinants of H2A.Z occupancy at specific nucleosomes and its relationship to transcription initiation remain unclear. To further elucidate the specificity of H2A.Z, we determined its genomic localization at single nucleosome resolution, as well as the localization of its chromatin remodelers Swr1 and Ino80. By analyzing H2A.Z occupancy in conjunction with RNA expression data that captures promoter-derived antisense initiation, we find that H2A.Z’s bimodal incorporation on either side of the NDR is not a general feature of TSS, but is specifically a marker for bidirectional transcription, such that the upstream flanking −1 H2A.Z-containing nucleosome is more appropriately considered as a +1 H2A.Z nucleosome for antisense transcription. The localization of H2A.Z almost exclusively at the +1 nucleosome suggests that a transcription-initiation dependent process could contribute to its specific incorporation.