The transcriptional interactome: gene expression in 3D

Chromatin Landscape Governing Murine Epidermal Differentiation

Chromatin landscape and regulatory networks are determinants in lineage specification and differentiation. To define the temporospatial differentiation axis in murine epidermal cells in vivo, we generated datasets profiling expression dynamics (RNA sequencing), chromatin accessibility (assay for transposase-accessible chromatin using sequencing), architecture (Hi-C), and histone modifications (chromatin immunoprecipitation followed by sequencing) in the epidermis. We show that many differentially regulated genes are suppressed during the differentiation process, with superenhancers controlling differentiation-specific epigenomic changes. Our data shows the relevance of the Dlx/Klf/Grhl combinatorial regulatory network in maintaining correct temporospatial gene expression during epidermal differentiation. We determined differential open compartments, topologically associating domain score, and looping in the basal cell and suprabasal cell epidermal fractions, with the evolutionarily conserved epidermal differentiation complex region showing distinct suprabasal cell-specific topologically associating domain and loop formation that coincided with superenhancer sites. Overall, our study provides a global genome-wide resource of chromatin dynamics that define unrecognized regulatory networks and the epigenetic control of Dlx3-bound superenhancer elements during epidermal differentiation.

Juvenile idiopathic arthritis-associated genetic loci exhibit spatially constrained gene regulatory effects across multiple tissues and immune cell types

Juvenile idiopathic arthritis (JIA) is an autoimmune, inflammatory joint disease with complex genetic etiology. Previous GWAS have found many genetic loci associated with JIA. However, the biological mechanism behind JIA remains unknown mainly because most risk loci are located in non-coding genetic regions. Interestingly, increasing evidence has found that regulatory elements in the non-coding regions can regulate the expression of distant target genes through spatial (physical) interactions. Here, we used information on the 3D genome organization (Hi-C data) to identify target genes that physically interact with SNPs within JIA risk loci. Subsequent analysis of these SNP-gene pairs using data from tissue and immune cell type-specific expression quantitative trait loci (eQTL) databases allowed the identification of risk loci that regulate the expression of their target genes. In total, we identified 59 JIA-risk loci that regulate the expression of 210 target genes across diverse tissues and immune cell types. Functional annotation of spatial eQTLs within JIA risk loci identified significant overlap with gene regulatory elements (i.e., enhancers and transcription factor binding sites). We found target genes involved in immune-related pathways such as antigen processing and presentation (e.g., ERAP2, HLA class I and II), the release of pro-inflammatory cytokines (e.g., LTBR, TYK2), proliferation and differentiation of specific immune cell types (e.g., AURKA in Th17 cells), and genes involved in physiological mechanisms related to pathological joint inflammation (e.g., LRG1 in arteries). Notably, many of the tissues where JIA-risk loci act as spatial eQTLs are not classically considered central to JIA pathology. Overall, our findings highlight the potential tissue and immune cell type-specific regulatory changes contributing to JIA pathogenesis. Future integration of our data with clinical studies can contribute to the development of improved JIA therapy.

Genome Tectonics: Linking Dynamic Genome Organization with Cellular Nutrients

BACKGROUND

Our daily intake of food provides nutrients for the maintenance of health, growth, and development. The field of nutrigenomics aims to link dietary intake/nutrients to changes in epigenetic status and gene expression.

SUMMARY

Although the relationship between our diet and our genes in under intense investigation, there is still a significant aspect of our genome that has received little attention with regard to this. In the past 15 years, the importance of genome organization has become increasingly evident, with research identifying small-scale local changes to large segments of the genome dynamically repositioning within the nucleus in response to/or mediating change in gene expression. The discovery of these dynamic processes and organization maybe as significant as dynamic plate tectonics is to geology, there is little information tying genome organization to specific nutrients or dietary intake.

KEY MESSAGES

Here, we detail key principles of genome organization and structure, with emphasis on genome folding and organization, and link how these contribute to our future understand of nutrigenomics.

Wheat Breeding, Transcription Factories, and Genetic Interactions: New Perspectives

Epistatic interactions and negative heterosis have been shown to be associated with interchromosomal interactions in wheat. Physical gene-gene interactions between co-regulated genes clustered in "transcription factories" have been documented, and a genome-wide atlas of functionally paired, interacting regulatory elements and genes of wheat recently produced. Integration of these new studies on gene and regulatory element interactions, co-regulation of gene expression in "transcription factories," and epigenetics generates new perspectives for wheat breeding and trait enhancement.

Parallel analysis of transcription, integration, and sequence of single HIV-1 proviruses

HIV-1-infected cells that persist despite antiretroviral therapy (ART) are frequently considered "transcriptionally silent," but active viral gene expression may occur in some cells, challenging the concept of viral latency. Applying an assay for profiling the transcriptional activity and the chromosomal locations of individual proviruses, we describe a global genomic and epigenetic map of transcriptionally active and silent proviral species and evaluate their longitudinal evolution in persons receiving suppressive ART. Using genome-wide epigenetic reference data, we show that proviral transcriptional activity is associated with activating epigenetic chromatin features in linear proximity of integration sites and in their inter- and intrachromosomal contact regions. Transcriptionally active proviruses were actively selected against during prolonged ART; however, this pattern was violated by large clones of virally infected cells that may outcompete negative selection forces through elevated intrinsic proliferative activity. Our results suggest that transcriptionally active proviruses are dynamically evolving under selection pressure by host factors.

A holistic view of mouse enhancer architectures reveals analogous pleiotropic effects and correlation with human disease

BACKGROUND

Efforts to elucidate the function of enhancers in vivo are underway but their vast numbers alongside differing enhancer architectures make it difficult to determine their impact on gene activity. By systematically annotating multiple mouse tissues with super- and typical-enhancers, we have explored their relationship with gene function and phenotype.

RESULTS

Though super-enhancers drive high total- and tissue-specific expression of their associated genes, we find that typical-enhancers also contribute heavily to the tissue-specific expression landscape on account of their large numbers in the genome. Unexpectedly, we demonstrate that both enhancer types are preferentially associated with relevant 'tissue-type' phenotypes and exhibit no difference in phenotype effect size or pleiotropy. Modelling regulatory data alongside molecular data, we built a predictive model to infer gene-phenotype associations and use this model to predict potentially novel disease-associated genes.

CONCLUSION

Overall our findings reveal that differing enhancer architectures have a similar impact on mammalian phenotypes whilst harbouring differing cellular and expression effects. Together, our results systematically characterise enhancers with predicted phenotypic traits endorsing the role for both types of enhancers in human disease and disorders.

The corrected gene proximity map for analyzing the 3D genome organization using Hi-C data

BACKGROUND

Genome-wide ligation-based assays such as Hi-C provide us with an unprecedented opportunity to investigate the spatial organization of the genome. Results of a typical Hi-C experiment are often summarized in a chromosomal contact map, a matrix whose elements reflect the co-location frequencies of genomic loci. To elucidate the complex structural and functional interactions between those genomic loci, networks offer a natural and powerful framework.

