Reference genome choice and filtering thresholds jointly influence phylogenomic analyses
Creators: name: Rick, Jessica
name: Brock, Chad
name: Lewanski, Alexander
name: Golcher-Benavides, Jimena
name: Wagner, Catherine
Publisher: name: Zenodo
Year: 2023
Subjects: subject: Phylogenetics
subject: imbalance
subject: Macroevolution
subject: Bioinformatics
subject: FOS: Computer and information sciences
subject: diversification rate
subject: Phylogenomics
subject: minor allele frequency
subject: reference genome
subject: RAD phylogenomics
subject: Genetics
subject: FOS: Biological sciences
subject: Ecology, Evolution, Behavior and Systematics

description: Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant-calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate to what extent the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find the greatest topological accuracy when filtering sites for minor allele count > 3–4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with minor allele count > 1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short read genomic data for phylogenetic inference.
description: Funding provided by: National Science FoundationCrossref Funder Registry ID: https://ror.org/021nxhr62Award Number: DEB-1556963
Funding provided by: National Institute of General Medical SciencesCrossref Funder Registry ID: https://ror.org/04q48ey07Award Number: 2P20GM103432
Funding provided by: National Aeronautics and Space AdministrationCrossref Funder Registry ID: https://ror.org/027ka1x80Award Number: NNX15AI08H
description: Data and supplementary material here are associated with the manuscript "Reference genome choice and filtering thresholds jointly influence phylogenetic analyses". Scripts can be found at https://github.com/jessicarick/refbias_scripts, and are archived on Zenodo at https://doi.org/10.5281/zenodo.5940690.