Noboru Jo Sakabe

Obesity-associated variants within FTO form long-range functional connections with IRX3

Genome-wide association studies (GWAS) have reproducibly associated variants within introns of FTO with increased risk for obesity and type 2 diabetes (T2D). Although the molecular mechanisms linking these noncoding variants with obesity are not immediately obvious, subsequent studies in mice demonstrated that FTO expression levels influence body mass and composition phenotypes. However, no direct connection between the obesity-associated variants and FTO expression or function has been made. Here we show that the obesity-associated noncoding sequences within FTO are functionally connected, at megabase distances, with the homeobox gene IRX3. The obesity-associated FTO region directly interacts with the promoters of IRX3 as well as FTO in the human, mouse and zebrafish genomes. Furthermore, long-range enhancers within this region recapitulate aspects of IRX3 expression, suggesting that the obesity-associated interval belongs to the regulatory landscape of IRX3. Consistent with this, obesity-associated single nucleotide polymorphisms are associated with expression of IRX3, but not FTO, in human brains. A direct link between IRX3 expression and regulation of body mass and composition is demonstrated by a reduction in body weight of 25 to 30% in Irx3-deficient mice, primarily through the loss of fat mass and increase in basal metabolic rate with browning of white adipose tissue. Finally, hypothalamic expression of a dominant-negative form of Irx3 reproduces the metabolic phenotypes of Irx3-deficient mice. Our data suggest that IRX3 is a functional long-range target of obesity-associated variants within FTO and represents a novel determinant of body mass and composition.

Genome-wide discovery of human heart enhancers

The various organogenic programs deployed during embryonic development rely on the precise expression of a multitude of genes in time and space. Identifying the cis-regulatory elements responsible for this tightly orchestrated regulation of gene expression is an essential step in understanding the genetic pathways involved in development. We describe a strategy to systematically identify tissue-specific cis-regulatory elements that share combinations of sequence motifs. Using heart development as an experimental framework, we employed a combination of Gibbs sampling and linear regression to build a classifier that identifies heart enhancers based on the presence and/or absence of various sequence features, including known and putative transcription factor (TF) binding specificities. In distinguishing heart enhancers from a large pool of random noncoding sequences, the performance of our classifier is vastly superior to four commonly used methods, with an accuracy reaching 92% in cross-validation. Furthermore, most of the binding specificities learned by our method resemble the specificities of TFs widely recognized as key players in heart development and differentiation, such as SRF, MEF2, ETS1, SMAD, and GATA. Using our classifier as a predictor, a genome-wide scan identified over 40,000 novel human heart enhancers. Although the classifier used no gene expression information, these novel enhancers are strongly associated with genes expressed in the heart. Finally, in vivo tests of our predictions in mouse and zebrafish achieved a validation rate of 62%, significantly higher than what is expected by chance. These results support the existence of underlying cis-regulatory codes dictating tissue-specific transcription in mammalian genomes and validate our enhancer classifier strategy as a method to uncover these regulatory codes.

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.

Tbx20 regulates a genetic program essential to adult mouse cardiomyocyte function

Human mutations in or variants of TBX20 are associated with congenital heart disease, cardiomyopathy, and arrhythmias. To investigate whether cardiac disease in patients with these conditions results from an embryonic or ongoing requirement for Tbx20 in myocardium, we ablated Tbx20 specifically in adult cardiomyocytes in mice. This ablation resulted in the onset of severe cardiomyopathy accompanied by arrhythmias, with death ensuing within 1 to 2 weeks of Tbx20 ablation. Accounting for this dramatic phenotype, we identified molecular signatures that posit Tbx20 as a central integrator of a genetic program that maintains cardiomyocyte function in the adult heart. Expression of a number of genes encoding critical transcription factors, ion channels, and cytoskeletal/myofibrillar proteins was downregulated consequent to loss of Tbx20. Genome-wide ChIP analysis of Tbx20-binding regions in the adult heart revealed that many of these genes were direct downstream targets of Tbx20 and uncovered a previously undescribed DNA-binding site for Tbx20. Bioinformatics and in vivo functional analyses revealed a cohort of transcription factors that, working with Tbx20, integrated multiple environmental signals to maintain ion channel gene expression in the adult heart. Our data provide insight into the mechanisms by which mutations in TBX20 cause adult heart disease in humans.

Dual transcriptional activator and repressor roles of TBX20 regulate adult cardiac structure and function