RESULTS

We propose a novel graph-theoretical framework, the Corrected Gene Proximity (CGP) map to study the effect of the 3D spatial organization of genes in transcriptional regulation. The starting point of the CGP map is a weighted network, the gene proximity map, whose weights are based on the contact frequencies between genes extracted from genome-wide Hi-C data. We derive a null model for the network based on the signal contributed by the 1D genomic distance and use it to "correct" the gene proximity for cell type 3D specific arrangements. The CGP map, therefore, provides a network framework for the 3D structure of the genome on a global scale. On human cell lines, we show that the CGP map can detect and quantify gene co-regulation and co-localization more effectively than the map obtained by raw contact frequencies. Analyzing the expression pattern of metabolic pathways of two hematopoietic cell lines, we find that the relative positioning of the genes, as captured and quantified by the CGP, is highly correlated with their expression change. We further show that the CGP map can be used to form an inter-chromosomal proximity map that allows large-scale abnormalities, such as chromosomal translocations, to be identified.

CONCLUSIONS

The Corrected Gene Proximity map is a map of the 3D structure of the genome on a global scale. It allows the simultaneous analysis of intra- and inter- chromosomal interactions and of gene co-regulation and co-localization more effectively than the map obtained by raw contact frequencies, thus revealing hidden associations between global spatial positioning and gene expression. The flexible graph-based formalism of the CGP map can be easily generalized to study any existing Hi-C datasets.

Wheat chromatin architecture is organized in genome territories and transcription factories

BACKGROUND

Polyploidy is ubiquitous in eukaryotic plant and fungal lineages, and it leads to the co-existence of several copies of similar or related genomes in one nucleus. In plants, polyploidy is considered a major factor in successful domestication. However, polyploidy challenges chromosome folding architecture in the nucleus to establish functional structures.

RESULTS

We examine the hexaploid wheat nuclear architecture by integrating RNA-seq, ChIP-seq, ATAC-seq, Hi-C, and Hi-ChIP data. Our results highlight the presence of three levels of large-scale spatial organization: the arrangement into genome territories, the diametrical separation between facultative and constitutive heterochromatin, and the organization of RNA polymerase II around transcription factories. We demonstrate the micro-compartmentalization of transcriptionally active genes determined by physical interactions between genes with specific euchromatic histone modifications. Both intra- and interchromosomal RNA polymerase-associated contacts involve multiple genes displaying similar expression levels.

CONCLUSIONS

Our results provide new insights into the physical chromosome organization of a polyploid genome, as well as on the relationship between epigenetic marks and chromosome conformation to determine a 3D spatial organization of gene expression, a key factor governing gene transcription in polyploids.

TADCompare: An R Package for Differential and Temporal Analysis of Topologically Associated Domains

Recent research using chromatin conformation capture technologies, such as Hi-C, has demonstrated the importance of topologically associated domains (TADs) and smaller chromatin loops, collectively referred hereafter as "interacting domains." Many such domains change during development or disease, and exhibit cell- and condition-specific differences. Quantification of the dynamic behavior of interacting domains will help to better understand genome regulation. Methods for comparing interacting domains between cells and conditions are highly limited. We developed TADCompare, a method for differential analysis of boundaries of interacting domains between two or more Hi-C datasets. TADCompare is based on a spectral clustering-derived measure called the eigenvector gap, which enables a loci-by-loci comparison of boundary differences. Using this measure, we introduce methods for identifying differential and consensus boundaries of interacting domains and tracking boundary changes over time. We further propose a novel framework for the systematic classification of boundary changes. Colocalization- and gene enrichment analysis of different types of boundary changes demonstrated distinct biological functionality associated with them. TADCompare is available on https://github.com/dozmorovlab/TADCompare and Bioconductor (submitted).

Topological structure analysis of chromatin interaction networks

BACKGROUND

Current Hi-C technologies for chromosome conformation capture allow to understand a broad spectrum of functional interactions between genome elements. Although significant progress has been made into analysis of Hi-C data to identify biologically significant features, many questions still remain open, in particular regarding potential biological significance of various topological features that are characteristic for chromatin interaction networks.

RESULTS

It has been previously observed that promoter capture Hi-C (PCHi-C) interaction networks tend to separate easily into well-defined connected components that can be related to certain biological functionality, however, such evidence was based on manual analysis and was limited. Here we present a novel method for analysis of chromatin interaction networks aimed towards identifying characteristic topological features of interaction graphs and confirming their potential significance in chromatin architecture. Our method automatically identifies all connected components with an assigned significance score above a given threshold. These components can be subjected afterwards to different assessment methods for their biological role and/or significance. The method was applied to the largest PCHi-C data set available to date that contains interactions for 17 haematopoietic cell types. The results demonstrate strong evidence of well-pronounced component structure of chromatin interaction networks and provide some characterisation of this component structure. We also performed an indicative assessment of potential biological significance of identified network components with the results confirming that the network components can be related to specific biological functionality.

CONCLUSIONS

The obtained results show that the topological structure of chromatin interaction networks can be well described in terms of isolated connected components of the network and that formation of these components can be often explained by biological features of functionally related gene modules. The presented method allows automatic identification of all such components and evaluation of their significance in PCHi-C dataset for 17 haematopoietic cell types. The method can be adapted for exploration of other chromatin interaction data sets that include information about sufficiently large number of different cell types, and, in principle, also for analysis of other kinds of cell type-specific networks.

Nucleome Dynamics during Retinal Development

More than 8,000 genes are turned on or off as progenitor cells produce the 7 classes of retinal cell types during development. Thousands of enhancers are also active in the developing retinae, many having features of cell- and developmental stage-specific activity. We studied dynamic changes in the 3D chromatin landscape important for precisely orchestrated changes in gene expression during retinal development by ultra-deep in situ Hi-C analysis on murine retinae. We identified developmental-stage-specific changes in chromatin compartments and enhancer-promoter interactions. We developed a machine learning-based algorithm to map euchromatin and heterochromatin domains genome-wide and overlaid it with chromatin compartments identified by Hi-C. Single-cell ATAC-seq and RNA-seq were integrated with our Hi-C and previous ChIP-seq data to identify cell- and developmental-stage-specific super-enhancers (SEs). We identified a bipolar neuron-specific core regulatory circuit SE upstream of Vsx2, whose deletion in mice led to the loss of bipolar neurons.

Discontinuous transcription

Numerous studies based on new single-cell and single-gene techniques show that individual genes can be transcribed in short bursts or pulses accompanied by changes in pulsing frequencies. Since so many examples of such discontinuous or fluctuating transcription have been found from prokaryotes to mammals, it now seems to be a common mode of gene expression. In this review we discuss the occurrence of the transcriptional fluctuations, the techniques used for their detection, their putative causes, kinetic characteristics, and probable physiological significance.

Genes responsive to rapamycin and serum deprivation are clustered on chromosomes and undergo reorganization within local chromatin environments

