Understanding Spatial Genome Organization: Methods and Insights

Changes in Circulating MicroRNA Levels as Potential Indicators of Training Adaptation in Professional Volleyball Players

The increasing demand placed on professional athletes to enhance their fitness and performance has prompted the search for new, more sensitive biomarkers of physiological ability. One such potential biomarker includes microRNA (miRNA) small regulatory RNA sequences. The study investigated the levels of the selected circulating miRNAs before and after a 10-week training cycle in 12 professional female volleyball players, as well as their association with cortisol, creatine kinase (CK), and interleukin 6 (IL-6), using the qPCR technique. Significant decreases in the miR-22 (0.40 ± 0.1 vs. 0.28 ± 0.12, p = 0.009), miR-17 (0.35 ± 0.13 vs. 0.23 ± 0.08; p = 0.039), miR-24 (0.09 ± 0.04 vs. 0.05 ± 0.02; p = 0.001), and miR-26a (0.11 ± 0.06 vs. 0.06 ± 0.04; p = 0.003) levels were observed after training, alongside reduced levels of cortisol and IL-6. The correlation analysis revealed associations between the miRNAs' relative quantity and the CK concentrations, highlighting their potential role in the muscle repair processes. The linear regression analysis indicated that miR-24 and miR-26a had the greatest impact on the CK levels. The study provides insights into the dynamic changes in the miRNA levels during training, suggesting their potential as biomarkers for monitoring the adaptive responses to exercise. Overall, the findings contribute to a better understanding of the physiological effects of exercise and the potential use of miRNAs, especially miR-24 and miR-26a, as biomarkers in sports science and medicine.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Epigenetic Analysis of the Dopamine Transporter Gene DAT1 with a Focus on Personality Traits in Athletes

Human phenotypes (traits) are determined by the selective use of a person's unique genotype (DNA sequence), following exposure to environmental stimuli, such as exercise. Inducing profound changes in epigenetics may be an underlying factor of the beneficial effects of exercise. This study aimed to investigate the association between methylation in the promoter region of the DAT1 gene and personality traits measured by the NEO-FFI questionnaire in a group of athletes. The study group included 163 athletes, and the control group consisted of 232 non-athletes. The obtained results show several significant differences between the studied groups of subjects. The Extraversion scale and the Conscientiousness scale results of the NEO-FFI are significantly higher in the group of athletes compared to controls. The total methylation and the number of methylated islands in the promoter region of the DAT1 gene are higher in the study group. Pearson's linear correlation between the total methylation, the number of methylated islands and the NEO-FFI shows significant results for the Extraversion and Agreeability scales. The total methylation and the number of methylated islands in the promoter region of the DAT1 gene are higher in the study group. Pearson's linear correlation between the total methylation, the number of methylated islands and the NEO-FFI shows significant results for the Extraversion and Agreeability scales. Our analysis of the methylation status of individual CpG sites revealed a new direction of research into the biological aspects of regulating dopamine release and personality traits in people practicing sports.

Characterization of network hierarchy reflects cell state specificity in genome organization

Dynamic chromatin structure acts as the regulator of transcription program in crucial processes including cancer and cell development, but a unified framework for characterizing chromatin structural evolution remains to be established. Here, we performed graph inferences on Hi-C data sets and derived the chromatin contact networks. We discovered significant decreases in information transmission efficiencies in chromatin of colorectal cancer (CRC) and T-cell acute lymphoblastic leukemia (T-ALL) compared to corresponding normal controls through graph statistics. Using network embedding in the Poincaré disk, the hierarchy depths of chromatin from CRC and T-ALL patients were found to be significantly shallower compared to their normal controls. A reverse trend of change in chromatin structure was observed during early embryo development. We found tissue-specific conservation of hierarchy order in chromatin contact networks. Our findings reveal the top-down hierarchy of chromatin organization, which is significantly attenuated in cancer.

Three-dimensional genome landscape comprehensively reveals patterns of spatial gene regulation in papillary and anaplastic thyroid cancers: a study using representative cell lines for each cancer type

BACKGROUND

Spatial chromatin structure is intricately linked with somatic aberrations, and somatic mutations of various cancer-related genes, termed co-mutations (CoMuts), occur in certain patterns during cancer initiation and progression. The functional mechanisms underlying these genetic events remain largely unclear in thyroid cancer (TC). With discrepant differentiation, papillary thyroid cancer (PTC) and anaplastic thyroid cancer (ATC) differ greatly in characteristics and prognosis. We aimed to reveal the spatial gene alterations and regulations between the two TC subtypes.

METHODS

We systematically investigated and compared the spatial co-mutations between ATC (8305C), PTC (BCPAP and TPC-1), and normal thyroid cells (Nthy-ori-3-1). We constructed a framework integrating whole-genome sequencing (WGS), high-throughput chromosome conformation capture (Hi-C), and transcriptome sequencing, to systematically detect the associations between the somatic co-mutations of cancer-related genes, structural variations (SVs), copy number variations (CNVs), and high-order chromatin conformation.

