Detection of gene loops by 3C in yeast

“Chromosome conformation capture” (3C) is a powerful method to detect physical interaction between any two genomic loci. 3C involves formaldehyde crosslinking to stabilize transient interactions, followed by restriction digestion, ligation and locus-specific PCR. Accordingly, 3C reveals complex three-dimensional interactions between distal genetic elements within intact cells at high resolution. Here, we describe a modified 3C protocol designed for detection of transient chromatin interactions in the yeast Saccharomyces cerevisiae. Using this protocol, we are able to detect juxtaposition of promoter and terminator regions of genes with ORFs as short as 1 kb in length. We anticipate that this method will be generally applicable to detect dynamic, short-range chromatin interactions and will facilitate the characterization of gene loops and their functional consequences.     

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.     

HiCORE: Hi-C Analysis for Identification of Core Chromatin Looping Regions with Higher Resolution

Genome-wide chromosome conformation capture (3C)-based high-throughput sequencing (Hi-C) has enabled identification of genome-wide chromatin loops. Because the Hi-C map with restriction fragment resolution is intrinsically associated with sparsity and stochastic noise, Hi-C data are usually binned at particular intervals; however, the binning method has limited reliability, especially at high resolution. Here, we describe a new method called HiCORE, which provides simple pipelines and algorithms to overcome the limitations of single-layered binning and predict core chromatin regions with three-dimensional physical interactions. In this approach, multiple layers of binning with slightly shifted genome coverage are generated, and interacting bins at each layer are integrated to infer narrower regions of chromatin interactions. HiCORE predicts chromatin looping regions with higher resolution, both in human and Arabidopsis genomes, and contributes to the identification of the precise positions of potential genomic elements in an unbiased manner.     

Invariant TAD Boundaries Constrain Cell-Type-Specific Looping Interactions between Promoters and Distal Elements around the CFTR Locus

Three-dimensional genome structure plays an important role in gene regulation. Globally, chromosomes are organized into active and inactive compartments while, at the gene level, looping interactions connect promoters to regulatory elements. Topologically associating domains (TADs), typically several hundred kilobases in size, form an intermediate level of organization. Major questions include how TADs are formed and how they are related to looping interactions between genes and regulatory elements. Here we performed a focused 5C analysis of a 2.8 Mb chromosome 7 region surrounding CFTR in a panel of cell types. We find that the same TAD boundaries are present in all cell types, indicating that TADs represent a universal chromosome architecture. Furthermore, we find that these TAD boundaries are present irrespective of the expression and looping of genes located between them. In contrast, looping interactions between promoters and regulatory elements are cell-type specific and occur mostly within TADs. This is exemplified by the CFTR promoter that in different cell types interacts with distinct sets of distal cell-type-specific regulatory elements that are all located within the same TAD. Finally, we find that long-range associations between loci located in different TADs are also detected, but these display much lower interaction frequencies than looping interactions within TADs. Interestingly, interactions between TADs are also highly cell-type-specific and often involve loci clustered around TAD boundaries. These data point to key roles of invariant TAD boundaries in constraining as well as mediating cell-type-specific long-range interactions and gene regulation.     

Microrheology for Hi-C Data Reveals the Spectrum of the Dynamic 3D Genome Organization

The one-dimensional information of genomic DNA is hierarchically packed inside the eukaryotic cell nucleus and organized in a three-dimensional (3D) space. Genome-wide chromosome conformation capture (Hi-C) methods have uncovered the 3D genome organization and revealed multiscale chromatin domains of compartments and topologically associating domains (TADs). Moreover, single-nucleosome live-cell imaging experiments have revealed the dynamic organization of chromatin domains caused by stochastic thermal fluctuations. However, the mechanism underlying the dynamic regulation of such hierarchical and structural chromatin units within the microscale thermal medium remains unclear. Microrheology is a way to measure dynamic viscoelastic properties coupling between thermal microenvironment and mechanical response. Here, we propose a new, to our knowledge, microrheology for Hi-C data to analyze the dynamic compliance property as a measure of rigidness and flexibility of genomic regions along with the time evolution. Our method allows the conversion of an Hi-C matrix into the spectrum of the dynamic rheological property along the genomic coordinate of a single chromosome. To demonstrate the power of the technique, we analyzed Hi-C data during the neural differentiation of mouse embryonic stem cells. We found that TAD boundaries behave as more rigid nodes than the intra-TAD regions. The spectrum clearly shows the dynamic viscoelasticity of chromatin domain formation at different timescales. Furthermore, we characterized the appearance of synchronous and liquid-like intercompartment interactions in differentiated cells. Together, our microrheology data derived from Hi-C data provide physical insights into the dynamics of the 3D genome organization.     

