PARNAS: Objectively Selecting the Most Representative Taxa on a Phylogeny
Author:
Markin AlexeyORCID, Wagle Sanket, Grover Siddhant, Vincent Baker Amy L.ORCID, Eulenstein OliverORCID, Anderson Tavis K.ORCID
Abstract
AbstractThe use of next-generation sequencing technology has enabled phylogenetic studies with hundreds of thousands of taxa. Such large-scale phylogenies have become a critical component in genomic epidemiology in pathogens such as SARS-CoV-2 and influenza A virus. However, detailed phenotypic characterization of pathogens or generating a computationally tractable dataset for detailed phylogenetic analyses requires bias free subsampling of taxa. To address this need, we propose parnas, an objective and flexible algorithm to sample and select taxa that best represent observed diversity by solving a generalized k-medoids problem on a phylogenetic tree. parnas solves this problem efficiently and exactly by novel optimizations and adapting algorithms from operations research. For more nuanced selections, taxa can be weighted with metadata or genetic sequence parameters, and the pool of potential representatives can be user-constrained. Motivated by influenza A virus genomic surveillance and vaccine design, parnas can be applied to identify representative taxa that optimally cover the diversity in a phylogeny within a specified distance radius. We demonstrated that parnas is more efficient and flexible than current approaches, and applied it to select representative influenza A virus in swine genes derived from over 5 years of genomic surveillance data. Our objective selection of 4 to 6 strains selected every two years from the 16 distinct genetic clades were sufficient to cover 80% of diversity circulating in US swine. We suggest that this method, through the objective selection of representatives in a phylogeny, provides criteria for rational multivalent vaccine design and for quantifying diversity. PARNAS is available at https://github.com/flu-crew/parnas.
Publisher
Cold Spring Harbor Laboratory
Reference35 articles.
1. Ranking viruses: measures of positional importance within networks define core viruses for rational polyvalent vaccine development 2. Population dynamics of cocirculating swine influenza A viruses in the United States from 2009 to 2012;Influenza and Other Respiratory Viruses,2013 3. Anderson, T. K. , Macken, C. A. , Lewis, N. S. , Scheuermann, R. H. , Van Reeth, K. , Brown, I. H. , Swenson, S. L. , Simon, G. , Saito, T. , Berhane, Y. , Ciacci-Zanella, J. , Pereda, A. , Davis, C. T. , Donis, R. O. , Webby, R. J. , and Vincent, A. L. (2016). A phylogeny-based global nomenclature system and automated annotation tool for H1 hemagglutinin genes from swine influenza A viruses. mSphere, 1(6). 4. Anderson, T. K. , Chang, J. , Arendsee, Z. W. , Venkatesh, D. , Souza, C. K. , Kimble, J. B. , Lewis, N. S. , Davis, C. T. , and Vincent, A. L. (2020). Swine influenza A viruses and the tangled relationship with humans. Cold Spring Harbor perspectives in medicine, page a038737. 5. octoflushow: an interactive tool describing spatial and temporal trends in the genetic diversity of influenza a virus in us swine;Microbiology Resource Announcements,2021
|
|