High-throughput mapping of the chromatin structure of human promoters

Translation of Epigenetics in Cell-Free DNA Liquid Biopsy Technology and Precision Oncology

Technological advancements in cell-free DNA (cfDNA) liquid biopsy have triggered exponential growth in numerous clinical applications. While cfDNA-based liquid biopsy has made significant strides in personalizing cancer treatment, the exploration and translation of epigenetics in liquid biopsy to clinical practice is still nascent. This comprehensive review seeks to provide a broad yet in-depth narrative of the present status of epigenetics in cfDNA liquid biopsy and its associated challenges. It highlights the potential of epigenetics in cfDNA liquid biopsy technologies with the hopes of enhancing its clinical translation. The momentum of cfDNA liquid biopsy technologies in recent years has propelled epigenetics to the forefront of molecular biology. We have only begun to reveal the true potential of epigenetics in both our understanding of disease and leveraging epigenetics in the diagnostic and therapeutic domains. Recent clinical applications of epigenetics-based cfDNA liquid biopsy revolve around DNA methylation in screening and early cancer detection, leading to the development of multi-cancer early detection tests and the capability to pinpoint tissues of origin. The clinical application of epigenetics in cfDNA liquid biopsy in minimal residual disease, monitoring, and surveillance are at their initial stages. A notable advancement in fragmentation patterns analysis has created a new avenue for epigenetic biomarkers. However, the widespread application of cfDNA liquid biopsy has many challenges, including biomarker sensitivity, specificity, logistics including infrastructure and personnel, data processing, handling, results interpretation, accessibility, and cost effectiveness. Exploring and translating epigenetics in cfDNA liquid biopsy technology can transform our understanding and perception of cancer prevention and management. cfDNA liquid biopsy has great potential in precision oncology to revolutionize conventional ways of early cancer detection, monitoring residual disease, treatment response, surveillance, and drug development. Adapting the implementation of liquid biopsy workflow to the local policy worldwide and developing point-of-care testing holds great potential to overcome global cancer disparity and improve cancer outcomes.

Effects of blood-processing protocols on cell-free DNA fragmentomics in plasma: Comparisons of one- and two-step centrifugations

BACKGROUND

Cell-free DNA (cfDNA) fragmentomic characteristics are promising analytes with abundant physiological signals for non-invasive disease diagnosis and monitoring. Previous studies on plasma cfDNA fragmentomics commonly employed a two-step centrifugation process for removing cell debris, involving a low-speed centrifugation followed by a high-speed centrifugation. However, the effects of centrifugation conditions on the analysis of cfDNA fragmentome remain uncertain.

METHODS

We collected blood samples from 10 healthy individuals and divided each sample into two aliquots for plasma preparation with one- and two-step centrifugation processes. We performed whole genome sequencing (WGS) of the plasma cfDNA in the two groups and comprehensively compared the cfDNA fragmentomic features. Additionally, we reanalyzed the fragmentomic features of cfDNA from 16 healthy individuals and 16 COVID-19 patients, processed through one- and two-step centrifugation in our previous study, to investigate the impact of centrifugation on disease signals.

RESULTS

Our results showed that there were no significant differences observed in the characteristics of nuclear cfDNA, including size, motif diversity score (MDS) of end motifs, and genome distribution, between plasma samples treated with one- and two-step centrifugation. The cfDNA size shortening in COVID-19 patients was observed in plasma samples with one- and two-step centrifugation methods. However, we observed a significantly higher relative abundance and longer size of cell-free mitochondrial DNA (mtDNA) in the one-step samples compared to the two-step samples. This difference in mtDNA caused by the one- and two-step centrifugation methods surpasses the pathological difference between COVID-19 patients and healthy individuals.

CONCLUSIONS

Our findings indicate that one-step low-speed centrifugation is a simple and potentially suitable method for analyzing nuclear cfDNA fragmentation characteristics. These results offer valuable guidance for cfDNA research in various clinical scenarios.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Endogenous aldehyde-induced DNA-protein crosslinks are resolved by transcription-coupled repair

DNA-protein crosslinks (DPCs) induced by aldehydes interfere with replication and transcription. Hereditary deficiencies in DPC repair and aldehyde clearance processes cause progeria, including Ruijs-Aalfs syndrome (RJALS) and AMeD syndrome (AMeDS) in humans. Although the elimination of DPC during replication has been well established, how cells overcome DPC lesions in transcription remains elusive. Here we show that endogenous aldehyde-induced DPC roadblocks are efficiently resolved by transcription-coupled repair (TCR). We develop a high-throughput sequencing technique to measure the genome-wide distribution of DPCs (DPC-seq). Using proteomics and DPC-seq, we demonstrate that the conventional TCR complex as well as VCP/p97 and the proteasome are required for the removal of formaldehyde-induced DPCs. TFIIS-dependent cleavage of RNAPII transcripts protects against transcription obstacles. Finally, a mouse model lacking both aldehyde clearance and TCR confirms endogenous DPC accumulation in actively transcribed regions. Collectively, our data provide evidence that transcription-coupled DPC repair (TC-DPCR) as well as aldehyde clearance are crucial for protecting against metabolic genotoxin, thus explaining the molecular pathogenesis of AMeDS and other disorders associated with defects in TCR, such as Cockayne syndrome.

Cold-induced deposition of bivalent H3K4me3-H3K27me3 modification and nucleosome depletion in Arabidopsis

Epigenetic regulation of gene expression plays a crucial role in plant development and environmental adaptation. The H3K4me3 and H3K27me3 have not only been discovered in the regulation of gene expression in multiple biological processes but also in responses to abiotic stresses in plants. However, evidence for the presence of both H3K4me3 and H3K27me3 on the same nucleosome is sporadic. Cold-induced deposition of bivalent H3K4me3-H3K27me3 modifications and nucleosome depletion over a considerable number of active genes is documented in potato tubers and provides clues on an additional role of the bivalent modifications. Limited by the available information of genes encoding PcG/TrxG proteins as well as their corresponding mutants in potatoes, the molecular mechanism underlying the cold-induced deposition of the bivalent mark remains elusive. In this study, we found a similar deposition of the bivalent H3K4me3-H3K27me3 mark over 2129 active genes in cold-treated Arabidopsis Col-0 seedlings. The expression levels of the bivalent mark-associated genes tend to be independent of bivalent modification levels. However, these genes were associated with greater chromatin accessibility, presumably to provide a distinct chromatin environment for gene expression. In mutants clf28 and lhp1, failure to deposit H3K27me3 in active genes upon cold treatment implies that the CLF is potentially involved in cold-induced deposition of H3K27me3, with assistance from LHP1. Failure to deposit H3K4me3 during cold treatment in atx1-2 suggests a regulatory role of ATX1 in the deposition of H3K4me3. In addition, we observed a cold-induced global reduction in nucleosome occupancy, which is potentially mediated by LHP1 in an H3K27me3-dependent manner.

