Matthew J. Blow

ChIP-seq accurately predicts tissue-specific activity of enhancers

A major yet unresolved quest in decoding the human genome is the identification of the regulatory sequences that control the spatial and temporal expression of genes. Distant-acting transcriptional enhancers are particularly challenging to uncover because they are scattered among the vast non-coding portion of the genome. Evolutionary sequence constraint can facilitate the discovery of enhancers, but fails to predict when and where they are active in vivo. Here we present the results of chromatin immunoprecipitation with the enhancer-associated protein p300 followed by massively parallel sequencing, and map several thousand in vivo binding sites of p300 in mouse embryonic forebrain, midbrain and limb tissue. We tested 86 of these sequences in a transgenic mouse assay, which in nearly all cases demonstrated reproducible enhancer activity in the tissues that were predicted by p300 binding. Our results indicate that in vivo mapping of p300 binding is a highly accurate means for identifying enhancers and their associated activities, and suggest that such data sets will be useful to study the role of tissue-specific enhancers in human biology and disease on a genome-wide scale.

The amphioxus genome illuminates vertebrate origins and cephalochordate biology

Cephalochordates, urochordates, and vertebrates evolved from a common ancestor over 520 million years ago. To improve our understanding of chordate evolution and the origin of vertebrates, we intensively searched for particular genes, gene families, and conserved noncoding elements in the sequenced genome of the cephalochordate Branchiostoma floridae, commonly called amphioxus or lancelets. Special attention was given to homeobox genes, opsin genes, genes involved in neural crest development, nuclear receptor genes, genes encoding components of the endocrine and immune systems, and conserved cis-regulatory enhancers. The amphioxus genome contains a basic set of chordate genes involved in development and cell signaling, including a fifteenth Hox gene. This set includes many genes that were co-opted in vertebrates for new roles in neural crest development and adaptive immunity. However, where amphioxus has a single gene, vertebrates often have two, three, or four paralogs derived from two whole-genome duplication events. In addition, several transcriptional enhancers are conserved between amphioxus and vertebrates--a very wide phylogenetic distance. In contrast, urochordate genomes have lost many genes, including a diversity of homeobox families and genes involved in steroid hormone function. The amphioxus genome also exhibits derived features, including duplications of opsins and genes proposed to function in innate immunity and endocrine systems. Our results indicate that the amphioxus genome is elemental to an understanding of the biology and evolution of nonchordate deuterostomes, invertebrate chordates, and vertebrates.

Targeted deletion of the 9p21 non-coding coronary artery disease risk interval in mice

Sequence polymorphisms in a 58-kilobase (kb) interval on chromosome 9p21 confer a markedly increased risk of coronary artery disease (CAD), the leading cause of death worldwide. The variants have a substantial effect on the epidemiology of CAD and other life-threatening vascular conditions because nearly one-quarter of Caucasians are homozygous for risk alleles. However, the risk interval is devoid of protein-coding genes and the mechanism linking the region to CAD risk has remained enigmatic. Here we show that deletion of the orthologous 70-kb non-coding interval on mouse chromosome 4 affects cardiac expression of neighbouring genes, as well as proliferation properties of vascular cells. Chr4Δ70kb/Δ70kb mice are viable, but show increased mortality both during development and as adults. Cardiac expression of two genes near the non-coding interval, Cdkn2a and Cdkn2b, is severely reduced in chr4Δ70kb/Δ70kb mice, indicating that distant-acting gene regulatory functions are located in the non-coding CAD risk interval. Allele-specific expression of Cdkn2b transcripts in heterozygous mice showed that the deletion affects expression through a cis-acting mechanism. Primary cultures of chr4Δ70kb/Δ70kb aortic smooth muscle cells exhibited excessive proliferation and diminished senescence, a cellular phenotype consistent with accelerated CAD pathogenesis. Taken together, our results provide direct evidence that the CAD risk interval has a pivotal role in regulation of cardiac Cdkn2a/b expression, and suggest that this region affects CAD progression by altering the dynamics of vascular cell proliferation.

