Abstract:Abstract Background Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high. Results We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin. Conclusion The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Enhancer Evolution across 20 Mammalian Species

Evolutionary innovations in conserved regulatory elements associate with developmental genes in mammals

Evolutionary innovation in conserved regulatory elements across the mammalian tree of life

Functional characterization of enhancer evolution in the primate lineage

Identification of H3K4me1-associated Proteins at Mammalian Enhancers

Prediction of gene regulatory enhancers across species reveals evolutionarily conserved sequence properties

Gapped-kmer sequence modeling robustly identifies regulatory vocabularies and distal enhancers conserved between evolutionarily distant mammals

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Evolution of transcriptional enhancers and animal diversity

The Evolution of Lineage-Specific Regulatory Activities in the Human Embryonic Limb

Coevolutionary Insights Between Promoters and Transcription Factors in the Plant and Animal Kingdoms

The Evolution of Bony Vertebrate Enhancers at Odds with Their Coding Sequence Landscape

Interactome Maps of Mouse Gene Regulatory Domains Reveal Basic Principles of Transcriptional Regulation.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Dynamic evolution of regulatory element ensembles in primate CD4+ T cells

LINE retrotransposons characterize mammalian tissue-specific and evolutionarily dynamic regulatory regions

Increased Enhancer—Promoter Interactions during Developmental Enhancer Activation in Mammals

Genome-wide analysis of mammalian promoter architecture and evolution

Increased enhancer–promoter interactions during developmental enhancer activation in mammals