ATAC-Seq Analysis of Accessible Chromatin: From Experimental Steps to Data Analysis

Accessible chromatin often represents gene regulatory elements, including promoters and enhancers, essential for gene expression. Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is one of the most popular techniques to investigate chromatin accessibility across the genome. Here we describe, step by step, a series of optimized experimental methods and bioinformatics pipelines for ATAC-seq analysis. As an example, we present an analysis of murine spermatogenic cells: a method to isolate germ cells, a reaction step using Tn5 transposase to insert sequencing adapters into accessible DNA, a library preparation method for high-throughput sequencing, and bioinformatics analysis of sequencing data. Overall, we introduce a framework of ATAC-seq analysis that can be applied to any cell population to identify cell-type-specific gene regulatory elements and their cis-regulatory networks.

Characterization of chromatin accessibility patterns in different mouse cell types using machine learning methods at single-cell resolution

Chromatin accessibility is a generic property of the eukaryotic genome, which refers to the degree of physical compaction of chromatin. Recent studies have shown that chromatin accessibility is cell type dependent, indicating chromatin heterogeneity across cell lines and tissues. The identification of markers used to distinguish cell types at the chromosome level is important to understand cell function and classify cell types. In the present study, we investigated transcriptionally active chromosome segments identified by sci-ATAC-seq at single-cell resolution, including 69,015 cells belonging to 77 different cell types. Each cell was represented by existence status on 20,783 genes that were obtained from 436,206 active chromosome segments. The gene features were deeply analyzed by Boruta, resulting in 3897 genes, which were ranked in a list by Monte Carlo feature selection. Such list was further analyzed by incremental feature selection (IFS) method, yielding essential genes, classification rules and an efficient random forest (RF) classifier. To improve the performance of the optimal RF classifier, its features were further processed by autoencoder, light gradient boosting machine and IFS method. The final RF classifier with MCC of 0.838 was constructed. Some marker genes such as H2-Dmb2, which are specifically expressed in antigen-presenting cells (e.g., dendritic cells or macrophages), and Tenm2, which are specifically expressed in T cells, were identified in this study. Our analysis revealed numerous potential epigenetic modification patterns that are unique to particular cell types, thereby advancing knowledge of the critical functions of chromatin accessibility in cell processes.

Beyond Sequence: Internucleosomal Interactions Dominate Array Assembly

The organization of the nucleosome array is a critical component of the chromatin assembly into higher order structure as well as its function. Here, we investigated the contributions of the DNA sequence and internucleosomal interactions on the organization of the nucleosomal arrays in compact structures using atomic force microscopy. We assembled nucleosomes on DNA substrates allowing for the formation of tetranucleosomes. We found that nucleosomes are capable of close positioning with no discernible space between them, even in the case of assembled dinucleosomes. This morphology of the array is in contrast with that observed for arrays assembled with repeats of the nucleosome positioning motifs separated by uniform spacers. Simulated assembly of tetranucleosomes by random placement along the substrates revealed that nucleosome array compaction is promoted by the interaction of the nucleosomes. We developed a theoretical model to account for the role of DNA sequence and internucleosomal interactions in the formation of the nucleosome structures. These findings suggest that, in the chromatin assembly, the affinity of the nucleosomes to the DNA sequence and the strengths of the internucleosomal interactions are the two major factors defining the compactness of the chromatin.

Systematic analysis of the effects of genetic variants on chromatin accessibility to decipher functional variants in non-coding regions

