Maximum Likelihood (ML) tree inference reconstructs phylogenies from Multiple Sequence Alignments (MSAs). Since MSAs are inherently noisy, ML tools may experience overfitting, whereby the inferred top...

Motivation: Phylogenetic trees describe the evolutionary history among biological species based on their genomic data. Maximum Likelihood (ML) based phylogenetic inference tools search for the tree and evolutionary model that best explain the observed genomic data. Given the independence of likelihood score calculations between different genomic sites, parallel computation is commonly deployed. This is followed by a parallel summation over the per-site scores to obtain the overall likelihood score of the tree. However, basic arithmetic operations on IEEE 754 floating-point numbers, such as addition and multiplication, inherently introduce rounding errors. Consequently, the order by which floating-point operations are executed affects the exact resulting likelihood value since these operations are not associative. Moreover, parallel reduction algorithms in numerical codes re-associate operations as a function of the core count and cluster network topology, inducing different round-off errors. These low-level deviations can cause heuristic searches to diverge and induce high-level result discrepancies (e.g., yield topologically distinct phylogenies). This effect has also been observed in multiple scientific fields, beyond phylogenetics. Results: We observe that varying the degree of parallelism results in diverging phylogenetic tree searches (high level results) for over 31 % out of 10 130 empirical datasets. More importantly, 8 % of these diverging datasets yield trees that are statistically significantly worse than the best known ML tree for the dataset (AU-test, p < 0.05). To alleviate this, we develop a variant of the widely used phylogenetic inference tool RAxML-NG, which does yield bit-reproducible results under varying core-counts, with a slowdown of only 0 to 12.7 % (median 0.8 %) on up to 768 cores. We further introduce the ReproRed reduction algorithm, which yields bit-identical results under varying core-counts, by maintaining a fixed operation order that is independent of the communication pattern. ReproRed is thus applicable to all associative reduction operations – in contrast to competitors, which are confined to summation. Our ReproRed reduction algorithm only exchanges the theoretical minimum number of messages, overlaps communication with computation, and utilizes fast base-cases for local reductions. ReproRed is able to all-reduce (via a subsequent broadcast) 4.1 · 106 operands across 48 to 768 cores in 19.7 to 48.61 µs, thereby exhibiting a slowdown of 13 to 93 % over a non-reproducible all-reduce algorithm. ReproRed outperforms the state-of-the-art reproducible all-reduction algorithm ReproBLAS (offers summation only) beyond 10 000 elements per core. In summary, we re-assess non-reproducibility in parallel phylogenetic inference, present the first bit-reproducible parallel phylogenetic inference tool, as well as introduce a general algorithm and open-source code for conducting reproducible associative parallel reduction operations. ### Competing Interest Statement The authors have declared no competing interest. European Research Council, https://ror.org/0472cxd90, 882500 European Union, https://ror.org/019w4f821, 101087081

Main new features Adaptive tree search heuristic (Togkousidis 2023) , (Togkousidis 2024) Multiple fast branch support measures Integrated model testing Full Changelog: 1.2.2...2.0-beta2

Taxonomic classification in biodiversity studies is the process of assigning the anonymous sequences of a marker gene (barcode) to a specific lineage using a reference database that contains named seq...

AbstractMotivation. Genotype datasets typically contain a large number of single nucleotide polymorphisms for a comparatively small number of individuals.

YouTube video by Alexandros Stamatakis

Παρότι αρκετές ελληνικές επιστημονικές ομάδες θα μπορούσαν να ηγηθούν ερευνητικών έργων λόγω του πλούτου των αρχαιολογικών και παλαιοντολογικών δειγμάτων της χώρας, συχνά περιορίζονται στον ρόλο του σ...

The need for effective and informed analysis of biological data is increasing with the explosive growth of genomic data. A phylogenetic framework is central to many molecular evolutionary approaches …

The need for effective and informed analysis of biological data is increasing with the explosive growth of genomic data. A phylogenetic framework is central to many molecular evolutionary approaches …

Abstract. Estimating the statistical robustness of the inferred tree(s) constitutes an integral part of most phylogenetic analyses. Commonly, one computes

The need for effective and informed analysis of biological data is increasing with the explosive growth of genomic data. A phylogenetic framework is central to many molecular evolutionary approaches …

The need for effective and informed analysis of biological data is increasing with the explosive growth of genomic data. A phylogenetic framework is central to many molecular evolutionary approaches …

The need for effective and informed analysis of biological data is increasing with the explosive growth of genomic data. A phylogenetic framework is central to many molecular evolutionary approaches …

The increasing availability of genomes and other “omes” for organisms across the tree of life has revolutionized the study of evolutionary processes. The comparative approach focused on different ele…

Share your videos with friends, family, and the world

Working with cognate data involves handling synonyms, that is, multiple words that describe the same concept in a language. In the early days of language phylogenetics it was recommended to select...

High-performance computing (HPC) constitutes an energy-hungry endeavor, and any efficiency gains via hardware and software advances are quickly (over-)compensated by increased consumption (rebound eff...