RESULTS

Spatial co-mutation hotspots were enriched around topologically associating domains (TADs) in TC. A common set of 227 boundaries were identified in both ATC and PTC, with significant overlaps between them. The spatial proximities of the co-mutated gene pairs in the two TC types were significantly greater than in the gene-level and overall backgrounds, and ATC cells had higher TAD contact frequency with CoMuts > 10 compared with PTC cells. Compared with normal thyroid cells, in ATC the number of the created novel three-dimensional chromatin structural domains increased by 10%, and the number of shifted TADs decreased by 7%. We found five TAD blocks with CoMut genes/events specific to ATC with certain mutations in genes including MAST-NSUN4, AM129B/TRUB2, COL5A1/PPP1R26, PPP1R26/GPSM1/CCDC183, and PRAC2/DLX4. For the majority of ATC and PTC cells, the HOXA10 and HIF2α signals close to the transcription start sites of CoMut genes within TADs were significantly stronger than those at the background. CNV breakpoints significantly overlapped with TAD boundaries in both TC subtypes. ATCs had more CNV losses overlapping with TAD boundaries, and noncoding SVs involved in intrachromosomal SVs, amplified inversions, and tandem duplication differed between ATC and PTC. TADs with short range were more abundant in ATC than PTC. More switches of A/B compartment types existed in ATC cells compared with PTC. Gene expression was significantly synchronized, and orchestrated by complex epigenetics and regulatory elements.

CONCLUSION

Chromatin interactions and gene alterations and regulations are largely heterogeneous in TC. CNVs and complex SVs may function in the TC genome by interplaying with TADs, and are largely different between ATC and PTC. Complexity of TC genomes, which are highly organized by 3D genome-wide interactions mediating mutational and structural variations and gene activation, may have been largely underappreciated. Our comprehensive analysis may provide key evidence and targets for more customized diagnosis and treatment of TC.

Double-strand break repair and mis-repair in 3D

DNA double-strand breaks (DSBs) are lesions that arise frequently from exposure to damaging agents as well as from ongoing physiological DNA transactions. Mis-repair of DSBs leads to rearrangements and structural variations in chromosomes, including insertions, deletions, and translocations implicated in disease. The DNA damage response (DDR) limits pathologic mutations and large-scale chromosome rearrangements. DSB repair initiates in 2D at DNA lesions with the stepwise recruitment of repair proteins and local chromatin remodeling which facilitates break accessibility. More complex structures are then formed via protein assembly into nanodomains and via genome folding into chromatin loops. Subsequently, 3D reorganization of DSBs is guided by clustering forces which drive the assembly of repair domains harboring multiple lesions. These domains are further stabilized and insulated into condensates via liquid-liquid phase-separation. Here, we discuss the benefits and risks associated with this 3D reorganization of the broken genome.

Sample-multiplexing approaches for single-cell sequencing

Single-cell sequencing is widely used in biological and medical studies. However, its application with multiple samples is hindered by inefficient sample processing, high experimental costs, ambiguous identification of true single cells, and technical batch effects. Here, we introduce sample-multiplexing approaches for single-cell sequencing in transcriptomics, epigenomics, genomics, and multiomics. In single-cell transcriptomics, sample multiplexing uses variants of native or artificial features as sample markers, enabling sample pooling and decoding. Such features include: (1) natural genetic variation, (2) nucleotide-barcode anchoring on cellular or nuclear membranes, (3) nucleotide-barcode internalization to the cytoplasm or nucleus, (4) vector-based barcode expression in cells, and (5) nucleotide-barcode incorporation during library construction. Other single-cell omics methods are based on similar concepts, particularly single-cell combinatorial indexing. These methods overcome current challenges, while enabling super-loading of single cells. Finally, selection guidelines are presented that can accelerate technological application.

The emerging landscape of spatial profiling technologies

Improved scale, multiplexing and resolution are establishing spatial nucleic acid and protein profiling methods as a major pillar for cellular atlas building of complex samples, from tissues to full organisms. Emerging methods yield omics measurements at resolutions covering the nano- to microscale, enabling the charting of cellular heterogeneity, complex tissue architectures and dynamic changes during development and disease. We present an overview of the developing landscape of in situ spatial genome, transcriptome and proteome technologies, exemplify their impact on cell biology and translational research, and discuss current challenges for their community-wide adoption. Among many transformative applications, we envision that spatial methods will map entire organs and enable next-generation pathology.

Multiple epigenetic factors co-localize with HMGN proteins in A-compartment chromatin

Background

Nucleosomal binding proteins, HMGN, is a family of chromatin architectural proteins that are expressed in all vertebrate nuclei. Although previous studies have discovered that HMGN proteins have important roles in gene regulation and chromatin accessibility, whether and how HMGN proteins affect higher order chromatin status remains unknown.

