Genome organization and long-range regulation of gene expression by enhancers

Unraveling the three-dimensional (3D) genome architecture in Neurodevelopmental Disorders (NDDs)

The human genome, comprising millions of pairs of bases, serves as the blueprint of life, encoding instructions for cellular processes. However, genomes are not merely linear sequences; rather, the complex of DNA and histones, known as chromatin, exhibits complex organization across various levels, which profoundly influence gene expression and cellular function. Central to understanding genome organization is the emerging field of three-dimensional (3D) genome studies. Utilizing advanced techniques such as Hi-C, researchers have unveiled non-random dispositions of genomic elements, highlighting their importance in transcriptional regulation and disease mechanisms. Topologically Associating Domains (TADs), that demarcate regions of chromatin with preferential internal interactions, play crucial roles in gene regulation and are increasingly implicated in various diseases such as cancer and schizophrenia. However, their role in Neurodevelopmental Disorders (NDDs) remains poorly understood. Here, we focus on TADs and 3D conservation across the evolution and between cell types in NDDs. The investigation into genome organization and its impact on disease has led to significant breakthroughs in understanding NDDs etiology such ASD (Autism Spectrum Disorder). By elucidating the wide spectrum of ASD manifestations, researchers aim to uncover the underlying genetic and epigenetic factors contributing to its heterogeneity. Moreover, studies linking TAD disruption to NDDs underscore the importance of spatial genome organization in maintaining proper brain development and function. In summary, this review highlights the intricate interplay between genome organization, transcriptional control, and disease pathology, shedding light on fundamental biological processes and offering insights into the mechanisms underlying NDDs like ASD.

The importance of functional genomics studies in precision rheumatology

Rheumatic diseases, those that affect the musculoskeletal system, cause significant morbidity. Among risk factors of these diseases is a significant genetic component. Recent advances in high-throughput omics techniques now allow a comprehensive profiling of patients at a genetic level through genome-wide association studies. Without functional interpretation of variants identified through these studies, clinical insight remains limited. Strategies include statistical fine-mapping that refine the list of variants in loci associated with disease, whilst colocalization techniques attempt to attribute function to variants that overlap a genetically active chromatin annotation. Functional validation using genome editing techniques can be used to further refine genetic signals and identify key pathways in cell types relevant to rheumatic disease biology. Insight gained from the combination of genetic studies and functional validation can be used to improve precision medicine in rheumatic diseases by allowing risk prediction and drug repositioning.

Dynamic RNA polymerase II occupancy drives differentiation of the intestine under the direction of HNF4

Terminal differentiation requires massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4,000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identify dynamic occupancy of RNA polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II occupancy and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of-function, chromatin immunoprecipitation sequencing (ChIP-seq), and immunoprecipitation (IP) mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II occupancy at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms that drive differentiation gene expression and find that pause-release of Pol II and post-transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in renewing adult tissue.

Regulus infers signed regulatory relations from few samples' information using discretization and likelihood constraints

MOTIVATION

Transcriptional regulation is performed by transcription factors (TF) binding to DNA in context-dependent regulatory regions and determines the activation or inhibition of gene expression. Current methods of transcriptional regulatory circuits inference, based on one or all of TF, regions and genes activity measurements require a large number of samples for ranking the candidate TF-gene regulation relations and rarely predict whether they are activations or inhibitions. We hypothesize that transcriptional regulatory circuits can be inferred from fewer samples by (1) fully integrating information on TF binding, gene expression and regulatory regions accessibility, (2) reducing data complexity and (3) using biology-based likelihood constraints to determine the global consistency between a candidate TF-gene relation and patterns of genes expressions and region activations, as well as qualify regulations as activations or inhibitions.

RESULTS

We introduce Regulus, a method which computes TF-gene relations from gene expressions, regulatory region activities and TF binding sites data, together with the genomic locations of all entities. After aggregating gene expressions and region activities into patterns, data are integrated into a RDF (Resource Description Framework) endpoint. A dedicated SPARQL (SPARQL Protocol and RDF Query Language) query retrieves all potential relations between expressed TF and genes involving active regulatory regions. These TF-region-gene relations are then filtered using biological likelihood constraints allowing to qualify them as activation or inhibition. Regulus provides signed relations consistent with public databases and, when applied to biological data, identifies both known and potential new regulators. Regulus is devoted to context-specific transcriptional circuits inference in human settings where samples are scarce and cell populations are closely related, using discretization into patterns and likelihood reasoning to decipher the most robust regulatory relations.

Dynamic RNA Polymerase II Recruitment Drives Differentiation of the Intestine under the direction of HNF4

Terminal differentiation requires a massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identified dynamic recruitment of RNA Polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II recruitment and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of- function, ChIP-seq and IP mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II recruitment at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms which drive differentiation gene expression and find pause-release of Pol II and post- transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in a renewing adult tissue.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

De novo assembly and annotation of the singing mouse genome

BACKGROUND

Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions.

RESULTS

Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types.

CONCLUSIONS

These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

Multiomic analysis of cohesin reveals that ZBTB transcription factors contribute to chromatin interactions

One bottleneck in understanding the principles of 3D chromatin structures is caused by the paucity of known regulators. Cohesin is essential for 3D chromatin organization, and its interacting partners are candidate regulators. Here, we performed proteomic profiling of the cohesin in chromatin and identified transcription factors, RNA-binding proteins and chromatin regulators associated with cohesin. Acute protein degradation followed by time-series genomic binding quantitation and BAT Hi-C analysis were conducted, and the results showed that the transcription factor ZBTB21 contributes to cohesin chromatin binding, 3D chromatin interactions and transcriptional repression. Strikingly, multiomic analyses revealed that the other four ZBTB factors interacted with cohesin, and double degradation of ZBTB21 and ZBTB7B led to a further decrease in cohesin chromatin occupancy. We propose that multiple ZBTB transcription factors orchestrate the chromatin binding of cohesin to regulate chromatin interactions, and we provide a catalog of many additional proteins associated with cohesin that warrant further investigation.

Recent advances in chromosome capture techniques unraveling 3D genome architecture in germ cells, health, and disease

