Genome-wide functional element detection using pairwise statistical alignment outperforms multiple genome...

Area of Science:

• Genomics
• Bioinformatics
• Computational Biology

Background:

• Comparative genomics utilizes sequence analysis to identify functional elements in silico.
• Tools like phastCons analyze multiple vertebrate genomes for conserved elements but overlook insertion/deletion events.
• Existing methods often rely on precomputed multiple sequence alignments, limiting their scope.

Purpose of the Study:

• To develop a novel method, GRAPeFoot, for simultaneous pairwise genome alignment and functional element annotation.
• To integrate the modeling of substitutions, insertions, and deletions (indels) for enhanced accuracy.
• To assess GRAPeFoot's performance against established tools using human and mouse genomes.

Main Methods:

• Developed GRAPeFoot for simultaneous alignment of two full genomes.
• Annotated conserved regions based on reduced rates of insertion, deletion, and substitution mutations.
• Compared GRAPeFoot predictions with phastCons results from the UCSC genome browser.

Main Results:

• GRAPeFoot identified constrained elements with rates comparable to phastCons, despite analyzing only two genomes.
• The method effectively integrates modeling of substitutions, indels, and purifying selection.
• Pairwise analysis demonstrated sensitivity similar to heuristic multi-genome analyses.

Conclusions:

• GRAPeFoot offers a sensitive approach for identifying functional genomic elements using pairwise comparisons.
• Integrated modeling of mutation types enhances the detection of conserved regions.
• The software and predictions are publicly available, facilitating further research.