Inference of gene regulatory interactions from deep sequencing data

This thesis describes our work on the inference of gene regulatory interactions from deep sequencing data. In the first part we start with an introduction to the problem of inferring transcription regulatory interactions (Chapter 2). Then we describe detailed methods of promoterome construction in human and mouse genomes from deep CAGE data (Chapter 3), and show how promoters together with gene expression data can be used to infer transcription regulatory interactions (Chapter 4). Finally, we show an application of this methodology to a human macrophage lineage undergoing differentiation accompanied by a detailed experimental validation of the predicted network structure (Chapter 5). Work presented in Chapter 5 comes from the FANTOM4 project. In the second part we focus on the regulatory functions of two small nucleolar RNAs (snoRNAs). In Chapter 8 we describe an atypical function of snoRNAs - regulation of the mRNA alternative splicing process by a particular class of mouse snoRNAs (MBII-52 variants). These are of great interest, as the locus encoding MBII-52 is linked to Prader-Willi Syndrome. Methods include in silico RNA hybridization screens and experimental confirmation of the predictions. In Chapter 9, we report a discovery of the first virus-encoded snoRNA in Epstein-Barr Virus (EBV). Again, we show that the function of this snoRNA (v-snoRNA1) is atypical: it is processed into small, 24 nt long fragments that can function as microRNAs

SwissRegulon, a database of genome-wide annotations of regulatory sites: recent updates