We previously demonstrated that genome reorganization, through chromosome territory repositioning, occurs concurrently with significant changes in gene expression in normal primary human fibroblasts treated with the drug rapamycin, or stimulated into quiescence. Although these events occurred concomitantly, it is unclear how specific changes in gene expression relate to reorganization of the genome at higher resolution. We used computational analyses, genome organization assays, and microscopy, to investigate the relationship between chromosome territory positioning and gene expression. We determined that despite relocation of chromosome territories, there was no substantial bias in the proportion of genes changing expression on any one chromosome, including chromosomes 10 and 18. Computational analyses identified that clusters of serum deprivation and rapamycin-responsive genes along the linear extent of chromosomes. Chromosome conformation capture (3C) analysis demonstrated the strengthening or loss of specific long-range chromatin interactions in response to rapamycin and quiescence induction, including a cluster of genes containing Interleukin-8 and several chemokine genes on chromosome 4. We further observed that the LIF gene, which is highly induced upon rapamycin treatment, strengthened interactions with up- and down-stream intergenic regions. Our findings indicate that the repositioning of chromosome territories in response to cell stimuli, this does not reflect gene expression changes occurring within physically clustered groups of genes.

Integrative transcriptome sequencing reveals extensive alternative trans-splicing and cis-backsplicing in human cells

Transcriptionally non-co-linear (NCL) transcripts can originate from trans-splicing (trans-spliced RNA; 'tsRNA') or cis-backsplicing (circular RNA; 'circRNA'). While numerous circRNAs have been detected in various species, tsRNAs remain largely uninvestigated. Here, we utilize integrative transcriptome sequencing of poly(A)- and non-poly(A)-selected RNA-seq data from diverse human cell lines to distinguish between tsRNAs and circRNAs. We identified 24,498 NCL events and found that a considerable proportion (20-35%) of them arise from both tsRNAs and circRNAs, representing extensive alternative trans-splicing and cis-backsplicing in human cells. We show that sequence generalities of exon circularization are also observed in tsRNAs. Recapitulation of NCL RNAs further shows that inverted Alu repeats can simultaneously promote the formation of tsRNAs and circRNAs. However, tsRNAs and circRNAs exhibit quite different, or even opposite, expression patterns, in terms of correlation with the expression of their co-linear counterparts, expression breadth/abundance, transcript stability, and subcellular localization preference. These results indicate that tsRNAs and circRNAs may play different regulatory roles and analysis of NCL events should take the joint effects of different NCL-splicing types and joint effects of multiple NCL events into consideration. This study describes the first transcriptome-wide analysis of trans-splicing and cis-backsplicing, expanding our understanding of the complexity of the human transcriptome.

Towards a Better Classification and Novel Therapies Based on the Genetics of Systemic Sclerosis

PURPOSE OF THE REVIEW

Nowadays, important advances have occurred in our understanding of the pathogenesis of systemic sclerosis (SSc), which is a rare immune-mediated inflammatory disease (IMID) characterized by vascular damage, immune imbalance, and fibrosis. Its etiology remains unknown; nevertheless, both environmental and genetic factors play a major role in the disease. This review will focus on the main advances made in the field of genetics of SSc.

RECENT FINDINGS

The assessment of how interindividual genetic variability affects disease onset and progression has enhanced our knowledge of disease biology, and this will eventually translate in the development of new diagnostic and therapeutic tools, which is the final goal of personalized medicine. We will provide an overview of the most relevant achievements in the genetics of SSc, its shared genetics among IMIDs with special attention on drug repurposing, current challenges for the functional characterization of risk variants, and future directions.

FOXO3 on the Road to Longevity: Lessons From SNPs and Chromatin Hubs

Health span is driven by a precise interplay between genes and the environment. Cell response to environmental cues is mediated by signaling cascades and genetic variants that affect gene expression by regulating chromatin plasticity. Indeed, they can promote the interaction of promoters with regulatory elements by forming active chromatin hubs.

FOXO3 encodes a transcription factor with a strong impact on aging and age-related phenotypes, as it regulates stress response, therefore affecting lifespan. A significant association has been shown between human longevity and several FOXO3 variants located in intron 2. This haplotype block forms a putative aging chromatin hub in which FOXO3 has a central role, as it modulates the physical connection and activity of neighboring genes involved in age-related processes.

Here we describe the role of FOXO3 and its single-nucleotide polymorphisms (SNPs) in healthy aging, with a focus on the enhancer region encompassing the SNP rs2802292, which upregulates FOXO3 expression and can promote the activity of the aging hub in response to different stress stimuli. FOXO3 protective effect on lifespan may be due to the accessibility of this region to transcription factors promoting its expression. This could in part explain the differences in FOXO3 association with longevity between genders, as its activity in females may be modulated by estrogens through estrogen receptor response elements located in the rs2802292-encompassing region. Altogether, the molecular mechanisms described here may help establish whether the rs2802292 SNP can be taken advantage of in predictive medicine and define the potential of targeting FOXO3 for age-related diseases.

The Third Intron of IRF8 Is a Cell-Type-Specific Chromatin Priming Element during Mouse Embryonal Stem Cell Differentiation

Interferon regulatory factor 8 (IRF8) is a nuclear transcription factor that plays a key role in the hierarchical differentiation of hematopoietic stem cells toward monocyte/dendritic cell lineages. Therefore, its expression is mainly limited to bone marrow-derived cells. The molecular mechanisms governing this cell-type-restricted expression have been described. However, the molecular mechanisms that are responsible for its silencing in non-hematopoietic cells are elusive. Recently, we demonstrated a role for IRF8 third intron in restricting its expression in non-hematopoietic cells. Furthermore, we showed that this intron alone is sufficient to promote repressed chromatin a cell-type-specific manner. Here we demonstrate the effect of the IRF8 third intron on chromatin conformation during murine embryonal stem cell differentiation. Using genome editing, we provide data showing that the third intron plays a key role in priming the chromatin state of the IRF8 locus during cell differentiation. It mediates dual regulatory effects in a cell-type-specific mode. It acts as a repressor element governing chromatin state of the IRF8 locus during embryonal stem cell differentiation to cardiomyocytes that are expression-restrictive cells. Conversely, it functions as an activator element that is essential for open chromatin structure during the differentiation of these cells to dendritic cells that are expression-permissive cells. Together, these results point to the role of IRF8 third intron as a cell-type-specific chromatin priming element during embryonal stem cell differentiation. These data add another layer to our understanding of the molecular mechanisms governing misexpression of a cell-type-specific gene such as IRF8.

NCLcomparator: systematically post-screening non-co-linear transcripts (circular, trans-spliced, or fusion RNAs) identified from various detectors

BACKGROUND

