Chromatin conformation analysis of primary patient tissue using a low input Hi-C method

The Transition from Quiescent to Activated States in Human Hematopoietic Stem Cells Is Governed by Dynamic 3D Genome Reorganization

Lifelong blood production requires long-term hematopoietic stem cells (LT-HSCs), marked by stemness states involving quiescence and self-renewal, to transition into activated short-term HSCs (ST-HSCs) with reduced stemness. As few transcriptional changes underlie this transition, we used single-cell and bulk assay for transposase-accessible chromatin sequencing (ATAC-seq) on human HSCs and hematopoietic stem and progenitor cell (HSPC) subsets to uncover chromatin accessibility signatures, one including LT-HSCs (LT/HSPC signature) and another excluding LT-HSCs (activated HSPC [Act/HSPC] signature). These signatures inversely correlated during early hematopoietic commitment and differentiation. The Act/HSPC signature contains CCCTC-binding factor (CTCF) binding sites mediating 351 chromatin interactions engaged in ST-HSCs, but not LT-HSCs, enclosing multiple stemness pathway genes active in LT-HSCs and repressed in ST-HSCs. CTCF silencing derepressed stemness genes, restraining quiescent LT-HSCs from transitioning to activated ST-HSCs. Hence, 3D chromatin interactions centrally mediated by CTCF endow a gatekeeper function that governs the earliest fate transitions HSCs make by coordinating disparate stemness pathways linked to quiescence and self-renewal.

Structural Variants Create New Topological-Associated Domains and Ectopic Retinal Enhancer-Gene Contact in Dominant Retinitis Pigmentosa

The cause of autosomal-dominant retinitis pigmentosa (adRP), which leads to loss of vision and blindness, was investigated in families lacking a molecular diagnosis. A refined locus for adRP on Chr17q22 (RP17) was delineated through genotyping and genome sequencing, leading to the identification of structural variants (SVs) that segregate with disease. Eight different complex SVs were characterized in 22 adRP-affected families with >300 affected individuals. All RP17 SVs had breakpoints within a genomic region spanning YPEL2 to LINC01476. To investigate the mechanism of disease, we reprogrammed fibroblasts from affected individuals and controls into induced pluripotent stem cells (iPSCs) and differentiated them into photoreceptor precursor cells (PPCs) or retinal organoids (ROs). Hi-C was performed on ROs, and differential expression of regional genes and a retinal enhancer RNA at this locus was assessed by qPCR. The epigenetic landscape of the region, and Hi-C RO data, showed that YPEL2 sits within its own topologically associating domain (TAD), rich in enhancers with binding sites for retinal transcription factors. The Hi-C map of RP17 ROs revealed creation of a neo-TAD with ectopic contacts between GDPD1 and retinal enhancers, and modeling of all RP17 SVs was consistent with neo-TADs leading to ectopic retinal-specific enhancer-GDPD1 accessibility. qPCR confirmed increased expression of GDPD1 and increased expression of the retinal enhancer that enters the neo-TAD. Altered TAD structure resulting in increased retinal expression of GDPD1 is the likely convergent mechanism of disease, consistent with a dominant gain of function. Our study highlights the importance of SVs as a genomic mechanism in unsolved Mendelian diseases.

CHESS enables quantitative comparison of chromatin contact data and automatic feature extraction

Dynamic changes in the three-dimensional (3D) organization of chromatin are associated with central biological processes, such as transcription, replication and development. Therefore, the comprehensive identification and quantification of these changes is fundamental to understanding of evolutionary and regulatory mechanisms. Here, we present Comparison of Hi-C Experiments using Structural Similarity (CHESS), an algorithm for the comparison of chromatin contact maps and automatic differential feature extraction. We demonstrate the robustness of CHESS to experimental variability and showcase its biological applications on (1) interspecies comparisons of syntenic regions in human and mouse models; (2) intraspecies identification of conformational changes in Zelda-depleted Drosophila embryos; (3) patient-specific aberrant chromatin conformation in a diffuse large B-cell lymphoma sample; and (4) the systematic identification of chromatin contact differences in high-resolution Capture-C data. In summary, CHESS is a computationally efficient method for the comparison and classification of changes in chromatin contact data.

Innovations in plant genetics adapting agriculture to climate change

Developing new genotypes of plants is one of the key options for adaptation of agriculture to climate change. Plants may be required to provide resilience in changed climates or support the migration of agriculture to new regions. Very different genotypes may be required to perform in the modified environments of protected agriculture. Consumers will continue to demand taste, convenience, healthy and safe food and sustainably and ethically produced food, despite the greater challenges of climate in the future. Improving the nutritional value of foods in response to climate change is a significant challenge. Genomic sequences of relevant germplasm and an understanding of the functional role of alleles controlling key traits will be an enabling platform for this innovation.