Results

We examined the roles that HMGN1 and HMGN2 proteins play in higher order chromatin structures in three different cell types. We interrogated data generated in situ, using several techniques, including Hi–C, Promoter Capture Hi–C, ChIP-seq, and ChIP–MS. Our results show that HMGN proteins occupy the A compartment in the 3D nucleus space. In particular, HMGN proteins occupy genomic regions involved in cell-type-specific long-range promoter–enhancer interactions. Interestingly, depletion of HMGN proteins in the three different cell types does not cause structural changes in higher order chromatin, i.e., in topologically associated domains (TADs) and in A/B compartment scores. Using ChIP-seq combined with mass spectrometry, we discovered protein partners that are directly associated with or neighbors of HMGNs on nucleosomes.

Conclusions

We determined how HMGN chromatin architectural proteins are positioned within a 3D nucleus space, including the identification of their binding partners in mononucleosomes. Our research indicates that HMGN proteins localize to active chromatin compartments but do not have major effects on 3D higher order chromatin structure and that their binding to chromatin is not dependent on specific protein partners.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00457-4.

A reappraisal of the form: function problem-theory and phenomenology

This paper is aimed at demonstrating that some geometrical and topological transformations and operations serve not only as promoters of many specific genetic and cellular events in multicellular living organisms, but also as initiators of the organization and regulation of their functions. Thus, changes in the form and structure of macromolecular and cellular systems must be directly associated to their functions. There are specific classes of enzymes that manipulate the geometry and topology of complex DNA-protein structures, and thereby they perform many important cellular processes, including segregation of daughter chromosomes, gene regulation, and DNA repair. We argue that form has an organizing power, hence a causal action, in the sense that it enables to induce functional events during different biological processes, at the supramolecular, cellular, and organismal levels of organization. Clearly, topological forms must be matched with specific kinetic and dynamical parameters to have a functional effectiveness in living systems. This effectiveness is remarkably apparent, to give an example, in the regulation of the genome functions and in cell activity. In more general terms, we try to show that the conformational plasticity of biological systems depends on different kinds of topological manipulations performed by specific families of enzymes. In doing so, they catalyze all those spatial and dynamical changes of biological structures that are suitable for the functions to be acted by the organism.

CRISPR-Cas-Mediated Tethering Recruits the Yeast HMR Mating-Type Locus to the Nuclear Periphery but Fails to Silence Gene Expression

To investigate the relationship between genome structure and function, we have developed a programmable CRISPR-Cas system for nuclear peripheral recruitment in yeast. We benchmarked this system at the HMR and GAL2 loci, both of which are well-characterized model systems for localization to the nuclear periphery. Using microscopy and gene silencing assays, we demonstrate that CRISPR-Cas-mediated tethering can recruit the HMR locus but does not detectably silence reporter gene expression. A previously reported Gal4-mediated tethering system does silence gene expression, and we demonstrate that the silencing effect has an unexpected dependence on the properties of the protein tether. The CRISPR-Cas system was unable to recruit GAL2 to the nuclear periphery. Our results reveal potential challenges for synthetic genome structure perturbations and suggest that distinct functional effects can arise from subtle structural differences in how genes are recruited to the periphery.

From DNA human sequence to the chromatin higher order organisation and its biological meaning: Using biomolecular interaction networks to understand the influence of structural variation on spatial genome organisation and its functional effect

The three-dimensional structure of the human genome has been proven to have a significant functional impact on gene expression. The high-order spatial chromatin is organised first by looping mediated by multiple protein factors, and then it is further formed into larger structures of topologically associated domains (TADs) or chromatin contact domains (CCDs), followed by A/B compartments and finally the chromosomal territories (CTs). The genetic variation observed in human population influences the multi-scale structures, posing a question regarding the functional impact of structural variants reflected by the variability of the genes expression patterns. The current methods of evaluating the functional effect include eQTLs analysis which uses statistical testing of influence of variants on spatially close genes. Rarely, non-coding DNA sequence changes are evaluated by their impact on the biomolecular interaction network (BIN) reflecting the cellular interactome that can be analysed by the classical graph-theoretic algorithms. Therefore, in the second part of the review, we introduce the concept of BIN, i.e. a meta-network model of the complete molecular interactome developed by integrating various biological networks. The BIN meta-network model includes DNA-protein binding by the plethora of protein factors as well as chromatin interactions, therefore allowing connection of genomics with the downstream biomolecular processes present in a cell. As an illustration, we scrutinise the chromatin interactions mediated by the CTCF protein detected in a ChIA-PET experiment in the human lymphoblastoid cell line GM12878. In the corresponding BIN meta-network the DNA spatial proximity is represented as a graph model, combined with the Proteins-Interaction Network (PIN) of human proteome using the Gene Association Network (GAN). Furthermore, we enriched the BIN with the signalling and metabolic pathways and Gene Ontology (GO) terms to assert its functional context. Finally, we mapped the Single Nucleotide Polymorphisms (SNPs) from the GWAS studies and identified the chromatin mutational hot-spots associated with a significant enrichment of SNPs related to autoimmune diseases. Afterwards, we mapped Structural Variants (SVs) from healthy individuals of 1000 Genomes Project and identified an interesting example of the missing protein complex associated with protein Q6GYQ0 due to a deletion on chromosome 14. Such an analysis using the meta-network BIN model is therefore helpful in evaluating the influence of genetic variation on spatial organisation of the genome and its functional effect in a cell.

