title: Probabilistic inference of viral quasispecies subject to recombination.
abstract: RNA viruses exist in their hosts as populations of different but related strains. The virus population, often called quasispecies, is shaped by a combination of genetic change and natural selection. Genetic change is due to both point mutations and recombination events. We present a jumping hidden Markov model that describes the generation of viral quasispecies and a method to infer its parameters from next-generation sequencing data. The model introduces position-specific probability tables over the sequence alphabet to explain the diversity that can be found in the population at each site. Recombination events are indicated by a change of state, allowing a single observed read to originate from multiple sequences. We present a specific implementation of the expectation maximization (EM) algorithm to find maximum a posteriori estimates of the model parameters and a method to estimate the distribution of viral strains in the quasispecies. The model is validated on simulated data, showing the advantage of explicitly taking the recombination process into account, and applied to reads obtained from a clinical HIV sample.
authors: Department of Biosystems Science and Engineering, ETH Zurich, Basel, Switzerland.
authors: Armin
authors: A
authors: Töpfer
authors: Osvaldo
authors: O
authors: Zagordi
authors: Sandhya
authors: S
authors: Prabhakaran
authors: Volker
authors: V
authors: Roth
authors: Eran
authors: E
authors: Halperin
authors: Niko
authors: N
authors: Beerenwinkel
mesh_terms: D000465:Algorithms; D014644:Genetic Variation; D016679:Genome, Viral; D006678:HIV / Q000145:classification / Q000235:genetics*; D015658:HIV Infections / Q000821:virology; D059014:High-Throughput Nucleotide Sequencing; D006801:Humans; D008390:Markov Chains; D010802:Phylogeny; D017354:Point Mutation; D011995:Recombination, Genetic; D012641:Selection, Genetic; D054299:env Gene Products, Human Immunodeficiency Virus / Q000235:genetics*
chemical_list: D054299:env Gene Products, Human Immunodeficiency Virus
references: IEEE/ACM Trans Comput Biol Bioinform. 2014 Jan-Feb;11(1):182-91
references: 26355517
references: BMC Bioinformatics. 2011;12 Suppl 6:S1
references: 21989211
references: BMC Bioinformatics. 2011 Apr 26;12:119
references: 21521499
references: Science. 1991 Nov 15;254(5034):963-9
references: 1683006
references: PLoS Comput Biol. 2008 May 09;4(4):e1000074
references: 18437230
references: Am J Hum Genet. 2006 Apr;78(4):629-44
references: 16532393
references: Curr Opin Virol. 2011 Nov;1(5):413-8
references: 22440844
references: PLoS Med. 2008 Jul 29;5(7):e158
references: 18666824
references: Naturwissenschaften. 1971 Oct;58(10):465-523
references: 4942363
references: BMC Bioinformatics. 2011 Jan 05;12:5
references: 21208435
references: Nucleic Acids Res. 2010 Nov;38(21):7400-9
references: 20671025
references: PLoS Comput Biol. 2010 Dec 16;6(12):e1001022
references: 21187908
references: J Comput Biol. 2010 Mar;17(3):417-28
references: 20377454
references: Nature. 2006 Jan 19;439(7074):344-8
references: 16327776
references: Proc Natl Acad Sci U S A. 2005 Jan 4;102(1):158-62
references: 15615859
references: Nat Rev Genet. 2010 Jan;11(1):31-46
references: 19997069
references: BMC Genomics. 2011 May 19;12:245
references: 21592414
references: BMC Bioinformatics. 2011 Jan 28;12:38
references: 21276213
references: J Comput Biol. 2002;9(5):747-60
references: 12487762
references: BMC Bioinformatics. 2006 May 22;7:265
references: 16716226