In eukaryotes, the genome does not emerge in a specific shape but rather as a hierarchial bundle within the nucleus. This multifaceted genome organization consists of multiresolution cellular structures, such as chromosome territories, compartments, and topologically associating domains, which are frequently defined by architecture, design proteins including CTCF and cohesin, and chromatin loops. This review briefly discusses the advances in understanding the basic rules of control, chromatin folding, and functional areas in early embryogenesis. With the use of chromosome capture techniques, the latest advancements in technologies for visualizing chromatin interactions come close to revealing 3D genome formation frameworks with incredible detail throughout all genomic levels, including at single-cell resolution. The possibility of detecting variations in chromatin architecture might open up new opportunities for disease diagnosis and prevention, infertility treatments, therapeutic approaches, desired exploration, and many other application scenarios.

Dynamic enhancer interactome promotes senescence and aging

Gene expression programs are regulated by enhancers which act in a context-specific manner, and can reside at great distances from their target genes. Extensive three-dimensional (3D) genome reorganization occurs in senescence, but how enhancer interactomes are reconfigured during this process is just beginning to be understood. Here we generated high-resolution contact maps of active enhancers and their target genes, assessed chromatin accessibility, and established one-dimensional maps of various histone modifications and transcription factors to comprehensively understand the regulation of enhancer configuration during senescence. Hyper-connected enhancer communities/cliques formed around genes that are highly expressed and within essential gene pathways in each cell state. In addition, motif analysis indicates the involvement of specific transcription factors in hyper-connected regulatory elements in each condition; importantly, MafK, a bZIP family transcription factor, was upregulated in senescence, and reduced expression of MafK ameliorated the senescence phenotypes. Because the accumulation of senescent cells is a key feature of aging, we further investigated enhancer connectomes in the liver of young and aged mice. Hyper-connected enhancer communities were identified during aging, which regulate essential genes that maintain cell differentiation and homeostasis. These findings reveal that hyper-connected enhancer communities correlate with high gene expression in senescence and aging and provide potential hotspots for therapeutic intervention in aging and age-associated diseases.

Chondrocyte De-Differentiation: Biophysical Cues to Nuclear Alterations

Autologous chondrocyte implantation (ACI) is a cell therapy to repair cartilage defects. In ACI a biopsy is taken from a non-load bearing area of the knee and expanded in-vitro. The expansion process provides the benefit of generating a large number of cells required for implantation; however, during the expansion these cells de-differentiate and lose their chondrocyte phenotype. In this review we focus on examining the de-differentiation phenotype from a mechanobiology and biophysical perspective, highlighting some of the nuclear mechanics and chromatin changes in chondrocytes seen during the expansion process and how this relates to the gene expression profile. We propose that manipulating chondrocyte nuclear architecture and chromatin organization will highlight mechanisms that will help to preserve the chondrocyte phenotype.

Insights into the cell fate decision-making processes from chromosome structural reorganizations

The cell fate decision-making process, which provides the capability of a cell transition to a new cell type, involves the reorganizations of 3D genome structures. Currently, the high temporal resolution picture of how the chromosome structural rearrangements occur and further influence the gene activities during the cell-state transition is still challenging to acquire. Here, we study the chromosome structural reorganizations during the cell-state transitions among the pluripotent embryonic stem cell, the terminally differentiated normal cell, and the cancer cell using a nonequilibrium landscape-switching model implemented in the molecular dynamics simulation. We quantify the chromosome (de)compaction pathways during the cell-state transitions and find that the two pathways having the same destinations can merge prior to reaching the final states. The chromosomes at the merging states have similar structural geometries but can differ in long-range compartment segregation and spatial distribution of the chromosomal loci and genes, leading to cell-type-specific transition mechanisms. We identify the irreversible pathways of chromosome structural rearrangements during the forward and reverse transitions connecting the same pair of cell states, underscoring the critical roles of nonequilibrium dynamics in the cell-state transitions. Our results contribute to the understanding of the cell fate decision-making processes from the chromosome structural perspective.

Deciphering aging at three-dimensional genomic resolution

Aging is characterized by progressive functional declines at the organismal, organic, and cellular levels and increased susceptibility to aging-related diseases. Epigenetic alteration is a hallmark of aging, senescent cells show epigenomic changes at multiple scales, such as 3D genome reorganization, alterations of histone modifications and chromatin accessibility, and DNA hypomethylation. Chromosome conformation capture (3C)-based technologies have enabled the generation of key information on genomic reorganizations during senescence. A comprehensive understanding of epigenomic alterations during aging will yield important insights into the underlying epigenetic mechanism for aging regulation, the identification of aging-related biomarkers, and the development of potential aging intervention targets.

Long-Distance Repression by Human Silencers: Chromatin Interactions and Phase Separation in Silencers

Three-dimensional genome organization represents an additional layer in the epigenetic regulation of gene expression. Active transcription controlled by enhancers or super-enhancers has been extensively studied. Enhancers or super-enhancers can recruit activators or co-activators to activate target gene expression through long-range chromatin interactions. Chromatin interactions and phase separation play important roles in terms of enhancer or super-enhancer functioning. Silencers are another major type of cis-regulatory element that can mediate gene regulation by turning off or reducing gene expression. However, compared to active transcription, silencer studies are still in their infancy. This review covers the current knowledge of human silencers, especially the roles of chromatin interactions and phase separation in silencers. This review also proposes future directions for human silencer studies.

Deciphering the Retinal Epigenome during Development, Disease and Reprogramming: Advancements, Challenges and Perspectives

Retinal neurogenesis is driven by concerted actions of transcription factors, some of which are expressed in a continuum and across several cell subtypes throughout development. While seemingly redundant, many factors diversify their regulatory outcome on gene expression, by coordinating variations in chromatin landscapes to drive divergent retinal specification programs. Recent studies have furthered the understanding of the epigenetic contribution to the progression of age-related macular degeneration, a leading cause of blindness in the elderly. The knowledge of the epigenomic mechanisms that control the acquisition and stabilization of retinal cell fates and are evoked upon damage, holds the potential for the treatment of retinal degeneration. Herein, this review presents the state-of-the-art approaches to investigate the retinal epigenome during development, disease, and reprogramming. A pipeline is then reviewed to functionally interrogate the epigenetic and transcriptional networks underlying cell fate specification, relying on a truly unbiased screening of open chromatin states. The related work proposes an inferential model to identify gene regulatory networks, features the first footprinting analysis and the first tentative, systematic query of candidate pioneer factors in the retina ever conducted in any model organism, leading to the identification of previously uncharacterized master regulators of retinal cell identity, such as the nuclear factor I, NFI. This pipeline is virtually applicable to the study of genetic programs and candidate pioneer factors in any developmental context. Finally, challenges and limitations intrinsic to the current next-generation sequencing techniques are discussed, as well as recent advances in super-resolution imaging, enabling spatio-temporal resolution of the genome.

