Genome-wide mapping and cryo-EM structural analyses of the overlapping tri-nucleosome composed of hexasome-hexasome-octasome moieties

The nucleosome is a fundamental unit of chromatin in which about 150 base pairs of DNA are wrapped around a histone octamer. The overlapping di-nucleosome has been proposed as a product of chromatin remodeling around the transcription start site, and previously found as a chromatin unit, in which about 250 base pairs of DNA continuously bind to the histone core composed of a hexamer and an octamer. In the present study, our genome-wide analysis of human cells suggests another higher nucleosome stacking structure, the overlapping tri-nucleosome, which wraps about 300-350 base-pairs of DNA in the region downstream of certain transcription start sites of actively transcribed genes. We determine the cryo-electron microscopy (cryo-EM) structure of the overlapping tri-nucleosome, in which three subnucleosome moieties, hexasome, hexasome, and octasome, are associated by short connecting DNA segments. Small angle X-ray scattering and coarse-grained molecular dynamics simulation analyses reveal that the cryo-EM structure of the overlapping tri-nucleosome may reflect its structure in solution. Our findings suggest that nucleosome stacking structures composed of hexasome and octasome moieties may be formed by nucleosome remodeling factors around transcription start sites for gene regulation.

Modeling population size independent tissue epigenomes by ChIL-seq with single thin sections

Recent advances in genome-wide technologies have enabled analyses using small cell numbers of even single cells. However, obtaining tissue epigenomes with cell-type resolution from large organs and tissues still remains challenging, especially when the available material is limited. Here, we present a ChIL-based approach for analyzing the diverse cellular dynamics at the tissue level using high-depth epigenomic data. "ChIL for tissues" allows the analysis of a single tissue section and can reproducibly generate epigenomic profiles from several tissue types, based on the distribution of target epigenomic states, tissue morphology, and number of cells. The proposed method enabled the independent evaluation of changes in cell populations and gene activation in cells from regenerating skeletal muscle tissues, using a statistical model of RNA polymerase II distribution on gene loci. Thus, the integrative analyses performed using ChIL can elucidate in vivo cell-type dynamics of tissues.

Genomic Profiling by ALaP-Seq Reveals Transcriptional Regulation by PML Bodies through DNMT3A Exclusion

The promyelocytic leukemia (PML) body is a phase-separated nuclear structure physically associated with chromatin, implying its crucial roles in genome functions. However, its role in transcriptional regulation is largely unknown. We developed APEX-mediated chromatin labeling and purification (ALaP) to identify the genomic regions proximal to PML bodies. We found that PML bodies associate with active regulatory regions across the genome and with ∼300 kb of the short arm of the Y chromosome (YS300) in mouse embryonic stem cells. The PML body association with YS300 is essential for the transcriptional activity of the neighboring Y-linked clustered genes. Mechanistically, PML bodies provide specific nuclear spaces that the de novo DNA methyltransferase DNMT3A cannot access, resulting in the steady maintenance of a hypo-methylated state at Y-linked gene promoters. Our study underscores a new mechanism for gene regulation in the 3D nuclear space and provides insights into the functional properties of nuclear structures for genome function.

Integrative comparison of the genomic and transcriptomic landscape between prostate cancer patients of predominantly African or European genetic ancestry

Men of predominantly African Ancestry (AA) have higher prostate cancer (CaP) incidence and worse survival than men of predominantly European Ancestry (EA). While socioeconomic factors drive this disparity, genomic factors may also contribute to differences in the incidence and mortality rates. To compare the prevalence of prostate tumor genomic alterations and transcriptomic profiles by patient genetic ancestry, we evaluated genomic profiles from The Cancer Genome Atlas (TCGA) CaP cohort (n = 498). Patient global and local genetic ancestry were estimated by computational algorithms using genotyping data; 414 (83.1%) were EA, 61 (12.2%) were AA, 11 (2.2%) were East Asian Ancestry (EAA), 10 (2.0%) were Native American (NA), and 2 (0.4%) were other ancestry. Genetic ancestry was highly concordant with self-identified race/ethnicity. Subsequent analyses were limited to 61 AA and 414 EA cases. Significant differences were observed by ancestry in the frequency of SPOP mutations (20.3% AA vs. 10.0% EA; p = 5.6×10-03), TMPRSS2-ERG fusions (29.3% AA vs. 39.6% EA; p = 4.4×10-02), and PTEN deletions/losses (11.5% AA vs. 30.2% EA; p = 3.5×10-03). Differentially expressed genes (DEGs) between AAs and EAs showed significant enrichment for prostate eQTL target genes (p = 8.09×10-48). Enrichment of highly expressed DEGs for immune-related pathways was observed in AAs, and for PTEN/PI3K signaling in EAs. Nearly one-third of DEGs (31.3%) were long non-coding RNAs (DE-lncRNAs). The proportion of DE-lncRNAs with higher expression in AAs greatly exceeded that with lower expression in AAs (p = 1.2×10-125). Both ChIP-seq and RNA-seq data suggested a stronger regulatory role for AR signaling pathways in DE-lncRNAs vs. non-DE-lncRNAs. CaP-related oncogenic lncRNAs, such as PVT1, PCAT1 and PCAT10/CTBP1-AS, were found to be more highly expressed in AAs. We report substantial heterogeneity in the prostate tumor genome and transcriptome between EA and AA. These differences may be biological contributors to racial disparities in CaP incidence and outcomes.

