An alignment-free model for comparison of regulatory sequences

Area of Science:

• Computational biology
• Genomics
• Bioinformatics

Background:

• Regulatory DNA sequences can maintain function despite significant sequence divergence over evolutionary time.
• Enhancers driving similar expression patterns pose challenges for traditional sequence alignment methods in silico.
• Alignment-free approaches are necessary for detecting functional conservation in divergent regulatory regions.

Purpose of the Study:

• To develop a novel computational framework, the Regulatory Region Scoring (RRS) model, for detecting functional conservation of regulatory sequences.
• To identify evolutionary and functional links between non-alignable enhancers using predicted transcription factor occupancy.
• To investigate the contribution of weak binding sites to the functional similarity of regulatory sequences.

Main Methods:

• Development of the Regulatory Region Scoring (RRS) model, an alignment-free computational framework.
• Prediction of transcription factor occupancy levels for regulatory sequences.
• Statistical analysis to determine functional and evolutionary links between enhancers based on binding site profiles (including presence, absence, and strength).

Main Results:

• The RRS model successfully detects functional and evolutionary links between non-alignable enhancers with statistical significance.
• Identified groups of enhancers likely to be under similar regulatory control.
• Demonstrated that the model captures sequence similarity through strong binding sites, weak binding sites, and statistically significant absence of sites.
• Provided evidence supporting the hypothesis that weak binding sites contribute to functional similarity.

Conclusions:

• The RRS model offers a new approach for identifying functional conservation in regulatory sequences, particularly for non-alignable regions.
• The model bridges the gap between detailed, data-intensive models and simpler, general models for regulatory sequence analysis.
• The findings highlight the importance of considering weak transcription factor binding sites for understanding enhancer function and evolution.