Modeling the Evolutionary Architectures of Transcribed Human Enhancer Sequences Reveals Distinct Origins,Functions, and Associations with Human Trait Variation

Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits.To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age (“simple” evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages (“complex” evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers.Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.     

Structural and functional conservation of the Drosophila doublesex splicing enhancer repeat elements.

We have compared the RNA sequences and secondary structures of the Drosophila melanogaster and Drosophila virilis doublesex (dsx) splicing enhancers. The sequences of the two splicing enhancers are highly divergent except for the presence of nearly identical 13-nt repeat elements (six in D. melanogaster and four in D. virilis) and a stretch of nucleotides at the 5' and 3' ends of the enhancers. In vitro RNA structure probing of the two enhancers revealed that the 13-nt repeats are predominantly single-stranded. Thus, both the primary sequences and single-stranded nature of the repeats are conserved between the two species. The significance of the primary sequence conservation was demonstrated by showing that the two enhancers are functionally interchangeable in Tra-/Tra2-dependent in vitro splicing. In addition, inhibition of splicing enhancer activity by antisense oligonucleotides complementary to the repeats demonstrated the importance of the conserved single-stranded structure of the repeats. In vitro binding studies revealed that Tra2 interacts with each of the D. melanogaster repeat elements, except for repeat 2, with affinities that are indistinguishable, whereas Tra binds nonspecifically to the enhancer. Taken together, these observations indicate that the organization of sequences within the dsx splicing enhancers of D. melanogaster and D. virilis results in a structure in which each of the repeat elements is single-stranded and therefore accessible for specific recognition by the RNA-binding domain of Tra2.     

Sequence-Specific and/or Stereospecific Constraints of the U3 Enhancer Elements of MCF 247-W Are Important for Pathogenicity

The oncogenic potential of many nonacute retroviruses is dependent on the duplication of the enhancer sequences present in the unique 3′ (U3) region of the long terminal repeat (LTR). In a molecular clone (MCF 247-W) of the murine leukemia virus MCF 247, a leukemogenic mink cell focus-inducing (MCF) virus, the U3 enhancer sequences are tandemly repeated in the LTR. We mutated the enhancer region of MCF 247-W to test the hypothesis that the duplicated enhancer sequences of this virus have a sequence-specific and/or a stereospecific role in enhancer function required for transformation. In one virus, we inserted 14 nucleotide bp into the novel sequence generated at the junction of the two enhancers to generate an MCF virus with an interrupted enhancer region. In the second virus, only one copy of the enhancer sequences was present. This second virus also lacked the junction sequence present between the two enhancers of MCF 247-W. Both viruses were less leukemogenic and had a longer mean latency period than MCF 247-W. These data indicate that the sequence generated at the junction of the two enhancers and/or the stereospecific arrangement of the two enhancer elements are required for the full oncogenic potential of MCF 247-W. We analyzed proviral LTRs within the c-myc locus in tumor DNAs from mice injected with the MCF virus with the interrupted enhancer region. Some of the proviral LTRs integrated upstream of c-myc contain enhancer regions that are larger than those of the injected virus. These results are consistent with the suggestion that the virus with an interrupted enhancer changes in vivo to perform its role in the transformation of T cells.     

Sequence elements essential for function of the Xenopus laevis ribosomal DNA enhancers.

The intergenic spacer region of the Xenopus laevis ribosomal DNA contains multiple elements which are either 60 or 81 base pairs long. Clusters of these elements have previously been shown to act as position- and distance-independent enhancers on an RNA polymerase I promoter when located in cis. By a combination of deletion and linker scanner mutagenesis we show that the sequences essential for enhancer function are located within a 56-base-pair region that is present in both the 60- and 81-base-pair repeats. Within the 56-base-pair region one linker scanner mutation was found to be relatively neutral, suggesting that each enhancer element may be composed of two smaller domains. Each 56-base-pair region appears to be an independent enhancer with multiple enhancers being additive in effect. We review the current evidence concerning the mechanism of action of these enhancers.     

A Short Sequence within Two Purine-Rich Enhancers Determines 5′ Splice Site Specificity

Purine-rich enhancers are exon sequences that promote inclusion of alternative exons, usually via activation of weak upstream 3′ splice sites. A recently described purine-rich enhancer from the caldesmon gene has an additional activity by which it directs selection of competing 5′ splice sites within an alternative exon. In this study, we have compared the caldesmon enhancer with another purine-rich enhancer from the chicken cardiac troponin T (cTNT) gene for the ability to regulate flanking splice sites. Although similar in sequence and length, the two enhancers demonstrated strikingly different specificities towards 5′ splice site choice when placed between competing 5′ splice sites in an internal exon. The 32-nucleotide caldesmon enhancer caused effective usage of the exon-internal 5′ splice site, whereas the 30-nucleotide cTNT enhancer caused effective usage of the exon-terminal 5′ splice site. Both enhancer-mediated splicing pathways represented modulation of the default pathway in which both 5′ splice sites were utilized. Each enhancer is multipartite, consisting of two purine-rich sequences of a simple (GAR)_n repeat interdigitated with two enhancer-specific sequences. The entire enhancer was necessary for maximal splice site selectivity; however, a 5- to 7-nucleotide region from the 3′ end of each enhancer dictated splice site selectivity. Mutations that interchanged this short region of the two enhancers switched specificity. The portion of the cTNT enhancer determinative for 5′ splice site selectivity was different than that shown to be maximally important for activation of a 3′ splice site, suggesting that enhancer environment can have a major impact on activity. These results are the first indication that individual purine-rich enhancers can differentiate between flanking splice sites. Furthermore, localization of the specificity of splice site choice to a short region within both enhancers indicates that subtle differences in enhancer sequence can have profound effects on the splicing pathway.     

Multiple Distinct Splicing Enhancers in the Protein-Coding Sequences of a Constitutively Spliced Pre-mRNA

We have identified multiple distinct splicing enhancer elements within protein-coding sequences of the constitutively spliced human β-globin pre-mRNA. Each of these highly conserved sequences is sufficient to activate the splicing of a heterologous enhancer-dependent pre-mRNA. One of these enhancers is activated by and binds to the SR protein SC35, whereas at least two others are activated by the SR protein SF2/ASF. A single base mutation within another enhancer element inactivates the enhancer but does not change the encoded amino acid. Thus, overlapping protein coding and RNA recognition elements may be coselected during evolution. These studies provide the first direct evidence that SR protein-specific splicing enhancers are located within the coding regions of constitutively spliced pre-mRNAs. We propose that these enhancers function as multisite splicing enhancers to specify 3′ splice-site selection.     

Prediction of transposable element derived enhancers using chromatin modification profiles

Experimentally characterized enhancer regions have previously been shown to display specific patterns of enrichment for several different histone modifications. We modelled these enhancer chromatin profiles in the human genome and used them to guide the search for novel enhancers derived from transposable element (TE) sequences. To do this, a computational approach was taken to analyze the genome-wide histone modification landscape characterized by the ENCODE project in two human hematopoietic cell types, GM12878 and K562. We predicted the locations of 2,107 and 1,448 TE-derived enhancers in the GM12878 and K562 cell lines respectively. A vast majority of these putative enhancers are unique to each cell line; only 3.5% of the TE-derived enhancers are shared between the two. We evaluated the functional effect of TE-derived enhancers by associating them with the cell-type specific expression of nearby genes, and found that the number of TE-derived enhancers is strongly positively correlated with the expression of nearby genes in each cell line. Furthermore, genes that are differentially expressed between the two cell lines also possess a divergent number of TE-derived enhancers in their vicinity. As such, genes that are up-regulated in the GM12878 cell line and down-regulated in K562 have significantly more TE-derived enhancers in their vicinity in the GM12878 cell line and vice versa. These data indicate that human TE-derived sequences are likely to be involved in regulating cell-type specific gene expression on a broad scale and suggest that the enhancer activity of TE-derived sequences is mediated by epigenetic regulatory mechanisms.     

Specific interaction of cellular factors with the B enhancer of polyoma virus.

Specific interactions between proteins from mouse 3T6 cells and the enhancer sequence of polyoma virus were detected using the method of band shifting on polyacrylamide gels. Proteins eluted from 3T6 nuclei using a buffer containing 0.55 M NaCl, formed a stable complex with the B enhancer of polyoma virus. At least two different factors are involved in this interaction. The contact sites which were mapped on the DNA sequence using DNase I footprinting correspond to a GC-rich palindrome surrounded by two sequences homologous respectively to the immunoglobulin and to the immunoglobulin and SV40 enhancers. Moreover Bal31 deletion analysis confirmed that similar sequences are required for the formation of the complex. In spite of a common function and partial sequence homology among some enhancers, neither the polyoma A enhancer, the mouse immunoglobulin heavy chain gene enhancer, nor the origin-promoter-enhancer region of SV40 efficiently competed with the polyoma B enhancer for the binding of these molecules.     

Novel conserved elements upstream of the H19 gene are transcribed and act as mesodermal enhancers

The reciprocally imprinted H19 and Igf2 genes form a co-ordinately regulated 130 kb unit in the mouse controlled by widely dispersed enhancers, epigenetically modified silencers and an imprinting control region (ICR). Comparative human and mouse genomic sequencing between H19 and Igf2 revealed two novel regions of strong homology upstream of the ICR termed H19 upstream conserved regions (HUCs). Mouse HUC1 and HUC2 act as potent enhancers capable of driving expression of an H19 reporter gene in a range of mesodermal tissues. Intriguingly, the HUC sequences are also transcribed bi-allelically in mouse and human, but their expression pattern in neural and endodermal tissues in day 13.5 embryos is distinct from their enhancer function. The location of the HUC mesodermal enhancers upstream of the ICR and H19, and their capacity for interaction with both H19 and Igf2 requires critical re-evaluation of the cis-regulation of imprinted gene expression of H19 and Igf2 in a range of mesodermal tissues. We propose that these novel sequences interact with the ICR at H19 and the epigenetically regulated silencer at differentially methylated region 1 (DMR1) of Igf2.     

Strong RNA splicing enhancers identified by a modified method of cycled selection interact with SR protein.

A modified method of cycled selection was used to characterize splicing enhancers for exon inclusion from a pool of beta-globin-based three exon/two intron pre-mRNAs with a variable number of random nucleotides incorporated in the internal exon. The pre-mRNAs generated by this method contained random sequences ranging from 0 to 18 nucleotides in length. This method was used to isolate particular splicing enhancer motifs from a previously enriched pool of extremely diverse enhancers. After four cycles of selection for mRNA containing the internal exon, a distinct enhancer motif (GACGAC...CAGCAG) was highly enriched. This motif served as strong splicing enhancers in a heterogeneous exon. We have shown here that the selected enhancer motif promotes exon inclusion through specific interaction with SRp30. We have also shown that although present in many of our selected splicing enhancers conforming to this motif, a typical purine-rich enhancer sequence is dispensable for either enhancer activity or binding with SRp30.