HPV-CCDC106 integration alters local chromosome architecture and hijacks an enhancer by three-dimensional genome structure remodeling in cervical cancer

Integration of human papillomavirus (HPV) DNA into the human genome is a reputed key driver of cervical cancer. However, the effects of HPV integration on chromatin structural organization and gene expression are largely unknown. We studied a cohort of 61 samples and identified an integration hot spot in the CCDC106 gene on chromosome 19. We then selected fresh cancer tissue that contained the unique integration loci at CCDC106 with no HPV episomal DNA and performed whole-genome, RNA, chromatin immunoprecipitation and high-throughput chromosome conformation capture (Hi-C) sequencing to identify the mechanisms of HPV integration in cervical carcinogenesis. Molecular analyses indicated that chromosome 19 exhibited significant genomic variation and differential expression densities, with correlation found between three-dimensional (3D) structural change and gene expression. Importantly, HPV integration divided one topologically associated domain (TAD) into two smaller TADs and hijacked an enhancer from PEG3 to CCDC106, with a decrease in PEG3 expression and an increase in CCDC106 expression. This expression dysregulation was further confirmed using 10 samples from our cohort, which exhibited the same HPV-CCDC106 integration. In summary, we found that HPV-CCDC106 integration altered local chromosome architecture and hijacked an enhancer via 3D genome structure remodeling. Thus, this study provides insight into the 3D structural mechanism underlying HPV integration in cervical carcinogenesis.     

Long-range heterochromatin association is mediated by silencing and double-strand DNA break repair proteins

The eukaryotic genome is highly organized in the nucleus, and this organization affects various nuclear processes. However, the molecular details of higher-order organization of chromatin remain obscure. In the present study, we show that the Saccharomyces cerevisiae silenced loci HML and HMR cluster in three-dimensional space throughout the cell cycle and independently of the telomeres. Long-range HML–HMR interactions require the homologous recombination (HR) repair pathway and phosphorylated H2A (γ-H2A). γ-H2A is constitutively present at silenced loci in unperturbed cells, its localization requires heterochromatin, and it is restricted to the silenced domain by the transfer DNA boundary element. SMC proteins and Scc2 localize to the silenced domain, and Scc2 binding requires the presence of γ-H2A. These findings illustrate a novel pathway for heterochromatin organization and suggest a role for HR repair proteins in genomic organization.     

Capturing native interactions: intrinsic methods to study chromatin conformation

The 3D organization of chromatin controls gene expression through spatial interactions between genomic loci. FISH and 3C‐based methods that are commonly used to study chromatin organization utilize chemical crosslinking, a step that may introduce biases in detectable chromatin interactions. In their recent study, Papantonis and colleagues (Brant et al, 2016 ) developed alternative new methods of detecting chromatin contacts without the use of chemical crosslinking agents. These tools increase the resolution and confidence at which interactions can be identified, and may be informative for chromatin interaction dynamics.     

An improved,chromosome-level genome of the giant panda (Ailuropoda melanoleuca)

BackgroundThe iconic giant panda (Ailuropoda melanoleuca), as both a flagship and umbrella species endemic to China, is a world famous symbol for wildlife conservation. The giant panda has several specific biological traits and holds a relatively small place in evolution. A high-quality genome of the giant panda is key to understanding the biology of this vulnerable species.FindingsWe generated a 2.48-Gb chromosome-level genome (GPv1) of the giant panda named “Jing Jing” with a contig N50 of 28.56 Mb and scaffold N50 of 134.17 Mb, respectively. The total length of chromosomes (n = 21) was 2.39-Gb, accounting for 96.4% of the whole genome. Compared with the previously published four genomes of the giant panda, our genome is characterized by the highest completeness and the correct sequence orientation. A gap-free and 850 kb length of immunoglobulin heavy-chain gene cluster was manually annotated in close proximity to the telomere of chromosome 14. Additionally, we developed an algorithm to predict the centromere position of each chromosome. We also constructed a complete chromatin structure for “Jing Jing”, which includes inter-chromosome interaction pattern, A/B compartment, topologically associated domain (TAD), TAD-clique and promoter-enhancer interaction (PEI).ConclusionsWe presented an improved chromosome-level genome and complete chromatin structure for the giant panda. This is a valuable resource for the future genetic and genomic studies on giant panda.     

Multiplexed single-cell profiling of chromatin states at genomic loci by expansion microscopy

