Characterizing the nuclease accessibility of DNA in human cells to map higher order structures of chromatin

In vivo assembly of complete eukaryotic nucleosomes and (H3-H4)-only non-canonical nucleosomal particles in the model bacterium Escherichia coli

As a fundamental unit for packaging genomic DNA into chromatin, the eukaryotic nucleosome core comprises a canonical octamer with two copies for each histone, H2A, H2B, H3, and H4, wrapped around with 147 base pairs of DNA. While H3 and H4 share structure-fold with archaeal histone-like proteins, the eukaryotic nucleosome core and the complete nucleosome (the core plus H1 histone) are unique to eukaryotes. To explore whether the eukaryotic nucleosome can assemble in prokaryotes and to reconstruct the possible route for its emergence in eukaryogenesis, we developed an in vivo system for assembly of nucleosomes in the model bacterium, Escherichia coli, and successfully reconstituted the core nucleosome, the complete nucleosome, and unexpectedly the non-canonical (H3-H4)4 octasome. The core and complete nucleosomes assembled in E. coli exhibited footprints typical of eukaryotic hosts after in situ micrococcal nuclease digestion. Additionally, they caused condensation of E. coli nucleoid. We also demonstrated the stable formation of non-canonical (H3-H4)2 tetrasome and (H3-H4)4 octasomes in vivo, which are suggested to be 'fossil complex' that marks the intermediate in the progressive development of eukaryotic nucleosome. The study presents a unique platform in a bacterium for in vivo assembly and studying the properties of non-canonical variants of nucleosome.

Dinochromosome Heterotermini with Telosomal Anchorages

Dinoflagellate birefringent chromosomes (BfCs) contain some of the largest known genomes, yet they lack typical nucleosomal micrococcal-nuclease protection patterns despite containing variant core histones. One BfC end interacts with extranuclear mitotic microtubules at the nuclear envelope (NE), which remains intact throughout the cell cycle. Ultrastructural studies, polarized light and fluorescence microscopy, and micrococcal nuclease-resistant profiles (MNRPs) revealed that NE-associated chromosome ends persisted post-mitosis. Histone H3K9me3 inhibition caused S-G2 delay in synchronous cells, without any effects at G1. Differential labeling and nuclear envelope swelling upon decompaction indicate an extension of the inner compartment into telosomal anchorages (TAs). Additionally, limited effects of low-concentration sirtinol on bulk BfCs, coupled with distinct mobility patterns in MNase-digested and psoralen-crosslinked nuclei observed on 2D gels, suggest that telomeric nucleosomes (TNs) are the primary histone structures. The absence of a nucleosomal ladder with cDNA probes, the presence of histone H2A and telomere-enriched H3.3 variants, along with the immuno-localization of H3 variants mainly at the NE further reinforce telomeric regions as the main nucleosomal domains. Cumulative biochemical and molecular analyses suggest that telomeric repeats constitute the major octameric MNRPs that provision chromosomal anchorage at the NE.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

nucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

Chromatin structure and dynamics: one nucleosome at a time

Eukaryotic genomes store information on many levels, including their linear DNA sequence, the posttranslational modifications of its constituents (epigenetic modifications), and its three-dimensional folding. Understanding how this information is stored and read requires multidisciplinary collaborations from many branches of science beyond biology, including physics, chemistry, and computer science. Concurrent recent developments in all these areas have enabled researchers to image the genome with unprecedented spatial and temporal resolution. In this review, we focus on what single-molecule imaging and tracking of individual proteins in live cells have taught us about chromatin structure and dynamics. Starting with the basics of single-molecule tracking (SMT), we describe some advantages over in situ imaging techniques and its current limitations. Next, we focus on single-nucleosome studies and what they have added to our current understanding of the relationship between chromatin dynamics and transcription. In celebration of Robert Feulgen's ground-breaking discovery that allowed us to start seeing the genome, we discuss current models of chromatin structure and future challenges ahead.

Interplay between the transcription preinitiation complex and the +1 nucleosome