Integrative analysis of the 3D genome structure reveals that CTCF maintains the properties of mouse female germline stem cells

The three-dimensional configuration of the genome ensures cell type-specific gene expression profiles by placing genes and regulatory elements in close spatial proximity. Here, we used in situ high-throughput chromosome conformation (in situ Hi-C), RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) to characterize the high-order chromatin structure signature of female germline stem cells (FGSCs) and identify its regulating key factor based on the data-driven of multiple omics data. By comparison with pluripotent stem cells (PSCs), adult stem cells (ASCs), and somatic cells at three major levels of chromatin architecture, A/B compartments, topologically associating domains, and chromatin loops, the chromatin architecture of FGSCs was most similar to that of other ASCs and largely different from that of PSCs and somatic cells. After integrative analysis of the three-dimensional chromatin structure, active compartment-associating loops (aCALs) were identified as a signature of high-order chromatin organization in FGSCs, which revealed that CCCTC-binding factor was a major factor to maintain the properties of FGSCs through regulation of aCALs. We found FGSCs belong to ASCs at chromatin structure level and characterized aCALs as the high-order chromatin structure signature of FGSCs. Furthermore, CTCF was identified to play a key role in regulating aCALS to maintain the biological functions of FGSCs. These data provide a valuable resource for future studies of the features of chromatin organization in mammalian stem cells and further understanding of the fundamental characteristics of FGSCs.

Exploring the role of eRNA in regulating gene expression

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.e., the combination of Type-I and Type-II). We show that with the increase of the E-P loop length, mRNA distribution can transition from unimodality to bimodality or vice versa in all the three regulation cases. However, in contrast to the other two regulations, Type-II regulation can lead to the highest mean mRNA level and the lowest mRNA noise, independent of the E-P loop length. These results would not only reveal the essential mechanism of how eRNA regulates gene expression, but also imply a new mechanism for phenotypic switching, namely the E-P loop can induce phenotypic switching.

Unraveling three-dimensional chromatin structural dynamics during spermatogonial differentiation

Spermatogonial stem cells (SSCs) are able to undergo both self-renewal and differentiation. Unlike self-renewal, which replenishes the SSC and progenitor pool, differentiation is an irreversible process committing cells to meiosis. Although the preparations for meiotic events in differentiating spermatogonia (Di-SG) are likely to be accompanied by alterations in chromatin structure, the three-dimensional chromatin architectural differences between SSCs and Di-SG, and the higher-order chromatin dynamics during spermatogonial differentiation, have not been systematically investigated. Here, we performed in situ high-throughput chromosome conformation capture, RNA-seq, and chromatin immunoprecipitation-sequencing analyses on porcine undifferentiated spermatogonia (which consist of SSCs and progenitors) and Di-SG. We identified that Di-SG exhibited less compact chromatin structural organization, weakened compartmentalization, and diminished topologically associating domains in comparison with undifferentiated spermatogonia, suggesting that diminished higher-order chromatin architecture in meiotic cells, as shown by recent reports, might be preprogrammed in Di-SG. Our data also revealed that A/B compartments, representing open or closed chromatin regions respectively, and topologically associating domains were related to dynamic gene expression during spermatogonial differentiation. Furthermore, we unraveled the contribution of promoter-enhancer interactions to premeiotic transcriptional regulation, which has not been accomplished in previous studies due to limited cell input and resolution. Together, our study uncovered the three-dimensional chromatin structure of SSCs/progenitors and Di-SG, as well as the interplay between higher-order chromatin architecture and dynamic gene expression during spermatogonial differentiation. These findings provide novel insights into the mechanisms for SSC self-renewal and differentiation and have implications for diagnosis and treatment of male sub-/infertility.

Understanding 3D Genome Organization and Its Effect on Transcriptional Gene Regulation Under Environmental Stress in Plant: A Chromatin Perspective

The genome of a eukaryotic organism is comprised of a supra-molecular complex of chromatin fibers and intricately folded three-dimensional (3D) structures. Chromosomal interactions and topological changes in response to the developmental and/or environmental stimuli affect gene expression. Chromatin architecture plays important roles in DNA replication, gene expression, and genome integrity. Higher-order chromatin organizations like chromosome territories (CTs), A/B compartments, topologically associating domains (TADs), and chromatin loops vary among cells, tissues, and species depending on the developmental stage and/or environmental conditions (4D genomics). Every chromosome occupies a separate territory in the interphase nucleus and forms the top layer of hierarchical structure (CTs) in most of the eukaryotes. While the A and B compartments are associated with active (euchromatic) and inactive (heterochromatic) chromatin, respectively, having well-defined genomic/epigenomic features, TADs are the structural units of chromatin. Chromatin architecture like TADs as well as the local interactions between promoter and regulatory elements correlates with the chromatin activity, which alters during environmental stresses due to relocalization of the architectural proteins. Moreover, chromatin looping brings the gene and regulatory elements in close proximity for interactions. The intricate relationship between nucleotide sequence and chromatin architecture requires a more comprehensive understanding to unravel the genome organization and genetic plasticity. During the last decade, advances in chromatin conformation capture techniques for unravelling 3D genome organizations have improved our understanding of genome biology. However, the recent advances, such as Hi-C and ChIA-PET, have substantially increased the resolution, throughput as well our interest in analysing genome organizations. The present review provides an overview of the historical and contemporary perspectives of chromosome conformation capture technologies, their applications in functional genomics, and the constraints in predicting 3D genome organization. We also discuss the future perspectives of understanding high-order chromatin organizations in deciphering transcriptional regulation of gene expression under environmental stress (4D genomics). These might help design the climate-smart crop to meet the ever-growing demands of food, feed, and fodder.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Offspring production of ovarian organoids derived from spermatogonial stem cells by defined factors with chromatin reorganization

