The Bioinformatics Core and The Garber Lab

Dr. Manuel Garber discusses the services and educational programs that are offered through the new Bioinformatics Core

A High-Throughput Chromatin Immunoprecipitation Approach Reveals Principles of Dynamic Gene Regulation in Mammals

Understanding the principles governing mammalian gene regulation has been hampered by the difficulty in measuring in vivo binding dynamics of large numbers of transcription factors (TF) to DNA. Here, we develop a high-throughput Chromatin ImmunoPrecipitation (HT-ChIP) method to systematically map protein-DNA interactions. HT-ChIP was applied to define the dynamics of DNA binding by 25 TFs and 4 chromatin marks at 4 time-points following pathogen stimulus of dendritic cells. Analyzing over 180,000 TF-DNA interactions we find that TFs vary substantially in their temporal binding landscapes. This data suggests a model for transcription regulation whereby TF networks are hierarchically organized into cell differentiation factors, factors that bind targets prior to stimulus to prime them for induction, and factors that regulate specific gene programs. Overlaying HT-ChIP data on gene-expression dynamics shows that many TF-DNA interactions are established prior to the stimuli, predominantly at immediate-early genes, and identified specific TF ensembles that coordinately regulate gene-induction

Evolutionary dynamics and tissue specificity of human long noncoding RNAs in six mammals