Eukaryotic transcription starts with the assembly of a preinitiation complex (PIC) on core promoters. Flanking this region is the +1 nucleosome, the first nucleosome downstream of the core promoter. While this nucleosome is rich in epigenetic marks and plays a key role in transcription regulation, how the +1 nucleosome interacts with the transcription machinery has been a long-standing question. Here, we summarize recent structural and functional studies of the +1 nucleosome in complex with the PIC. We specifically focus on how differently organized promoter-nucleosome templates affect the assembly of the PIC and PIC-Mediator on chromatin and result in distinct transcription initiation.

Is euchromatin really open in the cell?

Genomic DNA is wrapped around a core histone octamer and forms a nucleosome. In higher eukaryotic cells, strings of nucleosomes are irregularly folded as chromatin domains that act as functional genome units. According to a typical textbook model, chromatin can be categorized into two types, euchromatin and heterochromatin, based on its degree of compaction. Euchromatin is open, while heterochromatin is closed and condensed. However, is euchromatin really open in the cell? New evidence from genomics and advanced imaging studies has revealed that euchromatin consists of condensed liquid-like domains. Condensed chromatin seems to be the default chromatin state in higher eukaryotic cells. We discuss this novel view of euchromatin in the cell and how the revealed organization is relevant to genome functions.

The chromatin - triple helix connection

Mammalian genomes are extensively transcribed, producing a large number of coding and non-coding transcripts. A large fraction of the nuclear RNAs is physically associated with chromatin, functioning in gene activation and silencing, shaping higher-order genome organisation, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions. Different mechanisms allow the tethering of these chromatin-associated RNAs (caRNA) to chromosomes, including RNA binding proteins, the RNA polymerases and R-loops. In this review, we focus on the sequence-specific targeting of RNA to DNA by forming triple helical structures and describe its interplay with chromatin. It turns out that nucleosome positioning at triple helix target sites and the nucleosome itself are essential factors in determining the formation and stability of triple helices. The histone H3-tail plays a critical role in triple helix stabilisation, and the role of its epigenetic modifications in this process is discussed.

Changes in adenoviral chromatin organization precede early gene activation upon infection

Within the virion, adenovirus DNA associates with the virus-encoded, protamine-like structural protein pVII. Whether this association is organized, and how genome packaging changes during infection and subsequent transcriptional activation is currently unclear. Here, we combined RNA-seq, MNase-seq, ChIP-seq, and single genome imaging during early adenovirus infection to unveil the structure- and time-resolved dynamics of viral chromatin changes as well as their correlation with gene transcription. Our MNase mapping data indicates that the adenoviral genome is arranged in precisely positioned nucleoprotein particles with nucleosome-like characteristics, that we term adenosomes. We identified 238 adenosomes that are positioned by a DNA sequence code and protect about 60-70 bp of DNA. The incoming adenoviral genome is more accessible at early gene loci that undergo additional chromatin de-condensation upon infection. Histone H3.3 containing nucleosomes specifically replaces pVII at distinct genomic sites and at the transcription start sites of early genes. Acetylation of H3.3 is predominant at the transcription start sites and precedes transcriptional activation. Based on our results, we propose a central role for the viral pVII nucleoprotein architecture, which is required for the dynamic structural changes during early infection, including the regulation of nucleosome assembly prior to transcription initiation. Our study thus may aid the rational development of recombinant adenoviral vectors exhibiting sustained expression in gene therapy.

Structural basis of transcription reduction by a promoter-proximal +1 nucleosome

At active human genes, the +1 nucleosome is located downstream of the RNA polymerase II (RNA Pol II) pre-initiation complex (PIC). However, at inactive genes, the +1 nucleosome is found further upstream, at a promoter-proximal location. Here, we establish a model system to show that a promoter-proximal +1 nucleosome can reduce RNA synthesis in vivo and in vitro, and we analyze its structural basis. We find that the PIC assembles normally when the edge of the +1 nucleosome is located 18 base pairs (bp) downstream of the transcription start site (TSS). However, when the nucleosome edge is located further upstream, only 10 bp downstream of the TSS, the PIC adopts an inhibited state. The transcription factor IIH (TFIIH) shows a closed conformation and its subunit XPB contacts DNA with only one of its two ATPase lobes, inconsistent with DNA opening. These results provide a mechanism for nucleosome-dependent regulation of transcription initiation.