Genome-wide association study (GWAS) has identified thousands of single nucleotide polymorphisms (SNPs) associated with complex diseases and traits. However, deciphering the functions of these SNPs still faces challenges. Recent studies have shown that SNPs could alter chromatin accessibility and result in differences in tumor susceptibility between individuals. Therefore, systematically analyzing the effects of SNPs on chromatin accessibility could help decipher the functions of SNPs, especially those in non-coding regions. Using data from The Cancer Genome Atlas (TCGA), chromatin accessibility quantitative trait locus (caQTL) analysis was conducted to estimate the associations between genetic variants and chromatin accessibility. We analyzed caQTLs in 23 human cancer types and identified 9,478 caQTLs in breast carcinoma (BRCA). In BRCA, these caQTLs tend to alter the binding affinity of transcription factors, and open chromatin regions regulated by these caQTLs are enriched in regulatory elements. By integrating with eQTL data, we identified 141 caQTLs showing a strong signal for colocalization with eQTLs. We also identified 173 caQTLs in genome-wide association studies (GWAS) loci and inferred several possible target genes of these caQTLs. By performing survival analysis, we found that ~10% caQTLs potentially influence the prognosis of patients. To facilitate access to relevant data, we developed a user-friendly data portal, BCaQTL (http://gong_lab.hzau.edu.cn/caqtl_database), for data searching and downloading. Our work may facilitate fine-map regulatory mechanisms underlying risk loci of cancer and discover the biomarkers or therapeutic targets for cancer prognosis. The BCaQTL database will be an important resource for genetic and epigenetic studies.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

Background

A comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.

Results

Here we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.

Conclusions

Our data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00767-3.

Nutrient sensitive protein O-GlcNAcylation modulates the transcriptome through epigenetic mechanisms during embryonic neurogenesis

Protein O-GlcNAcylation is a dynamic, nutrient-sensitive mono-glycosylation deposited on numerous nucleo-cytoplasmic and mitochondrial proteins, including transcription factors, epigenetic regulators, and histones. However, the role of protein O-GlcNAcylation on epigenome regulation in response to nutrient perturbations during development is not well understood. Herein we recapitulated early human embryonic neurogenesis in cell culture and found that pharmacological up-regulation of O-GlcNAc levels during human embryonic stem cells' neuronal differentiation leads to up-regulation of key neurogenic transcription factor genes. This transcriptional de-repression is associated with reduced H3K27me3 and increased H3K4me3 levels on the promoters of these genes, perturbing promoter bivalency possibly through increased EZH2-Thr311 phosphorylation. Elevated O-GlcNAc levels also lead to increased Pol II-Ser5 phosphorylation and affect H2BS112O-GlcNAc and H2BK120Ub1 on promoters. Using an in vivo rat model of maternal hyperglycemia, we show similarly elevated O-GlcNAc levels and epigenetic dysregulations in the developing embryo brains because of hyperglycemia, whereas pharmacological inhibition of O-GlcNAc transferase (OGT) restored these molecular changes. Together, our results demonstrate O-GlcNAc mediated sensitivity of chromatin to nutrient status, and indicate how metabolic perturbations could affect gene expression during neurodevelopment.

DeNOPA: decoding nucleosome positions sensitively with sparse ATAC-seq data

As the basal bricks, the dynamics and arrangement of nucleosomes orchestrate the higher architecture of chromatin in a fundamental way, thereby affecting almost all nuclear biology processes. Thanks to its rather simple protocol, assay for transposase-accessible chromatin using sequencing (ATAC)-seq has been rapidly adopted as a major tool for chromatin-accessible profiling at both bulk and single-cell levels; however, to picture the arrangement of nucleosomes per se remains a challenge with ATAC-seq. In the present work, we introduce a novel ATAC-seq analysis toolkit, named decoding nucleosome organization profile based on ATAC-seq data (deNOPA), to predict nucleosome positions. Assessments showed that deNOPA outperformed state-of-the-art tools with ultra-sparse ATAC-seq data, e.g. no more than 0.5 fragment per base pair. The remarkable performance of deNOPA was fueled by the short fragment reads, which compose nearly half of sequenced reads in the ATAC-seq libraries and are commonly discarded by state-of-the-art nucleosome positioning tools. However, we found that the short fragment reads enrich information on nucleosome positions and that the linker regions were predicted by reads from both short and long fragments using Gaussian smoothing. Last, using deNOPA, we showed that the dynamics of nucleosome organization may not directly couple with chromatin accessibility in the cis-regulatory regions when human cells respond to heat shock stimulation. Our deNOPA provides a powerful tool with which to analyze the dynamics of chromatin at nucleosome position level with ultra-sparse ATAC-seq data.

Proximal-end bias from in-vitro reconstituted nucleosomes and the result on downstream data analysis

The most basic level of eukaryotic gene regulation is the presence or absence of nucleosomes on DNA regulatory elements. In an effort to elucidate in vivo nucleosome patterns, in vitro studies are frequently used. In vitro, short DNA fragments are more favorable for nucleosome formation, increasing the likelihood of nucleosome occupancy. This may in part result from the fact that nucleosomes prefer to form on the terminal ends of linear DNA. This phenomenon has the potential to bias in vitro reconstituted nucleosomes and skew results. If the ends of DNA fragments are known, the reads falling close to the ends are typically discarded. In this study we confirm the phenomenon of end bias of in vitro nucleosomes. We describe a method in which nearly identical libraries, with different known ends, are used to recover nucleosomes which form towards the terminal ends of fragmented DNA. Finally, we illustrate that although nucleosomes prefer to form on DNA ends, it does not appear to skew results or the interpretation thereof.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Comparative analysis and prediction of nucleosome positioning using integrative feature representation and machine learning algorithms

Background

Nucleosome plays an important role in the process of genome expression, DNA replication, DNA repair and transcription. Therefore, the research of nucleosome positioning has invariably received extensive attention. Considering the diversity of DNA sequence representation methods, we tried to integrate multiple features to analyze its effect in the process of nucleosome positioning analysis. This process can also deepen our understanding of the theoretical analysis of nucleosome positioning.

Results

Here, we not only used frequency chaos game representation (FCGR) to construct DNA sequence features, but also integrated it with other features and adopted the principal component analysis (PCA) algorithm. Simultaneously, support vector machine (SVM), extreme learning machine (ELM), extreme gradient boosting (XGBoost), multilayer perceptron (MLP) and convolutional neural networks (CNN) are used as predictors for nucleosome positioning prediction analysis, respectively. The integrated feature vector prediction quality is significantly superior to a single feature. After using principal component analysis (PCA) to reduce the feature dimension, the prediction quality of H. sapiens dataset has been significantly improved.

Conclusions

Comparative analysis and prediction on H. sapiens, C. elegans, D. melanogaster and S. cerevisiae datasets, demonstrate that the application of FCGR to nucleosome positioning is feasible, and we also found that integrative feature representation would be better.

The epigenetic basis of cellular heterogeneity

Single-cell sequencing-based methods for profiling gene transcript levels have revealed substantial heterogeneity in expression levels among morphologically indistinguishable cells. This variability has important functional implications for tissue biology and disease states such as cancer. Mapping of epigenomic information such as chromatin accessibility, nucleosome positioning, histone tail modifications and enhancer-promoter interactions in both bulk-cell and single-cell samples has shown that these characteristics of chromatin state contribute to expression or repression of associated genes. Advances in single-cell epigenomic profiling methods are enabling high-resolution mapping of chromatin states in individual cells. Recent studies using these techniques provide evidence that variations in different aspects of chromatin organization collectively define gene expression heterogeneity among otherwise highly similar cells.

Temporal gene regulation by p53 is associated with the rotational setting of its binding sites in nucleosomes

The tumor suppressor protein p53 is a DNA-binding transcription factor (TF) that, once activated, coordinates the expression of thousands of target genes. Increased p53 binding to gene promoters occurs shortly after p53 activation. Intriguingly, gene transcription exhibits differential kinetics with some genes being induced early (early genes) and others being induced late (late genes). To understand pre-binding factors contributing to the temporal gene regulation by p53, we performed time-course RNA sequencing experiments in human colon cancer cell line HCT116 treated with fluorouracil to identify early and late genes. Published p53 ChIP fragments co-localized with the early or late genes were used to uncover p53 binding sites (BS). We demonstrate that the BS associated with early genes are clustered around gene starts with decreased nucleosome occupancy. DNA analysis shows that these BS are likely exposed on nucleosomal surface if wrapped into nucleosomes, thereby facilitating stable interactions with and fast induction by p53. By contrast, p53 BS associated with late genes are distributed uniformly across the genes with increased nucleosome occupancy. Predicted rotational settings of these BS show limited accessibility. We therefore propose a hypothetical model in which the BS are fully, partially or not accessible to p53 in the nucleosomal context. The partial accessibility of the BS allows subunits of a p53 tetramer to bind, but the resulting p53-DNA complex may not be stable enough to recruit cofactors, which leads to delayed induction. Our work highlights the importance of DNA conformations of p53 BS in gene expression dynamics.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Chromosome Structural Mechanics Dictates the Local Spreading of Epigenetic Marks

We present a theoretical model that demonstrates the integral role chromosome organization and structural mechanics play in the spreading of histone modifications involved in epigenetic regulation. Our model shows that heterogeneous nucleosome positioning, and the resulting position-dependent mechanical properties, must be included to reproduce several qualitative features of experimental data of histone methylation spreading around an artificially induced "nucleation site." We show that our model recreates both the extent of spreading and the presence of a subdominant peak upstream of the transcription start site. Our model indicates that the spreading of epigenetic modifications is sensitive to heterogeneity in chromatin organization and the resulting variability in the chromatin's mechanical properties, suggesting that nucleosome spacing can directly control the conferral of epigenetic marks by modifying the structural mechanics of the chromosome. It further illustrates how the physical organization of the DNA polymer may play a significant role in re-establishing the epigenetic code upon cell division.

A Light-Inducible Strain for Genome-Wide Histone Turnover Profiling in Neurospora crassa

In chromatin, nucleosomes are composed of ∼146 bp of DNA wrapped around a histone octamer, and are highly dynamic structures subject to remodeling and exchange. Histone turnover has previously been implicated in various processes including the regulation of chromatin accessibility, segregation of chromatin domains, and dilution of histone marks. Histones in different chromatin environments may turnover at different rates, possibly with functional consequences. Neurospora crassa sports a chromatin environment that is more similar to that of higher eukaryotes than yeasts, which have been utilized in the past to explore histone exchange. We constructed a simple light-inducible system to profile histone exchange in N. crassa on a 3xFLAG-tagged histone H3 under the control of the rapidly inducible vvd promoter. After induction with blue light, incorporation of tagged H3 into chromatin occurred within 20 min. Previous studies of histone turnover involved considerably longer incubation periods and relied on a potentially disruptive change of medium for induction. We used this reporter to explore replication-independent histone turnover at genes and examine changes in histone turnover at heterochromatin domains in different heterochromatin mutant strains. In euchromatin, H3-3xFLAG patterns were almost indistinguishable from that observed in wild-type in all mutant backgrounds tested, suggesting that loss of heterochromatin machinery has little effect on histone turnover in euchromatin. However, turnover at heterochromatin domains increased with loss of trimethylation of lysine 9 of histone H3 or HP1, but did not depend on DNA methylation. Our reporter strain provides a simple yet powerful tool to assess histone exchange across multiple chromatin contexts.

Role of cell-type specific nucleosome positioning in inducible activation of mammalian promoters

The organization of nucleosomes across functional genomic elements represents a critical layer of control. Here, we present a strategy for high-resolution nucleosome profiling at selected genomic features, and use this to analyse dynamic nucleosome positioning at inducible and cell-type-specific mammalian promoters. We find that nucleosome patterning at inducible promoters frequently resembles that at active promoters, even before stimulus-driven activation. Accordingly, the nucleosome profile at many inactive inducible promoters is sufficient to predict cell-type-specific responsiveness. Induction of gene expression is generally not associated with major changes to nucleosome patterning, and a subset of inducible promoters can be activated without stable nucleosome depletion from their transcription start sites. These promoters are generally dependent on remodelling enzymes for their inducible activation, and exhibit transient nucleosome depletion only at alleles undergoing transcription initiation. Together, these data reveal how the responsiveness of inducible promoters to activating stimuli is linked to cell-type-specific nucleosome patterning.

Nucleosome organisation plays important roles in regulating functional genomic elements. Here, the authors use high-resolution profiling to analyse dynamic nucleosome positioning at inducible and cell-type-specific promoters, providing a global view of chromatin architecture at inducible promoters.

Centromere chromatin structure - Lessons from neocentromeres

Centromeres are highly specialized genomic loci that function during mitosis to maintain genome stability. Formed primarily on repetitive α-satellite DNA sequence characterisation of native centromeric chromatin structure has remained challenging. Fortuitously, neocentromeres are formed on a unique DNA sequence and represent an excellent model to interrogate centromeric chromatin structure. This review uncovers the specific findings from independent neocentromere studies that have advanced our understanding of canonical centromere chromatin structure.

An Efficient Exogenous Gene Insertion Site in CHO Cells with High Transcription Level to Enhance AID-Induced Mutation

Antibodies have been extensively used for the purpose of scientific research, clinical diagnosis, and therapy. Combination of in vitro somatic hypermutation and mammalian cell surface display has been an efficient technology for antibody or other proteins optimization, in which the efficiency of activation-induced cytidine deaminase (AID) mutations in genes is one of the most important key factors. Gene transcriptional level has been found to be positively proportional to AID-induced mutation frequency. Thus, construction of the cell clone bearing a gene of interest (GOI) with high transcription level can increase AID-induced mutations. In this study, a retargetable gene cassette is inserted onto predetermined chromosome site (ywhae gene site) which is among the genes with the highest as well as stable transcription, and is found that one subsite is suitable to be retargeted for efficient protein display in Chinese hamster ovary (CHO) cells. The resultant cell clone (T31) has higher and more stable transcription/expression than CHO-puro clone which was previously established through the strategy of random insertion followed by a high-throughput selection. It also possesses a significantly higher mutation frequency to GOI than CHO-puro cells; thus, it is a better clone for the in vitro improvement of antibody affinity, and probably other properties.

Born to run: control of transcription elongation by RNA polymerase II

The dynamic regulation of transcription elongation by RNA polymerase II (Pol II) is an integral part of the implementation of gene expression programmes during development. In most metazoans, the majority of transcribed genes exhibit transient pausing of Pol II at promoter-proximal regions, and the release of Pol II into gene bodies is controlled by many regulatory factors that respond to environmental and developmental cues. Misregulation of the elongation stage of transcription is implicated in cancer and other human diseases, suggesting that mechanistic understanding of transcription elongation control is therapeutically relevant. In this Review, we discuss the features, establishment and maintenance of Pol II pausing, the transition into productive elongation, the control of transcription elongation by enhancers and by factors of other cellular processes, such as topoisomerases and poly(ADP-ribose) polymerases (PARPs), and the potential of therapeutic targeting of the elongation stage of transcription by Pol II.

Transcriptional Landscape of PARs in Epithelial Malignancies

G protein-coupled receptors (GPCRs), the largest family of cell receptors, act as important regulators of diverse signaling pathways. Our understanding of the impact of GPCRs in tumors is emerging, yet there is no therapeutic platform based on GPCR driver genes. As cancer progresses, it disrupts normal epithelial organization and maintains the cells outside their normal niche. The dynamic and flexible microenvironment of a tumor contains both soluble and matrix-immobilized proteases that contribute to the process of cancer advancement. An example is the activation of cell surface protease-activated receptors (PARs). Mammalian PARs are a subgroup of GPCRs that form a family of four members, PAR_1–4, which are uniquely activated by proteases found in the microenvironment. PAR₁ and PAR₂ play central roles in tumor biology, and PAR₃ acts as a coreceptor. The significance of PAR₄ in neoplasia is just beginning to emerge. PAR₁ has been shown to be overexpressed in malignant epithelia, in direct correlation with tumor aggressiveness, but there is no expression in normal epithelium. In this review, the involvement of key transcription factors such as Egr1, p53, Twist, AP2, and Sp1 that control PAR₁ expression levels specifically, as well as hormone transcriptional regulation by both estrogen receptors (ER) and androgen receptors (AR) are discussed. The cloning of the human protease-activated receptor 2; Par2 (hPar2) promoter region and transcriptional regulation of estrogen (E₂) via binding of the E₂–ER complex to estrogen response elements (ERE) are shown. In addition, evidence that TEA domain 4 (TEAD₄) motifs are present within the hPar2 promoter is presented since the YAP oncogene, which plays a central part in tumor etiology, acts via the TEAD₄ transcription factor. As of now, no information is available on regulation of the hPar3 promoter. With regard to hPar4, only data showing CpG methylation promoter regulation is available. Characterization of the PAR transcriptional landscape may identify powerful targets for cancer therapies.

Chromatin Higher-Order Folding: A Perspective with Linker DNA Angles

The mechanism by which the "beads-on-a-string" nucleosome chain folds into various higher-order chromatin structures in eukaryotic cell nuclei is still poorly understood. The various models depicting higher-order chromatin as regular helical fibers and the very opposite "polymer melt" theory imply that interactions between nucleosome "beads" make the main contribution to the chromatin compaction. Other models in which the geometry of linker DNA "strings" entering and exiting the nucleosome define the three-dimensional structure predict that small changes in the linker DNA configuration may strongly affect nucleosome chain folding and chromatin higher-order structure. Among those studies, the cross-disciplinary approach pioneered by Jörg Langowski that combines computational modeling with biophysical and biochemical experiments was most instrumental for understanding chromatin higher-order structure in vitro. Strikingly, many recent studies, including genome-wide nucleosome interaction mapping and chromatin imaging, show an excellent agreement with the results of three-dimensional computational modeling based on the primary role of linker DNA geometry in chromatin compaction. This perspective relates nucleosome array models with experimental studies of nucleosome array folding in vitro and in situ. I argue that linker DNA configuration plays a key role in determining nucleosome chain flexibility, topology, and propensity for self-association, thus providing new implications for regulation of chromatin accessibility to DNA binding factors and RNA transcription machinery as well as long-range communications between distant genomic sites.

HIF-dependent and reversible nucleosome disassembly in hypoxia-inducible gene promoters

Hypoxia causes dramatic changes in gene expression profiles, and the mechanism of hypoxia-inducible transcription has been analyzed for use as a model system of stress-inducible gene regulation. In this study, changes in chromatin organization in promoters of hypoxia-inducible genes were investigated during hypoxia-reoxygenation conditions. Most of the hypoxia-inducible gene promoters were hypersensitive to DNase I under both normal and hypoxic conditions, and our data indicate an immediate recruitment of transcription factors under hypoxic conditions. In some of the hypoxia-inducible promoters, nucleosome-free DNA regions (NFRs) were established in parallel with hypoxia-induced transcription. We also show that the hypoxia-inducible formation of NFRs requires that hypoxia-inducible transcription factors (HIFs) bind to the promoters together with the transcriptional coactivator CBP. Within 1 h after the hypoxia exposure was ended (reoxygenation), HIF complexes were dissociated from the promoter regions. Within 24 h of reoxygenation, the hypoxia-induced transcription returned to basal levels and the nucleosome structure was reassembled in the hypoxia-inducible NFRs. Nucleosome reassembly required the function of the transcriptional coregulator SIN3A. Thus, reversible changes in nucleosome organization mediated by transcription factors are notable features of stress-inducible gene regulation.

Transcription start site profiling uncovers divergent transcription and enhancer-associated RNAs in Drosophila melanogaster

Background

High-resolution transcription start site (TSS) mapping in D. melanogaster embryos and cell lines has revealed a rich and detailed landscape of both cis- and trans-regulatory elements and factors. However, TSS profiling has not been investigated in an orthogonal in vivo setting. Here, we present a comprehensive dataset that links TSS dynamics with nucleosome occupancy and gene expression in the wandering third instar larva, a developmental stage characterized by large-scale shifts in transcriptional programs in preparation for metamorphosis.

Results

The data recapitulate major regulatory classes of TSSs, based on peak width, promoter-proximal polymerase pausing, and cis-regulatory element density. We confirm the paucity of divergent transcription units in D. melanogaster, but also identify notable exceptions. Furthermore, we identify thousands of novel initiation events occurring at unannotated TSSs that can be classified into functional categories by their local density of histone modifications. Interestingly, a sub-class of these unannotated TSSs overlaps with functionally validated enhancer elements, consistent with a regulatory role for “enhancer RNAs” (eRNAs) in defining developmental transcription programs.

Conclusions

High-depth TSS mapping is a powerful strategy for identifying and characterizing low-abundance and/or low-stability RNAs. Global analysis of transcription initiation patterns in a developing organism reveals a vast number of novel initiation events that identify potential eRNAs as well as other non-coding transcripts critical for animal development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4510-7) contains supplementary material, which is available to authorized users.