Introduction

Fate determination of germline stem cells remains poorly understood at the chromatin structure level.

Objectives

Our research hopes to develop successful offspring production of ovarian organoids derived from spermatogonial stem cells (SSCs) by defined factors.

Methods

The offspring production from oocytes transdifferentiated from mouse SSCs with tracking of transplanted SSCs in vivo, single cell whole exome sequencing, and in 3D cell culture reconstitution of the process of oogenesis derived from SSCs. The defined factors were screened with ovarian organoids. We uncovered extensive chromatin reorganization during SSC conversion into induced germline stem cells (iGSCs) using high throughput chromosome conformation.

Results

We demonstrate successful production of offspring from oocytes transdifferentiated from mouse spermatogonial stem cells (SSCs). Furthermore, we demonstrate direct induction of germline stem cells (iGSCs) differentiated into functional oocytes by transduction of H19, Stella, and Zfp57 and inactivation of Plzf in SSCs after screening with ovarian organoids. We uncovered extensive chromatin reorganization during SSC conversion into iGSCs, which was highly similar to female germline stem cells. We observed that although topologically associating domains were stable during SSC conversion, chromatin interactions changed in a striking manner, altering 35% of inactive and active chromosomal compartments throughout the genome.

Conclusion

We demonstrate successful offspring production of ovarian organoids derived from SSCs by defined factors with chromatin reorganization. These findings have important implications in various areas including mammalian gametogenesis, genetic and epigenetic reprogramming, biotechnology, and medicine.

The role of epigenetic modifications, long-range contacts, enhancers and topologically associating domains in the regulation of glioma grade-specific genes

Genome-wide studies have uncovered specific genetic alterations, transcriptomic patterns and epigenetic profiles associated with different glioma types. We have recently created a unique atlas encompassing genome-wide profiles of open chromatin, histone H3K27ac and H3Kme3 modifications, DNA methylation and transcriptomes of 33 glioma samples of different grades. Here, we intersected genome-wide atlas data with topologically associating domains (TADs) and demonstrated that the chromatin organization and epigenetic landscape of enhancers have a strong impact on genes differentially expressed in WHO low grade versus high grade gliomas. We identified TADs enriched in glioma grade-specific genes and/or epigenetic marks. We found the set of transcription factors, including REST, E2F1 and NFKB1, that are most likely to regulate gene expression in multiple TADs, containing specific glioma-related genes. Moreover, many genes associated with the cell-matrix adhesion Gene Ontology group, in particular 14 PROTOCADHERINs, were found to be regulated by long-range contacts with enhancers. Presented results demonstrate the existence of epigenetic differences associated with chromatin organization driving differential gene expression in gliomas of different malignancy.

Analysis of Major Depression Risk Genes Reveals Evolutionary Conservation, Shared Phenotypes, and Extensive Genetic Interactions

Major depressive disorder (MDD) affects around 15% of the population at some stage in their lifetime. It can be gravely disabling and it is associated with increased risk of suicide. Genetics play an important role; however, there are additional environmental contributions to the pathogenesis. A number of possible risk genes that increase liability for developing symptoms of MDD have been identified in genome-wide association studies (GWAS). The goal of this study was to characterize the MDD risk genes with respect to the degree of evolutionary conservation in simpler model organisms such as Caenorhabditis elegans and zebrafish, the phenotypes associated with variation in these genes and the extent of network connectivity. The MDD risk genes showed higher conservation in C. elegans and zebrafish than genome-to-genome comparisons. In addition, there were recurring themes among the phenotypes associated with variation of these risk genes in C. elegans. The phenotype analysis revealed enrichment for essential genes with pleiotropic effects. Moreover, the MDD risk genes participated in more interactions with each other than did randomly-selected genes from similar-sized gene sets. Syntenic blocks of risk genes with common functional activities were also identified. By characterizing evolutionarily-conserved counterparts to the MDD risk genes, we have gained new insights into pathogenetic processes relevant to the emergence of depressive symptoms in man.

Advantages of Heterotrophic Microalgae as a Host for Phytochemicals Production

Currently, most commercial recombinant technologies rely on host systems. However, each host has their own benefits and drawbacks, depending on the target products. Prokaryote host is lack of post-transcriptional and post-translational mechanisms, making them unsuitable for eukaryotic productions like phytochemicals. Even there are other eukaryote hosts (e.g., transgenic animals, mammalian cell, and transgenic plants), but those hosts have some limitations, such as low yield, high cost, time consuming, virus contamination, and so on. Thus, flexible platforms and efficient methods that can produced phytochemicals are required. The use of heterotrophic microalgae as a host system is interesting because it possibly overcome those obstacles. This paper presents a comprehensive review of heterotrophic microalgal expression host including advantages of heterotrophic microalgae as a host, genetic engineering of microalgae, genetic transformation of microalgae, microalgal engineering for phytochemicals production, challenges of microalgal hosts, key market trends, and future view. Finally, this review might be a directions of the alternative microalgae host for high-value phytochemicals production in the next few years.

Three-dimensional interactions between enhancers and promoters during intestinal differentiation depend upon HNF4

Cells in renewing tissues exhibit dramatic transcriptional changes as they differentiate. The contribution of chromatin looping to tissue renewal is incompletely understood. Enhancer-promoter interactions could be relatively stable as cells transition from progenitor to differentiated states; alternatively, chromatin looping could be as dynamic as the gene expression from their loci. The intestinal epithelium is the most rapidly renewing mammalian tissue. Proliferative cells in crypts of Lieberkühn sustain a stream of differentiated cells that are continually shed into the lumen. We apply chromosome conformation capture combined with chromatin immunoprecipitation (HiChIP) and sequencing to measure enhancer-promoter interactions in progenitor and differentiated cells of the intestinal epithelium. Despite dynamic gene regulation across the differentiation axis, we find that enhancer-promoter interactions are relatively stable. Functionally, we find HNF4 transcription factors are required for chromatin looping at target genes. Depletion of HNF4 disrupts local chromatin looping, histone modifications, and target gene expression. This study provides insights into transcriptional regulatory mechanisms governing homeostasis in renewing tissues.