Comprehensive computational analysis of epigenetic descriptors affecting CRISPR-Cas9 off-target activity

BACKGROUND

A common issue in CRISPR-Cas9 genome editing is off-target activity, which prevents the widespread use of CRISPR-Cas9 in medical applications. Among other factors, primary chromatin structure and epigenetics may influence off-target activity.

METHODS

In this work, we utilize crisprSQL, an off-target database, to analyze the effect of 19 epigenetic descriptors on CRISPR-Cas9 off-target activity. Termed as 19 epigenetic features/scores, they consist of 6 experimental epigenetic and 13 computed nucleosome organization-related features. In terms of novel features, 15 of the epigenetic scores are newly considered. The 15 newly considered scores consist of 13 freshly computed nucleosome occupancy/positioning scores and 2 experimental features (MNase and DRIP). The other 4 existing scores are experimental features (CTCF, DNase I, H3K4me3, RRBS) commonly used in deep learning models for off-target activity prediction. For data curation, MNase was aggregated from existing experimental nucleosome occupancy data. Based on the sequence context information available in crisprSQL, we also computed nucleosome occupancy/positioning scores for off-target sites.

RESULTS

To investigate the relationship between the 19 epigenetic features and off-target activity, we first conducted Spearman and Pearson correlation analysis. Such analysis shows that some computed scores derived from training-based models and training-free algorithms outperform all experimental epigenetic features. Next, we evaluated the contribution of all epigenetic features in two successful machine/deep learning models which predict off-target activity. We found that some computed scores, unlike all 6 experimental features, significantly contribute to the predictions of both models. As a practical research contribution, we make the off-target dataset containing all 19 epigenetic features available to the research community.

CONCLUSIONS

Our comprehensive computational analysis helps the CRISPR-Cas9 community better understand the relationship between epigenetic features and CRISPR-Cas9 off-target activity.

Nucleosome-Omics: A Perspective on the Epigenetic Code and 3D Genome Landscape

Genetic information is loaded on chromatin, which involves DNA sequence arrangement and the epigenetic landscape. The epigenetic information including DNA methylation, nucleosome positioning, histone modification, 3D chromatin conformation, and so on, has a crucial impact on gene transcriptional regulation. Out of them, nucleosomes, as basal chromatin structural units, play an important central role in epigenetic code. With the discovery of nucleosomes, various nucleosome-level technologies have been developed and applied, pushing epigenetics to a new climax. As the underlying methodology, next-generation sequencing technology has emerged and allowed scientists to understand the epigenetic landscape at a genome-wide level. Combining with NGS, nucleosome-omics (or nucleosomics) provides a fresh perspective on the epigenetic code and 3D genome landscape. Here, we summarized and discussed research progress in technology development and application of nucleosome-omics. We foresee the future directions of epigenetic development at the nucleosome level.

Phosphorylation of the FACT histone chaperone subunit SPT16 affects chromatin at RNA polymerase II transcriptional start sites in Arabidopsis

The heterodimeric histone chaperone FACT, consisting of SSRP1 and SPT16, contributes to dynamic nucleosome rearrangements during various DNA-dependent processes including transcription. In search of post-translational modifications that may regulate the activity of FACT, SSRP1 and SPT16 were isolated from Arabidopsis cells and analysed by mass spectrometry. Four acetylated lysine residues could be mapped within the basic C-terminal region of SSRP1, while three phosphorylated serine/threonine residues were identified in the acidic C-terminal region of SPT16. Mutational analysis of the SSRP1 acetylation sites revealed only mild effects. However, phosphorylation of SPT16 that is catalysed by protein kinase CK2, modulates histone interactions. A non-phosphorylatable version of SPT16 displayed reduced histone binding and proved inactive in complementing the growth and developmental phenotypes of spt16 mutant plants. In plants expressing the non-phosphorylatable SPT16 version we detected at a subset of genes enrichment of histone H3 directly upstream of RNA polymerase II transcriptional start sites (TSSs) in a region that usually is nucleosome-depleted. This suggests that some genes require phosphorylation of the SPT16 acidic region for establishing the correct nucleosome occupancy at the TSS of active genes.