Understanding three-dimensional chromatin organization in diploid genomes

The three-dimensional (3D) organization of chromatin in the nucleus of diploid eukaryotic organisms has fascinated biologists for many years. Despite major progress in chromatin conformation studies, current knowledge regarding the spatial organization of diploid (maternal and paternal) genomes is still limited. Recent advances in Hi-C technology and data processing approaches have enabled construction of diploid Hi-C contact maps. These maps greatly accelerated the pace of novel discoveries in haplotype-resolved 3D genome studies, revealing the role of allele biased chromatin conformation in transcriptional regulation. Here, we review emerging concepts and haplotype phasing strategies of Hi-C data in 3D diploid genome studies. We discuss new insights on homologous chromosomal organization and the interplay between allelic biased chromatin architecture and several nuclear functions, explaining how haplotype-resolved Hi-C technologies have been used to resolve important biological questions.

INT-Hi-C reveals distinct chromatin architecture in endosperm and leaf tissues of Arabidopsis

Higher-order chromatin structure undergoes striking changes in response to various developmental and environmental signals, causing distinct cell types to adopt specific chromatin organization. High throughput chromatin conformation capture (Hi-C) allows studying higher-order chromatin structure; however, this technique requires substantial amounts of starting material, which has limited the establishment of cell type-specific higher-order chromatin structure in plants. To overcome this limitation, we established a protocol that is applicable to a limited amount of nuclei by combining the INTACT (isolation of nuclei tagged in specific cell types) method and Hi-C (INT-Hi-C). Using this INT-Hi-C protocol, we generated Hi-C data from INTACT purified endosperm and leaf nuclei. Our INT-Hi-C data from leaf accurately reiterated chromatin interaction patterns derived from conventional leaf Hi-C data. We found that the higher-order chromatin organization of mixed leaf tissues and endosperm differs and that DNA methylation and repressive histone marks positively correlate with the chromatin compaction level. We furthermore found that self-looped interacting genes have increased expression in leaves and endosperm and that interacting intergenic regions negatively impact on gene expression in the endosperm. Last, we identified several imprinted genes involved in long-range and trans interactions exclusively in endosperm. Our study provides evidence that the endosperm adopts a distinct higher-order chromatin structure that differs from other cell types in plants and that chromatin interactions influence transcriptional activity.

MiRNAs: A Powerful Tool in Deciphering Gynecological Malignancies

Accumulated evidence on the clinical roles of microRNAs (miRNAs) in cancer prevention and control has revealed the emergence of new genetic techniques that have improved the understanding of the mechanisms essential for pathology induction and progression. Comprehension of the modifications and individual differences of miRNAs and their interactions in the pathogenesis of gynecological malignancies, together with an understanding of the phenotypic variations have considerably improved the management of the diagnosis and personalized treatment for different forms of cancer. In recent years, miRNAs have emerged as signaling molecules in biological pathways involved in different categories of cancer and it has been demonstrated that these molecules could regulate cancer-relevant processes, our focus being on malignancies of the gynecologic tract. The aim of this paper is to summarize novel research findings in the literature regarding the parts that miRNAs play in cancer-relevant processes, specifically regarding gynecological malignancy, while emphasizing their pivotal role in the disruption of cancer-related signaling pathways.

The Prospective Study of Epigenetic Regulatory Profiles in Sport and Exercise Monitored Through Chromosome Conformation Signatures

The integration of genetic and environmental factors that regulate the gene expression patterns associated with exercise adaptation is mediated by epigenetic mechanisms. The organisation of the human genome within three-dimensional space, known as chromosome conformation, has recently been shown as a dynamic epigenetic regulator of gene expression, facilitating the interaction of distal genomic regions due to tight and regulated packaging of chromosomes in the cell nucleus. Technological advances in the study of chromosome conformation mean a new class of biomarker-the chromosome conformation signature (CCS)-can identify chromosomal interactions across several genomic loci as a collective marker of an epigenomic state. Investigative use of CCSs in biological and medical research shows promise in identifying the likelihood that a disease state is present or absent, as well as an ability to prospectively stratify individuals according to their likely response to medical intervention. The association of CCSs with gene expression patterns suggests that there are likely to be CCSs that respond, or regulate the response, to exercise and related stimuli. The present review provides a contextual background to CCS research and a theoretical framework discussing the potential uses of this novel epigenomic biomarker within sport and exercise science and medicine.

