Genome-wide maps of transcription regulatory elements

Evaluation of Whole-Exome Enrichment Solutions: Lessons from the High-End of the Short-Read Sequencing Scale

Whole-exome sequencing has become a popular technique in research and clinical settings, assisting in disease diagnosis and increasing the understanding of disease pathogenesis. In this study, we aimed to compare common enrichment capture solutions available in the market. Peripheral blood-purified DNA samples were enriched with SureSelect^QXT V6 (Agilent) and various Illumina solutions: TruSeq DNA Nano, TruSeq DNA Exome, Nextera DNA Exome, and Illumina DNA Prep with Enrichment, and sequenced on a HiSeq 4000. We found that their percentage of duplicate reads was as much as 2 times higher than previously reported values for the previous HiSeq series. SureSelect^QXT and Illumina DNA Prep with Enrichment showed the best average on-target coverage, which improved when off-target regions were included. At high coverage levels and in shared bases, these two solutions and TruSeq DNA Exome provided three of the best performances. With respect to the number of small variants detected, SureSelect^QXT presented the lowest number of detected variants in target regions. When off-target regions were considered, its ability equalized to other solutions. Our results show SureSelect^QXT and Illumina DNA Prep with Enrichment to be the best enrichment capture solutions.

Initial high-resolution microscopic mapping of active and inactive regulatory sequences proves non-random 3D arrangements in chromatin domain clusters

BACKGROUND

The association of active transcription regulatory elements (TREs) with DNAse I hypersensitivity (DHS[+]) and an 'open' local chromatin configuration has long been known. However, the 3D topography of TREs within the nuclear landscape of individual cells in relation to their active or inactive status has remained elusive. Here, we explored the 3D nuclear topography of active and inactive TREs in the context of a recently proposed model for a functionally defined nuclear architecture, where an active and an inactive nuclear compartment (ANC-INC) form two spatially co-aligned and functionally interacting networks.

RESULTS

Using 3D structured illumination microscopy, we performed 3D FISH with differently labeled DNA probe sets targeting either sites with DHS[+], apparently active TREs, or DHS[-] sites harboring inactive TREs. Using an in-house image analysis tool, DNA targets were quantitatively mapped on chromatin compaction shaped 3D nuclear landscapes. Our analyses present evidence for a radial 3D organization of chromatin domain clusters (CDCs) with layers of increasing chromatin compaction from the periphery to the CDC core. Segments harboring active TREs are significantly enriched at the decondensed periphery of CDCs with loops penetrating into interchromatin compartment channels, constituting the ANC. In contrast, segments lacking active TREs (DHS[-]) are enriched toward the compacted interior of CDCs (INC).

CONCLUSIONS

Our results add further evidence in support of the ANC-INC network model. The different 3D topographies of DHS[+] and DHS[-] sites suggest positional changes of TREs between the ANC and INC depending on their functional state, which might provide additional protection against an inappropriate activation. Our finding of a structural organization of CDCs based on radially arranged layers of different chromatin compaction levels indicates a complex higher-order chromatin organization beyond a dichotomic classification of chromatin into an 'open,' active and 'closed,' inactive state.

Evidence of non-pancreatic beta cell-dependent roles of Tcf7l2 in the regulation of glucose metabolism in mice

Non-coding variation within TCF7L2 remains the strongest genetic determinant of type 2 diabetes risk in humans. A considerable effort has been placed in understanding the functional roles of TCF7L2 in pancreatic beta cells, despite evidence of TCF7L2 expression in various peripheral tissues important in glucose homeostasis. Here, we use a humanized mouse model overexpressing Tcf7l2, resulting in glucose intolerance, to infer the contribution of Tcf7l2 overexpression in beta cells and in other tissues to the metabolic phenotypes displayed by these mice. Restoring Tcf7l2 expression specifically in beta cells to endogenous levels, in face of its overexpression elsewhere, results in impaired insulin secretion, reduced beta cell number and islet area, corroborating data obtained in humans showing similar phenotypes as a result of manipulations leading to Tcf7l2 loss of function. Interestingly, the persistent overexpression of Tcf7l2 in non-pancreatic tissues results in a significant worsening in glucose tolerance in vivo, indicating that Tcf7l2 overexpression in beta cells does not account for the glucose intolerance in the Tcf7l2 overexpression mouse model. Collectively, these data posit that Tcf7l2 plays key roles in glucose metabolism through actions beyond pancreatic beta cells, and further points to functionally opposing cell-type specific effects for Tcf7l2 on the maintenance of balanced glucose metabolism, thereby urging a careful examination of its role in non-pancreatic tissues as well as its composite metabolic effects across distinct tissues. Uncovering these roles may lead to new therapeutic targets for type 2 diabetes.