ChIP-Seq identification of weakly conserved heart enhancers

Accurate control of tissue-specific gene expression plays a pivotal role in heart development, but few cardiac transcriptional enhancers have thus far been identified. Extreme noncoding-sequence conservation has successfully predicted enhancers that are active in many tissues but has failed to identify substantial numbers of heart-specific enhancers. Here, we used ChIP-Seq with the enhancer-associated protein p300 from mouse embryonic day 11.5 heart tissue to identify over 3,000 candidate heart enhancers genome wide. Compared to enhancers active in othertissues we studied at this time point, most candidate heart enhancers were less deeply conserved in vertebrate evolution. Nevertheless, transgenic mouse assays of 130 candidate regions revealed that most function reproducibly as enhancers active in the heart, irrespective of their degree of evolutionary constraint. These results provide evidence for a large population of poorly conserved heart enhancers and suggest that the evolutionary conservation of embryonic enhancers can vary depending on tissue type.

Large-scale discovery of enhancers from human heart tissue

Development and function of the human heart depend on the dynamic control of tissue-specific gene expression by distant-acting transcriptional enhancers. To generate an accurate genome-wide map of human heart enhancers, we used an epigenomic enhancer discovery approach and identified ∼6,200 candidate enhancer sequences directly from fetal and adult human heart tissue. Consistent with their predicted function, these elements were markedly enriched near genes implicated in heart development, function and disease. To further validate their in vivo enhancer activity, we tested 65 of these human sequences in a transgenic mouse enhancer assay and observed that 43 (66%) drove reproducible reporter gene expression in the heart. These results support the discovery of a genome-wide set of noncoding sequences highly enriched in human heart enhancers that is likely to facilitate downstream studies of the role of enhancers in development and pathological conditions of the heart.

Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads

BackgroundComprehensive annotation and quantification of transcriptomes are outstanding problems in functional genomics. While high throughput mRNA sequencing (RNA-Seq) has emerged as a powerful tool for addressing these problems, its success is dependent upon the availability and quality of reference genome sequences, thus limiting the organisms to which it can be applied.ResultsHere, we describe Rnnotator, an automated software pipeline that generates transcript models by de novo assembly of RNA-Seq data without the need for a reference genome. We have applied the Rnnotator assembly pipeline to two yeast transcriptomes and compared the results to the reference gene catalogs of these organisms. The contigs produced by Rnnotator are highly accurate (95%) and reconstruct full-length genes for the majority of the existing gene models (54.3%). Furthermore, our analyses revealed many novel transcribed regions that are absent from well annotated genomes, suggesting Rnnotator serves as a complementary approach to analysis based on a reference genome for comprehensive transcriptomics.ConclusionsThese results demonstrate that the Rnnotator pipeline is able to reconstruct full-length transcripts in the absence of a complete reference genome.

Fine Tuning of Craniofacial Morphology by Distant-Acting Enhancers