Sci-Hi-C: A single-cell Hi-C method for mapping 3D genome organization in large number of single cells

The highly dynamic nature of chromosome conformation and three-dimensional (3D) genome organization leads to cell-to-cell variability in chromatin interactions within a cell population, even if the cells of the population appear to be functionally homogeneous. Hence, although Hi-C is a powerful tool for mapping 3D genome organization, this heterogeneity of chromosome higher order structure among individual cells limits the interpretive power of population based bulk Hi-C assays. Moreover, single-cell studies have the potential to enable the identification and characterization of rare cell populations or cell subtypes in a heterogeneous population. However, it may require surveying relatively large numbers of single cells to achieve statistically meaningful observations in single-cell studies. By applying combinatorial cellular indexing to chromosome conformation capture, we developed single-cell combinatorial indexed Hi-C (sci-Hi-C), a high throughput method that enables mapping chromatin interactomes in large number of single cells. We demonstrated the use of sci-Hi-C data to separate cells by karytoypic and cell-cycle state differences and to identify cellular variability in mammalian chromosomal conformation. Here, we provide a detailed description of method design and step-by-step working protocols for sci-Hi-C.

Towards mapping the 3D genome through high speed single-molecule tracking of functional transcription factors in single living cells

How genomic DNA is organized in the nucleus is a long-standing question. We describe a single-molecule bioimaging method utilizing super-localization precision coupled to fully quantitative image analysis tools, towards determining snapshots of parts of the 3D genome architecture of model eukaryote budding yeast Saccharomyces cerevisiae with exceptional millisecond time resolution. We employ astigmatism imaging to enable robust extraction of 3D position data on genomically encoded fluorescent protein reporters that bind to DNA. Our relatively straightforward method enables snippets of 3D architectures of likely single genome conformations to be resolved captured via DNA-sequence specific binding proteins in single functional living cells.

Understanding the 3D genome: Emerging impacts on human disease

Recent burst of new technologies that allow for quantitatively delineating chromatin structure has greatly expanded our understanding of how the genome is organized in the three-dimensional (3D) space of the nucleus. It is now clear that the hierarchical organization of the eukaryotic genome critically impacts nuclear activities such as transcription, replication, as well as cellular and developmental events such as cell cycle, cell fate decision and embryonic development. In this review, we discuss new insights into how the structural features of the 3D genome hierarchy are established and maintained, how this hierarchy undergoes dynamic rearrangement during normal development and how its perturbation will lead to human disease, highlighting the accumulating evidence that links the diverse 3D genome architecture components to a multitude of human diseases and the emerging mechanisms by which 3D genome derangement causes disease phenotypes.

An Overview of Methods for Reconstructing 3-D Chromosome and Genome Structures from Hi-C Data

Over the past decade, methods for predicting three-dimensional (3-D) chromosome and genome structures have proliferated. This has been primarily due to the development of high-throughput, next-generation chromosome conformation capture (3C) technologies, which have provided next-generation sequencing data about chromosome conformations in order to map the 3-D genome structure. The introduction of the Hi-C technique-a variant of the 3C method-has allowed researchers to extract the interaction frequency (IF) for all loci of a genome at high-throughput and at a genome-wide scale. In this review we describe, categorize, and compare the various methods developed to map chromosome and genome structures from 3C data-particularly Hi-C data. We summarize the improvements introduced by these methods, describe the approach used for method evaluation, and discuss how these advancements shape the future of genome structure construction.

Mammalian Sex Chromosome Structure, Gene Content, and Function in Male Fertility

Mammalian sex chromosomes evolved from an ordinary pair of autosomes. The X chromosome is highly conserved, whereas the Y chromosome varies among species in size, structure, and gene content. Unlike autosomes that contain randomly mixed collections of genes, the sex chromosomes are enriched in testis-biased genes related to sexual development and reproduction, particularly in spermatogenesis and male fertility. This review focuses on how sex chromosome dosage compensation takes place and why meiotic sex chromosome inactivation occurs during spermatogenesis. Furthermore, the review also emphasizes how testis-biased genes are enriched on the sex chromosomes and their functions in male fertility. It is concluded that sex chromosomes are critical to sexual development and male fertility; however, our understanding of how sex chromosome genes direct sexual development and fertility has been hampered by the structural complexities of the sex chromosomes and by the multicopy nature of the testis gene families that also play a role in immunity, cancer development, and brain function.

Three-dimensional analysis reveals altered chromatin interaction by enhancer inhibitors harbors TCF7L2-regulated cancer gene signature

