Abstract:Cross-species comparison has emerged as a powerful paradigm for predicting cis-regulatory modules (CRMs) and understanding their evolution. The comparison requires reliable sequence alignment, which remains a challenging task for less conserved noncoding sequences. Furthermore, the existing models of DNA sequence evolution generally do not explicitly treat the special properties of CRM sequences. To address these limitations, we propose a model of CRM evolution that captures different modes of evolution of functional transcription factor binding sites (TFBSs) and the background sequences. A particularly novel aspect of our work is a probabilistic model of gains and losses of TFBSs, a process being recognized as an important part of regulatory sequence evolution. We present a computational framework that uses this model to solve the problems of CRM alignment and prediction. Our alignment method is similar to existing methods of statistical alignment but uses the conserved binding sites to improve alignment. Our CRM prediction method deals with the inherent uncertainties of binding site annotations and sequence alignment in a probabilistic framework. In simulated as well as real data, we demonstrate that our program is able to improve both alignment and prediction of CRM sequences over several state-of-the-art methods. Finally, we used alignments produced by our program to study binding site conservation in genome-wide binding data of key transcription factors in the Drosophila blastoderm, with two intriguing results: (i) the factor-bound sequences are under strong evolutionary constraints even if their neighboring genes are not expressed in the blastoderm and (ii) binding sites in distal bound sequences (relative to transcription start sites) tend to be more conserved than those in proximal regions. Our approach is implemented as software, EMMA (Evolutionary Model-based cis-regulatory Module Analysis), ready to be applied in a broad biological context.

DCAlign v1.0: Aligning biological sequences using co-evolution models and informed priors

Aligning biological sequences by exploiting residue conservation and coevolution

Small Coupling Expansion for Multiple Sequence Alignment

COATi: statistical pairwise alignment of protein-coding sequences

Enhancing coevolutionary signals in protein–protein interaction prediction through clade-wise alignment integration

SaAlign: Multiple DNA/RNA Sequence Alignment and Phylogenetic Tree Construction Tool for Ultra-Large Datasets and Ultra-Long Sequences Based on Suffix Array

Accurate Anchoring Alignment of Divergent Sequences

DNA^+Pro^: an Improved Progressive Multiple Sequence Alignment Algorithm for Evolutionary Analysis Using Combined DNA-Protein Sequences

Pairwise sequence alignment at arbitrarily large evolutionary distance

Multiple sequence alignment by ant colony optimization and divide-and-conquer

ReAligner: a program for refining DNA sequence multi-alignments

DiagAF: A More Accurate and Efficient Pre-Alignment Filter for Sequence Alignment.

CAUSA 2.0: Accurate and Consistent Evolutionary Analysis of Proteins Using Codon and Amino Acid Unified Sequence Alignments

GridWiseTM: The Benefits of a Transformed Energy System

Parametric Alignment of Drosophila Genomes

ProbAlign: a re-alignment method for long sequencing reads

Accelerating phylogeny-aware alignment with indel evolution using short time Fourier transform

A software tool for finding locally optimal alignments in protein and nucleic acid sequences

Alignment and Prediction of -Regulatory Modules Based on a Probabilistic Model of Evolution.

A Short Reads Alignment Algorithm Oriented To Massive Data