ClassyFlu: classification of influenza A viruses with Discriminatively trained profile-HMMs

• Published in: PloS one 2014-01, Vol.9 (1), p.e84558-e84558
• Main Authors: Van der Auwera, Sandra, Bulla, Ingo, Ziller, Mario, Pohlmann, Anne, Harder, Timm, Stanke, Mario
• Format: Article
• Language: English
• Subjects: Algorithms; Analysis; animal phylogenetics; Animals; Avian flu; Avian influenza; Avian influenza viruses; BASIC BIOLOGICAL SCIENCES; Bioinformatics; Biology; BLAST algorithm; Cladistic analysis; Cladistics; Classification; Computational Biology - methods; Computer Science; Databases, Factual; Deoxyribonucleic acid; Diagnostic systems; DNA; DNA sequencing; Epidemiology; Evolution; Exo-a-sialidase; Gene sequencing; Genomes; Genomics; Genotype; H5N1; Health aspects; Hemagglutinin Glycoproteins, Influenza Virus - genetics; Hemagglutinins; hidden Markov models; Humans; Immunology; Influenza; Influenza A; influenza A virus; Influenza A virus - classification; Influenza A virus - genetics; Internet; Laboratories; Lectins; Markov Chains; Markov processes; Mathematics; Medicine; Molecular biology; Molecular Sequence Annotation; Neuraminidase - genetics; Nucleotide sequence; Perl; phylogenetic analysis; Phylogenetics; Phylogeny; Retraining; sequence databases; Trends; Veterinary Science; viral pathogens; Viruses; Web Browser; Web services
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0084558

Accurate and rapid characterization of influenza A virus (IAV) hemagglutinin (HA) and neuraminidase (NA) sequences with respect to subtype and clade is at the basis of extended diagnostic services and implicit to molecular epidemiologic studies. ClassyFlu is a new tool and web service for the classification of IAV sequences of the HA and NA gene into subtypes and phylogenetic clades using discriminatively trained profile hidden Markov models (HMMs), one for each subtype or clade. ClassyFlu merely requires as input unaligned, full-length or partial HA or NA DNA sequences. It enables rapid and highly accurate assignment of HA sequences to subtypes H1-H17 but particularly focusses on the finer grained assignment of sequences of highly pathogenic avian influenza viruses of subtype H5N1 according to the cladistics proposed by the H5N1 Evolution Working Group. NA sequences are classified into subtypes N1-N10. ClassyFlu was compared to semiautomatic classification approaches using BLAST and phylogenetics and additionally for H5 sequences to the new "Highly Pathogenic H5N1 Clade Classification Tool" (IRD-CT) proposed by the Influenza Research Database. Our results show that both web tools (ClassyFlu and IRD-CT), although based on different methods, are nearly equivalent in performance and both are more accurate and faster than semiautomatic classification. A retraining of ClassyFlu to altered cladistics as well as an extension of ClassyFlu to other IAV genome segments or fragments thereof is undemanding. This is exemplified by unambiguous assignment to a distinct cluster within subtype H7 of sequences of H7N9 viruses which emerged in China early in 2013 and caused more than 130 human infections. http://bioinf.uni-greifswald.de/ClassyFlu is a free web service. For local execution, the ClassyFlu source code in PERL is freely available.