Distal regulatory elements influence the activity of gene promoters through chromatin looping. Chromosome conformation capture (3C) methods permit identification of chromatin contacts across different regions of the genome. However, due to limitations in the resolution of these methods, the detection of functional chromatin interactions remains a challenge. In the current study, we employ an integrated approach to define and characterize the functional chromatin contacts of human pancreatic cancer cells. We applied tethered chromatin capture to define classes of chromatin domains on a genome-wide scale. We identified three types of structural domains (topologically associated, boundary, and gap) and investigated the functional relationships of these domains with respect to chromatin state and gene expression. We uncovered six distinct sub-domains associated with epigenetic states. Interestingly, specific epigenetically active domains are sensitive to treatment with histone acetyltransferase (HAT) inhibitors and decrease in H3K27 acetylation levels. To examine whether the subdomains that change upon drug treatment are functionally linked to transcription factor regulation, we compared TCF7L2 chromatin binding and gene regulation to HAT inhibition. We identified a subset of coding RNA genes that together can stratify pancreatic cancer patients into distinct survival groups. Overall, this study describes a process to evaluate the functional features of chromosome architecture and reveals the impact of epigenetic inhibitors on chromosome architecture and identifies genes that may provide insight into disease outcome.

Chromatin interactions and expression quantitative trait loci reveal genetic drivers of multimorbidities

Clinical studies of non-communicable diseases identify multimorbidities that suggest a common set of predisposing factors. Despite the fact that humans have ~24,000 genes, we do not understand the genetic pathways that contribute to the development of multimorbid non-communicable disease. Here we create a multimorbidity atlas of traits based on pleiotropy of spatially regulated genes. Using chromatin interaction and expression Quantitative Trait Loci (eQTL) data, we analyse 20,782 variants (p < 5 × 10-6) associated with 1351 phenotypes to identify 16,248 putative spatial eQTL-eGene pairs that are involved in 76,013 short- and long-range regulatory interactions (FDR < 0.05) in different human tissues. Convex biclustering of spatial eGenes that are shared among phenotypes identifies complex interrelationships between nominally different phenotype-associated SNPs. Our approach enables the simultaneous elucidation of variant interactions with target genes that are drivers of multimorbidity, and those that contribute to unique phenotype associated characteristics.

Effective normalization for copy number variation in Hi-C data

BACKGROUND

Normalization is essential to ensure accurate analysis and proper interpretation of sequencing data, and chromosome conformation capture data such as Hi-C have particular challenges. Although several methods have been proposed, the most widely used type of normalization of Hi-C data usually casts estimation of unwanted effects as a matrix balancing problem, relying on the assumption that all genomic regions interact equally with each other.

RESULTS

In order to explore the effect of copy-number variations on Hi-C data normalization, we first propose a simulation model that predict the effects of large copy-number changes on a diploid Hi-C contact map. We then show that the standard approaches relying on equal visibility fail to correct for unwanted effects in the presence of copy-number variations. We thus propose a simple extension to matrix balancing methods that model these effects. Our approach can either retain the copy-number variation effects (LOIC) or remove them (CAIC). We show that this leads to better downstream analysis of the three-dimensional organization of rearranged genomes.

CONCLUSIONS

Taken together, our results highlight the importance of using dedicated methods for the analysis of Hi-C cancer data. Both CAIC and LOIC methods perform well on simulated and real Hi-C data sets, each fulfilling different needs.

Using DNase Hi-C techniques to map global and local three-dimensional genome architecture at high resolution

The folding and three-dimensional (3D) organization of chromatin in the nucleus critically impacts genome function. The past decade has witnessed rapid advances in genomic tools for delineating 3D genome architecture. Among them, chromosome conformation capture (3C)-based methods such as Hi-C are the most widely used techniques for mapping chromatin interactions. However, traditional Hi-C protocols rely on restriction enzymes (REs) to fragment chromatin and are therefore limited in resolution. We recently developed DNase Hi-C for mapping 3D genome organization, which uses DNase I for chromatin fragmentation. DNase Hi-C overcomes RE-related limitations associated with traditional Hi-C methods, leading to improved methodological resolution. Furthermore, combining this method with DNA capture technology provides a high-throughput approach (targeted DNase Hi-C) that allows for mapping fine-scale chromatin architecture at exceptionally high resolution. Hence, targeted DNase Hi-C will be valuable for delineating the physical landscapes of cis-regulatory networks that control gene expression and for characterizing phenotype-associated chromatin 3D signatures. Here, we provide a detailed description of method design and step-by-step working protocols for these two methods.

3C and 3C-based techniques: the powerful tools for spatial genome organization deciphering

It is well known that the chromosomes are organized in the nucleus and this spatial arrangement of genome play a crucial role in gene regulation and genome stability. Different techniques have been developed and applied to uncover the intrinsic mechanism of genome architecture, especially the chromosome conformation capture (3C) and 3C-derived methods. 3C and 3C-derived techniques provide us approaches to perform high-throughput chromatin architecture assays at the genome scale. However, the advantage and disadvantage of current methodologies of C-technologies have not been discussed extensively. In this review, we described and compared the methodologies of C-technologies used in genome organization studies with an emphasis on Hi-C method. We also discussed the crucial challenges facing current genome architecture studies based on 3C and 3C-derived technologies and the direction of future technologies to address currently outstanding questions in the field. These latest news contribute to our current understanding of genome structure, and provide a comprehensive reference for researchers to choose the appropriate method in future application. We consider that these constantly improving technologies will offer a finer and more accurate contact profiles of entire genome and ultimately reveal specific molecular machines govern its shape and function.