Use of the epigenetic toolbox to contextualize common variants associated with schizophrenia risk

Schizophrenia is a debilitating psychiatric disorder with a complex genetic architecture and limited understanding of its neuropathology, reflected by the lack of diagnostic measures and effective pharmacological treatments. Geneticists have recently identified more than 145 risk loci comprising hundreds of common variants of small effect sizes, most of which lie in noncoding genomic regions. This review will discuss how the epigenetic toolbox can be applied to contextualize genetic findings in schizophrenia. Progress in next-generation sequencing, along with increasing methodological complexity, has led to the compilation of genome-wide maps of DNA methylation, histone modifications, RNA expression, and more. Integration of chromatin conformation datasets is one of the latest efforts in deciphering schizophrenia risk, allowing the identification of genes in contact with regulatory variants across 100s of kilobases. Large-scale multiomics studies will facilitate the prioritization of putative causal risk variants and gene networks that contribute to schizophrenia etiology, informing clinical diagnostics and treatment downstream.
.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

Integrated paired-end enhancer profiling and whole-genome sequencing reveals recurrent CCNE1 and IGF2 enhancer hijacking in primary gastric adenocarcinoma

OBJECTIVE

Genomic structural variations (SVs) causing rewiring of cis-regulatory elements remain largely unexplored in gastric cancer (GC). To identify SVs affecting enhancer elements in GC (enhancer-based SVs), we integrated epigenomic enhancer profiles revealed by paired-end H3K27ac ChIP-sequencing from primary GCs with tumour whole-genome sequencing (WGS) data (PeNChIP-seq/WGS).

DESIGN

We applied PeNChIP-seq to 11 primary GCs and matched normal tissues combined with WGS profiles of >200 GCs. Epigenome profiles were analysed alongside matched RNA-seq data to identify tumour-associated enhancer-based SVs with altered cancer transcription. Functional validation of candidate enhancer-based SVs was performed using CRISPR/Cas9 genome editing, chromosome conformation capture assays (4C-seq, Capture-C) and Hi-C analysis of primary GCs.

RESULTS

PeNChIP-seq/WGS revealed ~150 enhancer-based SVs in GC. The majority (63%) of SVs linked to target gene deregulation were associated with increased tumour expression. Enhancer-based SVs targeting CCNE1, a key driver of therapy resistance, occurred in 8% of patients frequently juxtaposing diverse distal enhancers to CCNE1 proximal regions. CCNE1-rearranged GCs were associated with high CCNE1 expression, disrupted CCNE1 topologically associating domain (TAD) boundaries, and novel TAD interactions in CCNE1-rearranged primary tumours. We also observed IGF2 enhancer-based SVs, previously noted in colorectal cancer, highlighting a common non-coding genetic driver alteration in gastric and colorectal malignancies.

CONCLUSION

Integrated paired-end NanoChIP-seq and WGS of gastric tumours reveals tumour-associated regulatory SV in regions associated with both simple and complex genomic rearrangements. Genomic rearrangements may thus exploit enhancer-hijacking as a common mechanism to drive oncogene expression in GC.

Technologies to study spatial genome organization: beyond 3C

The way that chromatin is organized in three-dimensional nuclear space is now acknowledged as a factor critical for the major cell processes, like transcription, replication and cell division. Researchers have been armed with new molecular and imaging technologies to study this structure-to-function link of genomes, spearheaded by the introduction of the 'chromosome conformation capture' technology more than a decade ago. However, this technology is not without shortcomings, and novel variants and orthogonal approaches are being developed to overcome these. As a result, the field of nuclear organization is constantly fueled by methods of increasing resolution and/or throughput that strive to eliminate systematic biases and increase precision. In this review, we attempt to highlight the most recent advances in technology that promise to provide novel insights on how chromosomes fold and function.

Coolpup.py: versatile pile-up analysis of Hi-C data

MOTIVATION

Hi-C is currently the method of choice to investigate the global 3D organization of the genome. A major limitation of Hi-C is the sequencing depth required to robustly detect loops in the data. A popular approach used to mitigate this issue, even in single-cell Hi-C data, is genome-wide averaging (piling-up) of peaks, or other features, annotated in high-resolution datasets, to measure their prominence in less deeply sequenced data. However, current tools do not provide a computationally efficient and versatile implementation of this approach.

RESULTS

Here, we describe coolpup.py-a versatile tool to perform pile-up analysis on Hi-C data. We demonstrate its utility by replicating previously published findings regarding the role of cohesin and CTCF in 3D genome organization, as well as discovering novel details of Polycomb-driven interactions. We also present a novel variation of the pile-up approach that can aid the statistical analysis of looping interactions. We anticipate that coolpup.py will aid in Hi-C data analysis by allowing easy to use, versatile and efficient generation of pile-ups.