Chen et al. provide a survey of enhancer-promoter 3D looping in the intestinal epithelium by HiChIP, in vivo. They find that enhancer-promoter interactions are highly dependent upon the key intestinal transcription factor HNF4. Their findings provide insights into transcriptional regulatory mechanisms governing homeostasis in renewing tissues.

Promoter-interacting expression quantitative trait loci are enriched for functional genetic variants

Expression quantitative trait loci (eQTLs) studies provide associations of genetic variants with gene expression but fall short of pinpointing functionally important eQTLs. Here, using H3K27ac HiChIP assays, we mapped eQTLs overlapping active cis-regulatory elements that interact with their target gene promoters (promoter-interacting eQTLs, pieQTLs) in five common immune cell types (Database of Immune Cell Expression, Expression quantitative trait loci and Epigenomics (DICE) cis-interactome project). This approach allowed us to identify functionally important eQTLs and show mechanisms that explain their cell-type restriction. We also devised an approach to eQTL discovery that relies on HiChIP-based promoter interaction maps as a structural framework for deciding which SNPs to test for association with gene expression, and observe ultra-long-distance pieQTLs (>1 megabase away), including several disease-risk variants. We validated the functional role of pieQTLs using reporter assays, CRISPRi, dCas9-tiling guides and Cas9-mediated base-pair editing. In this article we present a method for functional eQTL discovery and provide insights into relevance of noncoding variants for cell-specific gene regulation and for disease association beyond conventional eQTL mapping.

Genomic retargeting of p53 and CTCF is associated with transcriptional changes during oncogenic HRas-induced transformation

Gene transcription is regulated by distant regulatory elements via combinatorial binding of transcription factors. It is increasingly recognized that alterations in chromatin state and transcription factor binding in these distant regulatory elements may have key roles in cancer development. Here we focused on the first stages of oncogene-induced carcinogenic transformation, and characterized the regulatory network underlying transcriptional changes associated with this process. Using Hi-C data, we observe spatial coupling between differentially expressed genes and their differentially accessible regulatory elements and reveal two candidate transcription factors, p53 and CTCF, as determinants of transcriptional alterations at the early stages of oncogenic HRas-induced transformation in human mammary epithelial cells. Strikingly, the malignant transcriptional reprograming is promoted by redistribution of chromatin binding of these factors without major variation in their expression level. Our results demonstrate that alterations in the regulatory landscape have a major role in driving oncogene-induced transcriptional reprogramming.

Microscopic Chromosomal Structural and Dynamical Origin of Cell Differentiation and Reprogramming

As an essential and fundamental process of life, cell development involves large-scale reorganization of the 3D genome architecture, which forms the basis of gene regulation. Here, a landscape-switching model is developed to explore the microscopic chromosomal structural origin of embryonic stem cell (ESC) differentiation and somatic cell reprogramming. It is shown that chromosome structure exhibits significant compartment-switching in the unit of topologically associating domain. It is found that the chromosome during differentiation undergoes monotonic compaction with spatial repositioning of active and inactive chromosomal loci toward the chromosome surface and interior, respectively. In contrast, an overexpanded chromosome, which exhibits universal localization of loci at the chromosomal surface with erasing the structural characteristics formed in the somatic cells, is observed during reprogramming. An early distinct differentiation pathway from the ESC to the terminally differentiated cell, giving rise to early bifurcation on the Waddington landscape for the ESC differentiation is suggested. The theoretical model herein including the non-equilibrium effects, draws a picture of the highly irreversible cell differentiation and reprogramming processes, in line with the experiments. The predictions provide a physical understanding of cell differentiation and reprogramming from the chromosomal structural and dynamical perspective and can be tested by future experiments.

Plant and animal chromatin three-dimensional organization: similar structures but different functions

Chromatin is the main carrier of genetic information and is non-randomly distributed within the nucleus. Next-generation sequence-based chromatin conformation capture technologies have enabled us to directly examine its three-dimensional organization at an unprecedented scale and resolution. In the best-studied mammalian models, chromatin folding can be broken down into three hierarchical levels, compartment, domains, and loops, which play important roles in transcriptional regulation. Although similar structures have now been identified in plants, they might not possess exactly the same functions as the mammalian ones. Here, we review recent Hi-C studies in plants, compare plant chromatin structures with their mammalian counterparts, and discuss the differences between plants with different genome sizes.

Exonuclease combinations reduce noises in 3D genomics technologies

Chromosome conformation-capture technologies are widely used in 3D genomics; however, experimentally, such methods have high-noise limitations and, therefore, require significant bioinformatics efforts to extract reliable distal interactions. Miscellaneous undesired linear DNAs, present during proximity-ligation, represent a main noise source, which needs to be minimized or eliminated. In this study, different exonuclease combinations were tested to remove linear DNA fragments from a circularized DNA preparation. This method efficiently removed linear DNAs, raised the proportion of annulation and increased the valid-pairs ratio from ∼40% to ∼80% for enhanced interaction detection in standard Hi-C. This strategy is applicable for development of various 3D genomics technologies, or optimization of Hi-C sequencing efficiency.

KLF3 Mediates Epidermal Differentiation through the Epigenomic Writer CBP

Impairments in the differentiation process can lead to skin diseases that can afflict ∼20% of the population. Thus, it is of utmost importance to understand the factors that promote the differentiation process. Here we identify the transcription factor KLF3 as a regulator of epidermal differentiation. Knockdown of KLF3 results in reduced differentiation gene expression and increased cell cycle gene expression. Over half of KLF3's genomic binding sites occur at active enhancers. KLF3 binds to active enhancers proximal to differentiation genes that are dependent upon KLF3 for expression. KLF3's genomic binding sites also highly overlaps with CBP, a histone acetyltransferase necessary for activating enhancers. Depletion of KLF3 causes reduced CBP localization at enhancers proximal to differentiation gene clusters, which leads to loss of enhancer activation but not priming. Our results suggest that KLF3 is necessary to recruit CBP to activate enhancers and drive epidermal differentiation gene expression.

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.