Three-dimensional genome architecture and emerging technologies: looping in disease

Genome compaction is a universal feature of cells and has emerged as a global regulator of gene expression. Compaction is maintained by a multitude of architectural proteins, long non-coding RNAs (lncRNAs), and regulatory DNA. Each component comprises interlinked regulatory circuits that organize the genome in three-dimensional (3D) space to manage gene expression. In this review, we update the current state of 3D genome catalogues and focus on how recent technological advances in 3D genomics are leading to an enhanced understanding of disease mechanisms. We highlight the use of genome-wide chromatin conformation capture (Hi-C) coupled with oligonucleotide capture technology (capture Hi-C) to map interactions between gene promoters and distal regulatory elements such as enhancers that are enriched for disease variants from genome-wide association studies (GWASs). We discuss how aberrations in architectural units are associated with various pathological outcomes, and explore how recent advances in genome and epigenome editing show great promise for a systematic understanding of complex genetic disorders. Our growing understanding of 3D genome architecture—coupled with the ability to engineer changes in it—may create novel therapeutic opportunities.

ARNT is a potential direct HIF-1 target gene in human Hep3B hepatocellular carcinoma cells

BACKGROUND

The transcription factor aryl hydrocarbon receptor nuclear translocator (ARNT) participates in the hypoxia-inducible factor (HIF) pathway which senses a decline in cellular oxygen tension. In hypoxia, HIF-1α and ARNT form the transcriptional active complex HIF-1 followed by the expression of target genes. ARNT is considered as constitutively expressed and unaffected by hypoxia. However, certain tumour cell lines derived from different entities are capable to elevate ARNT expression under hypoxic conditions which implies a survival benefit. It was demonstrated that high ARNT protein levels mediate radioresistance in tumour cells. Furthermore, a HIF-1α-driven feed-forward loop leading to augmented HIF signalling was discovered in Hep3B cells. Herein HIF-1α elevates the mRNA and protein expression of its binding partner ARNT in hypoxia. However, the detailed mechanism remained unclear. The objective of this study was to test whether HIF-1α might directly regulate ARNT expression by recruitment to the ARNT promoter.

METHODS

Chromatin immunoprecipitation (ChIP), CRISPR/Cas9 genome editing, Western blotting, quantitative RT-PCR and reporter gene assays were applied. The unpaired t test was used for statistical analysis.

RESULTS

ChIP assays revealed the binding of both HIF-1α and ARNT to the ARNT promoter in hypoxia. The relevance of this particular region for hypoxic ARNT induction was confirmed by CRISPR/Cas9 genome editing. ARNT normoxic basal expression and hypoxic inducibility was reduced in genome-edited Hep3B cells. This phenotype was accompanied with impaired HIF signalling and was rescued by ARNT overexpression.

CONCLUSIONS

The results indicate ARNT to be a putative HIF-1 target gene and a limiting factor in this model.

Using omics approaches to understand pulmonary diseases

Omics approaches are high-throughput unbiased technologies that provide snapshots of various aspects of biological systems and include: 1) genomics, the measure of DNA variation; 2) transcriptomics, the measure of RNA expression; 3) epigenomics, the measure of DNA alterations not involving sequence variation that influence RNA expression; 4) proteomics, the measure of protein expression or its chemical modifications; and 5) metabolomics, the measure of metabolite levels. Our understanding of pulmonary diseases has increased as a result of applying these omics approaches to characterize patients, uncover mechanisms underlying drug responsiveness, and identify effects of environmental exposures and interventions. As more tissue- and cell-specific omics data is analyzed and integrated for diverse patients under various conditions, there will be increased identification of key mechanisms that underlie pulmonary biological processes, disease endotypes, and novel therapeutics that are efficacious in select individuals. We provide a synopsis of how omics approaches have advanced our understanding of asthma, chronic obstructive pulmonary disease (COPD), acute respiratory distress syndrome (ARDS), idiopathic pulmonary fibrosis (IPF), and pulmonary arterial hypertension (PAH), and we highlight ongoing work that will facilitate pulmonary disease precision medicine.

Large-scale chromosome folding versus genomic DNA sequences: A discrete double Fourier transform technique

