Chromatin and epigenetic features of long-range gene regulation

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.

Learning the histone codes with large genomic windows and three-dimensional chromatin interactions using transformer

The quantitative characterization of the transcriptional control by histone modifications has been challenged by many computational studies, but most of them only focus on narrow and linear genomic regions around promoters, leaving a room for improvement. We present Chromoformer, a transformer-based, three-dimensional chromatin conformation-aware deep learning architecture that achieves the state-of-the-art performance in the quantitative deciphering of the histone codes in gene regulation. The core essence of Chromoformer architecture lies in the three variants of attention operation, each specialized to model individual hierarchy of transcriptional regulation involving from core promoters to distal elements in contact with promoters through three-dimensional chromatin interactions. In-depth interpretation of Chromoformer reveals that it adaptively utilizes the long-range dependencies between histone modifications associated with transcription initiation and elongation. We also show that the quantitative kinetics of transcription factories and Polycomb group bodies can be captured by Chromoformer. Together, our study highlights the great advantage of attention-based deep modeling of complex interactions in epigenomes.

The use of machine learning to discover regulatory networks controlling biological systems

Biological systems are composed of a vast web of multiscale molecular interactors and interactions. High-throughput technologies, both bulk and single cell, now allow for investigation of the properties and quantities of these interactors. Computational algorithms and machine learning methods then provide the tools to derive meaningful insights from the resulting data sets. One such approach is graphical network modeling, which provides a computational framework to explicitly model the molecular interactions within and between the cells comprising biological systems. These graphical networks aim to describe a putative chain of cause and effect between interacting molecules. This feature allows for determination of key molecules in a biological process, accelerated generation of mechanistic hypotheses, and simulation of experimental outcomes. We review the computational concepts and applications of graphical network models across molecular scales for both intracellular and intercellular regulatory biology, examples of successful applications, and the future directions needed to overcome current limitations.

MZF1 regulates α-globin gene transcription via long-range interactions in erythroid differentiation

Precise spatiotemporal gene expression regulation is crucial for human erythropoiesis. However, dramatic changes in the chromatin structure and transcriptome involved in α-globin gene expression during erythropoiesis still not fully understand. To identify candidate regulators for α-globin gene regulation, we carried out an integrated approach by integrating publicly available transcriptomic and epigenomic data. We computed active enhancers by overlapping enriched regions marked with H3K4me1 and H3K27ac and correlated their activity with mRNA expression. Next, we cataloged potential transcription factors via de novo motif analysis. We highlighted the discovery of potential novel transcription factor MZF1 of the α-globin gene in erythroid differentiation. To validate the role of MZF1, we quantified the expression level of MZF1 and α-globin gene in HSPCs, early erythroid progenitors and late erythroid precursors cells. Both the mRNA and protein expression patterns of MZF1 were consistent with the α-globin gene. Also, the qPCR result showed that the expression of the α-globin gene was significantly increased by the MZF1 overexpression. To further investigate the role of MZF1 regulating α-globin gene transcriptional activity during erythroid differentiation, we performed ChIP-qPCR at the α-globin locus. Our results showed that MZF1 recruitment both at 4 upstream HS sites and α-globin gene promoter in erythroid precursor cells. To determine the importance of the MZF1 to enhancer-promoter interaction at the α-globin locus, we compared interaction frequency before and after knockdown of MZF1 by chromosome conformation capture (3C) assay. Upon MZF1 depletion, both the expression of the α-globin gene and all 3C signals were significantly decreased. Taken together, MZF1 plays an important role in regulating α-globin gene expression by binding to long-region enhancers and α-globin gene promoter and facilitates the organization of specific 3D chromatin architecture in erythroid differentiation.

Profound alterations of the chromatin architecture at chromosome 11p15.5 in cells from Beckwith-Wiedemann and Silver-Russell syndromes patients

Beckwith-Wiedemann syndrome (BWS) and Silver-Russell syndrome (SRS) are imprinting-related disorders associated with genetic/epigenetic alterations of the 11p15.5 region, which harbours two clusters of imprinted genes (IGs). 11p15.5 IGs are regulated by the methylation status of imprinting control regions ICR1 and ICR2. 3D chromatin structure is thought to play a pivotal role in gene expression control; however, chromatin architecture models are still poorly defined in most cases, particularly for IGs. Our study aimed at elucidating 11p15.5 3D structure, via 3C and 3D FISH analyses of cell lines derived from healthy, BWS or SRS children. We found that, in healthy cells, IGF2/H19 and CDKN1C/KCNQ1OT1 domains fold in complex chromatin conformations, that facilitate the control of IGs mediated by distant enhancers. In patient-derived cell lines, we observed a profound impairment of such a chromatin architecture. Specifically, we identified a cross-talk between IGF2/H19 and CDKN1C/KCNQ1OT1 domains, consisting in in cis, monoallelic interactions, that are present in healthy cells but lost in patient cell lines: an inter-domain association that sees ICR2 move close to IGF2 on one allele, and to H19 on the other. Moreover, an intra-domain association within the CDKN1C/KCNQ1OT1 locus seems to be crucial for maintaining the 3D organization of the region.

ENdb: a manually curated database of experimentally supported enhancers for human and mouse

Enhancers are a class of cis-regulatory elements that can increase gene transcription by forming loops in intergenic regions, introns and exons. Enhancers, as well as their associated target genes, and transcription factors (TFs) that bind to them, are highly associated with human disease and biological processes. Although some enhancer databases have been published, most only focus on enhancers identified by high-throughput experimental techniques. Therefore, it is highly desirable to construct a comprehensive resource of manually curated enhancers and their related information based on low-throughput experimental evidences. Here, we established a comprehensive manually-curated enhancer database for human and mouse, which provides a resource for experimentally supported enhancers, and to annotate the detailed information of enhancers. The current release of ENdb documents 737 experimentally validated enhancers and their related information, including 384 target genes, 263 TFs, 110 diseases and 153 functions in human and mouse. Moreover, the enhancer-related information was supported by experimental evidences, such as RNAi, in vitro knockdown, western blotting, qRT-PCR, luciferase reporter assay, chromatin conformation capture (3C) and chromosome conformation capture-on-chip (4C) assays. ENdb provides a user-friendly interface to query, browse and visualize the detailed information of enhancers. The database is available at http://www.licpathway.net/ENdb.