AVAILABILITY AND IMPLEMENTATION

Coolpup.py is cross-platform, open-source and free (MIT licensed) software. Source code is available from https://github.com/Phlya/coolpuppy and it can be installed from the Python Packaging Index.

Human Endometrial Transcriptome and Progesterone Receptor Cistrome Reveal Important Pathways and Epithelial Regulators

CONTEXT

Poor uterine receptivity is one major factor leading to pregnancy loss and infertility. Understanding the molecular events governing successful implantation is hence critical in combating infertility.

OBJECTIVE

To define Progesterone Receptor (PGR)-regulated molecular mechanisms and epithelial roles in receptivity.

DESIGN

RNA-sequencing and PGR-ChIP-seq were conducted in parallel to identify PGR-regulated pathways during the Window of implantation (WOI) in endometrium of fertile women.

SETTING

Endometrial biopsies from the proliferative and mid-secretory phases were analyzed.

PATIENTS OR OTHER PARTICIPANTS

Participants were fertile, reproductive aged (18-37 years) women with normal cycle length, and without any history of dysmenorrhea, infertility, or irregular cycles. In total, 42 endometrial biopsies obtained from 42 women were analyzed in this study.

INTERVENTIONS

There were no interventions during this study.

MAIN OUTCOME MEASURES

Here we measured the alterations in gene expression and PGR occupancy in the genome during the WOI, based on the hypothesis that PGR binds uterine chromatin cycle dependently to regulate genes involved in uterine cell differentiation and function.

RESULTS

653 genes were identified with regulated PGR binding and differential expression during the WOI. These were involved in regulating inflammatory response, xenobiotic metabolism, epithelial mesenchymal transition, cell death, interleukin/Signal Transducer And Activator Of Transcription (STAT) signaling, estrogen response, and Mammalian target of rapamycin complex 1 (MTORC1) response. Transcriptome of the epithelium identified 3052 differentially expressed genes, of which 658 were uniquely regulated. Transcription factors Interferon Regulatory Factor 8 (IRF8) and Myocyte Enhancer Factor 2C (MEF2C) were found to be regulated in the epithelium during the WOI at the protein level, suggesting potentially important functions that are previously unrecognized.

CONCLUSION

PGR binds the genomic regions of genes regulating critical processes in uterine receptivity and function.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Cohesin Disrupts Polycomb-Dependent Chromosome Interactions in Embryonic Stem Cells

How chromosome organization is related to genome function remains poorly understood. Cohesin, loop extrusion, and CCCTC-binding factor (CTCF) have been proposed to create topologically associating domains (TADs) to regulate gene expression. Here, we examine chromosome conformation in embryonic stem cells lacking cohesin and find, as in other cell types, that cohesin is required to create TADs and regulate A/B compartmentalization. However, in the absence of cohesin, we identify a series of long-range chromosomal interactions that persist. These correspond to regions of the genome occupied by the polycomb repressive system and are dependent on PRC1. Importantly, we discover that cohesin counteracts these polycomb-dependent interactions, but not interactions between super-enhancers. This disruptive activity is independent of CTCF and insulation and appears to modulate gene repression by the polycomb system. Therefore, we discover that cohesin disrupts polycomb-dependent chromosome interactions to modulate gene expression in embryonic stem cells.

Cancer Is Associated with Alterations in the Three-Dimensional Organization of the Genome

The human genome is organized into topologically associating domains (TADs), which represent contiguous regions with a higher frequency of intra-interactions as opposed to inter-interactions. TADs contribute to gene expression regulation by restricting the interactions between their regulatory elements, and TAD disruption has been associated with cancer. Here, we provide a proof of principle that mutations within TADs can be used to predict the survival of cancer patients. Specifically, we constructed a set of 1467 consensus TADs representing the three-dimensional organization of the human genome and used Cox regression analysis to identify a total of 35 prognostic TADs in different cancer types. Interestingly, only 46% of the 35 prognostic TADs comprised genes with known clinical relevance. Moreover, in the vast majority of such cases, the prognostic value of the TAD was not directly related to the presence/absence of mutations in the gene(s), emphasizing the importance of regulatory mutations. In addition, we found that 34% of the prognostic TADs show strong structural perturbations in the cancer genome, consistent with the widespread, global epigenetic dysregulation often observed in cancer patients. In summary, this study elucidates the mechanisms through which non-coding variants may influence cancer progression and opens new avenues for personalized medicine.

The Many Faces of Gene Regulation in Cancer: A Computational Oncogenomics Outlook