Non-co-linear (NCL) transcripts consist of exonic sequences that are topologically inconsistent with the reference genome in an intragenic fashion (circular or intragenic trans-spliced RNAs) or in an intergenic fashion (fusion or intergenic trans-spliced RNAs). On the basis of RNA-seq data, numerous NCL event detectors have been developed and detected thousands of NCL events in diverse species. However, there are great discrepancies in the identification results among detectors, indicating a considerable proportion of false positives in the detected NCL events. Although several helpful guidelines for evaluating the performance of NCL event detectors have been provided, a systematic guideline for measurement of NCL events identified by existing tools has not been available.

RESULTS

We develop a software, NCLcomparator, for systematically post-screening the intragenic or intergenic NCL events identified by various NCL detectors. NCLcomparator first examine whether the input NCL events are potentially false positives derived from ambiguous alignments (i.e., the NCL events have an alternative co-linear explanation or multiple matches against the reference genome). To evaluate the reliability of the identified NCL events, we define the NCL score (NCLscore) based on the variation in the number of supporting NCL junction reads identified by the tools examined. Of the input NCL events, we show that the ambiguous alignment-derived events have relatively lower NCLscore values than the other events, indicating that an NCL event with a higher NCLscore has a higher level of reliability. To help selecting highly expressed NCL events, NCLcomparator also provides a series of useful measurements such as the expression levels of the detected NCL events and their corresponding host genes and the junction usage of the co-linear splice junctions at both NCL donor and acceptor sites.

CONCLUSION

NCLcomparator provides useful guidelines, with the input of identified NCL events from various detectors and the corresponding paired-end RNA-seq data only, to help users selecting potentially high-confidence NCL events for further functional investigation. The software thus helps to facilitate future studies into NCL events, shedding light on the fundamental biology of this important but understudied class of transcripts. NCLcomparator is freely accessible at https://github.com/TreesLab/NCLcomparator .

Trans-allelic mutational effects at the Peg3 imprinted locus

How one allele interacts with the other for the function of a gene is not well understood. In this study, we tested potential allelic interaction at the Peg3 imprinted locus with several mutant alleles targeting an Imprinting Control Region, the Peg3-DMR. According to the results, maternal deletion of the Peg3-DMR resulted in 2-fold up-regulation of two paternally expressed genes, Peg3 and Usp29. These trans-allelic mutational effects were observed consistently throughout various tissues with different developmental stages. These effects were also associated mainly with the genetic manipulation of the Peg3-DMR, but not with the other genomic changes within the Peg3 locus. The observed trans-allelic effects were unidirectional with the maternal influencing the paternal allele, but not with the opposite direction. Overall, the observed mutational effects suggest the presence of previously unrecognized trans-allelic regulation associated with the Peg3-DMR.

Trans-spliced long non-coding RNA: an emerging regulator of pluripotency

With dual capacities for unlimited self-renewal and pluripotent differentiation, pluripotent stem cells (PSCs) give rise to many cell types in our body and PSC culture systems provide an unparalleled opportunity to study early human development and disease. Accumulating evidence indicates that the molecular mechanisms underlying pluripotency maintenance in PSCs involve many factors. Among these regulators, recent studies have shown that long non-coding RNAs (lncRNAs) can affect the pluripotency circuitry by cooperating with master pluripotency-associated factors. Additionally, trans-spliced RNAs, which are generated by combining two or more pre-mRNA transcripts to produce a chimeric RNA, have been identified as regulators of various biological processes, including human pluripotency. In this review, we summarize and discuss current knowledge about the roles of lncRNAs, including trans-spliced lncRNAs, in controlling pluripotency.

Genomic meta-analysis of the interplay between 3D chromatin organization and gene expression programs under basal and stress conditions

Background

Our appreciation of the critical role of the genome’s 3D organization in gene regulation is steadily increasing. Recent 3C-based deep sequencing techniques elucidated a hierarchy of structures that underlie the spatial organization of the genome in the nucleus. At the top of this hierarchical organization are chromosomal territories and the megabase-scale A/B compartments that correlate with transcriptional activity within cells. Below them are the relatively cell-type-invariant topologically associated domains (TADs), characterized by high frequency of physical contacts between loci within the same TAD, and are assumed to function as regulatory units. Within TADs, chromatin loops bring enhancers and target promoters to close spatial proximity. Yet, we still have only rudimentary understanding how differences in chromatin organization between different cell types affect cell-type-specific gene expression programs that are executed under basal and challenged conditions.

Results

Here, we carried out a large-scale meta-analysis that integrated Hi–C data from thirteen different cell lines and dozens of ChIP-seq and RNA-seq datasets measured on these cells, either under basal conditions or after treatment. Pairwise comparisons between cell lines demonstrate a strong association between modulation of A/B compartmentalization, differential gene expression and transcription factor (TF) binding events. Furthermore, integrating the analysis of transcriptomes of different cell lines in response to various challenges, we show that A/B compartmentalization of cells under basal conditions significantly correlates not only with gene expression programs and TF binding profiles that are active under the basal condition but also with those induced in response to treatment. Yet, in pairwise comparisons between different cell lines, we find that a large portion of differential TF binding and gene induction events occur in genomic loci assigned to A compartment in both cell types, underscoring the role of additional critical factors in determining cell-type-specific transcriptional programs.

Conclusions

Our results further indicate the role of dynamic genome organization in regulation of differential gene expression between different cell types and the impact of intra-TAD enhancer–promoter interactions that are established under basal conditions on both the basal and treatment-induced gene expression programs.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0220-2) contains supplementary material, which is available to authorized users.

Storage, visualization, and navigation of 3D genomics data

The field of 3D genomics grew at increasing rates in the last decade. The volume and complexity of 2D and 3D data produced is progressively outpacing the capacities of the technology previously used for distributing genome sequences. The emergence of new technologies provides also novel opportunities for the development of innovative approaches. In this paper, we review the state-of-the-art computing technology, as well as the solutions adopted by the platforms currently available.

Heritable, Allele-Specific Chromosomal Looping between Tandem Promoters Specifies Promoter Usage of SHC1

One-half of the genes in the human genome contain alternative promoters, some of which generate products with opposing functions. Aberrant silencing or activation of such alternative promoters is associated with multiple diseases, including cancer, but little is known regarding the molecular mechanisms that control alternative promoter choice. The SHC1 gene encodes p46^Shc/p52^Shc and p66^Shc, proteins oppositely regulating anchorage-independent growth that are produced by transcription initiated from the upstream and downstream tandem promoters of SHC1, respectively. Here we demonstrate that activation of these promoters is mutually exclusive on separate alleles in single primary endothelial cells in a heritable fashion, ensuring expression of both transcripts by the cell. Peripheral blood lymphocytes that do not transcribe p66^Shc transcribed p52^Shc biallelically. This distinct monoallelic transcription pattern is established by allele-specific chromosomal looping between tandem promoters, which silences the upstream promoter. Our results reveal a new mechanism to control alternative promoter usage through higher-order chromatin structure.

Variable Extent of Lineage-Specificity and Developmental Stage-Specificity of Cohesin and CCCTC-Binding Factor Binding Within the Immunoglobulin and T Cell Receptor Loci