[Opposite Effects of Histone H1 and HMGN5 Protein on Distant Interactions in Chromatin]

Transcriptional enhancers in the cell nuclei typically interact with the target promoters in cis over long stretches of chromatin, but the mechanism of this communication remains unknown. Previously we have developed a defined in vitro system for quantitative analysis of the rate of distant enhancer-promoter communication (EPC) and have shown that the chromatin fibers maintain efficient distant EPC in cis. Here we investigate the roles of linker histone H1 and HMGN5 protein in EPC. A considerable negative effect of histone H1 on EPC depending on its C- and N-tails was shown. Protein HMGN5 that affects chromatin compaction and is associated with active chromatin counteracts EPC inhibition by H1. The data suggest that the efficiency of the interaction between the enhancer and the promoter depends on the structure and dynamics of the chromatin fiber localized between them and can be regulated by proteins associated with chromatin.

Integrating multiomics longitudinal data to reconstruct networks underlying lung development

A comprehensive understanding of the dynamic regulatory networks that govern postnatal alveolar lung development is still lacking. To construct such a model, we profiled mRNA, microRNA, DNA methylation, and proteomics of developing murine alveoli isolated by laser capture microdissection at 14 predetermined time points. We developed a detailed comprehensive and interactive model that provides information about the major expression trajectories, the regulators of specific key events, and the impact of epigenetic changes. Intersecting the model with single-cell RNA-Seq data led to the identification of active pathways in multiple or individual cell types. We then constructed a similar model for human lung development by profiling time-series human omics data sets. Several key pathways and regulators are shared between the reconstructed models. We experimentally validated the activity of a number of predicted regulators, leading to new insights about the regulation of innate immunity during lung development.

Dysregulated Transcriptional Control in Prostate Cancer

Recent advances in whole-genome and transcriptome sequencing of prostate cancer at different stages indicate that a large number of mutations found in tumors are present in non-protein coding regions of the genome and lead to dysregulated gene expression. Single nucleotide variations and small mutations affecting the recruitment of transcription factor complexes to DNA regulatory elements are observed in an increasing number of cases. Genomic rearrangements may position coding regions under the novel control of regulatory elements, as exemplified by the TMPRSS2-ERG fusion and the amplified enhancer identified upstream of the androgen receptor (AR) gene. Super-enhancers are increasingly found to play important roles in aberrant oncogenic transcription. Several players involved in these processes are currently being evaluated as drug targets and may represent new vulnerabilities that can be exploited for prostate cancer treatment. They include factors involved in enhancer and super-enhancer function such as bromodomain proteins and cyclin-dependent kinases. In addition, non-coding RNAs with an important gene regulatory role are being explored. The rapid progress made in understanding the influence of the non-coding part of the genome and of transcription dysregulation in prostate cancer could pave the way for the identification of novel treatment paradigms for the benefit of patients.

Building Transcription Factor Binding Site Models to Understand Gene Regulation in Plants

Transcription factors (TFs) are key cellular components that control gene expression. They recognize specific DNA sequences, the TF binding sites (TFBSs), and thus are targeted to specific regions of the genome where they can recruit transcriptional co-factors and/or chromatin regulators to fine-tune spatiotemporal gene regulation. Therefore, the identification of TFBSs in genomic sequences and their subsequent quantitative modeling is of crucial importance for understanding and predicting gene expression. Here, we review how TFBSs can be determined experimentally, how the TFBS models can be constructed in silico, and how they can be optimized by taking into account features such as position interdependence within TFBSs, DNA shape, and/or by introducing state-of-the-art computational algorithms such as deep learning methods. In addition, we discuss the integration of context variables into the TFBS modeling, including nucleosome positioning, chromatin states, methylation patterns, 3D genome architectures, and TF cooperative binding, in order to better predict TF binding under cellular contexts. Finally, we explore the possibilities of combining the optimized TFBS model with technological advances, such as targeted TFBS perturbation by CRISPR, to better understand gene regulation, evolution, and plant diversity.

An upstream enhancer regulates Gpihbp1 expression in a tissue-specific manner

Glycosylphosphatidylinositol-anchored high density lipoprotein-binding protein 1 (GPIHBP1), the protein that shuttles LPL to the capillary lumen, is essential for plasma triglyceride metabolism. When GPIHBP1 is absent, LPL remains stranded within the interstitial spaces and plasma triglyceride hydrolysis is impaired, resulting in severe hypertriglyceridemia. While the functions of GPIHBP1 in intravascular lipolysis are reasonably well understood, no one has yet identified DNA sequences regulating GPIHBP1 expression. In the current studies, we identified an enhancer element located ∼3.6 kb upstream from exon 1 of mouse Gpihbp1. To examine the importance of the enhancer, we used CRISPR/Cas9 genome editing to create mice lacking the enhancer (Gpihbp1Enh/Enh). Removing the enhancer reduced Gpihbp1 expression by >90% in the liver and by ∼50% in heart and brown adipose tissue. The reduced expression of GPIHBP1 was insufficient to prevent LPL from reaching the capillary lumen, and it did not lead to hypertriglyceridemia-even when mice were fed a high-fat diet. Compound heterozygotes (Gpihbp1Enh/- mice) displayed further reductions in Gpihbp1 expression and exhibited partial mislocalization of LPL (increased amounts of LPL within the interstitial spaces of the heart), but the plasma triglyceride levels were not perturbed. The enhancer element that we identified represents the first insight into DNA sequences controlling Gpihbp1 expression.

Identification and characterisation of hypomethylated DNA loci controlling quantitative resistance in Arabidopsis