Cell Type-Specific Survey of Epigenetic Modifications by Tandem Chromatin Immunoprecipitation Sequencing

The nervous system of higher eukaryotes is composed of numerous types of neurons and glia that together orchestrate complex neuronal responses. However, this complex pool of cells typically poses analytical challenges in investigating gene expression profiles and their epigenetic basis for specific cell types. Here, we developed a novel method that enables cell type-specific analyses of epigenetic modifications using tandem chromatin immunoprecipitation sequencing (tChIP-Seq). FLAG-tagged histone H2B, a constitutive chromatin component, was first expressed in Camk2a-positive pyramidal cortical neurons and used to purify chromatin in a cell type-specific manner. Subsequent chromatin immunoprecipitation using antibodies against H3K4me3-a chromatin modification mainly associated with active promoters-allowed us to survey the histone modifications in Camk2a-positive neurons. Indeed, tChIP-Seq identified hundreds of H3K4me3 modifications in promoter regions located upstream of genes associated with neuronal functions and genes with unknown functions in cortical neurons. tChIP-Seq provides a versatile approach to investigating the epigenetic modifications of particular cell types in vivo.

Regulation of chromatin folding by conformational variations of nucleosome linker DNA

Linker DNA conformational variability has been proposed to direct nucleosome array folding into more or less compact chromatin fibers but direct experimental evidence for such models are lacking. Here, we tested this hypothesis by designing nucleosome arrays with A-tracts at specific locations in the nucleosome linkers to induce inward (AT-IN) and outward (AT-OUT) bending of the linker DNA. Using electron microscopy and analytical centrifugation techniques, we observed spontaneous folding of AT-IN nucleosome arrays into highly compact structures, comparable to those induced by linker histone H1. In contrast, AT-OUT nucleosome arrays formed less compact structures with decreased nucleosome interactions similar to wild-type nucleosome arrays. Adding linker histone H1 further increased compaction of the A-tract arrays while maintaining structural differences between them. Furthermore, restriction nuclease digestion revealed a strongly reduced accessibility of nucleosome linkers in the compact AT-IN arrays. Electron microscopy analysis and 3D computational Monte Carlo simulations are consistent with a profound zigzag linker DNA configuration and closer nucleosome proximity in the AT-IN arrays due to inward linker DNA bending. We propose that the evolutionary preferred positioning of A-tracts in DNA linkers may control chromatin higher-order folding and thus influence cellular processes such as gene expression, transcription and DNA repair.