Cancer is a complex disease at many different levels. The molecular phenomenology of cancer is also quite rich. The mutational and genomic origins of cancer and their downstream effects on processes such as the reprogramming of the gene regulatory control and the molecular pathways depending on such control have been recognized as central to the characterization of the disease. More important though is the understanding of their causes, prognosis, and therapeutics. There is a multitude of factors associated with anomalous control of gene expression in cancer. Many of these factors are now amenable to be studied comprehensively by means of experiments based on diverse omic technologies. However, characterizing each dimension of the phenomenon individually has proven to fall short in presenting a clear picture of expression regulation as a whole. In this review article, we discuss some of the more relevant factors affecting gene expression control both, under normal conditions and in tumor settings. We describe the different omic approaches that we can use as well as the computational genomic analysis needed to track down these factors. Then we present theoretical and computational frameworks developed to integrate the amount of diverse information provided by such single-omic analyses. We contextualize this within a systems biology-based multi-omic regulation setting, aimed at better understanding the complex interplay of gene expression deregulation in cancer.

Structural Variants as a Basis for Targeted Therapies in Hematological Malignancies

Structural variants (SV) are changes in the genomic landscape that can alter gene expression levels and thus lead to disease development. The most common and best studied SVs in hematological malignancies are chromosomal translocations. Here, parts of two genes that are normally on different chromosomes come into close proximity due to a failure in DNA repair. As a consequence, fusion proteins which show a different function and/or cellular localization compared to the two original proteins are expressed, sometimes even at different levels. The identification of chromosomal translocations is often used to identify the specific disease a patient is suffering from. In addition, SVs such as deletions, duplications, inversions and single nucleotide polymorphisms (SNPs) can occur in hematopoietic cells and lead to their malignant transformations. Changes in the 3D genome structure have also recently been shown to impact disease development. In this review, we describe a variety of SVs occurring in different subtypes of hematological malignancies. Currently, most therapeutic approaches target fusion proteins which are the cellular product of chromosomal translocations. However, amplifications and SNPs also play a role in disease progression and can be targeted. We present some examples for different types of structural variants and how they are currently treated.

Visualising three-dimensional genome organisation in two dimensions

The three-dimensional organisation of the genome plays a crucial role in developmental gene regulation. In recent years, techniques to investigate this organisation have become more accessible to labs worldwide due to improvements in protocols and decreases in the cost of high-throughput sequencing. However, the resulting datasets are complex and can be challenging to analyse and interpret. Here, we provide a guide to visualisation approaches that can aid the interpretation of such datasets and the communication of biological results.

The Dynamic 3D Genome in Gametogenesis and Early Embryonic Development

During gametogenesis and early embryonic development, the chromatin architecture changes dramatically, and both the transcriptomic and epigenomic landscape are comprehensively reprogrammed. Understanding these processes is the holy grail in developmental biology and a key step towards evolution. The 3D conformation of chromatin plays a central role in the organization and function of nuclei. Recently, the dynamics of chromatin structures have been profiled in many model and non-model systems, from insects to mammals, resulting in an interesting comparison. In this review, we first introduce the research methods of 3D chromatin structure with low-input material suitable for embryonic study. Then, the dynamics of 3D chromatin architectures during gametogenesis and early embryonic development is summarized and compared between species. Finally, we discuss the possible mechanisms for triggering the formation of genome 3D conformation in early development.

A SY-Stematic approach towards understanding stem cell biology

The 2nd SY-Stem Symposium - a symposium for 'the next generation of stem cell researchers' - was held on the 21-23 March 2019 at the Vienna BioCenter in Austria. After the great success of the initial SY-Stem meeting in 2018, this year's event again focused on the work of young scientists. Here, we summarize the impressive amount of new research covering stem cell-related fields that was discussed at the meeting.

Androgen receptor enhancer usage and the chromatin regulatory landscape in human prostate cancers

The androgen receptor (AR) is commonly known as a key transcription factor in prostate cancer development, progression and therapy resistance. Genome-wide chromatin association studies revealed that transcriptional regulation by AR mainly depends on binding to distal regulatory enhancer elements that control gene expression through chromatin looping to gene promoters. Changes in the chromatin epigenetic landscape and DNA sequence can locally alter AR-DNA-binding capacity and consequently impact transcriptional output and disease outcome. The vast majority of reports describing AR chromatin interactions have been limited to cell lines, identifying numerous other factors and interacting transcription factors that impact AR chromatin interactions. Do these factors also impact AR cistromics - the genome-wide chromatin-binding landscape of AR - in vivo? Recent technological advances now enable researchers to identify AR chromatin-binding sites and their target genes in human specimens. In this review, we provide an overview of the different factors that influence AR chromatin binding in prostate cancer specimens, which is complemented with knowledge from cell line studies. Finally, we discuss novel perspectives on studying AR cistromics in clinical samples.