tagHi-C Reveals 3D Chromatin Architecture Dynamics during Mouse Hematopoiesis

Spatiotemporal chromatin reorganization during hematopoietic differentiation has not been comprehensively characterized, mainly because of the large numbers of starting cells required for current chromatin conformation capture approaches. Here, we introduce a low-input tagmentation-based Hi-C (tagHi-C) method to capture the chromatin structures of hundreds of cells. Using tagHi-C, we are able to map the spatiotemporal dynamics of chromatin structure in ten primary hematopoietic stem, progenitor, and differentiated cell populations from mouse bone marrow. Our results reveal that changes in compartment dynamics and the Rabl configuration occur during hematopoietic cell differentiation. We identify gene-body-associating domains (GADs) as general structures for highly expressed genes. Moreover, we extend the body of knowledge regarding genes influenced by genome-wide association study (GWAS) loci through spatial chromatin looping. Our study provides the tagHi-C method for studying the three-dimensional (3D) genome of a small number of cells and maps the comprehensive 3D chromatin landscape of bone marrow hematopoietic cells.

Chromosomal Mcm2-7 distribution and the genome replication program in species from yeast to humans

The spatio-temporal program of genome replication across eukaryotes is thought to be driven both by the uneven loading of pre-replication complexes (pre-RCs) across the genome at the onset of S-phase, and by differences in the timing of activation of these complexes during S phase. To determine the degree to which distribution of pre-RC loading alone could account for chromosomal replication patterns, we mapped the binding sites of the Mcm2-7 helicase complex (MCM) in budding yeast, fission yeast, mouse and humans. We observed similar individual MCM double-hexamer (DH) footprints across the species, but notable differences in their distribution: Footprints in budding yeast were more sharply focused compared to the other three organisms, consistent with the relative sequence specificity of replication origins in S. cerevisiae. Nonetheless, with some clear exceptions, most notably the inactive X-chromosome, much of the fluctuation in replication timing along the chromosomes in all four organisms reflected uneven chromosomal distribution of pre-replication complexes.

Gene-rich regions of the genome tend to replicate earlier in S phase than do repetitive and other non-genic regions. This may be an evolutionary consequence of the fact that replication later in S phase is associated with higher frequencies of mutation and genome rearrangement. Replication timing along the chromosome is determined by 1) events prior to S-phase that specify the locations where DNA replication can be initiated, referred to as origin licensing; and 2) the timing of activation of these licensed origins during S-phase, referred to as origin firing. To determine the relative importance of these two mechanisms, here we identify both the binding sites and the abundance of a key component of the origin licensing machinery in budding yeast, fission yeast, mice, and humans, namely the replicative helicase complex. We discovered that, with a few notable exceptions, which include the inactive X chromosome in mammals, the program of replication timing can be largely explained simply on the basis of origin licensing. Our results support a model for replication timing that emphasizes stochastic firing of origins that have been licensed before S phase begins.

Local states of chromatin compaction at transcription start sites control transcription levels