Decreased Enhancer-Promoter Proximity Accompanying Enhancer Activation

Enhancers can regulate the promoters of their target genes over very large genomic distances. It is widely assumed that mechanisms of enhancer action involve the reorganization of three-dimensional chromatin architecture, but this is poorly understood. The predominant model involves physical enhancer-promoter interaction by looping out the intervening chromatin. However, studying the enhancer-driven activation of the Sonic hedgehog gene (Shh), we have identified a change in chromosome conformation that is incompatible with this simple looping model. Using super-resolution 3D-FISH and chromosome conformation capture, we observe a decreased spatial proximity between Shh and its enhancers during the differentiation of embryonic stem cells to neural progenitors. We show that this can be recapitulated by synthetic enhancer activation, is impeded by chromatin-bound proteins located between the enhancer and the promoter, and appears to involve the catalytic activity of poly (ADP-ribose) polymerase. Our data suggest that models of enhancer-promoter communication need to encompass chromatin conformations other than looping.

•

Super-resolution microscopy reveals increased enhancer-promoter separation upon activation

•

Synthetic enhancer activation supports decreased enhancer-promoter proximity

•

Enhancer-promoter separation can be driven by poly(ADP-ribose) polymerase 1

Reprogramming of Meiotic Chromatin Architecture during Spermatogenesis

Chromatin organization undergoes drastic reconfiguration during gametogenesis. However, the molecular reprogramming of three-dimensional chromatin structure in this process remains poorly understood for mammals, including primates. Here, we examined three-dimensional chromatin architecture during spermatogenesis in rhesus monkey using low-input Hi-C. Interestingly, we found that topologically associating domains (TADs) undergo dissolution and reestablishment in spermatogenesis. Strikingly, pachytene spermatocytes, where synapsis occurs, are strongly depleted for TADs despite their active transcription state but uniquely show highly refined local compartments that alternate between transcribing and non-transcribing regions (refined-A/B). Importantly, such chromatin organization is conserved in mouse, where it remains largely intact upon transcription inhibition. Instead, it is attenuated in mutant spermatocytes, where the synaptonemal complex failed to be established. Intriguingly, this is accompanied by the restoration of TADs, suggesting that the synaptonemal complex may restrict TADs and promote local compartments. Thus, these data revealed extensive reprogramming of higher-order meiotic chromatin architecture during mammalian gametogenesis.

Single nucleosome imaging reveals loose genome chromatin networks via active RNA polymerase II

When a gene is activated, chromatin in the transcribed region is thought to be more open and dynamic. However, Nagashima et al. found that this is not necessarily the case—inhibition of transcription globally increases chromatin motion, revealing the existence of loose genome chromatin networks via transcriptional machinery.

Although chromatin organization and dynamics play a critical role in gene transcription, how they interplay remains unclear. To approach this issue, we investigated genome-wide chromatin behavior under various transcriptional conditions in living human cells using single-nucleosome imaging. While transcription by RNA polymerase II (RNAPII) is generally thought to need more open and dynamic chromatin, surprisingly, we found that active RNAPII globally constrains chromatin movements. RNAPII inhibition or its rapid depletion released the chromatin constraints and increased chromatin dynamics. Perturbation experiments of P-TEFb clusters, which are associated with active RNAPII, had similar results. Furthermore, chromatin mobility also increased in resting G0 cells and UV-irradiated cells, which are transcriptionally less active. Our results demonstrated that chromatin is globally stabilized by loose connections through active RNAPII, which is compatible with models of classical transcription factories or liquid droplet formation of transcription-related factors. Together with our computational modeling, we propose the existence of loose chromatin domain networks for various intra-/interchromosomal contacts via active RNAPII clusters/droplets.

Revealing Gene Function and Transcription Relationship by Reconstructing Gene-Level Chromatin Interaction

Chromatin is hierarchically organized in human interphase nuclei. Dynamic chromatin interactions are thought to influence gene transcription and cell fate determination. A consensus concept is that genes may form transcription factories within nucleus by spatially interaction. However, it is still not well known whether the function-related genes co-locate in three-dimensional (3D) space for co-transcription. Especially, there is a lack of visualization method that directly reflect the relationship between gene spatial interaction, gene function and co-transcription. In this study, we constructed three kinds of matrices based on gene ontology annotations, high-through chromosome conformation capture (Hi-C) data and RNA-seq data from twenty human tissues and cell lines. The comparative analysis for gene pairs revealed that 3D genome organization influences gene transcription predominantly at local scale. We found that the local genes within family clusters have similar transcription patterns. We also found that spatial reorganization of a histone gene cluster could control gene transcription. These observations suggest that function-related genes are close in space and activated or repressed together. Our work provided a framework for genome-wide studying the relationship between gene function, co-transcription and spatial interaction.

Deciphering Hi-C: from 3D genome to function

Hi-C is a commonly used technology in 3D genomics which can depict global chromatin interactions across eukaryotic genome. Integrating with different datasets, it can also be applied to studying various biological questions, such as nuclear organization, gene transcription regulation, spatiotemporal development, genome assembly, and cancer genomics. During the last decade, the development and application of Hi-C have dramatically changed the view of genome architecture, chromatin conformation, and gene interaction. So far, Hi-C-related studies remain vivacious and controversial; thus, a unified standard of library construction and bioinformatics analysis are urgently needed. In this review, we have summarized its history, development, methodologies, advances, applications, shortages, and future perspectives. We discuss a few limitations of the current Hi-C technologies and future directions for improvement and highlight how Hi-C can bridge 3D structure to gene function. This review will be helpful for scientists who want to engage in the 3D genomics field; it also shows some future tracks.

Surprising conservation of schizophrenia risk genes in lower organisms reflects their essential function and the evolution of genetic liability

