MCALIGN2: faster, accurate global pairwise alignment of non-coding DNA sequences based on explicit models of indel...

Area of Science:

• Genomics
• Bioinformatics
• Computational Biology

Background:

• Non-coding DNA is a vast genomic component across many species.
• Aligning non-coding DNA is crucial for understanding its function but is challenging.
• Existing alignment methods face difficulties with the complexity of non-coding sequences.

Purpose of the Study:

• To present an improved pair-hidden-Markov-Model (pair HMM) based method for global pairwise alignment of non-coding DNA.
• To enhance the accuracy and efficiency of non-coding DNA sequence alignment.

Main Methods:

• Developed MCALIGN2, a novel pair HMM based alignment method.
• Incorporated an explicit model for indel length frequency distribution.
• Utilized a deterministic global optimizer for maximum posterior probability alignment.
• Allowed for time-reversible models of nucleotide substitution.

Main Results:

• MCALIGN2 demonstrates excellent performance across various parameter values in simulations.
• The new method is up to ten times faster than the previous Monte Carlo based method (MCALIGN).
• MCALIGN2 shows higher accuracy in resolving indels compared to heuristic methods when using an accurate model.

Conclusions:

• MCALIGN2 generates superior alignment quality by integrating biological knowledge of indel distributions and nucleotide substitution models.
• The method outperforms other available sequence alignment tools for non-coding DNA alignment tasks.
• Explicit modeling improves alignment accuracy and provides a more robust approach for non-coding DNA analysis.