Evolution of lineage-specific functions in ancient cis-regulatory modules

Morphological evolution is driven both by coding sequence variation and by changes in regulatory sequences. However, how cis-regulatory modules (CRMs) evolve to generate entirely novel expression domains is largely unknown. Here, we reconstruct the evolutionary history of a lens enhancer located within a CRM that not only predates the lens, a vertebrate innovation, but bilaterian animals in general. Alignments of orthologous sequences from different deuterostomes sub-divide the CRM into a deeply conserved core and a more divergent flanking region. We demonstrate that all deuterostome flanking regions, including invertebrate sequences, activate gene expression in the zebrafish lens through the same ancient cluster of activator sites. However, levels of gene expression vary between species due to the presence of repressor motifs in flanking region and core. These repressor motifs are responsible for the relatively weak enhancer activity of tetrapod flanking regions. Ray-finned fish, however, have gained two additional lineage-specific activator motifs which in combination with the ancient cluster of activators and the core constitute a potent lens enhancer. The exploitation and modification of existing regulatory potential in flanking regions but not in the highly conserved core might represent a more general model for the emergence of novel regulatory functions in complex CRMs.     

Translation Enhancing ACA Motifs and Their Silencing by a Bacterial Small Regulatory RNA

GcvB is an archetypal multi-target small RNA regulator of genes involved in amino acid uptake or metabolism in enteric bacteria. Included in the GcvB regulon is the yifK locus, encoding a conserved putative amino acid transporter. GcvB inhibits yifK mRNA translation by pairing with a sequence immediately upstream from the Shine-Dalgarno motif. Surprisingly, we found that some target sequence mutations that disrupt pairing, and thus were expected to relieve repression, actually lower yifK expression and cause it not to respond to GcvB variants carrying the corresponding compensatory changes. Work prompted by these observations revealed that the GcvB target sequence in yifK mRNA includes elements that stimulate translation initiation. Replacing each base of an ACA trinucleotide near the center of the target sequence, by any other base, caused yifK expression to decrease. Effects were additive, with some triple replacements causing up to a 90% reduction. The enhancer activity did not require the ACA motif to be strictly positioned relative to the Shine-Dalgarno sequence, nor did it depend on a particular spacing between the latter and the initiating AUG. The dppA mRNA, another GcvB target, contains four ACA motifs at the target site. Quite strikingly, replacement of all four ACAs by random trinucleotide sequences yielded variants showing over 100-fold reduction in expression, virtually inactivating the gene. Altogether, these data identify the ACA motif as a translation-enhancing module and show that GcvB''s ability to antagonize the enhancer function in target mRNAs is quintessential to the regulatory effectiveness of this sRNA.     

Large‐scale analysis of Drosophila core promoter function using synthetic promoters

The core promoter plays a central role in setting metazoan gene expression levels, but how exactly it “computes” expression remains poorly understood. To dissect its function, we carried out a comprehensive structure–function analysis in Drosophila. First, we performed a genome‐wide bioinformatic analysis, providing an improved picture of the sequence motifs architecture. We then measured synthetic promoters’ activities of ~3,000 mutational variants with and without an external stimulus (hormonal activation), at large scale and with high accuracy using robotics and a dual luciferase reporter assay. We observed a strong impact on activity of the different types of mutations, including knockout of individual sequence motifs and motif combinations, variations of motif strength, nucleosome positioning, and flanking sequences. A linear combination of the individual motif features largely accounts for the combinatorial effects on core promoter activity. These findings shed new light on the quantitative assessment of gene expression in metazoans.     

MOPAT: a graph-based method to predict recurrent cis-regulatory modules from known motifs

The identification of cis-regulatory modules (CRMs) can greatly advance our understanding of eukaryotic regulatory mechanism. Current methods to predict CRMs from known motifs either depend on multiple alignments or can only deal with a small number of known motifs provided by users. These methods are problematic when binding sites are not well aligned in multiple alignments or when the number of input known motifs is large. We thus developed a new CRM identification method MOPAT (motif pair tree), which identifies CRMs through the identification of motif modules, groups of motifs co-occurring in multiple CRMs. It can identify 'orthologous' CRMs without multiple alignments. It can also find CRMs given a large number of known motifs. We have applied this method to mouse developmental genes, and have evaluated the predicted CRMs and motif modules by microarray expression data and known interacting motif pairs. We show that the expression profiles of the genes containing CRMs of the same motif module correlate significantly better than those of a random set of genes do. We also show that the known interacting motif pairs are significantly included in our predictions. Compared with several current methods, our method shows better performance in identifying meaningful CRMs.     

The differentiation-dependent inhibition of SV40 early promoter/enhancer activity in 3T3-L1 cells is mediated through promoter and not enhancer sequences