CCCTC-binding factor (CTCF) is largely responsible for the 3D architecture of the genome, in concert with the action of cohesin, through the creation of long-range chromatin loops. Cohesin is hypothesized to be the main driver of these long-range chromatin interactions by the process of loop extrusion. Here, we performed ChIP-seq for CTCF and cohesin in two stages each of T and B cell differentiation and examined the binding pattern in all six antigen receptor (AgR) loci in these lymphocyte progenitors and in mature T and B cells, ES cells, and fibroblasts. The four large AgR loci have many bound CTCF sites, most of which are only occupied in lymphocytes, while only the CTCF sites at the end of each locus near the enhancers or J genes tend to be bound in non-lymphoid cells also. However, despite the generalized lymphocyte restriction of CTCF binding in AgR loci, the Igκ locus is the only locus that also shows significant lineage-specificity (T vs. B cells) and developmental stage-specificity (pre-B vs. pro-B) in CTCF binding. We show that cohesin binding shows greater lineage- and stage-specificity than CTCF at most AgR loci, providing more specificity to the loops. We also show that the culture of pro-B cells in IL7, a common practice to expand the number of cells before ChIP-seq, results in a CTCF-binding pattern resembling pre-B cells, as well as other epigenetic and transcriptional characteristics of pre-B cells. Analysis of the orientation of the CTCF sites show that all sites within the large V portions of the Igh and TCRβ loci have the same orientation. This suggests either a lack of requirement for convergent CTCF sites creating loops, or indicates an absence of any loops between CTCF sites within the V region portion of those loci but only loops to the convergent sites at the D-J-enhancer end of each locus. The V region portions of the Igκ and TCRα/δ loci, by contrast, have CTCF sites in both orientations, providing many options for creating CTCF-mediated convergent loops throughout the loci. CTCF/cohesin loops, along with transcription factors, drives contraction of AgR loci to facilitate the creation of a diverse repertoire of antibodies and T cell receptors.

Adult hippocampal MeCP2 preserves the genomic responsiveness to learning required for long-term memory formation

MeCP2 is required both during postnatal neurodevelopment and throughout the adult life for brain function. Although it is well accepted that MeCP2 in the maturing nervous system is critical for establishing normal development, the functions of MeCP2 during adulthood are poorly understood. Particularly, the requirement of hippocampal MeCP2 for cognitive abilities in the adult is not studied. To characterize the role of MeCP2 in adult neuronal function and cognition, we used a temporal and region-specific disruption of MeCP2 expression in the hippocampus of adult male mice. We found that MeCP2 is required for long-term memory formation and that it controls the learning-induced transcriptional response of hippocampal neurons required for memory consolidation. Furthermore, we uncovered MeCP2 functions in the adult hippocampus that may underlie cognitive integrity. We showed that MeCP2 maintains the developmentally established chromatin configuration and epigenetic landscape of CA1 neurons throughout the adulthood, and that it regulates the expression of neuronal and immune-related genes in the adult hippocampus. Overall, our findings identify MeCP2 as a maintenance factor in the adult hippocampus that preserves signal responsiveness of the genome and allows for integrity of cognitive functions. This study provides new insight into how MeCP2 maintains adult brain functions, but also into the mechanisms underlying the cognitive impairments observed in RTT patients and highlights the understudied role of DNA methylation interpretation in adult cognitive processes.

Genome-Wide Association Study to Find Modifiers for Tetralogy of Fallot in the 22q11.2 Deletion Syndrome Identifies Variants in the GPR98 Locus on 5q14.3

BACKGROUND

The 22q11.2 deletion syndrome (22q11.2DS; DiGeorge syndrome/velocardiofacial syndrome) occurs in 1 of 4000 live births, and 60% to 70% of affected individuals have congenital heart disease, ranging from mild to severe. In our cohort of 1472 subjects with 22q11.2DS, a total of 62% (n=906) have congenital heart disease and 36% (n=326) of these have tetralogy of Fallot (TOF), comprising the largest subset of severe congenital heart disease in the cohort.

METHODS AND RESULTS

To identify common genetic variants associated with TOF in individuals with 22q11.2DS, we performed a genome-wide association study using Affymetrix 6.0 array and imputed genotype data. In our cohort, TOF was significantly associated with a genotyped single-nucleotide polymorphism (rs12519770, P=2.98×10-8) in an intron of the adhesion GPR98 (G-protein-coupled receptor V1) gene on chromosome 5q14.3. There was also suggestive evidence of association between TOF and several additional single-nucleotide polymorphisms in this region. Some genome-wide significant loci in introns or noncoding regions could affect regulation of genes nearby or at a distance. On the basis of this possibility, we examined existing Hi-C chromatin conformation data to identify genes that might be under shared transcriptional regulation within the region on 5q14.3. There are 6 genes in a topologically associated domain of chromatin with GPR98, including MEF2C (Myocyte-specific enhancer factor 2C). MEF2C is the only gene that is known to affect heart development in mammals and might be of interest with respect to 22q11.2DS.

CONCLUSIONS

In conclusion, common variants may contribute to TOF in 22q11.2DS and may function in cardiac outflow tract development.

What's Luck Got to Do with It: Single Cells, Multiple Fates, and Biological Nondeterminism

The field of single-cell biology has morphed from a philosophical digression at its inception, to a playground for quantitative biologists, to a major area of biomedical research. The last several years have witnessed an explosion of new technologies, allowing us to apply even more of the modern molecular biology toolkit to single cells. Conceptual progress, however, has been comparatively slow. Here, we provide a framework for classifying both the origins of the differences between individual cells and the consequences of those differences. We discuss how the concept of "random" differences is context dependent, and propose that rigorous definitions of inputs and outputs may bring clarity to the discussion. We also categorize ways in which probabilistic behavior may influence cellular function, highlighting studies that point to exciting future directions in the field.

The Molecular Revolution in Cutaneous Biology: Chromosomal Territories, Higher-Order Chromatin Remodeling, and the Control of Gene Expression in Keratinocytes

Three-dimensional organization of transcription in the nucleus and mechanisms controlling the global chromatin folding, including spatial interactions between the genes, noncoding genome elements, and epigenetic and transcription machinery, are essential for establishing lineage-specific gene expression programs during cell differentiation. Spatial chromatin interactions in the nucleus involving gene promoters and distal regulatory elements are currently considered major forces that drive cell differentiation and genome evolution in general, and such interactions are substantially reorganized during many pathological conditions. During terminal differentiation of the epidermal keratinocytes, the nucleus undergoes programmed transformation from highly active status, associated with execution of the genetic program of epidermal barrier formation, to a fully inactive condition and finally becomes a part of the keratinized cells of the cornified epidermal layer. This transition is accompanied by marked remodeling of the three-dimensional nuclear organization and microanatomy, including changes in the spatial arrangement of lineage-specific genes, nuclear bodies, and heterochromatin. This mini-review highlights the important landmarks in the accumulation of our current knowledge on three-dimensional organization of the nucleus, spatial arrangement of the genes, and their distal regulatory elements, and it provides an update on the mechanisms that control higher-order chromatin remodeling in the context of epidermal keratinocyte differentiation in the skin.