Variation in DNA methylation enables plants to inherit traits independently of changes to DNA sequence. Here, we have screened an Arabidopsis population of epigenetic recombinant inbred lines (epiRILs) for resistance against Hyaloperonospora arabidopsidis (Hpa). These lines share the same genetic background, but show variation in heritable patterns of DNA methylation. We identified four epigenetic quantitative trait loci (epiQTLs) that provide quantitative resistance without reducing plant growth or resistance to other (a)biotic stresses. Phenotypic characterisation and RNA-sequencing analysis revealed that Hpa-resistant epiRILs are primed to activate defence responses at the relatively early stages of infection. Collectively, our results show that hypomethylation at selected pericentromeric regions is sufficient to provide quantitative disease resistance, which is associated with genome-wide priming of defence-related genes. Based on comparisons of global gene expression and DNA methylation between the wild-type and resistant epiRILs, we discuss mechanisms by which the pericentromeric epiQTLs could regulate the defence-related transcriptome.

In plants, animals and microbes genetic information is encoded by DNA, which are made up of sequences of building blocks, called nucleotide bases. These sequences can be separated into sections known as genes that each encode specific traits. It was previously thought that only changes to the sequence of bases in a DNA molecule could alter the traits passed on to future generations. However, it has recently become clear that some traits can also be inherited through modifications to the DNA that do not alter its sequence.

One such modification is to attach a tag, known as a methyl group, to a nucleotide base known as cytosine. These methyl tags can be added to, or removed from, DNA to create different patterns of methylation. Previous studies have shown that plants whose DNA is less methylated than normal (‘hypo-methylated’) are more resistant to plant diseases. However, the location and identity of the hypo-methylated DNA regions controlling this resistance remained unknown. To address this problem, Furci, Jain et al. studied how DNA methylation in a small weed known as Arabidopsis thaliana affects how well the plants can resist a disease known as downy mildew.

Furci, Jain et al. studied a population of over 100 A. thaliana lines that have the same DNA sequences but different patterns of DNA methylation. The experiments identified four DNA locations that were less methylated in lines with enhanced resistance to downy mildew. Importantly, this form of resistance did not appear to reduce how well the plants grew, or make them less able to resist other diseases or environmental stresses.

The results of further experiments suggested that reduced methylation at the four DNA regions prime the plant’s immune system, enabling a faster and stronger activation of a multitude of defence genes across the genome after attack by downy mildew.

The next steps following on from this work are to investigate exactly how the four DNA regions with reduced methylation can prime so many different defence genes in the plant. Further research is also needed to determine whether it is possible to breed crop plants with lower levels of methylation at specific DNA locations to improve disease resistance, but without decreasing the amount and quality of food produced.

Energy required to pinch a DNA plectoneme

DNA supercoiling plays an important role from a biological point of view. One of its consequences at the supramolecular level is the formation of DNA superhelices named plectonemes. Normally separated by a distance on the order of 10 nm, the two opposite double strands of a DNA plectoneme must be brought closer if a protein or protein complex implicated in genetic regulation is to be bound simultaneously to both strands, as if the plectoneme was locally pinched. We propose an analytic calculation of the energetic barrier, of elastic nature, required to bring closer the two loci situated on the opposed double strands. We examine how this energy barrier scales with the DNA supercoiling. For physically relevant values of elastic parameters and of supercoiling density, we show that the energy barrier is in the k_{B}T range under physiological conditions, thus demonstrating that the limiting step to loci encounter is more likely the preceding plectoneme slithering bringing the two loci side by side.

Genetics in ankylosing spondylitis - Current state of the art and translation into clinical outcomes

Ankylosing spondylitis (AS) is the prototypic form of axial spondyloarthritis (axSpA). It is highly heritable, with studies conducted in twins and in unrelated cases and controls showing that the heritability for AS is much higher than those for inflammatory bowel disease or rheumatoid arthritis. To date, 116 loci have been identified, contributing to 28% of the genetic variation in the disease. These loci provide important clues into pathogenic pathways in the disease that have led to therapeutic advances such as the repositioning of IL-17 inhibitors in the disease. Much more research is currently required to determine the functional mechanisms by which the genetic associations operate, from which it is likely that novel therapeutic approaches will be developed.

The relationship between transgenerational acquired resistance and global DNA methylation in Arabidopsis

Progeny of heavily diseased plants develop transgenerational acquired resistance (TAR). In Arabidopsis, TAR can be transmitted over one stress-free generation. Although DNA methylation has been implicated in the regulation of TAR, the relationship between TAR and global DNA methylation remains unknown. Here, we characterised the methylome of TAR-expressing Arabidopsis at different generations after disease exposure. Global clustering of cytosine methylation revealed TAR-related patterns in the F3 generation, but not in the F1 generation. The majority of differentially methylated positions (DMPs) occurred at CG context in gene bodies. TAR in F3 progeny after one initial generation of disease, followed by two stress-free generations, was lower than TAR in F3 progeny after three successive generations of disease. This difference in TAR effectiveness was proportional to the intensity of differential methylation at a sub-set of cytosine positions. Comparison of TAR-related DMPs with previously characterised cytosine methylation in mutation accumulation lines revealed that ancestral disease stress preferentially acts on methylation-labile cytosine positions, but also extends to methylation-stable positions. Thus, the TAR-related impact of ancestral disease extends beyond stochastic variation in DNA methylation. Our study has shown that the Arabidopsis epigenome responds globally to disease in previous generations and we discuss its contribution to TAR.

Epigenetic and Cellular Diversity in the Brain through Allele-Specific Effects

The benefits of diploidy are considered to involve masking partially recessive mutations and increasing genetic diversity. Here, we review new studies showing evidence for diverse allele-specific expression and epigenetic states in mammalian brain cells, which suggest that diploidy expands the landscape of gene regulatory and expression programs in cells. Allele-specific expression has been thought to be restricted to a few specific classes of genes. However, new studies show novel genomic imprinting effects that are brain-region-, cell-type- and age-dependent. In addition, novel forms of random monoallelic expression that impact many autosomal genes have been described in vitro and in vivo. We discuss the implications for understanding the benefits of diploidy, and the mechanisms shaping brain development, function, and disease.

