Evidence of absence treated as absence of evidence: The effects of variation in the number and distribution of gaps treated as missing data on the results of standard maximum likelihood analysis

[Display omitted] •Likelihood scores vary negatively with gap opening cost and amount of missing data.•Negative relationship due to altered nucleotide homologies, not missing data per se.•Effect of alignment parameters and missing data on topology is unpredictable. Although numerous studies have dem...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Molecular phylogenetics and evolution 2021-01, Vol.154, p.106966-106966, Article 106966
Hauptverfasser: Jacob Machado, Denis, Castroviejo-Fisher, Santiago, Grant, Taran
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:[Display omitted] •Likelihood scores vary negatively with gap opening cost and amount of missing data.•Negative relationship due to altered nucleotide homologies, not missing data per se.•Effect of alignment parameters and missing data on topology is unpredictable. Although numerous studies have demonstrated the theoretical and empirical importance of treating gaps as insertion/deletion (indel) events in phylogenetic analyses, the standard approach to maximum likelihood (ML) analysis employed in the vast majority of empirical studies codes gaps as nucleotides of unknown identity (“missing data”). Therefore, it is imperative to understand the empirical consequences of different numbers and distributions of gaps treated as missing data. We evaluated the effects of variation in the number and distribution of gaps (i.e., no base, coded as IUPAC “.” or “–”) treated as missing data (i.e., any base, coded as “?” or IUPAC “N”) in standard ML analysis. We obtained alignments with variable numbers and arrangements of gaps by aligning seven diverse empirical datasets under different gap opening costs using MAFFT. We selected the optimal substitution model for each alignment using the corrected Akaike Information Criterion in jModelTest2 and searched for optimal trees using GARLI. We also employed a Monte Carlo approach to randomly replace nucleotides with gaps (treated as missing data) in an empirical dataset to understand more precisely the effects of varying their number and distribution. To compare alignments, we developed four new indices and used several existing measures to quantify the number and distribution of gaps in all alignments. Our most important finding is that ML scores correlate negatively with gap opening costs and the amount of missing data. However, this negative relationship is not due to the increase in missing data per se—which increases ML scores—but instead to the effect of gaps on nucleotide homology. These variables also cause significant but largely unpredictable effects on tree topology.
ISSN:1055-7903
1095-9513
DOI:10.1016/j.ympev.2020.106966