Evolutionary analysis of nucleosome positioning sequences based on New Symmetric Relative Entropy

New Symmetric Relative Entropy (NSRE) was applied innovatively to analyze the nucleosome sequences in S. cerevisiae, S. pombe and Drosophila. NSRE distributions could well reflect the characteristic differences of nucleosome sequences among three organisms, and the differences indicate a concerted evolution in the sequence usage of nucleosome. Further analysis about the nucleosomes around TSS shows that the constitutive property of +1/-1 nucleosomes in S. cerevisiae is different from that in S. pombe and Drosophila, which indicates that S. cerevisiae has a different transcription regulation mechanism based on nucleosome. However, in either case, the nucleosome dyad region is conserved and always has a higher NSRE. Base composition analysis shows that this conservative property in nucleosome dyad region is mainly determined by base A and T, and the dependence degrees on base A and T are consistent in three organisms.

Developmental transitions in Arabidopsis are regulated by antisense RNAs resulting from bidirectionally transcribed genes

Transcription terminators are DNA elements located at the 3' end of genes that ensure efficient cleavage of nascent RNA generating the 3' end of mRNA, as well as facilitating disengagement of elongating DNA-dependent RNA polymerase II. Surprisingly, terminators are also a potent source of antisense transcription. We have recently described an Arabidopsis antisense transcript originating from the 3' end of a master regulator of Arabidopsis thaliana seed dormancy DOG1. In this review, we discuss the broader implications of our discovery in light of recent developments in yeast and Arabidopsis. We show that, surprisingly, the key features of terminators that give rise to antisense transcription are preserved between Arabidopsis and yeast, suggesting a conserved mechanism. We also compare our discovery to known antisense-based regulatory mechanisms, highlighting the link between antisense-based gene expression regulation and major developmental transitions in plants.