A Pair of Maternal Chromosomes Derived from Meiotic Nondisjunction in Trisomy 21 Affects Nuclear Architecture and Transcriptional Regulation

Eukaryotic genomes are organised into complex higher-order structures within the nucleus, and the three-dimensional arrangement of chromosomes is functionally important for global gene regulation. The existence of supernumerary chromosome 21 in Down syndrome may perturb the nuclear architecture at different levels, which is normally optimised to maintain the physiological balance of gene expression. However, it has not been clearly elucidated whether and how aberrant configuration of chromosomes affects gene activities. To investigate the effects of trisomy 21 on nuclear organisation and gene expression, we performed three-dimensional fluorescent imaging analysis of chromosome-edited human induced pluripotent stem cells (iPSCs), which enabled identification of the parental origin of the three copies of chromosome 21. We found that two copies of maternal chromosomes resulting from meiotic nondisjunction had a higher tendency to form an adjacent pair and were located relatively distant from the nuclear membrane, suggesting the conserved interaction between these homologous chromosomes. Transcriptional profiling of parental-origin-specific corrected disomy 21 iPSC lines indicated upregulated expression of the maternal alleles for a group of genes, which was accompanied by a fluctuating expression pattern. These results suggest the unique effects of a pair of maternal chromosomes in trisomy 21, which may contribute to the pathological phenotype.

Identifying Causal Genes at the Multiple Sclerosis Associated Region 6q23 Using Capture Hi-C

Background

The chromosomal region 6q23 has been found to be associated with multiple sclerosis (MS) predisposition through genome wide association studies (GWAS). There are four independent single nucleotide polymorphisms (SNPs) associated with MS in this region, which spans around 2.5 Mb. Most GWAS variants associated with complex traits, including these four MS associated SNPs, are non-coding and their function is currently unknown. However, GWAS variants have been found to be enriched in enhancers and there is evidence that they may be involved in transcriptional regulation of their distant target genes through long range chromatin looping.

Aim

The aim of this work is to identify causal disease genes in the 6q23 locus by studying long range chromatin interactions, using the recently developed Capture Hi-C method in human T and B-cell lines. Interactions involving four independent associations unique to MS, tagged by rs11154801, rs17066096, rs7769192 and rs67297943 were analysed using Capture Hi-C Analysis of Genomic Organisation (CHiCAGO).

Results

We found that the pattern of chromatin looping interactions in the MS 6q23 associated region is complex. Interactions cluster in two regions, the first involving the rs11154801 region and a second containing the rs17066096, rs7769192 and rs67297943 SNPs. Firstly, SNPs located within the AHI1 gene, tagged by rs11154801, are correlated with expression of AHI1 and interact with its promoter. These SNPs also interact with other potential candidate genes such as SGK1 and BCLAF1. Secondly, the rs17066096, rs7769192 and rs67297943 SNPs interact with each other and with immune-related genes such as IL20RA, IL22RA2, IFNGR1 and TNFAIP3. Finally, the above-mentioned regions interact with each other and therefore, may co-regulate these target genes.

Conclusion

These results suggest that the four 6q23 variants, independently associated with MS, are involved in the regulation of several genes, including immune genes. These findings could help understand mechanisms of disease and suggest potential novel therapeutic targets.

Capture Hi-C identifies a novel causal gene, IL20RA, in the pan-autoimmune genetic susceptibility region 6q23

BACKGROUND

The identification of causal genes from genome-wide association studies (GWAS) is the next important step for the translation of genetic findings into biologically meaningful mechanisms of disease and potential therapeutic targets. Using novel chromatin interaction detection techniques and allele specific assays in T and B cell lines, we provide compelling evidence that redefines causal genes at the 6q23 locus, one of the most important loci that confers autoimmunity risk.

RESULTS

Although the function of disease-associated non-coding single nucleotide polymorphisms (SNPs) at 6q23 is unknown, the association is generally assigned to TNFAIP3, the closest gene. However, the DNA fragment containing the associated SNPs interacts through chromatin looping not only with TNFAIP3, but also with IL20RA, located 680 kb upstream. The risk allele of the most likely causal SNP, rs6927172, is correlated with both a higher frequency of interactions and increased expression of IL20RA, along with a stronger binding of both the NFκB transcription factor and chromatin marks characteristic of active enhancers in T-cells.

CONCLUSIONS

Our results highlight the importance of gene assignment for translating GWAS findings into biologically meaningful mechanisms of disease and potential therapeutic targets; indeed, monoclonal antibody therapy targeting IL-20 is effective in the treatment of rheumatoid arthritis and psoriasis, both with strong GWAS associations to this region.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

Expressed alleles of imprinted IGF2, DLK1 and MEG3 colocalize in 3D-preserved nuclei of porcine fetal cells

Background

To explore the relationship between spatial genome organization and gene expression in the interphase nucleus, we used a genomic imprinting model, which offers parental-specific gene expression. Using 3D FISH in porcine fetal liver cells, we compared the nuclear organization of the two parental alleles (expressed or not) of insulin-like growth factor 2 (IGF2), a paternally imprinted gene located on chromosome 2. We investigated whether its nuclear positioning favors specific locus associations. We also tested whether IGF2 is implicated in long-range chromatin trans-associations as previously shown in the mouse model species for its reciprocal imprinted gene H19.

Results

We focused on the 3D position of IGF2 alleles, with respect to their individual chromosome 2 territories. The paternally expressed allele was tagged with nascent RNA. There were no significant differences in the position of the two alleles (p = 0.06). To determine long-range chromatin trans-interactions, we chose 12 genes, some of which are known to be imprinted in mammalian model species and belong to a network of imprinted genes (i.e. SLC38A4, DLK1, MEG3, and ZAC1). We screened them and ABCG2, OSBP2, OSBPL1, RPL32, NF1, ZAR1, SEP15, GPC3 for associations with IGF2 in liver cells. All imprinted genes tested showed an association with IGF2. The DLK1/MEG3 locus showed the highest rate of colocalization. This gene association was confirmed by 3D FISH (in 20 % of the nuclei analyzed), revealing also the close proximity of chromosomes 2 and 7 (in 60 % of nuclei). Furthermore, our observations showed that the expressed paternal IGF2 allele is involved in this association. This IGF2-(DLK1/MEG3) association also occurred in a high percentage of fetal muscle cells (36 % of nuclei). Finally, we showed that nascent IGF2, DLK1 and MEG3 RNAs can associate in pairs or in a three-way combination.

Conclusion

Our results show that trans-associations occur between three imprinted genes IGF2, DLK1 and MEG3 both in fetal liver and muscle cells. All three expressed alleles associated in muscle cells. Our findings suggest that the 3D nuclear organization is linked to the transcriptional state of these genes.