Identification of genomic regulatory elements is essential for understanding the dynamics of cellular processes. This task has been substantially facilitated by the availability of genome sequences for many species and high-throughput data of transcripts and transcription factor (TF) binding. However, rigorous computational methods are necessary to derive accurate genome-wide annotations of regulatory sites from such data. SwissRegulon (http://swissregulon.unibas.ch) is a database containing genome-wide annotations of regulatory motifs, promoters and TF binding sites (TFBSs) in promoter regions across model organisms. Its binding site predictions were obtained with rigorous Bayesian probabilistic methods that operate on orthologous regions from related genomes, and use explicit evolutionary models to assess the evidence of purifying selection on each site. New in the current version of SwissRegulon is a curated collection of 190 mammalian regulatory motifs associated with ∼340 TFs, and TFBS annotations across a curated set of ∼35 000 promoters in both human and mouse. Predictions of TFBSs for Saccharomyces cerevisiae have also been significantly extended and now cover 158 of yeast's ∼180 TFs. All data are accessible through both an easily navigable genome browser with search functions, and as flat files that can be downloaded for further analysi

The Corepressor NCoR1 Antagonizes PGC-1α and Estrogen-Related Receptor α in the Regulation of Skeletal Muscle Function and Oxidative Metabolism

Skeletal muscle exhibits a high plasticity and accordingly can quickly adapt to different physiological and pathological stimuli by changing its phenotype largely through diverse epigenetic mechanisms. The nuclear receptor corepressor 1 (NCoR1) has the ability to mediate gene repression; however, its role in regulating biological programs in skeletal muscle is still poorly understood. We therefore studied the mechanistic and functional aspects of NCoR1 function in this tissue. NCoR1 muscle-specific knockout mice exhibited a 7.2% higher peak oxygen consumption (VO(2peak)), a 11% reduction in maximal isometric force, and increased ex vivo fatigue resistance during maximal stimulation. Interestingly, global gene expression analysis revealed a high overlap between the effects of NCoR1 deletion and peroxisome proliferator-activated receptor gamma (PPARγ) coactivator 1α (PGC-1α) overexpression on oxidative metabolism in muscle. Importantly, PPARβ/δ and estrogen-related receptor α (ERRα) were identified as common targets of NCoR1 and PGC-1α with opposing effects on the transcriptional activity of these nuclear receptors. In fact, the repressive effect of NCoR1 on oxidative phosphorylation gene expression specifically antagonizes PGC-1α-mediated coactivation of ERRα. We therefore delineated the molecular mechanism by which a transcriptional network controlled by corepressor and coactivator proteins determines the metabolic properties of skeletal muscle, thus representing a potential therapeutic target for metabolic diseases

Embryonic stem cell-specific microRNAs contribute to pluripotency by inhibiting regulators of multiple differentiation pathways

The findings that microRNAs (miRNAs) are essential for early development in many species and that embryonic miRNAs can reprogram somatic cells into induced pluripotent stem cells suggest that these miRNAs act directly on transcriptional and chromatin regulators of pluripotency. To elucidate the transcription regulatory networks immediately downstream of embryonic miRNAs, we extended the motif activity response analysis approach that infers the regulatory impact of both transcription factors (TFs) and miRNAs from genome-wide expression states. Applying this approach to multiple experimental data sets generated from mouse embryonic stem cells (ESCs) that did or did not express miRNAs of the ESC-specific miR-290-295 cluster, we identified multiple TFs that are direct miRNA targets, some of which are known to be active during cell differentiation. Our results provide new insights into the transcription regulatory network downstream of ESC-specific miRNAs, indicating that these miRNAs act on cell cycle and chromatin regulators at several levels and downregulate TFs that are involved in the innate immune respons

Methods for analyzing deep sequencing expression data: constructing the human and mouse promoterome with deepCAGE data

A set of methods is presented for normalization, quantification of noise and co-expression analysis for gene expression studies using deep sequencing

The snoRNA MBII-52 (SNORD 115) is processed into smaller RNAs and regulates alternative splicing

The loss of HBII-52 and related C/D box small nucleolar RNA (snoRNA) expression units have been implicated as a cause for the Prader-Willi syndrome (PWS). We recently found that the C/D box snoRNA HBII-52 changes the alternative splicing of the serotonin receptor 2C pre-mRNA, which is different from the traditional C/D box snoRNA function in non-mRNA methylation. Using bioinformatic predictions and experimental verification, we identified five pre-mRNAs (DPM2, TAF1, RALGPS1, PBRM1 and CRHR1) containing alternative exons that are regulated by MBII-52, the mouse homolog of HBII-52. Analysis of a single member of the MBII-52 cluster of snoRNAs by RNase protection and northern blot analysis shows that the MBII-52 expressing unit generates shorter RNAs that originate from the full-length MBII-52 snoRNA through additional processing steps. These novel RNAs associate with hnRNPs and not with proteins associated with canonical C/D box snoRNAs. Our data indicate that not a traditional C/D box snoRNA MBII-52, but a processed version lacking the snoRNA stem is the predominant MBII-52 RNA missing in PWS. This processed snoRNA functions in alternative splice-site selection. Its substitution could be a therapeutic principle for PW

Identification of downstream effectors of retinoic acid specifying the zebrafish pancreas by integrative genomics.

Retinoic acid (RA) is a key signal for the specification of the pancreas. Still, the gene regulatory cascade triggered by RA in the endoderm remains poorly characterized. In this study, we investigated this regulatory network in zebrafish by combining RNA-seq, RAR ChIP-seq and ATAC-seq assays. By analysing the effect of RA and of the RA receptor (RAR) inverse-agonist BMS493 on the transcriptome and on the chromatin accessibility of endodermal cells, we identified a large set of genes and regulatory regions regulated by RA signalling. RAR ChIP-seq further defined the direct RAR target genes in zebrafish, including hox genes as well as several pancreatic regulators like mnx1, insm1b, hnf1ba and gata6. Comparison of zebrafish and murine RAR ChIP-seq data highlighted the conserved direct target genes and revealed that some RAR sites are under strong evolutionary constraints. Among them, a novel highly conserved RAR-induced enhancer was identified downstream of the HoxB locus and driving expression in the nervous system and in the gut in a RA-dependent manner. Finally, ATAC-seq data unveiled the role of the RAR-direct targets Hnf1ba and Gata6 in opening chromatin at many regulatory loci upon RA treatment

FANTOM4 EdgeExpressDB: an integrated database of promoters, genes, microRNAs, expression dynamics and regulatory interactions

EdgeExpressDB is a novel database and set of interfaces for interpreting biological networks and comparing large high-throughput expression datasets

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations