Fully automated sequence alignment methods are comparable to, and much faster than, traditional methods in large data sets: an example with hepatitis B virus.

Autor: Catanach TA; Ornithology Department, Academy of Natural Sciences of Drexel University, Philadelphia, PA, United States of America.; Illinois Natural History Survey, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America.; Department of Wildlife and Fisheries Sciences, Texas A&M University, College Station, TX, United States of America., Sweet AD; Illinois Natural History Survey, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America.; Department of Entomology, Purdue University, West Lafayette, IN, United States of America., Nguyen ND; Computer Science and Engineering, University of San Diego, California, La Jolla, CA, United States of America., Peery RM; Department of Biology, University of Alberta, Edmonton, Alberta, Canada.; Department of Plant Biology, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America., Debevec AH; School of Integrative Biology, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America., Thomer AK; School of Information, University of Michigan-Ann Arbor, Ann Arbor, MI, United States of America., Owings AC; Program in Ecology, Evolution, and Conservation Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America., Boyd BM; Illinois Natural History Survey, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America.; Department of Entomology, University of Georga, Athens, GA, United States of America., Katz AD; Illinois Natural History Survey, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America.; Department of Entomology, University of Illinois at Urbana-Champaign, Champaign, IL, United States of America., Soto-Adames FN; Florida State Collection of Arthropods, Florida Department of Agriculture and Consumer Services, Gainesville, FL, United States of America.; Department of Entomology and Nematology, University of Florida, Gainesville, FL, United States of America., Allen JM; Biology Department, University of Nevada, Reno, Reno, NV, United States of America.
Jazyk: angličtina
Zdroj: PeerJ [PeerJ] 2019 Jan 03; Vol. 7, pp. e6142. Date of Electronic Publication: 2019 Jan 03 (Print Publication: 2019).
DOI: 10.7717/peerj.6142
Abstrakt: Aligning sequences for phylogenetic analysis (multiple sequence alignment; MSA) is an important, but increasingly computationally expensive step with the recent surge in DNA sequence data. Much of this sequence data is publicly available, but can be extremely fragmentary (i.e., a combination of full genomes and genomic fragments), which can compound the computational issues related to MSA. Traditionally, alignments are produced with automated algorithms and then checked and/or corrected "by eye" prior to phylogenetic inference. However, this manual curation is inefficient at the data scales required of modern phylogenetics and results in alignments that are not reproducible. Recently, methods have been developed for fully automating alignments of large data sets, but it is unclear if these methods produce alignments that result in compatible phylogenies when compared to more traditional alignment approaches that combined automated and manual methods. Here we use approximately 33,000 publicly available sequences from the hepatitis B virus (HBV), a globally distributed and rapidly evolving virus, to compare different alignment approaches. Using one data set comprised exclusively of whole genomes and a second that also included sequence fragments, we compared three MSA methods: (1) a purely automated approach using traditional software, (2) an automated approach including by eye manual editing, and (3) more recent fully automated approaches. To understand how these methods affect phylogenetic results, we compared resulting tree topologies based on these different alignment methods using multiple metrics. We further determined if the monophyly of existing HBV genotypes was supported in phylogenies estimated from each alignment type and under different statistical support thresholds. Traditional and fully automated alignments produced similar HBV phylogenies. Although there was variability between branch support thresholds, allowing lower support thresholds tended to result in more differences among trees. Therefore, differences between the trees could be best explained by phylogenetic uncertainty unrelated to the MSA method used. Nevertheless, automated alignment approaches did not require human intervention and were therefore considerably less time-intensive than traditional approaches. Because of this, we conclude that fully automated algorithms for MSA are fully compatible with older methods even in extremely difficult to align data sets. Additionally, we found that most HBV diagnostic genotypes did not correspond to evolutionarily-sound groups, regardless of alignment type and support threshold. This suggests there may be errors in genotype classification in the database or that HBV genotypes may need a revision.
Competing Interests: The authors declare there are no competing interests.
Databáze: MEDLINE