Introduction The shape of the face is one of the most distinctive features among humans, and differences in facial morphology have substantial implications in areas such as social interaction, psychology, forensics, and clinical genetics. Craniofacial shape is highly heritable, including the normal spectrum of morphological variation as well as susceptibility to major craniofacial birth defects. In this study, we explored the role of transcriptional enhancers in the development of the craniofacial complex. Our study is based on the rationale that such enhancers, which can be hundreds of kilobases away from their target genes, regulate the spatial patterns, levels, and timing of gene expression in normal development. Craniofacial developmental enhancers contribute to craniofacial morphology. We identified distant-acting transcriptional enhancers active in the developing craniofacial complex and studied their activity patterns in detail in transgenic mice (left). Selected enhancers were deleted from the genome in mice in order to examine their role in modulating craniofacial morphology, which revealed subtle but significant effects of enhancers on the shape of the face and skull (right). Methods To identify distant-acting enhancers active during craniofacial development, we used chromatin immunoprecipitation on embryonic mouse face tissue followed by sequencing to identify noncoding genome regions bound by the enhancer-associated p300 protein. We used LacZ reporter assays in transgenic mice and optical projection tomography (OPT) to determine three-dimensional expression patterns of a subset of these candidate enhancers. Last, we deleted three of the craniofacial enhancers from the mouse genome to assess their effect on gene expression and craniofacial morphology during development. Results We identified more than 4000 candidate enhancer sequences predicted to be active in the developing craniofacial complex. The majority of these sequences are at least partially conserved between humans and mice, and many are located in chromosomal regions associated with normal facial morphology or craniofacial birth defects. Characterization of more than 200 candidate enhancer sequences in transgenic mice revealed a remarkable spatial complexity of in vivo expression patterns. Targeted deletions of three craniofacial enhancers near genes with known roles in craniofacial development resulted in changes of expression of those genes as well as quantitatively subtle but definable alterations of craniofacial shape. Discussion Our analysis identifies enhancers that fine tune expression of genes during craniofacial development in mice. These results support that variation in the sequence or copy number of craniofacial enhancers may contribute to the spectrum of facial variation we find in human populations. Because many craniofacial enhancers are located in genome regions associated with craniofacial birth defects, such as clefts of the lip and palate, our results also offer a starting point for exploring the contribution of noncoding sequences to these disorders. No Two Faces Are Alike Gene disruptions can cause severe dysmorphologies like cleft palate, but what causes the subtle shifts in facial morphology that make each face unique? Studying mice, Attanasio et al. (1241006) identified over 4000 candidate genetic enhancers around genes driving craniofacial development. To avoid the challenge of recognizing individual mouse faces, optical projection tomography was used to link changes in facial morphology with alterations in the function of specific enhancers. Targeted deletion of individual craniofacial enhancers from the mouse genome sculpts facial shapes. The shape of the human face and skull is largely genetically determined. However, the genomic basis of craniofacial morphology is incompletely understood and hypothesized to involve protein-coding genes, as well as gene regulatory sequences. We used a combination of epigenomic profiling, in vivo characterization of candidate enhancer sequences in transgenic mice, and targeted deletion experiments to examine the role of distant-acting enhancers in craniofacial development. We identified complex regulatory landscapes consisting of enhancers that drive spatially complex developmental expression patterns. Analysis of mouse lines in which individual craniofacial enhancers had been deleted revealed significant alterations of craniofacial shape, demonstrating the functional importance of enhancers in defining face and skull morphology. These results demonstrate that enhancers are involved in craniofacial development and suggest that enhancer sequence variation contributes to the diversity of human facial morphology.

A High-Resolution Enhancer Atlas of the Developing Telencephalon

The mammalian telencephalon plays critical roles in cognition, motor function, and emotion. Though many of the genes required for its development have been identified, the distant-acting regulatory sequences orchestrating their in vivo expression are mostly unknown. Here, we describe a digital atlas of in vivo enhancers active in subregions of the developing telencephalon. We identified more than 4,600 candidate embryonic forebrain enhancers and studied the in vivo activity of 329 of these sequences in transgenic mouse embryos. We generated serial sets of histological brain sections for 145 reproducible forebrain enhancers, resulting in a publicly accessible web-based data collection comprising more than 32,000 sections. We also used epigenomic analysis of human and mouse cortex tissue to directly compare the genome-wide enhancer architecture in these species. These data provide a primary resource for investigating gene regulatory mechanisms of telencephalon development and enable studies of the role of distant-acting enhancers in neurodevelopmental disorders.

Mutations in the DLG3 Gene Cause Nonsyndromic X-Linked Mental Retardation