The ‘open’ and ‘compact’ regions of chromatin are considered to be regions of active and silent transcription, respectively. However, individual genes produce transcripts at different levels, suggesting that transcription output does not depend on the simple open-compact conversion of chromatin, but on structural variations in chromatin itself, which so far have remained elusive. In this study, weakly crosslinked chromatin was subjected to sedimentation velocity centrifugation, which fractionated the chromatin according to its degree of compaction. Open chromatin remained in upper fractions, while compact chromatin sedimented to lower fractions depending on the level of nucleosome assembly. Although nucleosomes were evenly detected in all fractions, histone H1 was more highly enriched in the lower fractions. H1 was found to self-associate and crosslinked to histone H3, suggesting that H1 bound to H3 interacts with another H1 in an adjacent nucleosome to form compact chromatin. Genome-wide analyses revealed that nearly the entire genome consists of compact chromatin without differences in compaction between repeat and non-repeat sequences; however, active transcription start sites (TSSs) were rarely found in compact chromatin. Considering the inverse correlation between chromatin compaction and RNA polymerase binding at TSSs, it appears that local states of chromatin compaction determine transcription levels.

Conserved noncoding sequences provide insights into regulatory sequence and loss of gene expression in maize

Thousands of species will be sequenced in the next few years; however, understanding how their genomes work, without an unlimited budget, requires both molecular and novel evolutionary approaches. We developed a sensitive sequence alignment pipeline to identify conserved noncoding sequences (CNSs) in the Andropogoneae tribe (multiple crop species descended from a common ancestor ∼18 million years ago). The Andropogoneae share similar physiology while being tremendously genomically diverse, harboring a broad range of ploidy levels, structural variation, and transposons. These contribute to the potential of Andropogoneae as a powerful system for studying CNSs and are factors we leverage to understand the function of maize CNSs. We found that 86% of CNSs were comprised of annotated features, including introns, UTRs, putative cis-regulatory elements, chromatin loop anchors, noncoding RNA (ncRNA) genes, and several transposable element superfamilies. CNSs were enriched in active regions of DNA replication in the early S phase of the mitotic cell cycle and showed different DNA methylation ratios compared to the genome-wide background. More than half of putative cis-regulatory sequences (identified via other methods) overlapped with CNSs detected in this study. Variants in CNSs were associated with gene expression levels, and CNS absence contributed to loss of gene expression. Furthermore, the evolution of CNSs was associated with the functional diversification of duplicated genes in the context of maize subgenomes. Our results provide a quantitative understanding of the molecular processes governing the evolution of CNSs in maize.

Peculiarities of Plasmodium falciparum Gene Regulation and Chromatin Structure

The highly complex life cycle of the human malaria parasite, Plasmodium falciparum, is based on an orchestrated and tightly regulated gene expression program. In general, eukaryotic transcription regulation is determined by a combination of sequence-specific transcription factors binding to regulatory DNA elements and the packaging of DNA into chromatin as an additional layer. The accessibility of regulatory DNA elements is controlled by the nucleosome occupancy and changes of their positions by an active process called nucleosome remodeling. These epigenetic mechanisms are poorly explored in P. falciparum. The parasite genome is characterized by an extraordinarily high AT-content and the distinct architecture of functional elements, and chromatin-related proteins also exhibit high sequence divergence compared to other eukaryotes. Together with the distinct biochemical properties of nucleosomes, these features suggest substantial differences in chromatin-dependent regulation. Here, we highlight the peculiarities of epigenetic mechanisms in P. falciparum, addressing chromatin structure and dynamics with respect to their impact on transcriptional control. We focus on the specialized chromatin remodeling enzymes and discuss their essential function in P. falciparum gene regulation.

Molecular and computational approaches to map regulatory elements in 3D chromatin structure

Epigenetic marks do not change the sequence of DNA but affect gene expression in a cell-type specific manner by altering the activities of regulatory elements. Development of new molecular biology assays, sequencing technologies, and computational approaches enables us to profile the human epigenome in three-dimensional structure genome-wide. Here we describe various molecular biology techniques and bioinformatic tools that have been developed to measure the activities of regulatory elements and their chromatin interactions. Moreover, we list currently available three-dimensional epigenomic data sets that are generated in various human cell types and tissues to assist in the design and analysis of research projects.

InVivo Cancer-Based Functional Genomics