Bleomycin analogues preferentially cleave at the transcription start sites of actively transcribed genes in human cells

Bleomycin (BLM) is a cancer chemotherapeutic agent that is used in the treatment of several types of tumours. The cytotoxicity of three BLM analogues, BLM Z, 6'-deoxy-BLM Z and zorbamycin (ZBM), was determined in human HeLa cells in comparison with BLM. It was found that the IC₅₀ values were 2.9μM for 6'-deoxy-BLM Z, 3.2μM for BLM Z, 4.4μM for BLM and 7.9μM for ZBM in HeLa cells. Using next-generation Illumina DNA sequencing techniques, the genome-wide cleavage of DNA by the BLM analogues was determined in human HeLa cells and compared with BLM. It was ascertained that BLM, 6'-deoxy-BLM Z and ZBM preferentially cleaved at the transcription start sites of actively transcribed genes in human cells. The degree of preferential cleavage at the transcription start sites was quantified and an inverse correlation with the IC₅₀ values was observed. This indicated that the degree of preferential cleavage at transcription start sites is an important component in determining the cytotoxicity of BLM analogues.

Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

Sequences that influence nucleosome positioning in promoter regions, and their relation to gene regulation, have been the topic of much research over the last decade. In yeast, significant nucleosome-depleted regions are found, which facilitate transcription. With the arrival of nucleosome positioning maps for the human genome, it was discovered that in our genome, unlike in that of yeast, promoters encode for high nucleosome occupancy. In this work, we look at the genomes of a range of different organisms, to provide a catalog of nucleosome positioning signals in promoters across the tree of life. We utilize a computational model of the nucleosome, based on crystallographic analyses of the structure and elasticity of the nucleosome, to predict the nucleosome positioning signals in promoter regions. To be able to apply our model to large genomic datasets, we introduce an approximative scheme that makes use of the limited range of correlations in nucleosomal sequence preferences to create a computationally efficient approximation of the full biophysical model. Our predictions show that a clear distinction between unicellular and multicellular life is visible in the intrinsically encoded nucleosome affinity. Furthermore, the strength of the nucleosome positioning signals correlates with the complexity of the organism. We conclude that encoding for high nucleosome occupancy, as in the human genome, is in fact a universal feature of multicellular life.

Mathematical Modeling of Avidity Distribution and Estimating General Binding Properties of Transcription Factors from Genome-Wide Binding Profiles

The shape of the experimental frequency distributions (EFD) of diverse molecular interaction events quantifying genome-wide binding is often skewed to the rare but abundant quantities. Such distributions are systematically deviated from standard power-law functions proposed by scale-free network models suggesting that more explanatory and predictive probabilistic model(s) are needed. Identification of the mechanism-based data-driven statistical distributions that provide an estimation and prediction of binding properties of transcription factors from genome-wide binding profiles is the goal of this analytical survey. Here, we review and develop an analytical framework for modeling, analysis, and prediction of transcription factor (TF) DNA binding properties detected at the genome scale. We introduce a mixture probabilistic model of binding avidity function that includes nonspecific and specific binding events. A method for decomposition of specific and nonspecific TF-DNA binding events is proposed. We show that the Kolmogorov-Waring (KW) probability function (PF), modeling the steady state TF binding-dissociation stochastic process, fits well with the EFD for diverse TF-DNA binding datasets. Furthermore, this distribution predicts total number of TF-DNA binding sites (BSs), estimating specificity and sensitivity as well as other basic statistical features of DNA-TF binding when the experimental datasets are noise-rich and essentially incomplete. The KW distribution fits equally well to TF-DNA binding activity for different TFs including ERE, CREB, STAT1, Nanog, and Oct4. Our analysis reveals that the KW distribution and its generalized form provides the family of power-law-like distributions given in terms of hypergeometric series functions, including standard and generalized Pareto and Waring distributions, providing flexible and common skewed forms of the transcription factor binding site (TFBS) avidity distribution function. We suggest that the skewed binding events may be due to a wide range of evolutionary processes of creating weak avidity TFBS associated with random mutations, while the rare high-avidity binding sites (i.e., high-avidity evolutionarily conserved canonical e-boxes) rarely occurred. These, however, may be positively selected in microevolution.

The eukaryotic promoter database in its 30th year: focus on non-vertebrate organisms

We present an update of the Eukaryotic Promoter Database EPD (http://epd.vital-it.ch), more specifically on the EPDnew division, which contains comprehensive organisms-specific transcription start site (TSS) collections automatically derived from next generation sequencing (NGS) data. Thanks to the abundant release of new high-throughput transcript mapping data (CAGE, TSS-seq, GRO-cap) the database could be extended to plant and fungal species. We further report on the expansion of the mass genome annotation (MGA) repository containing promoter-relevant chromatin profiling data and on improvements for the EPD entry viewers. Finally, we present a new data access tool, ChIP-Extract, which enables computational biologists to extract diverse types of promoter-associated data in numerical table formats that are readily imported into statistical analysis platforms such as R.

Genome-wide positioning of bivalent mononucleosomes

Background

Bivalent chromatin refers to overlapping regions containing activating histone H3 Lys4 trimethylation (H3K4me3) and inactivating H3K27me3 marks. Existence of such bivalent marks on the same nucleosome has only recently been suggested. Previous genome-wide efforts to characterize bivalent chromatin have focused primarily on individual marks to define overlapping zones of bivalency rather than mapping positions of truly bivalent mononucleosomes.

Results

Here, we developed an efficacious sequential ChIP technique for examining global positioning of individual bivalent nucleosomes. Using next generation sequencing approaches we show that although individual H3K4me3 and H3K27me3 marks overlap in broad zones, bivalent nucleosomes are focally enriched in the vicinity of the transcription start site (TSS). These seem to occupy the H2A.Z nucleosome positions previously described as salt-labile nucleosomes, and are correlated with low gene expression. Although the enrichment profiles of bivalent nucleosomes show a clear dependency on CpG island content, they demonstrate a stark anti-correlation with methylation status.

Conclusions

We show that regional overlap of H3K4me3 and H3K27me3 chromatin tend to be upstream to the TSS, while bivalent nucleosomes with both marks are mainly promoter proximal near the TSS of CpG island-containing genes with poised/low expression. We discuss the implications of the focal enrichment of bivalent nucleosomes around the TSS on the poised chromatin state of promoters in stem cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0221-6) contains supplementary material, which is available to authorized users.

Fast, Accurate and Automatic Ancient Nucleosome and Methylation Maps with epiPALEOMIX

The first epigenomes from archaic hominins (AH) and ancient anatomically modern humans (AMH) have recently been characterized, based, however, on a limited number of samples. The extent to which ancient genome-wide epigenetic landscapes can be reconstructed thus remains contentious. Here, we present epiPALEOMIX, an open-source and user-friendly pipeline that exploits post-mortem DNA degradation patterns to reconstruct ancient methylomes and nucleosome maps from shotgun and/or capture-enrichment data. Applying epiPALEOMIX to the sequence data underlying 35 ancient genomes including AMH, AH, equids and aurochs, we investigate the temporal, geographical and preservation range of ancient epigenetic signatures. We first assess the quality of inferred ancient epigenetic signatures within well-characterized genomic regions. We find that tissue-specific methylation signatures can be obtained across a wider range of DNA preparation types than previously thought, including when no particular experimental procedures have been used to remove deaminated cytosines prior to sequencing. We identify a large subset of samples for which DNA associated with nucleosomes is protected from post-mortem degradation, and nucleosome positioning patterns can be reconstructed. Finally, we describe parameters and conditions such as DNA damage levels and sequencing depth that limit the preservation of epigenetic signatures in ancient samples. When such conditions are met, we propose that epigenetic profiles of CTCF binding regions can be used to help data authentication. Our work, including epiPALEOMIX, opens for further investigations of ancient epigenomes through time especially aimed at tracking possible epigenetic changes during major evolutionary, environmental, socioeconomic, and cultural shifts.

Evidence of selection for an accessible nucleosomal array in human

BACKGROUND

Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.

RESULTS

We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.

CONCLUSION

The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.

Long-range regulators of the lncRNA HOTAIR enhance its prognostic potential in breast cancer

Predicting response to endocrine therapy and survival in oestrogen receptor positive breast cancer is a significant clinical challenge and novel prognostic biomarkers are needed. Long-range regulators of gene expression are emerging as promising biomarkers and therapeutic targets for human diseases, so we have explored the potential of distal enhancer elements of non-coding RNAs in the prognostication of breast cancer survival. HOTAIR is a long non-coding RNA that is overexpressed, promotes metastasis and is predictive of decreased survival. Here, we describe a long-range transcriptional enhancer of the HOTAIR gene that binds several hormone receptors and associated transcription factors, interacts with the HOTAIR promoter and augments transcription. This enhancer is dependent on Forkhead-Box transcription factors and functionally interacts with a novel alternate HOTAIR promoter. HOTAIR expression is negatively regulated by oestrogen, positively regulated by FOXA1 and FOXM1, and is inversely correlated with oestrogen receptor and directly correlated with FOXM1 in breast tumours. The combination of HOTAIR and FOXM1 enables greater discrimination of endocrine therapy responders and non-responders in patients with oestrogen receptor positive breast cancer. Consistent with this, HOTAIR expression is increased in cell-line models of endocrine resistance. Analysis of breast cancer gene expression data indicates that HOTAIR is co-expressed with FOXA1 and FOXM1 in HER2-enriched tumours, and these factors enhance the prognostic power of HOTAIR in aggressive HER2+ breast tumours. Our study elucidates the transcriptional regulation of HOTAIR, identifies HOTAIR and its regulators as novel biomarkers of patient response to endocrine therapy and corroborates the importance of transcriptional enhancers in cancer.

Identifying and retargeting transcriptional hot spots in the human genome