Electronic supplementary material

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

Nuclear organization and 3D chromatin architecture in cognition and neuropsychiatric disorders

The current view of neuroplasticity depicts the changes in the strength and number of synaptic connections as the main physical substrate for behavioral adaptation to new experiences in a changing environment. Although transcriptional regulation is known to play a role in these synaptic changes, the specific contribution of activity-induced changes to both the structure of the nucleus and the organization of the genome remains insufficiently characterized. Increasing evidence indicates that plasticity-related genes may work in coordination and share architectural and transcriptional machinery within discrete genomic foci. Here we review the molecular and cellular mechanisms through which neuronal nuclei structurally adapt to stimuli and discuss how the perturbation of these mechanisms can trigger behavioral malfunction.

Post-transcriptional and translational regulation modulates gene co-expression behavior in more synchronized pace to carry out molecular function in the cell

MOTIVATION

Biological processes involve much complex interplay between cellular molecules at different molecular levels, and this interplay may exhibit various co-expression patterns explicitly representing the cellular inner regulation mechanism. Whereas, coexpression patterns cannot be necessarily conserved across the different molecular levels for complex regulation processes involved even after transcripts being produced. Investigation of co-expression propagation from transcript level to protein level will reflect inner regulation effects in function states of cells.

RESULTS

In this study, we perform a comparative analysis of gene coexpression patterns in Plasmodium falciparum. We investigate coexpression patterns propagation from transcript level to protein level to reveal the underlying biological meaning of post-transcriptional and translational mechanism. Our systems-level approach shows after posttranscriptional and translational regulation gene co-expression pace at protein level is mechanistically adjusted to higher synchronicity. Moreover, co-expression patterns at protein level are more linked to function categories, such as co-expression at the same time point is more related with binding categories, and co-expression delayed by several time points is more related with activity categories. Therefore, posttranscriptional and translational regulation modulates co-expression relationships between molecules for meeting the function demands.

The Chromatin Signature of Pluripotency: Establishment and Maintenance

The revolutionary discovery that somatic cells can be reprogrammed by a defined set transcription factors to induced pluripotent stem cells (iPSCs) changed dramatically the way we perceive cell fate determination. Importantly, iPSCs, similar to embryo-derived stem cells (ESCs), are characterized by a remarkable developmental plasticity and the capacity to self-renew "indefinitely" under appropriate culture conditions, opening new avenues for personalized therapy and disease modeling. Elucidating the molecular mechanisms that maintain, induce, or alter stem cell identity is crucial for a deeper understanding of cell fate determination and potential translational applications. Intense research over the last 10 years exploiting technological advances in epigenomics and genome editing has unraveled many of the mysteries of pluripotent identity enabling novel and efficient ways to manipulate it for biomedical purposes. In this review, we focus on the chromatin and epigenetic characteristics that distinguish stem cells from somatic cells and their dynamic changes during differentiation and reprogramming.

Mining 3D genome structure populations identifies major factors governing the stability of regulatory communities

Three-dimensional (3D) genome structures vary from cell to cell even in an isogenic sample. Unlike protein structures, genome structures are highly plastic, posing a significant challenge for structure-function mapping. Here we report an approach to comprehensively identify 3D chromatin clusters that each occurs frequently across a population of genome structures, either deconvoluted from ensemble-averaged Hi-C data or from a collection of single-cell Hi-C data. Applying our method to a population of genome structures (at the macrodomain resolution) of lymphoblastoid cells, we identify an atlas of stable inter-chromosomal chromatin clusters. A large number of these clusters are enriched in binding of specific regulatory factors and are therefore defined as ‘Regulatory Communities.' We reveal two major factors, centromere clustering and transcription factor binding, which significantly stabilize such communities. Finally, we show that the regulatory communities differ substantially from cell to cell, indicating that expression variability could be impacted by genome structures.

3D genome structures are plastic and vary from cell to cell even in an isogenic sample. Here, the authors present an approach to identify frequent 3D chromatin clusters across a population of genome structures, either deconvoluted from ensemble-averaged Hi-C data or from a collection of single-cell Hi-C data.

Systems Biology Approaches for Understanding Genome Architecture

The linear and three-dimensional arrangement and composition of chromatin in eukaryotic genomes underlies the mechanisms directing gene regulation. Understanding this organization requires the integration of many data types and experimental results. Here we describe the approach of integrating genome-wide protein-DNA binding data to determine chromatin states. To investigate spatial aspects of genome organization, we present a detailed description of how to run stochastic simulations of protein movements within a simulated nucleus in 3D. This systems level approach enables the development of novel questions aimed at understanding the basic mechanisms that regulate genome dynamics.

Minor Loops in Major Folds: Enhancer-Promoter Looping, Chromatin Restructuring, and Their Association with Transcriptional Regulation and Disease

The organization and folding of chromatin within the nucleus can determine the outcome of gene expression. Recent technological advancements have enabled us to study chromatin interactions in a genome-wide manner at high resolution. These studies have increased our understanding of the hierarchy and dynamics of chromatin domains that facilitate cognate enhancer–promoter looping, defining the transcriptional program of different cell types. In this review, we focus on vertebrate chromatin long-range interactions as they relate to transcriptional regulation. In addition, we describe how the alteration of boundaries that mark discrete regions in the genome with high interaction frequencies within them, called topological associated domains (TADs), could lead to various phenotypes, including human diseases, which we term as “TADopathies.”

Capture Hi-C reveals novel candidate genes and complex long-range interactions with related autoimmune risk loci

Genome-wide association studies have been tremendously successful in identifying genetic variants associated with complex diseases. The majority of association signals are intergenic and evidence is accumulating that a high proportion of signals lie in enhancer regions. We use Capture Hi-C to investigate, for the first time, the interactions between associated variants for four autoimmune diseases and their functional targets in B- and T-cell lines. Here we report numerous looping interactions and provide evidence that only a minority of interactions are common to both B- and T-cell lines, suggesting interactions may be highly cell-type specific; some disease-associated SNPs do not interact with the nearest gene but with more compelling candidate genes (for example, FOXO1, AZI2) often situated several megabases away; and finally, regions associated with different autoimmune diseases interact with each other and the same promoter suggesting common autoimmune gene targets (for example, PTPRC, DEXI and ZFP36L1).

There is evidence that a proportion of the polymorphisms identified by genome-wide association studies lie in enchancer regions. Here the authors use Capture Hi-C to investigate the interaction with targets in autoimmune disease, showing interactions can be long range and cell-type specific.

NCLscan: accurate identification of non-co-linear transcripts (fusion, trans-splicing and circular RNA) with a good balance between sensitivity and precision