Pinpointing the underlying mechanisms that drive tumorigenesis in human patients is a prerequisite for identifying suitable therapeutic targets for precision medicine. In contrast to cell culture systems, mouse models are highly favored for evaluating tumor progression and therapeutic response in a more realistic in vivo context. The past decade has witnessed a dramatic increase in the number of functional genomic studies using diverse mouse models, including in vivo clustered regularly interspaced short palindromic repeats (CRISPR) and RNA interference (RNAi) screens, and these have provided a wealth of knowledge addressing multiple essential questions in translational cancer research. We compare the multiple mouse systems and genomic tools that are commonly used for in vivo screens to illustrate their strengths and limitations. Crucial components of screen design and data analysis are also discussed.

Identification of determinants of differential chromatin accessibility through a massively parallel genome-integrated reporter assay

A key mechanism in cellular regulation is the ability of the transcriptional machinery to physically access DNA. Transcription factors interact with DNA to alter the accessibility of chromatin, which enables changes to gene expression during development or disease or as a response to environmental stimuli. However, the regulation of DNA accessibility via the recruitment of transcription factors is difficult to study in the context of the native genome because every genomic site is distinct in multiple ways. Here we introduce the multiplexed integrated accessibility assay (MIAA), an assay that measures chromatin accessibility of synthetic oligonucleotide sequence libraries integrated into a controlled genomic context with low native accessibility. We apply MIAA to measure the effects of sequence motifs on cell type-specific accessibility between mouse embryonic stem cells and embryonic stem cell-derived definitive endoderm cells, screening 7905 distinct DNA sequences. MIAA recapitulates differential accessibility patterns of 100-nt sequences derived from natively differential genomic regions, identifying E-box motifs common to epithelial-mesenchymal transition driver transcription factors in stem cell-specific accessible regions that become repressed in endoderm. We show that a single binding motif for a key regulatory transcription factor is sufficient to open chromatin, and classify sets of stem cell-specific, endoderm-specific, and shared accessibility-modifying transcription factor motifs. We also show that overexpression of two definitive endoderm transcription factors, T and Foxa2, results in changes to accessibility in DNA sequences containing their respective DNA-binding motifs and identify preferential motif arrangements that influence accessibility.

Enhanced Detection of DNA Viruses in the Cerebrospinal Fluid of Encephalitis Patients Using Metagenomic Next-Generation Sequencing

The long and expanding list of viral pathogens associated with causing encephalitis confounds current diagnostic procedures, and in up to 50% of cases, the etiology remains undetermined. Sequence-agnostic metagenomic next-generation sequencing (mNGS) obviates the need to specify targets in advance and thus has great potential in encephalitis diagnostics. However, the low relative abundance of viral nucleic acids in clinical specimens poses a significant challenge. Our protocol employs two novel techniques to selectively remove human material at two stages, significantly increasing the representation of viral material. Our bioinformatic workflow using open source protein- and nucleotide sequence-matching software balances sensitivity and specificity in diagnosing and characterizing any DNA viruses present. A panel of 12 cerebrospinal fluid (CSFs) from encephalitis cases was retrospectively interrogated by mNGS, with concordant results in seven of nine samples with a definitive DNA virus diagnosis, and a different herpesvirus was identified in the other two. In two samples with an inconclusive diagnosis, DNA viruses were detected and in a virus-negative sample, no viruses were detected. This assay has the potential to detect DNA virus infections in cases of encephalitis of unknown etiology and to improve the current screening tests by identifying new and emerging agents.

Interplay between TERT promoter mutations and methylation culminates in chromatin accessibility and TERT expression

