Abstract:BackgroundSequence alignments form part of many investigations in molecular biology, including the determination of phylogenetic relationships, the prediction of protein structure and function, and the measurement of evolutionary rates. However, to obtain meaningful results, a significant degree of sequence similarity is required to ensure that the alignments are accurate and the inferences correct. Limitations arise when sequence similarity is low, which is particularly problematic when working with fast-evolving genes, evolutionary distant taxa, genomes with nucleotide biases, and cases of convergent evolution.ResultsA novel approach was conceptualized to address the "low sequence similarity" alignment problem. We developed an alignment algorithm termed FIRE (F unctional I nference using the R ates of E volution), which aligns sequences using the evolutionary rate at codon sites, as measured by the dN/dS ratio, rather than nucleotide or amino acid residues. FIRE was used to test the hypotheses that evolutionary rates can be used to align sequences and that the alignments may be used to infer protein domain function. Using a range of test data, we found that aligning domains based on evolutionary rates was possible even when sequence similarity was very low (for example, antibody variable regions). Furthermore, the alignment has the potential to infer protein domain function, indicating that domains with similar functions are subject to similar evolutionary constraints. These data suggest that an evolutionary rate-based approach to sequence analysis (particularly when combined with structural data) may be used to study cases of convergent evolution or when sequences have very low similarity. However, when aligning homologous gene sets with sequence similarity, FIRE did not perform as well as the best traditional alignment algorithms indicating that the conventional approach of aligning residues as opposed to evolutionary rates remains the method of choice in these cases.ConclusionsFIRE provides proof of concept that it is possible to align sequences and infer domain function by using evolutionary rates rather than residue similarity. This represents a new approach to sequence analysis with a wide range of potential applications in molecular biology.

Sensitive protein alignments at tree-of-life scale using DIAMOND

Fast and sensitive protein alignment using DIAMOND

Sequence alignment using large protein structure alphabets improves sensitivity to remote homologs

A comparative study of available software for high-accuracy homology modeling: from sequence alignments to structural models

DiGAlign: Versatile and Interactive Visualization of Sequence Alignment for Comparative Genomics

A Sequence-Based Evolutionary Distance Method for Phylogenetic Analysis of Highly Divergent Proteins.

Aligning biological sequences by exploiting residue conservation and coevolution

Revealing divergent evolution, identifying circular permutations and detecting active-sites by protein structure comparison

Lapidary: Identifying and reporting amino acid sequences in metagenomes using sequence reads and Diamond

Map of Life: Measuring and Visualizing Species' Relatedness with "Molecular Distance Maps"

SHARK enables sensitive detection of evolutionary homologs and functional analogs in unalignable and disordered sequences

DecentTree: scalable Neighbour-Joining for the genomic era

ML-DSP: Machine Learning with Digital Signal Processing for ultrafast, accurate, and scalable genome classification at all taxonomic levels

DAMBE5: A Comprehensive Software Package for Data Analysis in Molecular Biology and Evolution

Exploring large protein sequence space through homology- and representation-based hierarchical clustering

Phylogenetic Profiles as a Unified Framework for Measuring Protein Structure, Function and Evolution

DIAMOND2GO: A rapid Gene Ontology assignment and enrichment tool for functional genomics

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor

Alignment-free sequence comparison based on next-generation sequencing reads.

Evolutionary rates at codon sites may be used to align sequences and infer protein domain function

Disentangling Cobionts and Contamination in Long-Read Genomic Data using Sequence Composition