TPC1 has two variant isoforms, and their removal has different effects on endo-lysosomal functions compared to loss of TPC2

Organelle ion homeostasis within the endo-lysosomal system is critical for physiological functions. Two-pore channels (TPCs) are cation channels that reside in endo-lysosomal organelles, and overexpression results in endo-lysosomal trafficking defects. However, the impact of a lack of TPC expression on endo-lysosomal trafficking is unknown. Here, we characterize Tpcn1 expression in two transgenic mouse lines (Tpcn1^XG716 and Tpcn1^T159) and show expression of a novel evolutionarily conserved Tpcn1B transcript from an alternative promoter, raising important questions regarding the status of Tpcn1 expression in mice recently described to be Tpcn1 knockouts. We show that the transgenic Tpcn1^T159 line lacks expression of both Tpcn1 isoforms in all tissues analyzed. Using mouse embryonic fibroblasts (MEFs) from Tpcn1^−/− and Tpcn2^−/− animals, we show that a lack of Tpcn1 or Tpcn2 expression has no significant impact on resting endo-lysosomal pH or morphology. However, differential effects in endo-lysosomal function were observed upon the loss of Tpcn1 or Tpcn2 expression; thus, while Tpcn1^−/− MEFs have impaired trafficking of cholera toxin from the plasma membrane to the Golgi apparatus, Tpcn2^−/− MEFs show slower kinetics of ligand-induced platelet-derived growth factor receptor β (PDGFRβ) degradation, which is dependent on trafficking to lysosomes. Our findings indicate that TPC1 and TPC2 have important but distinct roles in the endo-lysosomal pathway.

The human imprintome: regulatory mechanisms, methods of ascertainment, and roles in disease susceptibility

Imprinted genes form a special subset of the genome, exhibiting monoallelic expression in a parent-of-origin-dependent fashion. This monoallelic expression is controlled by parental-specific epigenetic marks, which are established in gametogenesis and early embryonic development and are persistent in all somatic cells throughout life. We define this specific set of cis-acting epigenetic regulatory elements as the imprintome, a distinct and specially tasked subset of the epigenome. Imprintome elements contain DNA methylation and histone modifications that regulate monoallelic expression by affecting promoter accessibility, chromatin structure, and chromatin configuration. Understanding their regulation is critical because a significant proportion of human imprinted genes are implicated in complex diseases. Significant species variation in the repertoire of imprinted genes and their epigenetic regulation, however, will not allow model organisms solely to be used for this crucial purpose. Ultimately, only the human will suffice to accurately define the human imprintome.

Beyond the ENCODE project: using genomics and epigenomics strategies to study enhancer evolution

The complex expression patterns observed for many genes are often regulated by distal transcription enhancers. Changes in the nucleotide sequences of enhancers may therefore lead to changes in gene expression, representing a central mechanism by which organisms evolve. With the development of the experimental technique of chromatin immunoprecipitation (ChIP), in which discrete regions of the genome bound by specific proteins can be identified, it is now possible to identify transcription factor binding events (putative cis-regulatory elements) in entire genomes. Comparing protein-DNA binding maps allows us, for the first time, to attempt to identify regulatory differences and infer global patterns of change in gene expression across species. Here, we review studies that used genome-wide ChIP to study the evolution of enhancers. The trend is one of high divergence of cis-regulatory elements between species, possibly compensated by extensive creation and loss of regulatory elements and rewiring of their target genes. We speculate on the meaning of the differences observed and discuss that although ChIP experiments identify the biochemical event of protein-DNA interaction, it cannot determine whether the event results in a biological function, and therefore more studies are required to establish the effect of divergence of binding events on species-specific gene expression.

Detection of regulatory SNPs in human genome using ChIP-seq ENCODE data