The telomerase reverse transcriptase (TERT) gene is responsible for telomere maintenance in germline and stem cells, and is re-expressed in 90% of human cancers. CpG methylation in the TERT promoter (TERTp) was correlated with TERT mRNA expression. Furthermore, two hotspot mutations in TERTp, dubbed C228T and C250T, have been revealed to facilitate binding of transcription factor ETS/TCF and subsequent TERT expression. This study aimed to elucidate the combined contribution of epigenetic (promoter methylation and chromatin accessibility) and genetic (promoter mutations) mechanisms in regulating TERT gene expression in healthy skin samples and in melanoma cell lines (n = 61). We unexpectedly observed that the methylation of TERTp was as high in a subset of healthy skin cells, mainly keratinocytes, as in cutaneous melanoma cell lines. In spite of the high promoter methylation fraction in wild-type (WT) samples, TERT mRNA was only expressed in the melanoma cell lines with either high methylation or intermediate methylation in combination with TERT mutations. TERTp methylation was positively correlated with chromatin accessibility and TERT mRNA expression in 8 melanoma cell lines. Cooperation between epigenetic and genetic mechanisms were best observed in heterozygous mutant cell lines as chromosome accessibility preferentially concerned the mutant allele. Combined, these results suggest a complex model in which TERT expression requires either a widely open chromatin state in TERTp-WT samples due to high methylation throughout the promoter or a combination of moderate methylation fraction/chromatin accessibility in the presence of the C228T or C250T mutations.

Protect-seq: genome-wide profiling of nuclease inaccessible domains reveals physical properties of chromatin

In metazoan cell nuclei, heterochromatin constitutes large chromatin domains that are in close contact with the nuclear lamina. These heterochromatin/lamina-associated domains (LADs) domains are difficult to profile and warrants a simpler and direct method. Here we report a new method, Protect-seq, aimed at identifying regions of heterochromatin via resistance to nuclease degradation followed by next-generation sequencing (NGS). We performed Protect-seq on the human colon cancer cell line HCT-116 and observed overlap with previously curated LADs. We provide evidence that these protected regions are enriched for and can distinguish between the repressive histone modification H3K9me3, H3K9me2 and H3K27me3. Moreover, in human cells the loss of H3K9me3 leads to an increase in chromatin accessibility and loss of Protect-seq signal. For further validation, we performed Protect-seq in the fibrosarcoma cell line HT1080 and found a similar correlation with previously curated LADs and repressive histone modifications. In sum, Protect-seq is an efficient technique that allows rapid identification of nuclease resistant chromatin, which correlate with heterochromatin and radial positioning.

Biased visibility in Hi-C datasets marks dynamically regulated condensed and decondensed chromatin states genome-wide

Background

Proximity ligation based techniques, like Hi-C, involve restriction digestion followed by ligation of formaldehyde cross-linked chromatin. Distinct chromatin states can impact the restriction digestion, and hence the visibility in the contact maps, of engaged loci. Yet, the extent and the potential impact of digestion bias remain obscure and under-appreciated in the literature.

Results

Through analysis of 45 Hi-C datasets, lamina-associated domains (LADs), inactive X-chromosome in mammals, and polytene bands in fly, we first established that the DNA in condensed chromatin had lesser accessibility to restriction endonucleases used in Hi-C as compared to that in decondensed chromatin. The observed bias was independent of known systematic biases, was not appropriately corrected by existing computational methods, and needed an additional optimization step. We then repurposed this bias to identify novel condensed domains outside LADs, which were bordered by insulators and were dynamically associated with the polycomb mediated epigenetic and transcriptional states during development.

Conclusions

Our observations suggest that the corrected one-dimensional read counts of existing Hi-C datasets can be reliably repurposed to study the gene-regulatory dynamics associated with chromatin condensation and decondensation, and that the existing Hi-C datasets should be interpreted with cautions.

Physical and data structure of 3D genome

With the textbook view of chromatin folding based on the 30-nm fiber being challenged, it has been proposed that interphase DNA has an irregular 10-nm nucleosome polymer structure whose folding philosophy is unknown. Nevertheless, experimental advances suggest that this irregular packing is associated with many nontrivial physical properties that are puzzling from a polymer physics point of view. Here, we show that the reconciliation of these exotic properties necessitates modularizing three-dimensional genome into tree data structures on top of, and in striking contrast to, the linear topology of DNA double helix. These functional modules need to be connected and isolated by an open backbone that results in porous and heterogeneous packing in a quasi-self-similar manner, as revealed by our electron and optical imaging. Our multiscale theoretical and experimental results suggest the existence of higher-order universal folding principles for a disordered chromatin fiber to avoid entanglement and fulfill its biological functions.

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.