The ongoing requirement in adult heart for transcription factors with key roles in cardiac development is not well understood. We recently demonstrated that TBX20, a transcriptional regulator required for cardiac development, has key roles in the maintenance of functional and structural phenotypes in adult mouse heart. Conditional ablation of Tbx20 in adult cardiomyocytes leads to a rapid onset and progression of heart failure, with prominent conduction and contractility phenotypes that lead to death. Here we describe a more comprehensive molecular characterization of the functions of TBX20 in adult mouse heart. Coupling genome-wide chromatin immunoprecipitation and transcriptome analyses (RNA-Seq), we identified a subset of genes that change expression in Tbx20 adult cardiomyocyte-specific knockout hearts which are direct downstream targets of TBX20. This analysis revealed a dual role for TBX20 as both a transcriptional activator and a repressor, and that each of these functions regulates genes with very specialized and distinct molecular roles. We also show how TBX20 binds to its targets genome-wide in a context-dependent manner, using various cohorts of co-factors to either promote or repress distinct genetic programs within adult heart. Our integrative approach has uncovered several novel aspects of TBX20 and T-box protein function within adult heart. Sequencing data accession number (http://www.ncbi.nlm.nih.gov/geo): GSE30943.

HIF1α Represses Cell Stress Pathways to Allow Proliferation of Hypoxic Fetal Cardiomyocytes.

Transcriptional mediators of cell stress pathways, including HIF1α, ATF4, and p53, are key to normal development and play critical roles in disease, including ischemia and cancer. Despite their importance, mechanisms by which pathways mediated by these transcription factors interact with one another are not fully understood. In addressing the controversial role of HIF1α in cardiomyocytes (CMs) during heart development, we discovered a mid-gestational requirement for HIF1α for proliferation of hypoxic CMs, involving metabolic switching and a complex interplay among HIF1α, ATF4, and p53. Loss of HIF1α resulted in activation of ATF4 and p53, the latter inhibiting CM proliferation. Bioinformatic and biochemical analyses revealed unexpected mechanisms by which HIF1α intersects with ATF4 and p53 pathways. Our results highlight previously undescribed roles of HIF1α and interactions among major cell stress pathways that could be targeted to enhance proliferation of CMs in ischemia and may have relevance to other diseases, including cancer.

Genome‐wide maps of transcription regulatory elements

Expression of eukaryotic genes with complex spatial‐temporal regulation during development requires finer regulation than that of genes with simpler expression patterns. Given the high degree of conservation of the developmental gene set across distantly related phylogenetic taxa, it is argued that evolutionary variation has occurred by tweaking regulation of expression of developmental genes, rather than by changes in genes themselves. Complex regulation is often achieved through the coordinated action of transcription regulatory elements spread across the genome up to tens of kilobases from the promoters of their target genes. Disruption of regulatory elements has been implicated in several diseases and studies showing associations between disease traits and nonprotein coding variation hint for a role of regulatory elements as cause of diseases. Therefore, the identification and mapping of regulatory elements in genome scale is crucial to understand how gene expression is regulated, how organisms evolve, and to identify sequence variation causing diseases. Previously developed experimental techniques have been adapted to identify regulatory elements in genome scale and high‐throughput, allowing a global view of their biological roles. We review methods as chromatin immunoprecipitation, DNase I hypersensitivity, and computational approaches and how they have been employed to generate maps of histone modifications, open chromatin, nucleosome positioning, and transcription factor binding regions in whole mammalian genomes. Given the importance of non‐promoter elements in gene regulation and the recent explosion in the number of studies devoted to them, we focus on these elements and discuss the insights on gene regulation being obtained by these studies. Copyright

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.

Probing chromatin landscape reveals roles of endocardial TBX20 in septation

Mutations in the T-box transcription factor TBX20 are associated with multiple forms of congenital heart defects, including cardiac septal abnormalities, but our understanding of the contributions of endocardial TBX20 to heart development remains incomplete. Here, we investigated how TBX20 interacts with endocardial gene networks to drive the mesenchymal and myocardial movements that are essential for outflow tract and atrioventricular septation. Selective ablation of Tbx20 in murine endocardial lineages reduced the expression of extracellular matrix and cell migration genes that are critical for septation. Using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), we identified accessible chromatin within endocardial lineages and intersected these data with TBX20 ChIP-seq and chromatin loop maps to determine that TBX20 binds a conserved long-range enhancer to regulate versican (Vcan) expression. We also observed reduced Vcan expression in Tbx20-deficient mice, supporting a direct role for TBX20 in Vcan regulation. Further, we show that the Vcan enhancer drove reporter gene expression in endocardial lineages in a TBX20-binding site-dependent manner. This work illuminates gene networks that interact with TBX20 to orchestrate cardiac septation and provides insight into the chromatin landscape of endocardial lineages during septation.

Reducing mitochondrial reads in ATAC-seq using CRISPR/Cas9

ATAC-seq is a high-throughput sequencing technique that identifies open chromatin. Depending on the cell type, ATAC-seq samples may contain ~20–80% of mitochondrial sequencing reads. As the regions of open chromatin of interest are usually located in the nuclear genome, mitochondrial reads are typically discarded from the analysis. We tested two approaches to decrease wasted sequencing in ATAC-seq libraries generated from lymphoblastoid cell lines: targeted cleavage of mitochondrial DNA fragments using CRISPR technology and removal of detergent from the cell lysis buffer. We analyzed the effects of these treatments on the number of usable (unique, non-mitochondrial) reads and the number and quality of peaks called, including peaks identified in enhancers and transcription start sites. Both treatments resulted in considerable reduction of mitochondrial reads (1.7 and 3-fold, respectively). The removal of detergent, however, resulted in increased background and fewer peaks. The highest number of peaks and highest quality data was obtained by preparing samples with the original ATAC-seq protocol (using detergent) and treating them with CRISPR. This strategy reduced the amount of sequencing required to call a high number of peaks, which could lead to cost reduction when performing ATAC-seq on large numbers of samples and in cell types that contain a large amount of mitochondria.