A vast amount of SNPs derived from genome-wide association studies are represented by non-coding ones, therefore exacerbating the need for effective identification of regulatory SNPs (rSNPs) among them. However, this task remains challenging since the regulatory part of the human genome is annotated much poorly as opposed to coding regions. Here we describe an approach aggregating the whole set of ENCODE ChIP-seq data in order to search for rSNPs, and provide the experimental evidence of its efficiency. Its algorithm is based on the assumption that the enrichment of a genomic region with transcription factor binding loci (ChIP-seq peaks) indicates its regulatory function, and thereby SNPs located in this region are more likely to influence transcription regulation. To ensure that the approach preferably selects functionally meaningful SNPs, we performed enrichment analysis of several human SNP datasets associated with phenotypic manifestations. It was shown that all samples are significantly enriched with SNPs falling into the regions of multiple ChIP-seq peaks as compared with the randomly selected SNPs. For experimental verification, 40 SNPs falling into overlapping regions of at least 7 TF binding loci were selected from OMIM. The effect of SNPs on the binding of the DNA fragments containing them to the nuclear proteins from four human cell lines (HepG2, HeLaS3, HCT-116, and K562) has been tested by EMSA. A radical change in the binding pattern has been observed for 29 SNPs, besides, 6 more SNPs also demonstrated less pronounced changes. Taken together, the results demonstrate the effective way to search for potential rSNPs with the aid of ChIP-seq data provided by ENCODE project.

Breakpoint regions of ETO gene involved in (8; 21) leukemic translocations are enriched in acetylated histone H3

One of the most frequent chromosomal translocation found in patients with acute myeloid leukemia (AML) is the t(8;21). This translocation involves the RUNX1 and ETO genes. The breakpoints regions for t(8;21) are located at intron 5 and intron 1 of the RUNX1 and ETO gene respectively. To date, no homologous sequences have been found in these regions to explain their recombination. The breakpoint regions of RUNX1 gene are characterized by the presence of DNasaI hypersensitive sites and topoisomerase II cleavage sites, but no information exists about complementary regions of ETO gene. Here, we report analysis of chromatin structure of ETO breakpoint regions. Chromatin immunoprecipitation (ChIP) were performed with antibodies specific to acetylated histone H3, H4, and total histone H1. Nucleosomal distribution at the ETO locus was evaluated by determining total levels of histone H3. Our data show that in myeloid cells, the breakpoint regions at the ETO gene are enriched in hyperacetylated histone H3 compared to a control region of similar size where no translocations have been described. Moreover, acetylated H4 associates with both the whole ETO breakpoint regions as well as the control intron. Interestingly, we observed no H1 association either at the breakpoint regions or the control region of the ETO gene. Our data indicate that a common chromatin structure enriched in acetylated histones is present in breakpoint regions involved in formation of (8;21) leukemic translocation.

Sequence signatures extracted from proximal promoters can be used to predict distal enhancers

BACKGROUND

Gene expression is controlled by proximal promoters and distal regulatory elements such as enhancers. While the activity of some promoters can be invariant across tissues, enhancers tend to be highly tissue-specific.

RESULTS

We compiled sets of tissue-specific promoters based on gene expression profiles of 79 human tissues and cell types. Putative transcription factor binding sites within each set of sequences were used to train a support vector machine classifier capable of distinguishing tissue-specific promoters from control sequences. We obtained reliable classifiers for 92% of the tissues, with an area under the receiver operating characteristic curve between 60% (for subthalamic nucleus promoters) and 98% (for heart promoters). We next used these classifiers to identify tissue-specific enhancers, scanning distal non-coding sequences in the loci of the 200 most highly and lowly expressed genes. Thirty percent of reliable classifiers produced consistent enhancer predictions, with significantly higher densities in the loci of the most highly expressed compared to lowly expressed genes. Liver enhancer predictions were assessed in vivo using the hydrodynamic tail vein injection assay. Fifty-eight percent of the predictions yielded significant enhancer activity in the mouse liver, whereas a control set of five sequences was completely negative.

CONCLUSIONS

We conclude that promoters of tissue-specific genes often contain unambiguous tissue-specific signatures that can be learned and used for the de novo prediction of enhancers.

Capturing the regulatory interactions of eukaryote genomes