We have identified truncating mutations in the human DLG3 (neuroendocrine dlg) gene in 4 of 329 families with moderate to severe X-linked mental retardation. DLG3 encodes synapse-associated protein 102 (SAP102), a member of the membrane-associated guanylate kinase protein family. Neuronal SAP102 is expressed during early brain development and is localized to the postsynaptic density of excitatory synapses. It is composed of three amino-terminal PDZ domains, an src homology domain, and a carboxyl-terminal guanylate kinase domain. The PDZ domains interact directly with the NR2 subunits of the NMDA glutamate receptor and with other proteins responsible for NMDA receptor localization, immobilization, and signaling. The mutations identified in this study all introduce premature stop codons within or before the third PDZ domain, and it is likely that this impairs the ability of SAP102 to interact with the NMDA receptor and/or other proteins involved in downstream NMDA receptor signaling pathways. NMDA receptors have been implicated in the induction of certain forms of synaptic plasticity, such as long-term potentiation and long-term depression, and these changes in synaptic efficacy have been proposed as neural mechanisms underlying memory and learning. The disruption of NMDA receptor targeting or signaling, as a result of the loss of SAP102, may lead to altered synaptic plasticity and may explain the intellectual impairment observed in individuals with DLG3 mutations.

Rapid Quantification of Mutant Fitness in Diverse Bacteria by Sequencing Randomly Bar-Coded Transposons

ABSTRACT Transposon mutagenesis with next-generation sequencing (TnSeq) is a powerful approach to annotate gene function in bacteria, but existing protocols for TnSeq require laborious preparation of every sample before sequencing. Thus, the existing protocols are not amenable to the throughput necessary to identify phenotypes and functions for the majority of genes in diverse bacteria. Here, we present a method, random bar code transposon-site sequencing (RB-TnSeq), which increases the throughput of mutant fitness profiling by incorporating random DNA bar codes into Tn5 and mariner transposons and by using bar code sequencing (BarSeq) to assay mutant fitness. RB-TnSeq can be used with any transposon, and TnSeq is performed once per organism instead of once per sample. Each BarSeq assay requires only a simple PCR, and 48 to 96 samples can be sequenced on one lane of an Illumina HiSeq system. We demonstrate the reproducibility and biological significance of RB-TnSeq with Escherichia coli, Phaeobacter inhibens, Pseudomonas stutzeri, Shewanella amazonensis, and Shewanella oneidensis. To demonstrate the increased throughput of RB-TnSeq, we performed 387 successful genome-wide mutant fitness assays representing 130 different bacterium-carbon source combinations and identified 5,196 genes with significant phenotypes across the five bacteria. In P. inhibens, we used our mutant fitness data to identify genes important for the utilization of diverse carbon substrates, including a putative d-mannose isomerase that is required for mannitol catabolism. RB-TnSeq will enable the cost-effective functional annotation of diverse bacteria using mutant fitness profiling. IMPORTANCE A large challenge in microbiology is the functional assessment of the millions of uncharacterized genes identified by genome sequencing. Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach to assign phenotypes and functions to genes. However, the current strategies for TnSeq are too laborious to be applied to hundreds of experimental conditions across multiple bacteria. Here, we describe an approach, random bar code transposon-site sequencing (RB-TnSeq), which greatly simplifies the measurement of gene fitness by using bar code sequencing (BarSeq) to monitor the abundance of mutants. We performed 387 genome-wide fitness assays across five bacteria and identified phenotypes for over 5,000 genes. RB-TnSeq can be applied to diverse bacteria and is a powerful tool to annotate uncharacterized genes using phenotype data. A large challenge in microbiology is the functional assessment of the millions of uncharacterized genes identified by genome sequencing. Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach to assign phenotypes and functions to genes. However, the current strategies for TnSeq are too laborious to be applied to hundreds of experimental conditions across multiple bacteria. Here, we describe an approach, random bar code transposon-site sequencing (RB-TnSeq), which greatly simplifies the measurement of gene fitness by using bar code sequencing (BarSeq) to monitor the abundance of mutants. We performed 387 genome-wide fitness assays across five bacteria and identified phenotypes for over 5,000 genes. RB-TnSeq can be applied to diverse bacteria and is a powerful tool to annotate uncharacterized genes using phenotype data.