Schizophrenia is a devastating psychiatric illness that affects approximately 1% of the population. Genetic variation in multiple genes causes elevated risk for the disorder, but the molecular basis is inadequately understood and it is not clear how risk genes have evolved and persisted in the genome. To address these issues, we have identified orthologs/homologs of 344 schizophrenia risk genes (from the Psychiatric Genomics Consortium dataset) in lower organisms, including C. elegans, Drosophila and zebrafish, along with phenotypes produced by genetic disruption in C. elegans. Schizophrenia risk genes were evolutionarily conserved at significantly higher rates in C. elegans (81%) and zebrafish (88%) than genes in general for these two species (40-70%). The risk-gene equivalents were highly (~3-fold) enriched for essential genes consistent with polygenic mutation threshold models, which propose that genetic susceptibility results from the inevitable expression of harmful combinations of risk variants in the population. Most notably, numerous examples of cross-species synteny revealed how blocks of risk genes geared toward a shared biological purpose coalesced into proximity during evolution. We obtained initial evidence that schizophrenia risk genes affected different stages of development, potentially allowing differential modulation by the environment. Taken together, studies of the conservation of schizophrenia risk genes in simple model organisms provided novel insights into the molecular basis for genetic susceptibility to a complex human psychiatric disorder.

Chromatin as a nuclear spring

The nucleus in eukaryotic cells is the site for genomic functions such as RNA transcription, DNA replication, and DNA repair/recombination. However, the nucleus is subjected to various mechanical forces associated with diverse cellular activities, including contraction, migration, and adhesion. Although it has long been assumed that the lamina structure, underlying filamentous mesh-work of the nuclear envelope, plays an important role in resisting mechanical forces, the involvement of compact chromatin in mechanical resistance has also recently been suggested. However, it is still unclear how chromatin functions to cope with the stresses. To address this issue, we studied the mechanical responses of human cell nuclei by combining a force measurement microscopy setup with controlled biochemical manipulation of chromatin. We found that nuclei with condensed chromatin possess significant elastic rigidity, whereas the nuclei with a decondensed chromatin are considerably soft. Further analyses revealed that the linker DNA and nucleosome-nucleosome interactions via histone tails in the chromatin act together to generate a spring-like restoring force that resists nuclear deformation. The elastic restoring force is likely to be generated by condensed chromatin domains, consisting of interdigitated or "melted" 10-nm nucleosome fibers. Together with other recent studies, it is suggested that chromatin functions not only as a "memory device" to store, replicate, and express the genetic information for various cellular functions but also as a "nuclear spring" to resist and respond to mechanical forces.

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.

Comparing and Contrasting the Effects of Drosophila Condensin II Subunit dCAP-D3 Overexpression and Depletion in Vivo

The Condensin II complex plays important, conserved roles in genome organization throughout the cell cycle and in the regulation of gene expression. Previous studies have linked decreased Condensin II subunit expression with a variety of diseases. Here, we show that elevated levels of Condensin II subunits are detected in somatic cancers. To evaluate potential biological effects of elevated Condensin II levels, we overexpressed the Condensin II subunit, dCAP-D3 in Drosophila melanogaster larval tissues and examined the effects on the mitotic- and interphase-specific functions of Condensin II. Interestingly, while ubiquitous overexpression resulted in pupal lethality, tissue specific overexpression of dCAP-D3 caused formation of nucleoplasmic protein aggregates which slowed mitotic prophase progression, mimicking results observed when dCAP-D3 levels are depleted. Surprisingly, dCAP-D3 aggregate formation resulted in faster transitions from metaphase to anaphase. Overexpressed dCAP-D3 protein failed to precipitate other Condensin II subunits in nondividing tissues, but did cause changes to gene expression which occurred in a manner opposite of what was observed when dCAP-D3 levels were depleted in both dividing and nondividing tissues. Our findings show that altering dCAP-D3 levels in either direction has detrimental effects on mitotic timing, the regulation of gene expression, and organism development. Taken together, these data suggest that the different roles for Condensin II throughout the cell cycle may be independent of each other and/or that dCAP-D3 may possess functions that are separate from those involving its association with the Condensin II complex. If conserved, these findings could have implications for tumors harboring elevated CAP-D3 levels.

Gene regulation in the 3D genome

The spatial organization of the genome is essential for the precise control of gene expression. Recent advances in sequencing and imaging technologies allow us to explore the 3D genome and its relationship to gene regulation at an unprecedented scale. In this review, we provide an overview of lessons learned from studying the chromatin structure and their implications in communications between gene promoters and distal cis-regulatory elements, such as enhancers. We first review the current knowledge of general genome organization, followed by the importance of chromatin folding in gene regulation. Next, we proceed to a brief survey of the recently developed chromosome conformation capture technologies, as well as most widely adopted read-outs from such data. We then introduce two emerging models that offer explanations regarding how distal enhancers achieve transcriptional control of target genes in the 3D genome. Last, we discuss the promising prospects of leveraging knowledge in chromatin spatial organization for studying complex diseases and traits.

REPA: Applying Pathway Analysis to Genome-Wide Transcription Factor Binding Data

Pathway analysis has been extensively applied to aid in the interpretation of the results of genome-wide transcription profiling studies, and has been shown to successfully find associations between the biological phenomena under study and biological pathways. There are two widely used approaches of pathway analysis: over-representation analysis, and gene set analysis. Recently genome-wide transcription factor binding data has become widely available allowing for the application of pathway analysis to this type of data. In this work, we developed regulatory enrichment pathway analysis (REPA) to apply gene set analysis to genome-wide transcription factor binding data to infer associations between transcription factors and biological pathways. We used the transcription factor binding data generated by the ENCODE project, and gene sets from the Molecular Signatures and KEGG databases. Our results showed that 54 percent of the predictions examined have literature support and that REPA's recall is roughly 54 percent. This level of precision is promising as several of REPA's predictions are expected to be novel and can be used to guide new research avenues. In addition, the results of our case studies showed that REPA enhances the interpretation of genome-wide transcription profiling studies by suggesting putative regulators behind the observed transcriptional responses.

Genetic Dissection of a Supergene Implicates Tfap2a in Craniofacial Evolution of Threespine Sticklebacks