A key finding from early genomics research is the remarkable consistency in the number of protein-coding regions across diverse species. This has led many researchers to look to the cis-regulatory elements of genes as the fundamental influence behind evolving gene function and subsequent species diversification. Historically, since these elements are often located in vast intergenic and intronic regions of the genome, their identification has been recalcitrant. Now, with the deluge of whole-genome data from representatives of numerous eukaryotic lineages, various approaches have enabled us to begin to recognize features that characterize regulatory regions of the genome. Here we endeavour to collate these approaches in order to give an overview of the complexities involved in extrapolating regulatory signatures. The resource provided by the escalating richness of whole-genome datasets enables more sophisticated modelling of these regulatory signatures yet at the same time introduces increasing potential for noise. While we are only at the advent of making these discoveries, the next decade promises to be a very exciting and rewarding time for genome researchers.

Dynamics of enhancer chromatin signatures mark the transition from pluripotency to cell specification during embryogenesis

The generation of distinctive cell types that form different tissues and organs requires precise, temporal and spatial control of gene expression. This depends on specific cis-regulatory elements distributed in the noncoding DNA surrounding their target genes. Studies performed on mammalian embryonic stem cells and Drosophila embryos suggest that active enhancers form part of a defined chromatin landscape marked by histone H3 lysine 4 mono-methylation (H3K4me1) and histone H3 lysine 27 acetylation (H3K27ac). Nevertheless, little is known about the dynamics and the potential roles of these marks during vertebrate embryogenesis. Here, we provide genomic maps of H3K4me1/me3 and H3K27ac at four developmental time-points of zebrafish embryogenesis and analyze embryonic enhancer activity. We find that (1) changes in H3K27ac enrichment at enhancers accompany the shift from pluripotency to tissue-specific gene expression, (2) in early embryos, the peaks of H3K27ac enrichment are bound by pluripotent factors such as Nanog, and (3) the degree of evolutionary conservation is higher for enhancers that become marked by H3K27ac at the end of gastrulation, suggesting their implication in the establishment of the most conserved (phylotypic) transcriptome that is known to occur later at the pharyngula stage.

Comparative studies of gene expression and the evolution of gene regulation

The hypothesis that differences in gene regulation have an important role in speciation and adaptation is more than 40 years old. With the advent of new sequencing technologies, we are able to characterize and study gene expression levels and associated regulatory mechanisms in a large number of individuals and species at an unprecedented resolution and scale. We have thus gained new insights into the evolutionary pressures that shape gene expression levels and have developed an appreciation for the relative importance of evolutionary changes in different regulatory genetic and epigenetic mechanisms. The current challenge is to link gene regulatory changes to adaptive evolution of complex phenotypes. Here we mainly focus on comparative studies in primates and how they are complemented by studies in model organisms.

Transcriptional enhancers in development and disease

Distal transcription enhancers are cis-regulatory elements that promote gene expression, enabling spatiotemporal control of genetic programs such as those required in metazoan developmental processes. Because of their importance, their disruption can lead to disease.

Methods to detect selection on noncoding DNA

Vast tracts of noncoding DNA contain elements that regulate gene expression in higher eukaryotes. Describing these regulatory elements and understanding how they evolve represent major challenges for biologists. Advances in the ability to survey genome-scale DNA sequence data are providing unprecedented opportunities to use evolutionary models and computational tools to identify functionally important elements and the mode of selection acting on them in multiple species. This chapter reviews some of the current methods that have been developed and applied on noncoding DNA, what they have shown us, and how they are limited. Results of several recent studies reveal that a significantly larger fraction of noncoding DNA in eukaryotic organisms is likely to be functional than previously believed, implying that the functional annotation of most noncoding DNA in these organisms is largely incomplete. In Drosophila, recent studies have further suggested that a large fraction of noncoding DNA divergence observed between species may be the product of recurrent adaptive substitution. Similar studies in humans have revealed a more complex pattern, with signatures of recurrent positive selection being largely concentrated in conserved noncoding DNA elements. Understanding these patterns and the extent to which they generalize to other organisms awaits the analysis of forthcoming genome-scale polymorphism and divergence data from more species.

Distal interleukin-1β (IL-1β) response element of human matrix metalloproteinase-13 (MMP-13) binds activator protein 1 (AP-1) transcription factors and regulates gene expression