Mammalian cell line development requires streamlined methodologies that will reduce both the cost and time to identify candidate cell lines. Improvements in site-specific genomic editing techniques can result in flexible, predictable, and robust cell line engineering. However, an outstanding question in the field is the specific site of integration. Here, we seek to identify productive loci within the human genome that will result in stable, high expression of heterologous DNA. Using an unbiased, random integration approach and a green fluorescent reporter construct, we identify ten single-integrant, recombinant human cell lines that exhibit stable, high-level expression. From these cell lines, eight unique corresponding integration loci were identified. These loci are concentrated in non-protein coding regions or intronic regions of protein coding genes. Expression mapping of the surrounding genes reveals minimal disruption of endogenous gene expression. Finally, we demonstrate that targeted de novo integration at one of the identified loci, the 12(th) exon-intron region of the GRIK1 gene on chromosome 21, results in superior expression and stability compared to the standard, illegitimate integration approach at levels approaching 4-fold. The information identified here along with recent advances in site-specific genomic editing techniques can lead to expedited cell line development.

Nucleosome positioning changes during human embryonic stem cell differentiation

Nucleosomes are the basic unit of chromatin. Nucleosome positioning (NP) plays a key role in transcriptional regulation and other biological processes. To better understand NP we used MNase-seq to investigate changes that occur as human embryonic stem cells (hESCs) transition to nascent mesoderm and then to smooth muscle cells (SMCs). Compared to differentiated cell derivatives, nucleosome occupancy at promoters and other notable genic sites, such as exon/intron junctions and adjacent regions, in hESCs shows a stronger correlation with transcript abundance and is less influenced by sequence content. Upon hESC differentiation, genes being silenced, but not genes being activated, display a substantial change in nucleosome occupancy at their promoters. Genome-wide, we detected a shift of NP to regions of higher G+C content as hESCs differentiate to SMCs. Notably, genomic regions with higher nucleosome occupancy harbor twice as many G↔C changes but fewer than half A↔T changes, compared to regions with lower nucleosome occupancy. Finally, our analysis indicates that the hESC genome is not rearranged and has a sequence mutation rate resembling normal human genomes. Our study reveals another unique feature of hESC chromatin, and sheds light on the relationship between nucleosome occupancy and sequence G+C content.

Whole genome nucleosome sequencing identifies novel types of forensic markers in degraded DNA samples

In the case of mass disasters, missing persons and forensic caseworks, highly degraded biological samples are often encountered. It can be a challenge to analyze and interpret the DNA profiles from these samples. Here we provide a new strategy to solve the problem by taking advantage of the intrinsic structural properties of DNA. We have assessed the in vivo positions of more than 35 million putative nucleosome cores in human leukocytes using high-throughput whole genome sequencing, and identified 2,462 single nucleotide variations (SNVs), 128 insertion-deletion polymorphisms (indels). After comparing the sequence reads with 44 STR loci commonly used in forensics, five STRs (TH01, TPOX, D18S51, DYS391, and D10S1248)were matched. We compared these “nucleosome protected STRs” (NPSTRs) with five other non-NPSTRs using mini-STR primer design, real-time PCR, and capillary gel electrophoresis on artificially degraded DNA. Moreover, genotyping performance of the five NPSTRs and five non-NPSTRs was also tested with real casework samples. All results show that loci located in nucleosomes are more likely to be successfully genotyped in degraded samples. In conclusion, after further strict validation, these markers could be incorporated into future forensic and paleontology identification kits, resulting in higher discriminatory power for certain degraded sample types.

Sequence features accurately predict genome-wide MeCP2 binding in vivo

Methyl-CpG binding protein 2 (MeCP2) is critical for proper brain development and expressed at near-histone levels in neurons, but the mechanism of its genomic localization remains poorly understood. Using high-resolution MeCP2-binding data, we show that DNA sequence features alone can predict binding with 88% accuracy. Integrating MeCP2 binding and DNA methylation in a probabilistic graphical model, we demonstrate that previously reported genome-wide association with methylation is in part due to MeCP2's affinity to GC-rich chromatin, a result replicated using published data. Furthermore, MeCP2 co-localizes with nucleosomes. Finally, MeCP2 binding downstream of promoters correlates with increased expression in Mecp2-deficient neurons.

MeCP2 is critical for proper brain development, and mutations in the gene encoding MeCP2 are responsible for several neurological disorders. Here, the authors show that the previously reported genome-wide preference of MeCP2 to methylated CpGs is in part due to MeCP2's affinity to GC-rich chromatin.

Altered primary chromatin structures and their implications in cancer development

BACKGROUND

Cancer development is a complex process involving both genetic and epigenetic changes. Genetic changes in oncogenes and tumor-suppressor genes are generally considered as primary causes, since these genes may directly regulate cellular growth. In addition, it has been found that changes in epigenetic factors, through mutation or altered gene expression, may contribute to cancer development. In the nucleus of eukaryotic cells DNA and histone proteins form a structure called chromatin which consists of nucleosomes that, like beads on a string, are aligned along the DNA strand. Modifications in chromatin structure are essential for cell type-specific activation or repression of gene transcription, as well as other processes such as DNA repair, DNA replication and chromosome segregation. Alterations in epigenetic factors involved in chromatin dynamics may accelerate cell cycle progression and, ultimately, result in malignant transformation. Abnormal expression of remodeler and modifier enzymes, as well as histone variants, may confer to cancer cells the ability to reprogram their genomes and to yield, maintain or exacerbate malignant hallmarks. At the end, genetic and epigenetic alterations that are encountered in cancer cells may culminate in chromatin changes that may, by altering the quantity and quality of gene expression, promote cancer development.

METHODS

During the last decade a vast number of studies has uncovered epigenetic abnormalities that are associated with the (anomalous) packaging and remodeling of chromatin in cancer genomes. In this review I will focus on recently published work dealing with alterations in the primary structure of chromatin resulting from imprecise arrangements of nucleosomes along DNA, and its functional implications for cancer development.

CONCLUSIONS

The primary chromatin structure is regulated by a variety of epigenetic mechanisms that may be deregulated through gene mutations and/or gene expression alterations. In recent years, it has become evident that changes in chromatin structure may coincide with the occurrence of cancer hallmarks. The functional interrelationships between such epigenetic alterations and cancer development are just becoming manifest and, therefore, the oncology community should continue to explore the molecular mechanisms governing the primary chromatin structure, both in normal and in cancer cells, in order to improve future approaches for cancer detection, prevention and therapy, as also for circumventing drug resistance.