In nature, multiple adaptive phenotypes often coevolve and can be controlled by tightly linked genetic loci known as supergenes. Dissecting the genetic basis of these linked phenotypes is a major challenge in evolutionary genetics. Multiple freshwater populations of threespine stickleback fish (Gasterosteus aculeatus) have convergently evolved two constructive craniofacial traits, longer branchial bones and increased pharyngeal tooth number, likely as adaptations to dietary differences between marine and freshwater environments. Prior QTL mapping showed that both traits are partially controlled by overlapping genomic regions on chromosome 21 and that a regulatory change in Bmp6 likely underlies the tooth number QTL. Here, we mapped the branchial bone length QTL to a 155 kb, eight-gene interval tightly linked to, but excluding the coding regions of Bmp6 and containing the candidate gene Tfap2a Further recombinant mapping revealed this bone length QTL is separable into at least two loci. During embryonic and larval development, Tfap2a was expressed in the branchial bone primordia, where allele specific expression assays revealed the freshwater allele of Tfap2a was expressed at lower levels relative to the marine allele in hybrid fish. Induced loss-of-function mutations in Tfap2a revealed an essential role in stickleback craniofacial development and show that bone length is sensitive to Tfap2a dosage in heterozygotes. Combined, these results suggest that closely linked but genetically separable changes in Bmp6 and Tfap2a contribute to a supergene underlying evolved skeletal gain in multiple freshwater stickleback populations.

Insights about genome function from spatial organization of the genome

Over the last 15 years, development of chromosome conformation capture (3C) and its subsequent high-throughput variants in conjunction with the fast development of sequencing technology has allowed investigators to generate large volumes of data giving insights into the spatial three-dimensional (3D) architecture of the genome. This huge data has been analyzed and validated using various statistical, mathematical, genomics, and biophysical tools in order to examine the chromosomal interaction patterns, understand the organization of the chromosome, and find out functional implications of the interactions. This review summarizes the data generated by several large-scale high-throughput chromosome conformation capture studies and the functional implications obtained from the data analyses. We also discuss emerging results on factors (both CCCTC binding factor (CTCF) related and CTCF independent) that could contribute to looping interactions.

The biology and polymer physics underlying large-scale chromosome organization

Chromosome large-scale organization is a beautiful example of the interplay between physics and biology. DNA molecules are polymers and thus belong to the class of molecules for which physicists have developed models and formulated testable hypotheses to understand their arrangement and dynamic properties in solution, based on the principles of polymer physics. Biologists documented and discovered the biochemical basis for the structure, function and dynamic spatial organization of chromosomes in cells. The underlying principles of chromosome organization have recently been revealed in unprecedented detail using high-resolution chromosome capture technology that can simultaneously detect chromosome contact sites throughout the genome. These independent lines of investigation have now converged on a model in which DNA loops, generated by the loop extrusion mechanism, are the basic organizational and functional units of the chromosome.

Graphene oxide regulates cox2 in human embryonic kidney 293T cells via epigenetic mechanisms: dynamic chromosomal interactions

To extend the applications of engineered nanomaterials, such as graphene oxide (GO), it is necessary to minimize cytotoxicity. However, the mechanisms underlying this cytotoxicity are unclear. Dynamic chromosomal interactions have been used to illustrate the molecular bases of gene expression, which offers a more sensitive and cutting-edge technology to elucidate complex biological processes associated with epigenetic regulations. In this study, the role of GO-triggered chromatin interactions in the activation of cox2, a hallmark of inflammation, was investigated in normal human cells. Using chromosome conformation capture technology, we showed that GO triggers physical interactions between the downstream enhancer and the cox2 promoter in human embryonic kidney 293T (293T) via p65 and p300 complex-mediated dynamic chromatin looping, which was required for high cox2 expression. Moreover, tumor necrosis factor-α (TNF-α), located upstream of the p65 signaling pathway, contributed to the regulation of cox2 activation through dynamic chromatin architecture. Compared with pristine GO and aminated GO (GO-NH₂), poly (acrylic acid)-functionalized GO (GO-PAA) induced a weaker inflammatory response and a weaker effect on chromatin architecture. Our results mechanistically link GO-mediated chromatin interactions with the regulation of cox2 and suggest that GO derivatives may minimize toxicity in practical applications.

Genes uniquely expressed in human growth plate chondrocytes uncover a distinct regulatory network

Background

Chondrogenesis is the earliest stage of skeletal development and is a highly dynamic process, integrating the activities and functions of transcription factors, cell signaling molecules and extracellular matrix proteins. The molecular mechanisms underlying chondrogenesis have been extensively studied and multiple key regulators of this process have been identified. However, a genome-wide overview of the gene regulatory network in chondrogenesis has not been achieved.

Results

In this study, employing RNA sequencing, we identified 332 protein coding genes and 34 long non-coding RNA (lncRNA) genes that are highly selectively expressed in human fetal growth plate chondrocytes. Among the protein coding genes, 32 genes were associated with 62 distinct human skeletal disorders and 153 genes were associated with skeletal defects in knockout mice, confirming their essential roles in skeletal formation. These gene products formed a comprehensive physical interaction network and participated in multiple cellular processes regulating skeletal development. The data also revealed 34 transcription factors and 11,334 distal enhancers that were uniquely active in chondrocytes, functioning as transcriptional regulators for the cartilage-selective genes.

Conclusions

Our findings revealed a complex gene regulatory network controlling skeletal development whereby transcription factors, enhancers and lncRNAs participate in chondrogenesis by transcriptional regulation of key genes. Additionally, the cartilage-selective genes represent candidate genes for unsolved human skeletal disorders.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4378-y) contains supplementary material, which is available to authorized users.

Cis-regulatory landscapes of four cell types of the retina

The retina is composed of ∼50 cell-types with specific functions for the process of vision. Identification of the cis-regulatory elements active in retinal cell-types is key to elucidate the networks controlling this diversity. Here, we combined transcriptome and epigenome profiling to map the regulatory landscape of four cell-types isolated from mouse retinas including rod and cone photoreceptors as well as rare inter-neuron populations such as horizontal and starburst amacrine cells. Integration of this information reveals sequence determinants and candidate transcription factors for controlling cellular specialization. Additionally, we refined parallel reporter assays to enable studying the transcriptional activity of large collection of sequences in individual cell-types isolated from a tissue. We provide proof of concept for this approach and its scalability by characterizing the transcriptional capacity of several hundred putative regulatory sequences within individual retinal cell-types. This generates a catalogue of cis-regulatory regions active in retinal cell types and we further demonstrate their utility as potential resource for cellular tagging and manipulation.