The collagenase matrix metalloproteinase-13 (MMP-13) plays an important role in the destruction of cartilage in arthritic joints. MMP-13 expression is strongly up-regulated in arthritis, largely because of stimulation by inflammatory cytokines such as IL-1β. Treatment of chondrocytes with IL-1β induces transcription of MMP-13 in vitro. IL-1β signaling converges upon the activator protein-1 transcription factors, which have been shown to be required for IL-1β-induced MMP-13 gene expression. Using chromatin immunoprecipitation (ChIP), we detected activator protein-1 binding within an evolutionarily conserved DNA sequence ∼20 kb 5' relative to the MMP-13 transcription start site (TSS). Also using ChIP, we detected histone modifications and binding of RNA polymerase II within this conserved region, all of which are consistent with transcriptional activation. Chromosome conformation capture indicates that chromosome looping brings this region in close proximity with the MMP-13 TSS. Finally, a luciferase reporter construct driven by a component of the conserved region demonstrated an expression pattern similar to that of endogenous MMP-13. These data suggest that a conserved region at 20 kb upstream from the MMP-13 TSS includes a distal transcriptional response element of MMP-13, which contributes to MMP-13 gene expression.

Global mapping of cell type-specific open chromatin by FAIRE-seq reveals the regulatory role of the NFI family in adipocyte differentiation

Identification of regulatory elements within the genome is crucial for understanding the mechanisms that govern cell type-specific gene expression. We generated genome-wide maps of open chromatin sites in 3T3-L1 adipocytes (on day 0 and day 8 of differentiation) and NIH-3T3 fibroblasts using formaldehyde-assisted isolation of regulatory elements coupled with high-throughput sequencing (FAIRE-seq). FAIRE peaks at the promoter were associated with active transcription and histone modifications of H3K4me3 and H3K27ac. Non-promoter FAIRE peaks were characterized by H3K4me1+/me3-, the signature of enhancers, and were largely located in distal regions. The non-promoter FAIRE peaks showed dynamic change during differentiation, while the promoter FAIRE peaks were relatively constant. Functionally, the adipocyte- and preadipocyte-specific non-promoter FAIRE peaks were, respectively, associated with genes up-regulated and down-regulated by differentiation. Genes highly up-regulated during differentiation were associated with multiple clustered adipocyte-specific FAIRE peaks. Among the adipocyte-specific FAIRE peaks, 45.3% and 11.7% overlapped binding sites for, respectively, PPARγ and C/EBPα, the master regulators of adipocyte differentiation. Computational motif analyses of the adipocyte-specific FAIRE peaks revealed enrichment of a binding motif for nuclear family I (NFI) transcription factors. Indeed, ChIP assay showed that NFI occupy the adipocyte-specific FAIRE peaks and/or the PPARγ binding sites near PPARγ, C/EBPα, and aP2 genes. Overexpression of NFIA in 3T3-L1 cells resulted in robust induction of these genes and lipid droplet formation without differentiation stimulus. Overexpression of dominant-negative NFIA or siRNA-mediated knockdown of NFIA or NFIB significantly suppressed both induction of genes and lipid accumulation during differentiation, suggesting a physiological function of these factors in the adipogenic program. Together, our study demonstrates the utility of FAIRE-seq in providing a global view of cell type-specific regulatory elements in the genome and in identifying transcriptional regulators of adipocyte differentiation.

Appendage expression driven by the Hoxd Global Control Region is an ancient gnathostome feature

The evolutionary transition of the fins of fish into tetrapod limbs involved genetic changes to developmental systems that resulted in novel skeletal patterns and functions. Approaches to understanding this issue have entailed the search for antecedents of limb structure in fossils, genes, and embryos. Comparative genetic analyses have produced ambiguous results: although studies of posterior Hox genes from homology group 13 (Hoxa-13 and Hoxd-13) reveal similarities in gene expression between the distal segments of fins and limbs, this functional homology has not been supported by genomic comparisons of the activity of their cis-regulatory elements, namely the Hoxd Global Control Region. Here, we show that cis-regulatory elements driving Hoxd gene expression in distal limbs are present in fish. Using an interspecies transgenesis approach, we find functional conservation between gnathostome Hoxd enhancers, demonstrating that orthologous sequences from tetrapods, zebrafish and skate can drive reporter gene expression in mouse limbs and zebrafish fins. Our results support the notion that some of the novelties associated with tetrapod limbs arose by modification of deeply conserved cis- and trans-acting mechanisms of Hox regulation in gnathostomes.