Analysis of RNA-seq data often detects numerous ‘non-co-linear’ (NCL) transcripts, which comprised sequence segments that are topologically inconsistent with their corresponding DNA sequences in the reference genome. However, detection of NCL transcripts involves two major challenges: removal of false positives arising from alignment artifacts and discrimination between different types of NCL transcripts (trans-spliced, circular or fusion transcripts). Here, we developed a new NCL-transcript-detecting method (‘NCLscan’), which utilized a stepwise alignment strategy to almost completely eliminate false calls (>98% precision) without sacrificing true positives, enabling NCLscan outperform 18 other publicly-available tools (including fusion- and circular-RNA-detecting tools) in terms of sensitivity and precision, regardless of the generation strategy of simulated dataset, type of intragenic or intergenic NCL event, read depth of coverage, read length or expression level of NCL transcript. With the high accuracy, NCLscan was applied to distinguishing between trans-spliced, circular and fusion transcripts on the basis of poly(A)- and nonpoly(A)-selected RNA-seq data. We showed that circular RNAs were expressed more ubiquitously, more abundantly and less cell type-specifically than trans-spliced and fusion transcripts. Our study thus describes a robust pipeline for the discovery of NCL transcripts, and sheds light on the fundamental biology of these non-canonical RNA events in human transcriptome.

Long-Range Chromatin Interactions

To accommodate genomes in the limited space of the cell nucleus and ensure the correct execution of gene expression programs, genomes are packaged in complex fashion in the three-dimensional cell nucleus. As a consequence of the extensive higher-order organization of chromosomes, distantly located genomic regions on the same or distinct chromosomes undergo long-range interactions. This article discusses the nature of long interactions, mechanisms of their formation, and their emerging functional roles in gene regulation and genome maintenance.

Nucleolar tethering mediates pairing between the IgH and Myc loci

Gene loci on different chromosomes can preferentially colocalize in the cell nucleus. However, many of the mechanisms mediating this spatial proximity remain to be elucidated. The IgH locus on Chromosome 12 and the Myc locus on Chromosome 15 are a well-studied model for gene colocalization in murine B cells, where the two loci are positioned in close proximity at a higher than expected frequency. These gene loci are also partners in the chromosomal translocation that causes murine plasmacytoma and Burkitt’s lymphoma. Because both Chromosome 12 and Chromosome 15 carry nucleolar organizer regions (NORs) in the most commonly studied mouse strains, we hypothesized that NOR-mediated tethering of the IgH and Myc loci to shared nucleoli could serve as a mechanism to drive IgH:Myc colocalization. Using mouse strains that naturally carry nucleolar organizer regions (NORs) on different sets of chromosomes, we establish that IgH and Myc are positioned proximal to nucleoli in a NOR dependent manner and show that their joint association with nucleoli significantly increases the frequency of IgH and Myc pairing. Thus we demonstrate that simple nucleolar tethering can increase the colocalization frequency of genes on NOR-bearing chromosomes.

Mango: a bias-correcting ChIA-PET analysis pipeline

MOTIVATION

Chromatin Interaction Analysis by Paired-End Tag sequencing (ChIA-PET) is an established method for detecting genome-wide looping interactions at high resolution. Current ChIA-PET analysis software packages either fail to correct for non-specific interactions due to genomic proximity or only address a fraction of the steps required for data processing. We present Mango, a complete ChIA-PET data analysis pipeline that provides statistical confidence estimates for interactions and corrects for major sources of bias including differential peak enrichment and genomic proximity.

RESULTS

Comparison to the existing software packages, ChIA-PET Tool and ChiaSig revealed that Mango interactions exhibit much better agreement with high-resolution Hi-C data. Importantly, Mango executes all steps required for processing ChIA-PET datasets, whereas ChiaSig only completes 20% of the required steps. Application of Mango to multiple available ChIA-PET datasets permitted the independent rediscovery of known trends in chromatin loops including enrichment of CTCF, RAD21, SMC3 and ZNF143 at the anchor regions of interactions and strong bias for convergent CTCF motifs.

AVAILABILITY AND IMPLEMENTATION

Mango is open source and distributed through github at https://github.com/dphansti/mango.

CONTACT

mpsnyder@standford.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Nucleosome positioning and composition modulate in silico chromatin flexibility

The dynamic organization of chromatin plays an essential role in the regulation of gene expression and in other fundamental cellular processes. The underlying physical basis of these activities lies in the sequential positioning, chemical composition, and intermolecular interactions of the nucleosomes-the familiar assemblies of ∼150 DNA base pairs and eight histone proteins-found on chromatin fibers. Here we introduce a mesoscale model of short nucleosomal arrays and a computational framework that make it possible to incorporate detailed structural features of DNA and histones in simulations of short chromatin constructs. We explore the effects of nucleosome positioning and the presence or absence of cationic N-terminal histone tails on the 'local' inter-nucleosomal interactions and the global deformations of the simulated chains. The correspondence between the predicted and observed effects of nucleosome composition and numbers on the long-range communication between the ends of designed nucleosome arrays lends credence to the model and to the molecular insights gleaned from the simulated structures. We also extract effective nucleosome-nucleosome potentials from the simulations and implement the potentials in a larger-scale computational treatment of regularly repeating chromatin fibers. Our results reveal a remarkable effect of nucleosome spacing on chromatin flexibility, with small changes in DNA linker length significantly altering the interactions of nucleosomes and the dimensions of the fiber as a whole. In addition, we find that these changes in nucleosome positioning influence the statistical properties of long chromatin constructs. That is, simulated chromatin fibers with the same number of nucleosomes exhibit polymeric behaviors ranging from Gaussian to worm-like, depending upon nucleosome spacing. These findings suggest that the physical and mechanical properties of chromatin can span a wide range of behaviors, depending on nucleosome positioning, and that care must be taken in the choice of models used to interpret the experimental properties of long chromatin fibers.

Rbs1, a new protein implicated in RNA polymerase III biogenesis in yeast Saccharomyces cerevisiae

Little is known about the RNA polymerase III (Pol III) complex assembly and its transport to the nucleus. We demonstrate that a missense cold-sensitive mutation, rpc128-1007, in the sequence encoding the C-terminal part of the second largest Pol III subunit, C128, affects the assembly and stability of the enzyme. The cellular levels and nuclear concentration of selected Pol III subunits were decreased in rpc128-1007 cells, and the association between Pol III subunits as evaluated by coimmunoprecipitation was also reduced. To identify the proteins involved in Pol III assembly, we performed a genetic screen for suppressors of the rpc128-1007 mutation and selected the Rbs1 gene, whose overexpression enhanced de novo tRNA transcription in rpc128-1007 cells, which correlated with increased stability, nuclear concentration, and interaction of Pol III subunits. The rpc128-1007 rbs1Δ double mutant shows a synthetic growth defect, indicating that rpc128-1007 and rbs1Δ function in parallel ways to negatively regulate Pol III assembly. Rbs1 physically interacts with a subset of Pol III subunits, AC19, AC40, and ABC27/Rpb5. Additionally, Rbs1 interacts with the Crm1 exportin and shuttles between the cytoplasm and nucleus. We postulate that Rbs1 binds to the Pol III complex or subcomplex and facilitates its translocation to the nucleus.