Long intergenic noncoding RNAs (lincRNAs) play diverse regulatory roles in human development and disease, but little is known about their evolutionary history and constraint. Here, we characterize human lincRNA expression patterns in nine tissues across six mammalian species and multiple individuals. Of the 1898 human lincRNAs expressed in these tissues, we find orthologous transcripts for 80% in chimpanzee, 63% in rhesus, 39% in cow, 38% in mouse, and 35% in rat. Mammalian-expressed lincRNAs show remarkably strong conservation of tissue specificity, suggesting that it is selectively maintained. In contrast, abundant splice-site turnover suggests that exact splice sites are not critical. Relative to evolutionarily young lincRNAs, mammalian-expressed lincRNAs show higher primary sequence conservation in their promoters and exons, increased proximity to protein-coding genes enriched for tissue-specific functions, fewer repeat elements, and more frequent single-exon transcripts. Remarkably, we find that ∼20% of human lincRNAs are not expressed beyond chimpanzee and are undetectable even in rhesus. These hominid-specific lincRNAs are more tissue specific, enriched for testis, and faster evolving within the human lineage.National Institutes of Health (U.S.) (U54-HG004555)National Institutes of Health (U.S.) (R01-HG004037)National Science Foundation (U.S.) (CAREER 0644282

Single Cell Transcriptomics Reveals Dysregulated Cellular and Molecular Networks in a Fragile X Syndrome model [preprint]

Despite advances in understanding the pathophysiology of Fragile X syndrome (FXS), its molecular bases are still poorly understood. Whole brain tissue expression profiles have proved surprisingly uninformative. We applied single cell RNA sequencing to profile a FXS mouse model. We found that FXS results in a highly cell type specific effect and it is strongest among different neuronal types. We detected a downregulation of mRNAs bound by FMRP and this effect is prominent in neurons. Metabolic pathways including translation are significantly upregulated across all cell types with the notable exception of excitatory neurons. These effects point to a potential difference in the activity of mTOR pathways, and together with other dysregulated pathways suggest an excitatory-inhibitory imbalance in the FXS cortex which is exacerbated by astrocytes. Our data demonstrate the cell-type specific complexity of FXS and provide a resource for interrogating the biological basis of this disorder

Computational methods for transcriptome annotation and quantification using RNA-seq

High-throughput RNA sequencing (RNA-seq) promises a comprehensive picture of the transcriptome, allowing for the complete annotation and quantification of all genes and their isoforms across samples. Realizing this promise requires increasingly complex computational methods. These computational challenges fall into three main categories: (i) read mapping, (ii) transcriptome reconstruction and (iii) expression quantification. Here we explain the major conceptual and practical challenges, and the general classes of solutions for each category. Finally, we highlight the interdependence between these categories and discuss the benefits for different biological applications

A High-Throughput Chromatin Immunoprecipitation Approach Reveals Principles of Dynamic Gene Regulation in Mammals

Understanding the principles governing mammalian gene regulation has been hampered by the difficulty in measuring in vivo binding dynamics of large numbers of transcription factors (TF) to DNA. Here, we develop a high-throughput Chromatin ImmunoPrecipitation (HT-ChIP) method to systematically map protein-DNA interactions. HT-ChIP was applied to define the dynamics of DNA binding by 25 TFs and 4 chromatin marks at 4 time-points following pathogen stimulus of dendritic cells. Analyzing over 180,000 TF-DNA interactions we find that TFs vary substantially in their temporal binding landscapes. This data suggests a model for transcription regulation whereby TF networks are hierarchically organized into cell differentiation factors, factors that bind targets prior to stimulus to prime them for induction, and factors that regulate specific gene programs. Overlaying HT-ChIP data on gene-expression dynamics shows that many TF-DNA interactions are established prior to the stimuli, predominantly at immediate-early genes, and identified specific TF ensembles that coordinately regulate gene-induction

DolphinNext: a distributed data processing platform for high throughput genomics

BACKGROUND: The emergence of high throughput technologies that produce vast amounts of genomic data, such as next-generation sequencing (NGS) is transforming biological research. The dramatic increase in the volume of data, the variety and continuous change of data processing tools, algorithms and databases make analysis the main bottleneck for scientific discovery. The processing of high throughput datasets typically involves many different computational programs, each of which performs a specific step in a pipeline. Given the wide range of applications and organizational infrastructures, there is a great need for highly parallel, flexible, portable, and reproducible data processing frameworks. Several platforms currently exist for the design and execution of complex pipelines. Unfortunately, current platforms lack the necessary combination of parallelism, portability, flexibility and/or reproducibility that are required by the current research environment. To address these shortcomings, workflow frameworks that provide a platform to develop and share portable pipelines have recently arisen. We complement these new platforms by providing a graphical user interface to create, maintain, and execute complex pipelines. Such a platform will simplify robust and reproducible workflow creation for non-technical users as well as provide a robust platform to maintain pipelines for large organizations. RESULTS: To simplify development, maintenance, and execution of complex pipelines we created DolphinNext. DolphinNext facilitates building and deployment of complex pipelines using a modular approach implemented in a graphical interface that relies on the powerful Nextflow workflow framework by providing 1. A drag and drop user interface that visualizes pipelines and allows users to create pipelines without familiarity in underlying programming languages. 2. Modules to execute and monitor pipelines in distributed computing environments such as high-performance clusters and/or cloud 3. Reproducible pipelines with version tracking and stand-alone versions that can be run independently. 4. Modular process design with process revisioning support to increase reusability and pipeline development efficiency. 5. Pipeline sharing with GitHub and automated testing 6. Extensive reports with R-markdown and shiny support for interactive data visualization and analysis. CONCLUSION: DolphinNext is a flexible, intuitive, web-based data processing and analysis platform that enables creating, deploying, sharing, and executing complex Nextflow pipelines with extensive revisioning and interactive reporting to enhance reproducible results

Genome-wide assessment of post-transcriptional control in the fly brain

Post-transcriptional control of gene expression has central importance during development and adulthood and in physiology in general. However, little is known about the extent of post-transcriptional control of gene expression in the brain. Most post-transcriptional regulatory effectors (e.g., miRNAs) destabilize target mRNAs by shortening their polyA tails. Hence, the fraction of a given mRNA that it is fully polyadenylated should correlate with its stability and serves as a good measure of post-transcriptional control. Here, we compared RNA-seq datasets from fly brains that were generated either from total (rRNA-depleted) or polyA-selected RNA. By doing this comparison we were able to compute a coefficient that measures the extent of post-transcriptional control for each brain-expressed mRNA. In agreement with current knowledge, we found that mRNAs encoding ribosomal proteins, metabolic enzymes, and housekeeping genes are among the transcripts with least post-transcriptional control, whereas mRNAs that are known to be highly unstable, like circadian mRNAs and mRNAs expressing synaptic proteins and proteins with neuronal functions, are under strong post-transcriptional control. Surprisingly, the latter group included many specific groups of genes relevant to brain function and behavior. In order to determine the importance of miRNAs in this regulation, we profiled miRNAs from fly brains using oligonucleotide microarrays. Surprisingly, we did not find a strong correlation between the expression levels of miRNAs in the brain and the stability of their target mRNAs; however, genes identified as highly regulated post-transcriptionally were strongly enriched for miRNA targets. This demonstrates a central role of miRNAs for modulating the levels and turnover of brain-specific mRNAs in the fly

A comparative genomics multitool for scientific discovery and conservation

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity

The Tec kinase ITK differentially optimizes NFAT, NF-κB, and MAPK signaling during early T cell activation to regulate graded gene induction [preprint]

The strength of peptide:MHC interactions with the T cell receptor (TCR) is correlated with the time to first cell division, the relative scale of the effector cell response, and the graded expression of activation-associated proteins like IRF4. To regulate T cell activation programming, the TCR and the TCR proximal kinase ITK simultaneously trigger many biochemically separate TCR signaling cascades. T cells lacking ITK exhibit selective impairments in effector T cell responses after activation, but under the strongest signaling conditions ITK activity is dispensable. To gain insight into whether TCR signal strength and ITK activity tune observed graded gene expression through unequal activation of disparate signaling pathways, we examined Erk1/2 activation and NFAT, NF-κB translocation in naive OT-I CD8+ cell nuclei. We observed consistent digital activation of NFAT1 and Erk-MAPK, but NF-κB displayed dynamic, graded activation in response to variation in TCR signal strength and was tunable by treatment with an ITK inhibitor. Inhibitor-treated cells showed dampened induction of AP-1 factors Fos and Fosb, NF-κB response gene transcripts, and survival factor Il2 transcripts. ATAC-seq analysis also revealed genomic regions most sensitive to ITK inhibition were enriched for NF-κB and AP-1 motifs. Specific inhibition of NF-κB during peptide stimulation tuned expression of early gene products like c-Fos. Together, these data indicate a key role for ITK in orchestrating optimal activation of separate TCR downstream pathways, specifically aiding NF-κB activation. More broadly, we revealed a mechanism by which variation in TCR signal strength can produce patterns of graded gene expression in activated T cells