Epigenetic aspects of floral homeotic genes in relation to sexual dimorphism in the dioecious plant Mercurialis annua

In plants, dioecy characterizes species that carry male and female flowers on separate plants and it occurs in about 6% of angiosperms; however, the molecular mechanisms that underlie dioecy are essentially unknown. The ability for sex-reversal by hormone application raises the hypothesis that the genes required for the expression of both sexes are potentially functional but are regulated by epigenetic means. In this study, proteomic analysis of nuclear proteins isolated from flower buds of females, males, and feminized males of the dioecious plant Mercurialis annua revealed differential expression of nuclear proteins that are implicated in chromatin structure and function, including floral homeotic proteins. Focusing on floral genes, we found that class B genes were mainly expressed in male flowers, while class D genes, as well as SUPERMAN-like genes, were mainly expressed in female flowers. Cytokinin-induced feminization of male plants was associated with down-regulation of male-specific genes concomitantly with up-regulation of female-specific genes. No correlation was found between the expression of class B and D genes and the changes in DNA methylation or chromatin conformation of these genes. Thus, we could not confirm DNA methylation or chromatin conformation of floral genes to be the major determinant regulating sexual dimorphisms. Instead, determination of sex in M. annua might be controlled upstream of floral genes by one or more sex-specific factors that affect hormonal homeostasis. A comprehensive model is proposed for sex-determination in M. annua.

Quantitative MNase-seq accurately maps nucleosome occupancy levels

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, its aggressive endo-/exo-nuclease activities make MNase-seq unreliable for determining nucleosome occupancies, because cleavages within linker regions produce oligo- and mono-nucleosomes, whereas cleavages within nucleosomes destroy them. Here, we introduce a theoretical framework for predicting nucleosome occupancies and an experimental protocol with appropriate spike-in normalization that confirms our theory and provides accurate occupancy levels over an MNase digestion time course. As with human cells, we observe no overall differences in nucleosome occupancies between Drosophila euchromatin and heterochromatin, which implies that heterochromatic compaction does not reduce MNase accessibility of linker DNA.

Accessibility of promoter DNA is not the primary determinant of chromatin-mediated gene regulation

DNA accessibility is thought to be of major importance in regulating gene expression. We test this hypothesis using a restriction enzyme as a probe of chromatin structure and as a proxy for transcription factors. We measured the digestion rate and the fraction of accessible DNA at almost all genomic AluI sites in budding yeast and mouse liver nuclei. Hepatocyte DNA is more accessible than yeast DNA, consistent with longer linkers between nucleosomes, suggesting that nucleosome spacing is a major determinant of accessibility. DNA accessibility varies from cell to cell, such that essentially no sites are accessible or inaccessible in every cell. AluI sites in inactive mouse promoters are accessible in some cells, implying that transcription factors could bind without activating the gene. Euchromatin and heterochromatin have very similar accessibilities, suggesting that transcription factors can penetrate heterochromatin. Thus, DNA accessibility is not likely to be the primary determinant of gene regulation.

H3K18Ac as a Marker of Cancer Progression and Potential Target of Anti-Cancer Therapy

Acetylation and deacetylation are posttranslational modifications (PTMs) which affect the regulation of chromatin structure and its remodeling. Acetylation of histone 3 at lysine placed on position 18 (H3K18Ac) plays an important role in driving progression of many types of cancer, including breast, colon, lung, hepatocellular, pancreatic, prostate, and thyroid cancer. The aim of this review is to analyze and discuss the newest findings regarding the role of H3K18Ac and acetylation of other histones in carcinogenesis. We summarize the level of H3K18Ac in different cancer cell lines and analyze its association with patients’ outcomes, including overall survival (OS), progression-free survival (PFS), and disease-free survival (DFS). Finally, we describe future perspectives of cancer therapeutic strategies based on H3K18 modifications.