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.

Genome-wide comparative analysis reveals human-mouse regulatory landscape and evolution

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

Mouse regulatory DNA landscapes reveal global principles of cis-regulatory evolution

A Map of the Cis-Regulatory Sequences in the Mouse Genome

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

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

Comparative Epigenomic Annotation of Regulatory DNA

Eukaryotic Regulatory Element Conservation Analysis and Identification Using Comparative Genomics

Reorganization of 3D Genome Structure May Contribute to Gene Regulatory Evolution in Primates

A comparative genome analysis of gene expression reveals different regulatory mechanisms between mouse and human embryo pre-implantation development

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

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

A census of human transcription factors: function, expression and evolution

Topologically associating domains are disrupted by evolutionary genome rearrangements forming species-specific enhancer connections in mice and humans

The Transcription Factor Binding Landscape of Mouse Development

Comparative studies of gene expression and the evolution of gene regulation

Stem cell transcriptional profiles from mouse subspecies reveal cis -regulatory evolution at translation genes

Stem cell transcriptional profiles from mouse subspecies reveal cis-regulatory evolution at translation genes

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

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

Characterization And Comparison Of The Tissue-Related Modules In Human And Mouse