Using state-of-the-art techniques combining imaging methods and high-throughput genomic mapping tools leaded to the significant progress in detailing chromosome architecture of various organisms. However, a gap still remains between the rapidly growing structural data on the chromosome folding and the large-scale genome organization. Could a part of information on the chromosome folding be obtained directly from underlying genomic DNA sequences abundantly stored in the databanks? To answer this question, we developed an original discrete double Fourier transform (DDFT). DDFT serves for the detection of large-scale genome regularities associated with domains/units at the different levels of hierarchical chromosome folding. The method is versatile and can be applied to both genomic DNA sequences and corresponding physico-chemical parameters such as base-pairing free energy. The latter characteristic is closely related to the replication and transcription and can also be used for the assessment of temperature or supercoiling effects on the chromosome folding. We tested the method on the genome of E. coli K-12 and found good correspondence with the annotated domains/units established experimentally. As a brief illustration of further abilities of DDFT, the study of large-scale genome organization for bacteriophage PHIX174 and bacterium Caulobacter crescentus was also added. The combined experimental, modeling, and bioinformatic DDFT analysis should yield more complete knowledge on the chromosome architecture and genome organization.

Review: Regulation of the cancer epigenome by long non-coding RNAs

Long non-coding RNAs have emerged as highly versatile players in the regulation of gene expression in development and human disease, particularly cancer. Hundreds of lncRNAs become dysregulated across tumor types, and multiple lncRNAs have demonstrated functions as tumor-suppressors or oncogenes. Furthermore, studies have demonstrated that dysregulation of lncRNAs results in alterations of the epigenome in cancer cells, potentially providing a novel mechanism for the massive epigenomic alterations observed in many tumors. Here, we highlight and provide some illustrious examples of lncRNAs in various epigenetic regulatory processes, including coordination of chromatin dynamics, regulation of DNA methylation, modulation of other non-coding RNAs and mRNA stability, and control of epigenetic substrate availability through altered tumor metabolism. In light of all these known and emerging functions in epigenetic regulation of tumorigenesis and cancer progression, lncRNAs represent attractive targets for future therapeutic strategies in cancer.

Network Reconstruction Reveals that Valproic Acid Activates Neurogenic Transcriptional Programs in Adult Brain Following Traumatic Injury

Objectives

To determine the mechanism of action of valproic acid (VPA) in the adult central nervous system (CNS) following traumatic brain injury (TBI) and hemorrhagic shock (HS).

Methods

Data were analyzed from different sources, including experiments in a porcine model, data from postmortem human brain, published studies, public and commercial databases.

Results

The transcriptional program in the CNS following TBI, HS, and VPA treatment includes activation of regulatory pathways that enhance neurogenesis and suppress gliogenesis. Genes which encode the transcription factors (TFs) that specify neuronal cell fate, including MEF2D, MYT1L, NEUROD1, PAX6 and TBR1, and their target genes, are induced by VPA. VPA represses genes responsible for oligodendrogenesis, maintenance of white matter, T-cell activation, angiogenesis, and endothelial cell proliferation, adhesion and chemotaxis. NEUROD1 has regulatory interactions with 38% of the genes regulated by VPA in a swine model of TBI and HS in adult brain. Hi-C spatial mapping of a VPA pharmacogenomic SNP in the GRIN2B gene shows it is part of a transcriptional hub that contacts 12 genes that mediate chromatin-mediated neurogenesis and neuroplasticity.

Conclusions

Following TBI and HS, this study shows that VPA administration acts in the adult brain through differential activation of TFs responsible for neurogenesis, genes responsible for neuroplasticity, and repression of TFs that specify oligodendrocyte cell fate, endothelial cell chemotaxis and angiogenesis.

Short title: Mechanism of action of valproic acid in traumatic brain injury

Electronic supplementary material

The online version of this article (doi:10.1007/s11095-017-2130-6) contains supplementary material, which is available to authorized users.

Massively multiplex single-cell Hi-C

We present single-cell combinatorial indexed Hi-C (sciHi-C), a method that applies combinatorial cellular indexing to chromosome conformation capture. In this proof of concept, we generate and sequence six sciHi-C libraries comprising a total of 10,696 single cells. We use sciHi-C data to separate cells by karyotypic and cell-cycle state differences and identify cell-to-cell heterogeneity in mammalian chromosomal conformation. Our results demonstrate that combinatorial indexing is a generalizable strategy for single-cell genomics.

Cajal body function in genome organization and transcriptome diversity

Nuclear bodies contribute to non-random organization of the human genome and nuclear function. Using a major prototypical nuclear body, the Cajal body, as an example, we suggest that these structures assemble at specific gene loci located across the genome as a result of high transcriptional activity. Subsequently, target genes are physically clustered in close proximity in Cajal body-containing cells. However, Cajal bodies are observed in only a limited number of human cell types, including neuronal and cancer cells. Ultimately, Cajal body depletion perturbs splicing kinetics by reducing target small nuclear RNA (snRNA) transcription and limiting the levels of spliceosomal snRNPs, including their modification and turnover following each round of RNA splicing. As such, Cajal bodies are capable of shaping the chromatin interaction landscape and the transcriptome by influencing spliceosome kinetics. Future studies should concentrate on characterizing the direct influence of Cajal bodies upon snRNA gene transcriptional dynamics. Also see the video abstract here.