An algorithm for mapping positively selected members of quasispecies-type viruses.
BACKGROUND
Many RNA viruses do not have a single, representative genome but instead form a set of related variants that has been called a quasispecies. The sequence variability of such viruses presents a significant bioinformatics challenge. In order for the sequence information to be understood, the complete mutational spectrum needs to be distilled to a biologically relevant and analyzable representation.


RESULTS
Here, we develop a "selection mapping" algorithm--QUASI--that identifies the positively selected variants of viral proteins. The key to the selection mapping algorithm is the identification of particular replacement mutations that are overabundant relative to silent mutations at each codon (e.g., threonine at hemagglutinin position 262). Selection mapping identifies such replacement mutations as positively selected. Conversely, selection mapping recognizes negatively selected variants as mutational "noise" (e.g., serine at hemagglutinin position 262).


CONCLUSIONS
Selection mapping is a fundamental improvement over earlier methods (e.g., dN/dS) that identify positive selection at codons but do not identify which amino acids at these codons confer selective advantage. Using QUASI's selection maps, we characterize the selected mutational landscapes of influenza A H3 hemagglutinin, HIV-1 reverse transcriptase, and HIV-1 gp120.
Stewart
J J
JJ
BioGenetic Ventures, 1100 Olive Way, Suite 300, Seattle, WA 98101 jjs@alumni.princeton.edu
Watts
P
P
Litwin
S
S
D000465:Algorithms; D000595:Amino Acid Sequence; D003062:Codon / Q000235:genetics; D019295:Computational Biology / Q000379:methods; D018604:Epitope Mapping / Q000379:methods; D014644:Genetic Variation; D015699:HIV Envelope Protein gp120 / Q000235:genetics; D015497:HIV-1 / Q000235:genetics*; D006389:Hemagglutinins, Viral / Q000235:genetics; D009980:Influenza A virus / Q000235:genetics*; D008969:Molecular Sequence Data; D012984:Software; D013045:Species Specificity
D003062:Codon; D015699:HIV Envelope Protein gp120; D006389:Hemagglutinins, Viral