Novel long-range regulatory mechanisms controlling PKD2 gene expression

Background

Cis-regulatory elements control gene expression over large distances through the formation of chromatin loops, which allow contact between enhancers and gene promoters. Alterations in cis-acting regulatory systems could be linked to human genetic diseases. Here, we analyse the spatial organization of a large region spanning the polycystic kidney disease 2 (PKD2) gene, one of the genes responsible of autosomal dominant polycystic kidney disease (ADPKD).

Results

By using chromosome conformation capture carbon copy (5C) technology in primary human renal cyst epithelial cells, we identify novel contacts of the PKD2 promoter with chromatin regions, which display characteristics of regulatory elements. In parallel, by using functional analysis with a reporter assay, we demonstrate that three DNAse I hypersensitive sites regions are involved in the regulation of PKD2 gene expression.

Conclusions

Finally, through alignment of CCCTC-binding factor (CTCF) sites, we suggest that these novel enhancer elements are brought to the PKD2 promoter by chromatin looping via the recruitment of CTCF.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4892-6) contains supplementary material, which is available to authorized users.

Relationship of the Interaction Between Two Quantitative Trait Loci with γ-Globin Expression in β-Thalassemia Intermedia Patients

Globin switching is a significant factor on blood hemoglobin (Hb) level but its molecular mechanisms have not yet been identified, however, several quantitative trait loci (QTL) and polymorphisms involved regions on chromosomes 2p, 6q, 8q and X account for variation in the γ-globin expression level. We studied the effect of interaction between a region on intron six of the TOX gene, chromosome 8q (chr8q) and XmnI locus on the γ-globin promoter, chr11p on γ-globin expression in 150 β-thalassemia intermedia (β-TI) patients, evaluated by statistical interaction analysis. Our results showed a significant interaction between one QTL on intron six of the TOX gene (rs9693712) and XmnI locus that effect γ-globin expression. Interchromosomal interaction mediates through transcriptional machanisms to preserve true genome architectural features, chromosomes localization and DNA bending. This interaction can be a part of the unknown molecular mechanism of globin switching and regulation of gene expression.

Integrating ChIP-seq with other functional genomics data

Transcription is regulated by transcription factor (TF) binding at promoters and distal regulatory elements and histone modifications that control the accessibility of these elements. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) has become the standard assay for identifying genome-wide protein-DNA interactions in vitro and in vivo. As large-scale ChIP-seq data sets have been collected for different TFs and histone modifications, their potential to predict gene expression can be used to test hypotheses about the mechanisms of gene regulation. In addition, complementary functional genomics assays provide a global view of chromatin accessibility and long-range cis-regulatory interactions that are being combined with TF binding and histone remodeling to study the regulation of gene expression. Thus, ChIP-seq analysis is now widely integrated with other functional genomics assays to better understand gene regulatory mechanisms. In this review, we discuss advances and challenges in integrating ChIP-seq data to identify context-specific chromatin states associated with gene activity. We describe the overall computational design of integrating ChIP-seq data with other functional genomics assays. We also discuss the challenges of extending these methods to low-input ChIP-seq assays and related single-cell assays.

From Genotype to Phenotype: A Primer on the Functional Follow-up of Genome-Wide Association Studies in Cardiovascular Disease

Genome-wide association studies (GWASs) have implicated many human genomic loci in the development of complex traits. The loci identified by these studies are potentially involved in novel pathways that contribute to disease pathophysiology. However, eventual therapeutic targeting of these pathways relies on bridging the gap between genetic association and function, a task that first requires validation of causal genetic variants, casual genes, and directionality of effect. Executing this task requires basic knowledge of interpreting GWAS results and prioritizing candidates for further study, in addition to understanding the experimental methods available for evaluating candidate variants. Here we review the basic genetic principles of genome-wide association studies, the computational and experimental tools used for identifying causal variants and genes, and salient illustrative examples of how cardiovascular loci have undergone functional investigation.

Signaling by Steroid Hormones in the 3D Nuclear Space

Initial studies showed that ligand-activated hormone receptors act by binding to the proximal promoters of individual target genes. Genome-wide studies have now revealed that regulation of transcription by steroid hormones mainly depends on binding of the receptors to distal regulatory elements. Those distal elements, either enhancers or silencers, act on the regulation of target genes by chromatin looping to the gene promoters. In the nucleus, this level of chromatin folding is integrated within dynamic higher orders of genome structures, which are organized in a non-random fashion. Terminally differentiated cells exhibit a tissue-specific three-dimensional (3D) organization of the genome that favors or restrains the activity of transcription factors and modulates the function of steroid hormone receptors, which are transiently activated upon hormone exposure. Conversely, integration of the hormones signal may require modifications of the 3D organization to allow appropriate transcriptional outcomes. In this review, we summarize the main levels of organization of the genome, review how they can modulate the response to steroids in a cell specific manner and discuss the role of receptors in shaping and rewiring the structure in response to hormone. Taking into account the dynamics of 3D genome organization will contribute to a better understanding of the pleiotropic effects of steroid hormones in normal and cancer cells.

The emerging landscape of in vitro and in vivo epigenetic allelic effects

Epigenetic mechanisms that cause maternally and paternally inherited alleles to be expressed differently in offspring have the potential to radically change our understanding of the mechanisms that shape disease susceptibility, phenotypic variation, cell fate, and gene expression. However, the nature and prevalence of these effects in vivo have been unclear and are debated. Here, I consider major new studies of epigenetic allelic effects in cell lines and primary cells and in vivo. The emerging picture is that these effects take on diverse forms, and this review attempts to clarify the nature of the different forms that have been uncovered for genomic imprinting and random monoallelic expression (RME). I also discuss apparent discrepancies between in vitro and in vivo studies. Importantly, multiple studies suggest that allelic effects are prevalent and can be developmental stage- and cell type-specific. I propose some possible functions and consider roles for allelic effects within the broader context of gene regulatory networks, cellular diversity, and plasticity. Overall, the field is ripe for discovery and is in need of mechanistic and functional studies.

