TopHap: Rapid inference of key phylogenetic structures from common haplotypes in large genome collections with limited diversity

Motivation Building reliable phylogenies from very large collections of sequences with a limited number of phylogenetically informative sites is challenging because sequencing errors and recurrent/backward mutations interfere with the phylogenetic signal, confounding true evolutionary relationships. Massive global efforts of sequencing genomes and reconstructing the phylogeny of SARS-CoV-2 strains exemplify these difficulties since there are only hundreds of phylogenetically informative sites and millions of genomes. For such datasets, we set out to develop a method for building the phylogenetic tree of genomic haplotypes consisting of positions harboring common variants to improve the signal-to-noise ratio for more accurate phylogenetic inference of resolvable phylogenetic features. Results We present the TopHap approach that determines spatiotemporally common haplotypes of common variants and builds their phylogeny at a fraction of the computational time of traditional methods. To assess topological robustness, we develop a bootstrap resampling strategy that resamples genomes spatiotemporally. The application of TopHap to build a phylogeny of 68,057 genomes (68KG) produced an evolutionary tree of major SARS-CoV-2 haplotypes. This phylogeny is concordant with the mutation tree inferred using the co-occurrence pattern of mutations and recovers key phylogenetic relationships from more traditional analyses. We also evaluated alternative roots of the SARS-CoV-2 phylogeny and found that the earliest sampled genomes in 2019 likely evolved by four mutations of the most recent common ancestor of all SARS-CoV-2 genomes. An application of TopHap to more than 1 million genomes reconstructed the most comprehensive evolutionary relationships of major variants, which confirmed the 68KG phylogeny and provided evolutionary origins of major variants of concern. Availability TopHap is available on the web at https://github.com/SayakaMiura/TopHap. Contact s.kumar@temple.edu

[1]  J. Bloom Recovery of Deleted Deep Sequencing Data Sheds More Light on the Early Wuhan SARS-CoV-2 Epidemic , 2021, bioRxiv.

[2]  Russell B. Corbett-Detig,et al.  Ultrafast Sample placement on Existing tRees (UShER) enables real-time phylogenetics for the SARS-CoV-2 pandemic , 2021, Nature Genetics.

[3]  Sudhir Kumar,et al.  MEGA11: Molecular Evolutionary Genetics Analysis Version 11 , 2021, Molecular biology and evolution.

[4]  David Haussler,et al.  Ultrafast Sample Placement on Existing Trees (UShER) Empowers Real-Time Phylogenetics for the SARS-CoV-2 Pandemic , 2020, bioRxiv.

[5]  Maxwell D. Sanderford,et al.  An evolutionary portrait of the progenitor SARS-CoV-2 and its dominant offshoots in COVID-19 pandemic , 2020, bioRxiv.

[6]  S. Salzberg,et al.  Rapid detection of inter-clade recombination in SARS-CoV-2 with Bolotie , 2020, bioRxiv.

[7]  Benoit Morel,et al.  Phylogenetic Analysis of SARS-CoV-2 Data Is Difficult , 2020, bioRxiv.

[8]  Xingguang Li,et al.  Phylogenetic and phylodynamic analyses of SARS-CoV-2 , 2020, Virus Research.

[9]  Edward C. Holmes,et al.  A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology , 2020, Nature Microbiology.

[10]  R. Nielsen,et al.  Assessing Uncertainty in the Rooting of the SARS-CoV-2 Phylogeny , 2020, bioRxiv.

[11]  F. Balloux,et al.  Emergence of genomic diversity and recurrent mutations in SARS-CoV-2 , 2020, Infection, Genetics and Evolution.

[12]  X. Wan,et al.  Are pangolins the intermediate host of the 2019 novel coronavirus (SARS-CoV-2)? , 2020, PLoS pathogens.

[13]  E. Holmes,et al.  The proximal origin of SARS-CoV-2 , 2020, Nature Medicine.

[14]  X. Wan,et al.  Are pangolins the intermediate host of the 2019 novel coronavirus (2019-nCoV) ? , 2020, bioRxiv.

[15]  Kai Zhao,et al.  A pneumonia outbreak associated with a new coronavirus of probable bat origin , 2020, Nature.

[16]  Alexey M. Kozlov,et al.  RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference , 2018, bioRxiv.

[17]  Nicola De Maio,et al.  Reversible polymorphism-aware phylogenetic models and their application to tree inference. , 2016, Journal of theoretical biology.

[18]  N. Beerenwinkel,et al.  Tree inference for single-cell data , 2016, Genome Biology.

[19]  N. Navin,et al.  The first five years of single-cell cancer genomics and beyond , 2015, Genome research.

[20]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[21]  Denis Krompass,et al.  Performance, Accuracy, and Web Server for Evolutionary Placement of Short Sequence Reads under Maximum Likelihood , 2011, Systematic biology.

[22]  Remco R. Bouckaert,et al.  DensiTree: making sense of sets of phylogenetic trees , 2010, Bioinform..

[23]  E. Virginia Armbrust,et al.  pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree , 2010, BMC Bioinformatics.

[24]  Douglas E. Soltis,et al.  Applying the Bootstrap in Phylogeny Reconstruction , 2003 .

[25]  P. Lewis A likelihood approach to estimating phylogeny from discrete morphological character data. , 2001, Systematic biology.

[26]  J. Felsenstein CONFIDENCE LIMITS ON PHYLOGENIES: AN APPROACH USING THE BOOTSTRAP , 1985, Evolution; international journal of organic evolution.

[27]  Toshiyuki Yamamoto,et al.  CONFLICT OF INTEREST: None declared. , 2013 .

[28]  L. O.,et al.  A Likelihood Approach to Estimating Phylogeny from Discrete Morphological Character Data , 2002 .