Using a plasmid bearing chloramphenicol acetyltransferase (CAT) gene controlled by Simian virus 40 (SV40) early promoter/enhancer complex (pA0cat), we analyzed functional enhancer motifs in 3T3-L1 fibroblast and adipocyte cells. Deletion mutant series of pA0 at the enhancer complex showed that gene expression both in fibroblast and adipocyte cells was dependent on a similar set of enhancer motifs. When pA0 was introduced into 3T3-L1 fibroblasts and the cells were induced to differentiate into adipocytes, CAT activity expressed in fibroblasts was suppressed. Experiments with the deletion mutants at the enhancer complex showed that the suppression was not related to any enhancer motif, and CAT activity was observed with a plasmid having only the promoter sequence. When pA0cat was co-transfected with excess of promoter sequence, the suppression in adipocytes was counteracted. This suggested that negativetrans-acting factors of the promoter sequence were responsible for the suppression in adipocytes.Abbreviations CAT chloramphenicol acetyltransferase - CAT the gene encoding CAT - SV40 Simian virus 40 - Asc-P ascorbic acid phosphate     

TELS: A Novel Computational Framework for Identifying Motif Signatures of Transcribed Enhancers

In mammalian cells, transcribed enhancers (TrEns) play important roles in the initiation of gene expression and maintenance of gene expression levels in a spatiotemporal manner. One of the most challenging questions is how the genomic characteristics of enhancers relate to enhancer activities. To date, only a limited number of enhancer sequence characteristics have been investigated, leaving space for exploring the enhancers’ DNA code in a more systematic way. To address this problem, we developed a novel computational framework, Transcribed Enhancer Landscape Search (TELS), aimed at identifying predictive cell type/tissue-specific motif signatures of TrEns. As a case study, we used TELS to compile a comprehensive catalog of motif signatures for all known TrEns identified by the FANTOM5 consortium across 112 human primary cells and tissues. Our results confirm that combinations of different short motifs characterize in an optimized manner cell type/tissue-specific TrEns. Our study is the first to report combinations of motifs that maximize classification performance of TrEns exclusively transcribed in one cell type/tissue from TrEns exclusively transcribed in different cell types/tissues. Moreover, we also report 31 motif signatures predictive of enhancers’ broad activity. TELS codes and material are publicly available at http://www.cbrc.kaust.edu.sa/TELS.     

Single nucleotide polymorphisms in the upstream regulatory region alter the expression of myostatin

The expression of the gene encoding myostatin (MSTN), the product of which is a negative regulator of skeletal muscle growth and development in mammals, is regulated by many cis-regulatory elements, including enhancer box (E-box) motifs. While E-box motif mutants of MSTN exhibit altered expression of myostatin in many animal models, the phenotypes of these mutations in chicken are not investigated. In this study, we cloned and sequenced the full encoded DNA sequence of MSTN gene and its upstream promoter region in Wenshang Luhua chicken breed. After analysis of the sequence, 13 E-box motifs were identified in the MSTN promoter region, which were denoted by E1 to E13 according to their positions in the region. Although many single nucleotide polymorphisms (SNPs) were revealed in the MSTN promoter region, only two SNPs were in the E-boxes, i.e., the first nucleotide of the E3 and the fifth nucleotide of E4. The effects of these two polymorphisms on the expression of MSTN gene were explored both with MSTN-GFP reporter constructs in vitro and real-time PCR in vivo. The results suggested that the E-boxes in the chicken MSTN promoter region are involved in the regulation of myostatin expression and the polymorphisms in E3 and E4 altered the expression of myostatin.     

A systematic approach to identify functional motifs within vertebrate developmental enhancers

Uncovering the cis-regulatory logic of developmental enhancers is critical to understanding the role of non-coding DNA in development. However, it is cumbersome to identify functional motifs within enhancers, and thus few vertebrate enhancers have their core functional motifs revealed. Here we report a combined experimental and computational approach for discovering regulatory motifs in developmental enhancers. Making use of the zebrafish gene expression database, we computationally identified conserved non-coding elements (CNEs) likely to have a desired tissue-specificity based on the expression of nearby genes. Through a high throughput and robust enhancer assay, we tested the activity of ∼ 100 such CNEs and efficiently uncovered developmental enhancers with desired spatial and temporal expression patterns in the zebrafish brain. Application of de novo motif prediction algorithms on a group of forebrain enhancers identified five top-ranked motifs, all of which were experimentally validated as critical for forebrain enhancer activity. These results demonstrate a systematic approach to discover important regulatory motifs in vertebrate developmental enhancers. Moreover, this dataset provides a useful resource for further dissection of vertebrate brain development and function.