Genome organization: connecting the developmental origins of disease and genetic variation

An adverse early life environment can increase the risk of metabolic and other disorders later in life. Genetic variation can modify an individual's susceptibility to these environmental challenges. These gene by environment interactions are important, but difficult, to dissect. The nucleus is the primary organelle where environmental responses impact directly on the genetic variants within the genome, resulting in changes to the biology of the genome and ultimately the phenotype. Understanding genome biology requires the integration of the linear DNA sequence, epigenetic modifications and nuclear proteins that are present within the nucleus. The interactions between these layers of information may be captured in the emergent spatial genome organization. As such genome organization represents a key research area for decoding the role of genetic variation in the Developmental Origins of Health and Disease.

Nucleosome-free DNA regions differentially affect distant communication in chromatin

Communication between distantly spaced genomic regions is one of the key features of gene regulation in eukaryotes. Chromatin per se can stimulate efficient enhancer-promoter communication (EPC); however, the role of chromatin structure and dynamics in this process remains poorly understood. Here we show that nucleosome spacing and the presence of nucleosome-free DNA regions can modulate chromatin structure/dynamics and, in turn, affect the rate of EPC in vitro and in silico. Increasing the length of internucleosomal linker DNA from 25 to 60 bp results in more efficient EPC. The presence of longer nucleosome-free DNA regions can positively or negatively affect the rate of EPC, depending upon the length and location of the DNA region within the chromatin fiber. Thus the presence of histone-free DNA regions can differentially affect the efficiency of EPC, suggesting that gene regulation over a distance could be modulated by changes in the length of internucleosomal DNA spacers.

Long-range control of gene expression via RNA-directed DNA methylation

RNA-mediated transcriptional silencing, in plants known as RNA-directed DNA methylation (RdDM), is a conserved process where small interfering RNA (siRNA) and long non-coding RNA (lncRNA) help establish repressive chromatin modifications. This process represses transposons and affects the expression of protein-coding genes. We found that in Arabidopsis thaliana AGO4 binding sites are often located distant from genes differentially expressed in ago4. Using Hi-C to compare interactions between genotypes, we show that RdDM-targeted loci have the potential to engage in chromosomal interactions, but these interactions are inhibited in wild-type conditions. In mutants defective in RdDM, the frequency of chromosomal interactions at RdDM targets is increased. This includes increased frequency of interactions between Pol V methylated sites and distal genes that are repressed by RdDM. We propose a model, where RdDM prevents the formation of chromosomal interactions between genes and their distant regulatory elements.

Prenatal Exposure to a Maternal High-Fat Diet Affects Histone Modification of Cardiometabolic Genes in Newborn Rats

Infants born to women with diabetes or obesity are exposed to excess circulating fuels during fetal heart development and are at higher risk of cardiac diseases. We have previously shown that late-gestation diabetes, especially in conjunction with a maternal high-fat (HF) diet, impairs cardiac functions in rat-offspring. This study investigated changes in genome-wide histone modifications in newborn hearts from rat-pups exposed to maternal diabetes and HF-diet. Chromatin-immunoprecipitation-sequencing revealed a differential peak distribution on gene promoters in exposed pups with respect to acetylation of lysines 9 and 14 and to trimethylation of lysines 4 and 27 in histone H3 (all, false discovery rate, FDR < 0.1). In the HF-diet exposed offspring, 54% of the annotated genes showed the gene-activating mark trimethylated lysine 4. Many of these genes (1) are associated with the "metabolic process" in general and particularly with "positive regulation of cholesterol biosynthesis" (FDR = 0.03); (2) overlap with 455 quantitative trait loci for blood pressure, body weight, serum cholesterol (all, FDR < 0.1); and (3) are linked to cardiac disease susceptibility/progression, based on disease ontology analyses and scientific literature. These results indicate that maternal HF-diet changes the cardiac histone signature in offspring suggesting a fuel-mediated epigenetic reprogramming of cardiac tissue in utero.

Deep homology in the age of next-generation sequencing

The principle of homology is central to conceptualizing the comparative aspects of morphological evolution. The distinctions between homologous or non-homologous structures have become blurred, however, as modern evolutionary developmental biology (evo-devo) has shown that novel features often result from modification of pre-existing developmental modules, rather than arising completely de novo. With this realization in mind, the term 'deep homology' was coined, in recognition of the remarkably conserved gene expression during the development of certain animal structures that would not be considered homologous by previous strict definitions. At its core, it can help to formulate an understanding of deeper layers of ontogenetic conservation for anatomical features that lack any clear phylogenetic continuity. Here, we review deep homology and related concepts in the context of a gene expression-based homology discussion. We then focus on how these conceptual frameworks have profited from the recent rise of high-throughput next-generation sequencing. These techniques have greatly expanded the range of organisms amenable to such studies. Moreover, they helped to elevate the traditional gene-by-gene comparison to a transcriptome-wide level. We will end with an outlook on the next challenges in the field and how technological advances might provide exciting new strategies to tackle these questions.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Sex-specific chromatin landscapes in an ultra-compact chordate genome

Background

In multicellular organisms, epigenome dynamics are associated with transitions in the cell cycle, development, germline specification, gametogenesis and inheritance. Evolutionarily, regulatory space has increased in complex metazoans to accommodate these functions. In tunicates, the sister lineage to vertebrates, we examine epigenome adaptations to strong secondary genome compaction, sex chromosome evolution and cell cycle modes.

Results

Across the 70 MB Oikopleura dioica genome, we profiled 19 histone modifications, and RNA polymerase II, CTCF and p300 occupancies, to define chromatin states within two homogeneous tissues with distinct cell cycle modes: ovarian endocycling nurse nuclei and mitotically proliferating germ nuclei in testes. Nurse nuclei had active chromatin states similar to other metazoan epigenomes, with large domains of operon-associated transcription, a general lack of heterochromatin, and a possible role of Polycomb PRC2 in dosage compensation. Testis chromatin states reflected transcriptional activity linked to spermatogenesis and epigenetic marks that have been associated with establishment of transgenerational inheritance in other organisms. We also uncovered an unusual chromatin state specific to the Y-chromosome, which combined active and heterochromatic histone modifications on specific transposable elements classes, perhaps involved in regulating their activity.

Conclusions

Compacted regulatory space in this tunicate genome is accompanied by reduced heterochromatin and chromatin state domain widths. Enhancers, promoters and protein-coding genes have conserved epigenomic features, with adaptations to the organization of a proportion of genes in operon units. We further identified features specific to sex chromosomes, cell cycle modes, germline identity and dosage compensation, and unusual combinations of histone PTMs with opposing consensus functions.

Electronic supplementary material

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

Inferring a role for methylation of intergenic DNA in the regulation of genes aberrantly expressed in precursor B-cell acute lymphoblastic leukemia

A complete understanding of the mechanisms involved in the development of pre-B ALL is lacking. In this study, we integrated DNA methylation data and gene expression data to elucidate the impact of aberrant intergenic DNA methylation on gene expression in pre-B ALL. We found a subset of differentially methylated intergenic loci that were associated with altered gene expression in pre-B ALL patients. Notably, 84% of these regions were also bound by transcription factors (TF) known to play roles in differentiation and B-cell development in a lymphoblastoid cell line. Further, an overall downregulation of eRNA transcripts was observed in pre-B ALL patients and these transcripts were associated with the downregulation of putative target genes involved in B-cell migration, proliferation, and apoptosis. The identification of novel putative regulatory regions highlights the significance of intergenic DNA sequences and may contribute to the identification of new therapeutic targets for the treatment of pre-B ALL.

Extraoral bitter taste receptors in health and disease

Bitter taste receptors (TAS2Rs or T2Rs) belong to the superfamily of seven-transmembrane G protein-coupled receptors, which are the targets of >50% of drugs currently on the market. Canonically, T2Rs are located in taste buds of the tongue, where they initiate bitter taste perception. However, accumulating evidence indicates that T2Rs are widely expressed throughout the body and mediate diverse nontasting roles through various specialized mechanisms. It has also become apparent that T2Rs and their polymorphisms are associated with human disorders. In this review, we summarize the physiological and pathophysiological roles that extraoral T2Rs play in processes as diverse as innate immunity and reproduction, and the major challenges in this emerging field.

Genetics of Thyroid-Stimulating Hormone Receptor-Relevance for Autoimmune Thyroid Disease

Production of thyroid-stimulating hormone receptor (TSHR) antibodies represents the hallmark of Graves' disease (GD) pathogenesis. Thus, for more than two decades the TSHR gene has been at the center of studies intended to elucidate its contribution to disease pathology. The advent of genome-wide association technology allowed to establish a strong association of the TSHR gene with GD. Subsequent fine-mapping studies narrowed the disease-susceptibility region to a 40 kb sequence in intron 1, where at least five GD-associated SNPs in tight linkage disequilibrium were identified. The current challenge is to understand the functional mechanisms by which these polymorphisms modify physiological processes and trigger disease. The aim of this review is to summarize the current knowledge on the role of the TSHR gene in GD pathogenesis, which has been gained through linkage and association studies, as well as to discuss the emerging mechanisms underlying biological implications of TSHR variants in the development of GD.

The Promoter and Multiple Enhancers of the pou4f3 Gene Regulate Expression in Inner Ear Hair Cells

Few enhancers that target gene expression to inner ear hair cells (HCs) have been identified. Using transgenic analysis of enhanced green fluorescent protein (eGFP) reporter constructs and bioinformatics, we evaluated the control of pou4f3 gene expression, since it is expressed only in HCs within the inner ear and continues to be expressed throughout life. An 8.5-kb genomic DNA fragment 5' to the start codon, containing three regions of high cross-species homology, drove expression in all embryonic and neonatal HCs, and adult vestibular and inner HCs, but not adult outer HCs. Transgenes with 0.4, 0.8, 2.5, or 6.5 kb of 5' DNA did not produce HC expression. However, addition of the region from 6.5 to 7.2 kb produced expression in vestibular HCs and neonatal basal turn outer HCs, which also implicated the region from 7.2 to 8.5 kb in inner and apical outer HC expression. Deletion of the region from 0.4 to 5.5 kb 5' from the 8.5-kb construct did not affect HC expression, further indicating lack of HC regulatory elements. When the region from 1 to 0.4 kb was replaced with the minimal promoter of the Ela1 gene, HC expression was maintained but at a drastically reduced level. Bioinformatics identified regions of highly conserved sequence outside of the 8.5 kb, which contained POU4F3-, GFI1-, and LHX3-binding sites. These regions may be involved in maintaining POU4F3 expression in adult outer HCs. Our results identify separate enhancers at various locations that direct expression to different HC types at different ages and determine that 0.4 kb of upstream sequence determines expression level. These data will assist in the identification of mutations in noncoding, regulatory regions of this deafness gene.

Combinatorial control of plant gene expression

Combinatorial gene regulation provides a mechanism by which relatively small numbers of transcription factors can control the expression of a much larger number of genes with finely tuned temporal and spatial patterns. This is achieved by transcription factors assembling into complexes in a combinatorial fashion, exponentially increasing the number of genes that they can target. Such an arrangement also increases the specificity and affinity for the cis-regulatory sequences required for accurate target gene expression. Superimposed on this transcription factor combinatorial arrangement is the increasing realization that histone modification marks expand the regulatory information, which is interpreted by histone readers and writers that are part of the regulatory apparatus. Here, we review the progress in these areas from the perspective of plant combinatorial gene regulation, providing examples of different regulatory solutions and comparing them to other metazoans. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Towards genome-wide prediction and characterization of enhancers in plants

Enhancers are important cis-regulatory DNA elements that regulate transcription programs by recruiting transcription factors and directing them to the promoters of target genes in a cell-type/tissue-specific manner. The expression of a gene can be regulated by one or multiple enhancers at different developmental stages and/or in different tissues. Enhancers are difficult to identify because of their unpredictable positions relative to their cognate promoters. Remarkably, only a handful of enhancers have been identified in plant species largely due to the lack of general approaches for enhancer identification. Extensive genomic and epigenomic research in mammalian species has revealed that the genomic locations of enhancers can be predicted based on the binding sites of transcriptional co-factors and several distinct features associated with open chromatin. Here we review the methodologies used in enhancer prediction in mammalian species. We also review the recent applications of these methodologies in Arabidopsis thaliana and discuss the future directions of enhancer identification in plants. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Transcriptional, epigenetic and retroviral signatures identify regulatory regions involved in hematopoietic lineage commitment

Genome-wide approaches allow investigating the molecular circuitry wiring the genetic and epigenetic programs of human somatic stem cells. Hematopoietic stem/progenitor cells (HSPC) give rise to the different blood cell types; however, the molecular basis of human hematopoietic lineage commitment is poorly characterized. Here, we define the transcriptional and epigenetic profile of human HSPC and early myeloid and erythroid progenitors by a combination of Cap Analysis of Gene Expression (CAGE), ChIP-seq and Moloney leukemia virus (MLV) integration site mapping. Most promoters and transcripts were shared by HSPC and committed progenitors, while enhancers and super-enhancers consistently changed upon differentiation, indicating that lineage commitment is essentially regulated by enhancer elements. A significant fraction of CAGE promoters differentially expressed upon commitment were novel, harbored a chromatin enhancer signature, and may identify promoters and transcribed enhancers driving cell commitment. MLV-targeted genomic regions co-mapped with cell-specific active enhancers and super-enhancers. Expression analyses, together with an enhancer functional assay, indicate that MLV integration can be used to identify bona fide developmentally regulated enhancers. Overall, this study provides an overview of transcriptional and epigenetic changes associated to HSPC lineage commitment, and a novel signature for regulatory elements involved in cell identity.

The analysis of novel distal Cebpa enhancers and silencers using a transcriptional model reveals the complex regulatory logic of hematopoietic lineage specification

C/EBPα plays an instructive role in the macrophage-neutrophil cell-fate decision and its expression is necessary for neutrophil development. How Cebpa itself is regulated in the myeloid lineage is not known. We decoded the cis-regulatory logic of Cebpa, and two other myeloid transcription factors, Egr1 and Egr2, using a combined experimental-computational approach. With a reporter design capable of detecting both distal enhancers and silencers, we analyzed 46 putative cis-regulatory modules (CRMs) in cells representing myeloid progenitors, and derived early macrophages or neutrophils. In addition to novel enhancers, this analysis revealed a surprisingly large number of silencers. We determined the regulatory roles of 15 potential transcriptional regulators by testing 32,768 alternative sequence-based transcriptional models against CRM activity data. This comprehensive analysis allowed us to infer the cis-regulatory logic for most of the CRMs. Silencer-mediated repression of Cebpa was found to be effected mainly by TFs expressed in non-myeloid lineages, highlighting a previously unappreciated contribution of long-distance silencing to hematopoietic lineage resolution. The repression of Cebpa by multiple factors expressed in alternative lineages suggests that hematopoietic genes are organized into densely interconnected repressive networks instead of hierarchies of mutually repressive pairs of pivotal TFs. More generally, our results demonstrate that de novo cis-regulatory dissection is feasible on a large scale with the aid of transcriptional modeling. Current address: Department of Biology, University of North Dakota, 10 Cornell Street, Stop 9019, Grand Forks, ND 58202-9019, USA.

An ankylosing spondylitis-associated genetic variant in the IL23R-IL12RB2 intergenic region modulates enhancer activity and is associated with increased Th1-cell differentiation

Objectives

To explore the functional basis for the association between ankylosing spondylitis (AS) and single-nucleotide polymorphisms (SNPs) in the IL23R-IL12RB2 intergenic region.

Methods

We performed conditional analysis on genetic association data and used epigenetic data on chromatin remodelling and transcription factor (TF) binding to identify the primary AS-associated IL23R-IL12RB2 intergenic SNP. Functional effects were tested in luciferase reporter assays in HEK293T cells and allele-specific TF binding was investigated by electrophoretic mobility gel shift assays. IL23R and IL12RB2 mRNA levels in CD4+ T cells were compared between cases homozygous for the AS-risk ‘A’ allele and the protective ‘G’ allele. The proportions of interleukin (IL)-17A+ and interferon (IFN)-γ+ CD4+ T-cells were measured by fluorescence-activated cell sorting and compared between these AS-risk and protective genotypes.

Results

Conditional analysis identified rs11209032 as the probable causal SNP within a 1.14 kb putative enhancer between IL23R and IL12RB2. Reduced luciferase activity was seen for the risk allele (p<0.001) and reduced H3K4me1 methylation observed in CD4+ T-cells from ‘A/A’ homozygotes (p=0.02). The binding of nuclear extract to the risk allele was decreased ∼3.5-fold compared with the protective allele (p<0.001). The proportion of IFN-γ+ CD4+ T-cells was increased in ‘A/A’ homozygotes (p=0.004), but neither IL23R nor IL12RB2 mRNA was affected.

Conclusions

The rs11209032 SNP downstream of IL23R forms part of an enhancer, allelic variation of which may influence Th1-cell numbers. Homozygosity for the risk ‘A’ allele is associated with more IFN-γ-secreting (Th1) cells. Further work is necessary to explain the mechanisms for these important observations.

Long-range looping of a locus control region drives tissue-specific chromatin packing within a multigene cluster