Genome-wide analysis of the spatiotemporal regulation of firing and dormant replication origins in human cells

In metazoan cells, only a limited number of mini chromosome maintenance (MCM) complexes are fired during S phase, while the majority remain dormant. Several methods have been used to map replication origins, but such methods cannot identify dormant origins. Herein, we determined MCM7-binding sites in human cells using ChIP-Seq, classified them into firing and dormant origins using origin data and analysed their association with various chromatin signatures. Firing origins, but not dormant origins, were well correlated with open chromatin regions and were enriched upstream of transcription start sites (TSSs) of transcribed genes. Aggregation plots of MCM7 signals revealed minimal difference in the efficacy of MCM loading between firing and dormant origins. We also analysed common fragile sites (CFSs) and found a low density of origins at these sites. Nevertheless, firing origins were enriched upstream of the TSSs. Based on the results, we propose a model in which excessive MCMs are actively loaded in a genome-wide manner, irrespective of chromatin status, but only a fraction are passively fired in chromatin areas with an accessible open structure, such as regions upstream of TSSs of transcribed genes. This plasticity in the specification of replication origins may minimize collisions between replication and transcription.

MobiSeq: De novo SNP discovery in model and non-model species through sequencing the flanking region of transposable elements

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions.

Chd2 regulates chromatin for proper gene expression toward differentiation in mouse embryonic stem cells

Chromatin reorganization is necessary for pluripotent stem cells, including embryonic stem cells (ESCs), to acquire lineage potential. However, it remains unclear how ESCs maintain their characteristic chromatin state for appropriate gene expression upon differentiation. Here, we demonstrate that chromodomain helicase DNA-binding domain 2 (Chd2) is required to maintain the differentiation potential of mouse ESCs. Chd2-depleted ESCs showed suppressed expression of developmentally regulated genes upon differentiation and subsequent differentiation defects without affecting gene expression in the undifferentiated state. Furthermore, chromatin immunoprecipitation followed by sequencing revealed alterations in the nucleosome occupancy of the histone variant H3.3 for developmentally regulated genes in Chd2-depleted ESCs, which in turn led to elevated trimethylation of the histone H3 lysine 27. These results suggest that Chd2 is essential in preventing suppressive chromatin formation for developmentally regulated genes and determines subsequent effects on developmental processes in the undifferentiated state.

Deciphering the genomic targets of alkylating polyamide conjugates using high-throughput sequencing

Chemically engineered small molecules targeting specific genomic sequences play an important role in drug development research. Pyrrole-imidazole polyamides (PIPs) are a group of molecules that can bind to the DNA minor-groove and can be engineered to target specific sequences. Their biological effects rely primarily on their selective DNA binding. However, the binding mechanism of PIPs at the chromatinized genome level is poorly understood. Herein, we report a method using high-throughput sequencing to identify the DNA-alkylating sites of PIP-indole-seco-CBI conjugates. High-throughput sequencing analysis of conjugate 2: showed highly similar DNA-alkylating sites on synthetic oligos (histone-free DNA) and on human genomes (chromatinized DNA context). To our knowledge, this is the first report identifying alkylation sites across genomic DNA by alkylating PIP conjugates using high-throughput sequencing.

Exploration of nucleosome positioning patterns in transcription factor function

The binding of transcription factors (TFs) triggers activation of specific chromatin regions through the recruitment and activation of RNA polymerase. Unique nucleosome positioning (NP) occurs during gene expression and has been suggested to be involved in various other chromatin functions. However, the diversity of NP that can occur for each function has not been clarified. Here we used MNase-Seq data to evaluate NP around 258 cis-regulatory elements in the mouse genome. Principal component analysis of the 258 elements revealed that NP consisted of five major patterns. Furthermore, the five NP patterns had predictive power for the level of gene expression. We also demonstrated that selective NP patterns appeared around TF binding sites. These results suggest that the NP patterns are correlated to specific functions on chromatin.

Chd5 Regulates MuERV-L/MERVL Expression in Mouse Embryonic Stem Cells Via H3K27me3 Modification and Histone H3.1/H3.2

Chd5 is an essential factor for neuronal differentiation and spermatogenesis and is a known tumor suppressor. H3K27me3 and H3K4un are modifications recognized by Chd5; however, it remains unclear how Chd5 remodels chromatin structure. We completely disrupted the Chd5 locus using the CRISPR-Cas9 system to generate a 52 kbp long deletion and analyzed Chd5 function in mouse embryonic stem cells. Our findings show that Chd5 represses murine endogenous retrovirus-L (MuERV-L/MERVL), an endogenous retrovirus-derived retrotransposon, by regulating H3K27me3 and H3.1/H3.2 function.