Proper regulation of genome architecture and activity is essential for the development and function of multicellular organisms. Histone modifications, acting in combination, specify these activity states at individual genomic loci. However, the methods used to study these modifications often require either a large number of cells or are limited to targeting one histone mark at a time. Here, we developed a new method called Single Cell Evaluation of Post-TRanslational Epigenetic Encoding (SCEPTRE) that uses Expansion Microscopy (ExM) to visualize and quantify multiple histone modifications at non-repetitive genomic regions in single cells at a spatial resolution of ∼75 nm. Using SCEPTRE, we distinguished multiple histone modifications at a single housekeeping gene, quantified histone modification levels at multiple developmentally-regulated genes in individual cells, and evaluated the relationship between histone modifications and RNA polymerase II loading at individual loci. We find extensive variability in epigenetic states between individual gene loci hidden from current population-averaged measurements. These findings establish SCEPTRE as a new technique for multiplexed detection of combinatorial chromatin states at single genomic loci in single cells.     

Predicting chromatin organization using histone marks

Genome-wide mapping of three dimensional chromatin organization is an important yet technically challenging task. To aid experimental effort and to understand the determinants of long-range chromatin interactions, we have developed a computational model integrating Hi-C and histone mark ChIP-seq data to predict two important features of chromatin organization: chromatin interaction hubs and topologically associated domain (TAD) boundaries. Our model accurately and robustly predicts these features across datasets and cell types. Cell-type specific histone mark information is required for prediction of chromatin interaction hubs but not for TAD boundaries. Our predictions provide a useful guide for the exploration of chromatin organization.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0740-z) contains supplementary material, which is available to authorized users.     

Identifying multi-locus chromatin contacts in human cells using tethered multiple 3C

Background

Several recently developed experimental methods, each an extension of the chromatin conformation capture (3C) assay, have enabled the genome-wide profiling of chromatin contacts between pairs of genomic loci in 3D. Especially in complex eukaryotes, data generated by these methods, coupled with other genome-wide datasets, demonstrated that non-random chromatin folding correlates strongly with cellular processes such as gene expression and DNA replication.

Results

We describe a genome architecture assay, tethered multiple 3C (TM3C), that maps genome-wide chromatin contacts via a simple protocol of restriction enzyme digestion and religation of fragments upon agarose gel beads followed by paired-end sequencing. In addition to identifying contacts between pairs of loci, TM3C enables identification of contacts among more than two loci simultaneously. We use TM3C to assay the genome architectures of two human cell lines: KBM7, a near-haploid chronic leukemia cell line, and NHEK, a normal diploid human epidermal keratinocyte cell line. We confirm that the contact frequency maps produced by TM3C exhibit features characteristic of existing genome architecture datasets, including the expected scaling of contact probabilities with genomic distance, megabase scale chromosomal compartments and sub-megabase scale topological domains. We also confirm that TM3C captures several known cell type-specific contacts, ploidy shifts and translocations, such as Philadelphia chromosome formation (Ph+) in KBM7. We confirm a subset of the triple contacts involving the IGF2-H19 imprinting control region (ICR) using PCR analysis for KBM7 cells. Our genome-wide analysis of pairwise and triple contacts demonstrates their preference for linking open chromatin regions to each other and for linking regions with higher numbers of DNase hypersensitive sites (DHSs) to each other. For near-haploid KBM7 cells, we infer whole genome 3D models that exhibit clustering of small chromosomes with each other and large chromosomes with each other, consistent with previous studies of the genome architectures of other human cell lines.

Conclusion

TM3C is a simple protocol for ascertaining genome architecture and can be used to identify simultaneous contacts among three or four loci. Application of TM3C to a near-haploid human cell line revealed large-scale features of chromosomal organization and multi-way chromatin contacts that preferentially link regions of open chromatin.

Electronic supplementary material

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

r3Cseq: an R/Bioconductor package for the discovery of long-range genomic interactions from chromosome conformation capture and next-generation sequencing data

The coupling of chromosome conformation capture (3C) with next-generation sequencing technologies enables the high-throughput detection of long-range genomic interactions, via the generation of ligation products between DNA sequences, which are closely juxtaposed in vivo. These interactions involve promoter regions, enhancers and other regulatory and structural elements of chromosomes and can reveal key details of the regulation of gene expression. 3C-seq is a variant of the method for the detection of interactions between one chosen genomic element (viewpoint) and the rest of the genome. We present r3Cseq, an R/Bioconductor package designed to perform 3C-seq data analysis in a number of different experimental designs. The package reads a common aligned read input format, provides data normalization, allows the visualization of candidate interaction regions and detects statistically significant chromatin interactions, thus greatly facilitating hypothesis generation and the interpretation of experimental results. We further demonstrate its use on a series of real-world applications.