The relationships of higher order chromatin organization to mammalian gene expression remain incompletely defined. The human Growth Hormone (hGH) multigene cluster contains five gene paralogs. These genes are selectively activated in either the pituitary or the placenta by distinct components of a remote locus control region (LCR). Prior studies have revealed that appropriate activation of the placental genes is dependent not only on the actions of the LCR, but also on the multigene composition of the cluster itself. Here, we demonstrate that the hGH LCR 'loops' over a distance of 28 kb in primary placental nuclei to make specific contacts with the promoters of the two GH genes in the cluster. This long-range interaction sequesters the GH genes from the three hCS genes which co-assemble into a tightly packed 'hCS chromatin hub'. Elimination of the long-range looping, via specific deletion of the placental LCR components, triggers a dramatic disruption of the hCS chromatin hub. These data reveal a higher-order structural pathway by which long-range looping from an LCR impacts on local chromatin architecture that is linked to tissue-specific gene regulation within a multigene cluster.

Identification of a plant isoflavonoid that causes biliary atresia

Biliary atresia (BA) is a rapidly progressive and destructive fibrotic disorder of unknown etiology affecting the extrahepatic biliary tree of neonates. Epidemiological studies suggest that an environmental factor, such as a virus or toxin, is the cause of the disease, although none have been definitively established. Several naturally occurring outbreaks of BA in Australian livestock have been associated with the ingestion of unusual plants by pregnant animals during drought conditions. We used a biliary secretion assay in zebrafish to isolate a previously undescribed isoflavonoid, biliatresone, from Dysphania species implicated in a recent BA outbreak. This compound caused selective destruction of the extrahepatic, but not intrahepatic, biliary system of larval zebrafish. A mutation that enhanced biliatresone toxicity mapped to a region of the zebrafish genome that has conserved synteny with an established human BA susceptibility locus. The toxin also caused loss of cilia in neonatal mouse extrahepatic cholangiocytes in culture and disrupted cell polarity and monolayer integrity in cholangiocyte spheroids. Together, these findings provide direct evidence that BA could be initiated by perinatal exposure to an environmental toxin.

Exploring the mechanisms of genome-wide long-range interactions: interpreting chromosome organization

Developments in chromosome conformation capture (3C) technologies have revealed that the three-dimensional organization of a genome leads widely separated functional elements to reside in close proximity. However, the mechanisms responsible for mediating long-range interactions are still not completely known. In this review, we firstly evaluate and compare the current seven 3C-based methods, summarize their advantages and discuss their limitations to our current understanding of genome structure. Then, software packages available to perform the analysis of 3C-based data are described. Moreover, we review the insights into the two main mechanisms of long-range interactions, which regulate gene expression by bringing together promoters and distal regulatory elements and by creating structural domains that contain functionally related genes with similar expression landscape. At last, we summarize what is known about the mediating factors involved in stimulation/repression of long-range interactions, such as transcription factors and noncoding RNAs.

Epigenetics/Programming in the HPA Axis

The hypothalamic-pituitary-adrenal axis provides physiological adaptations to various environmental stimuli in mammals. These stimuli including maternal care, diet, immune challenge, stress, and others have the potential to stably modify or program the functioning of the HPA axis when experienced early in life or at later critical stages of development. Epigenetic mechanisms mediate the biological embedding of environmental stimuli or conditions. These changes are influenced by the genotype and both, environment and genotype contribute to the development of a specific phenotype with regard to the stress response that might be more susceptible or resilient to the development of mental conditions. The effects of stress might be a result of cumulative stress or a mismatch between the environments experienced early in life versus the conditions much later. These effects including the associated epigenetic modifications are potentially reversible.

GenomicInteractions: An R/Bioconductor package for manipulating and investigating chromatin interaction data

Background

Precise quantitative and spatiotemporal control of gene expression is necessary to ensure proper cellular differentiation and the maintenance of homeostasis. The relationship between gene expression and the spatial organisation of chromatin is highly complex, interdependent and not completely understood. The development of experimental techniques to interrogate both the higher-order structure of chromatin and the interactions between regulatory elements has recently lead to important insights on how gene expression is controlled. The ability to gain these and future insights is critically dependent on computational tools for the analysis and visualisation of data produced by these techniques.

Results and conclusion

We have developed GenomicInteractions, a freely available R/Bioconductor package designed for processing, analysis and visualisation of data generated from various types of chromosome conformation capture experiments. The package allows the easy annotation and summarisation of large genome-wide datasets at both the level of individual interactions and sets of genomic features, and provides several different methods for interrogating and visualising this type of data. We demonstrate this package’s utility by showing example analyses performed on interaction datasets generated using Hi-C and ChIA-PET.

Electronic supplementary material

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

Functionally conserved enhancers with divergent sequences in distant vertebrates

BACKGROUND

To examine the contributions of sequence and function conservation in the evolution of enhancers, we systematically identified enhancers whose sequences are not conserved among distant groups of vertebrate species, but have homologous function and are likely to be derived from a common ancestral sequence. Our approach combined comparative genomics and epigenomics to identify potential enhancer sequences in the genomes of three groups of distantly related vertebrate species.

RESULTS

We searched for sequences that were conserved within groups of closely related species but not between groups of more distant species, and were associated with an epigenetic mark of enhancer activity. To facilitate inferring orthology between non-conserved sequences, we limited our search to introns whose orthology could be unambiguously established by mapping the bracketing exons. We show that a subset of these non-conserved but syntenic sequences from the mouse and zebrafish genomes have homologous functions in a zebrafish transgenic enhancer assay. The conserved expression patterns driven by these enhancers are probably associated with short transcription factor-binding motifs present in the divergent sequences.

CONCLUSIONS

We have identified numerous potential enhancers with divergent sequences but a conserved function. These results indicate that selection on function, rather than sequence, may be a common mode of enhancer evolution; evidence for selection at the sequence level is not a necessary criterion to define a gene regulatory element.