- Research article
- Open Access
Detection of viral sequence fragments of HIV-1 subfamilies yet unknown
© Unterthiner et al; licensee BioMed Central Ltd. 2011
- Received: 28 September 2010
- Accepted: 11 April 2011
- Published: 11 April 2011
Methods of determining whether or not any particular HIV-1 sequence stems - completely or in part - from some unknown HIV-1 subtype are important for the design of vaccines and molecular detection systems, as well as for epidemiological monitoring. Nevertheless, a single algorithm only, the Branching Index (BI), has been developed for this task so far. Moving along the genome of a query sequence in a sliding window, the BI computes a ratio quantifying how closely the query sequence clusters with a subtype clade. In its current version, however, the BI does not provide predicted boundaries of unknown fragments.
We have developed Unknown Subtype Finder (USF), an algorithm based on a probabilistic model, which automatically determines which parts of an input sequence originate from a subtype yet unknown. The underlying model is based on a simple profile hidden Markov model (pHMM) for each known subtype and an additional pHMM for an unknown subtype. The emission probabilities of the latter are estimated using the emission frequencies of the known subtypes by means of a (position-wise) probabilistic model for the emergence of new subtypes. We have applied USF to SIV and HIV-1 sequences formerly classified as having emerged from an unknown subtype. Moreover, we have evaluated its performance on artificial HIV-1 recombinants and non-recombinant HIV-1 sequences. The results have been compared with the corresponding results of the BI.
Our results demonstrate that USF is suitable for detecting segments in HIV-1 sequences stemming from yet unknown subtypes. Comparing USF with the BI shows that our algorithm performs as good as the BI or better.
- Query Sequence
- Emission Probability
- Nucleotide Frequency
- Profile Hide Markov Model
- Branching Index
An accurate and reliable classification of viral sequences data for human immunodeficiency virus-1 (HIV-1) and some other viruses of interest is important for epidemiological studies. It facilitates the understanding of the influence of genetic diversity on host immune response and provides therapeutic decision support [1–3]. As HIV-1 is, however, one of the genetically most variable viruses and genomic recombinations are frequent in HIV-1 , the task of classifying corresponding viral sequence data is a challenging one.
HIV-1 is classified into three main phylogenetic groups (M, N, and O), introduced into humans by separate zoonotic events (all stemming from simian immunodeficiency viruses (SIVs) in chimpanzees . The M group is responsible for the HIV pandemic, and it is divided into nine subtypes, with subtype A and F being subdivided into subsubtypes . Inter-subtype recombination occurs very frequently among HIV-1 subtypes : So far, 48 circulating recombinant forms have been reported .
Up to now, about fifty tools for classification of HIV genomes, recognition of recombinants, and breakpoint detection have been developed. Examples are the REGA HIV-1 Subtyping Tool , the Recombinant Identification Program (RIP) , the jumping profile Hidden Markov Model (jpHMM) [11, 12], the Recco , and the oligonucleotide-based method introduced in . Nevertheless, to our knowledge, however, only one algorithm, called the Branching Index (BI), has been developed for deciding whether an HIV-1 sequence in question stems - completely or in part - from a subtype still unknown [15, 16]. Notice that it is impossible to deduce unknown sequence segments using an existing subtype classification method, based on a probabilistic model such as jpHMM, and to identify regions of low a posteriori probabilities for all of the well known subfamilies (see paragraph 'Discussion and conclusions - Miscellaneous').
In view of the large and rapidly growing quantity of sequence data, the need for a fully automatic tool for pinning down boundaries of unknown fragments is increasing. Since the BI is based on a sliding window approach, it only provides a visualization of the breakpoint positions, but no report of their exact position. We have addressed this problem by developing a model-based algorithm, which automatically detects those boundaries by taking a multiple sequence alignment (MSA) grouped into subfamilies as a basis.
A comparison of our algorithm with the BI, regarding scope and performance, is carried out in the section 'Results - Comparison'.
The main input into our algorithm consists of i) an MSA representing the known sequences, with its sequences grouped into subfamilies, ii) a query sequence, iii) a classification of the query sequence with respect to the subfamilies (i.e. each position of the query sequence has to be assigned to a subfamily from the MSA) as main input.
Before the core algorithm is carried out, we take a preparatory step, allowing for the input of unaligned query sequences. More precisely, we align the query sequence to the given alignment with ClustalW , and remove columns which i) constitute insertions in the alignment by the query sequence, or ii) contain too many gaps in the alignment (we use a threshold of 10% gaps).
The main idea of our algorithm (in the following referred to as Unknown Subtype Finder or USF) is that of constructing two simple pHMMs (allowing neither deletions nor insertions): The first pHMM models the sequence of predicted subtypes F i for each position (in the following pHMM K) and the second pHMM models an unknown subtype (in the following pHMM U). That is, for the example given in Figure 1, the first ten positions of pHMM K are modelled on the basis of the nucleotide frequencies of Subtype B at those positions and the last eight positions on the basis of the frequencies of Subtype A. In addition to the transitions within these two pHMMs, we allow for jumps between them.
pHMMs K and U
pHMMs are widely used for modelling nucleotide and protein sequence families for the purpose of database searching (see [19, 20]). In particular, they are used to model the position-wise nucleotide distribution in an MSA. Standard pHMMs also allow for the modelling of insertions and deletions in the query sequence. But we do not use insertion or deletion states, as the sequences are already aligned (The high conservation of HIV-1 sequences allows for this approach). Hence, except for the initial and final states, our pHMMs are composed of so-called match states only. For decoding the most probable path through our model, we use the Viterbi algorithm .
We model pHMM K in the conventional way: For each position i in the alignment, we model the emission probabilities of the i-th state of the pHMM K on the basis of the nucleotide frequencies of F i . To this end, choosing a Bayesian approach to model the emission frequencies, we assume that the a priori distribution of is a Dirichlet distribution (see ), with parameter (estimated in ). The parameter may be interpreted as pseudo counts which are added to the nucleotide frequencies. The emission probabilities then are the corresponding relative frequencies of these modified nucleotide frequencies.
For pHMM U, we have to choose another approach, as the empirical nucleotide frequencies of an unknown subtype are not available. Hence, we try to deduce reasonable emission probabilities of an unknown subtype on the basis of the nucleotide frequencies of the known subtypes. For more details, see the paragraph 'Emission probabilities of pHMM U' in this subsection.
Jumps between pHMMs
Emission probabilities of pHMM U
In order to model the emission probabilities of pHMM U, we rely on the observation that for almost all sites for HIV-1 at least some of the subtypes share the same emission probabilities. In fact, for the majority of sites, it would be most plausible to assign equal emission probabilities to all subtypes. Neglecting the trivial case of all subtypes having the same emission probability assigned to, the phenomenon that some but not all of the subtypes show equal emission probabilities could be explained biologically as follows: If a site allows for more than one nucleotide to be present (i.e., if at least two alleles are observed), there are very few, discrete characteristics of the virus which determine the fitness of the virus, depending on the nucleotide present at the respective site. As the characteristics at a particular site are small in number and discrete, the number of corresponding nucleotide distributions is also small. To clarify that, let us assume that for a site i the dependence of the virus fitness on the nucleotide at site i is determined by a binary characteristics (values 0 and 1) of the virus. Then i) for the value 0, it might be that the virus can only survive if adenine is present at site i (leading to a nucleotide distribution where adenine has a probability very near to one), ii) for the value 1 the virus can survive if cytosine is present, with a significant disadvantage with respect to its fitness (leading to a nucleotide distribution where adenine has a probability of, say, about 90% and cytosine one about 10%). In the following we will call the different nucleotide distributions (resp. emission probabilities) at the site "sources". In the example just given there are two sources.
Examples of calculation of emission probabilities
To construct the left-hand part of the model, we use a Bayesian approach to determine position-wise an optimal number of sources and how the subtypes should be assigned to the sources. For each source the emission probabilities are estimated on the basis of the emission frequencies of the subtypes assigned to the source. The probability, with which a source is chosen, is proportional to the number of subtypes assigned to it. The right-hand part is modelled by a Dirichlet distribution with the same value for the parameter as in paragraph 'pHMMs K and U' of this subsection. We denote the a priori probability of a source involved, but yet unknown, by η ≪ 1. The estimates of the emission probabilities of the unknown subtype are obtained by averaging over both parts of the model, i.e. we use the expectations corresponding to the emission probabilities under the model. The details of the estimation procedure are given in the next subsection.
Details of pHMM U
Let 1,..., S be the subtype indices. If sources have been assigned to all subtypes, we speak of a source combination. The individual sources in a combination of r sources are indexed by 1,..., r. The space of all source combinations is denoted by Q, the source of subtype i by q i . For each source j of a source combination , we denote the subtypes assigned to source j by . That is, if S = 4 and the subtypes 1, 2, and 4 are assigned to Source 1, and the Subtype 3 to Source 2, we have m1 = 3, m2 = 1, , and (Notice that r and the are defined for a particular source combination, but that for the sake of readability we do not identify that source by an additional index, in case several sources are considered). The number of nucleotides, generally denoted by N, is equal to 5 in this case (We treat gaps as ordinary nucleotides). The nucleotide frequencies of subtype i at a fixed position of the genome are denoted by .
Prior probability of number of sources
Using (1) and the AIC (Akaike Information Criterion ), we deduce the most plausible source combination for each site and with that the most plausible number of sources. Estimating the ρ j as the empirical frequencies of the number of sources (considering all eligible sites), we obtain the values (ρ j )j = 1,2,3= (0.85, 0.09, 0.06). For the sake of computational efficiency, we restrict the number of sources to values lower or equal to 3. Notice that the number of sources to which one can restrict the algorithm depends on the scale of the intersubtype variation of the virus genome at the informative sites of the genome.
Estimation of emission probabilities
In this section, we present the results of i) the calibration of USF on a) artificial HIV-1 recombinants and b) non-recombinant HIV-1 sequences designated as having emerged from a known subtype, ii) the application of USF to a) SIV sequences and b) sequences designated as unknown in the LANL HIV database (in the following called "Subtype U" sequences), and iii) the comparison of USF and BI.
More precisely, we generate the following two sets of test sequences: (i) A set T sens for measuring the sensitivity with respect to the ability of the algorithm to detect genome segments stemming from an unknown subtype, and (ii) a set T spec for measuring the specificity. The set T sens is composed of 229 sequences generated by taking a sequence from subtypes A-D and F-G and replacing a segment of this sequence by a segment of a sequence from some other subtype. We call the subtype of the major part of the genome the 'base subtype' and the subtype of the inserted part of the genome the 'insertion subtype'. A preliminary analysis shows that in case the subtypes H, J, or K have been assigned to the query sequence (or a part of it), USF is not suitable for a reliable detection of unclassifiable genome parts. Hence, for the role of a base subtype, those subtypes are excluded from our analysis. Nevertheless, segments of them may play the role of insertion subtypes. Segments of the subtypes B and D may not be combined, due to the small phylogenetic distance of those subtypes. Moreover, the replaced segments have a length of 1000 positions and their position has been chosen randomly. T spec is composed of 265 sequences sampled from the genome-length sequences being classified as subtype A-D or F-G in the LANL HIV database (50 for all subtypes except for the subtypes F and G, for which only 35 and 30, respectively, sequences were available). For T spec the sequences were assigned to the subtype they stem from according to their LANL HIV database designation. Therefore, if classified correctly, the complete sequence is classified as known. Any detected unknown regions are counted as false positives. For T sens we determine the subtype classification using the jpHMM, excluding the subtypes H, J, and K from the assignable subtypes.
To facilitate the testing technically, we restrict our analysis to the positions 808 to 8781 with respect to HXB . Covering this part of the genome, we analyse the performance of USF in relatively conserved regions, as well as highly variable ones and we do not have to cope with the low number of sequences available for covering the boundary parts of the genome.
Theoretical determination of η
We have tried also to determine η by a theoretical approach. More precisely, we have simulated unknown subtypes by excluding a subtype from the data based on which the emission probabilities of pHMM U were estimated. We then have chosen η such that the emission frequencies of the excluded subtype is estimated best (with respect to maximum likelihood). Unfortunately, this approach has filed to values of η smaller by an entire order of magnitude than the values found by means of the calibration described above. Consequently, we refrain from using this theoretical approach.
SIV sequences and Subtype U sequences
In order to check whether USF correctly classifies very divergent sequences, we have applied it to five full-length SIV genomes (AF103818, DQ373063, EF394356, U42720, X52154) from different parts of the SIV clade. As before, we did not allow for assigning subtypes H, J, and K in the subtype classification. In the same way we have tested the 8 full-length Subtype U genomes (AF286236, AF457101, AY046058, EF029066, EF029067, EF029068, EF029069, FJ388921). Except the Subtype U sequence AY046058, all sequences have been correctly identified as completely unknown (about 8% of the genome have not been classified as unknown).
Comparison with the BI
The BI is a method based on distance and phylogeny. It determines which parts of a query sequence should be classified among known sequences. Moving along the genome of a query sequence with a sliding window, the BI computes a ratio quantifying how closely the query sequence clusters with a subtype clade. On the basis of this quantity, it determines whether the respective part of a query sequence is unclassifiable with respect to the known subtypes.
We apply the BI to a subset of T spec , as well as the SIV and Subtype U sequences used in the evaluation described in the subsection 'Results - SIV sequences and Subtype U sequences'. As we had to carry out the testing manually, using the web interface of the BI , we had to confine ourselves to a limited number of sequences from T spec and could not test the BI on T sens at all. (For the purpose of the latter, it would have been necessary to reestimate the parameters of the BI after having removed a subtype from the training data. That, however, the web interface available does not allow.)
Application of the BI to the 5 SIV sequences and the 8 Subtype U sequences from the subsection 'Results - SIV sequences and Subtype U sequences' yields valid results in 3 and 4 cases, respectively. Out of these 7 sequences, all but one Subtype U sequence (AY046058) are classified correctly as completely unknown, with about 6% of the genome of AY046058 being misclassified.
Using the two-sided Wilcoxon signed-rank test, the version of USF without postprocessing performs significantly better (with respect to our position-wise measure) than the BI for the subtypes A and B. For the Subtype F, the BI is significantly better than USF (p = 0.05). For the other subtypes, this test does not yield significant results. If USF is used in the version equipped with postprocessing, it yields significantly better results than the BI for the subtypes A and B, with the differences on the other subtypes being highly insignificant.
We have presented USF, a tool for detection of unclassifiable segments in viral sequences. Using a probabilistic, model-driven approach, the tool is suitable in principle for all species (or other taxa) which are subdivided into subfamilies i) without too many indels separating the subfamilies and ii) where the phylogenetic distances between the subfamilies are not too inhomogeneous.
We have applied USF to i) artificial recombinants of two subtypes (excluding one subtype from the training data to simulate an unknown subtype), ii) sequences designated (in the LANL HIV database) as originating from a pure subtype, iii) SIV sequences, and iv) Subtype U sequences. As far as feasible, we have compared our results with the only other tool available with the same scope, the Branching Index (BI).
Performance of USF
Analyzing the performance of USF by subtype, one can see that it performs considerably better (with respect to specificity) on the subtypes A-C than on D, F, and G, whereas it does not yield acceptable results for the subtypes H, J, and K. Its unsatisfactory performance on the last three subtypes does not come unexpectedly: The subtypes H, J, and K are composed of only 2 or 3 complete genome sequences, and that does not allow for a realistic modelling of the emission probabilities of a pHMM without using an information sharing protocol (see ). The weaker performance for subtypes D, F, and G might also be explicable by this effect, with the situation being obfuscated for the Subtype F by the fact that this subtype is divided into two subsubtypes.
The results of the application of USF to artificial recombinants can be explained in part also by the size of the involved subtypes: The poorest results are achieved when subtypes G or J, which both belong to the subtypes of smaller size, act as base subtype. Obviously, the size of the insertion subtype should not have any impact on the performance of USF (and the results also do not suggest that). Astonishingly, there does not seem to be a correlation between the phylogenetic distance of a pair of base and insertion subtypes and the performance of USF on the respective pair: Testing T sens involves 46 pairs of subtypes. Considering the 13 pairs with the lowest phylogenetic distance, none of them is among the 3 poorest performing pairs and 3 are among the 7 poorest performing. As we have observed a poor performance of USF when the subtypes B and D are the base and insertion subtypes, we may conclude that, if the phylogenetic distance of the subtype pair is above a certain threshold, the performance of USF does not seem to depend on how remotely the subtypes are related exactly.
Specificity of USF & BI
Comparing USF (employing the removal of very short segments in the outcome) with the BI with respect to specificity, USF, roughly speaking, performs better on some of the large size subtypes (A and B), whereas there are no significant differences on the large size Subtype C and the smaller size subtypes D, F, and G.
Sensitivity of USF & BI
For a comparison of the sensitivities we had to restrict ourselves to the SIV and Subtype U sequences. In spite of the importance of the sensitivity to assess the performance of USF and BI, the analysis of this characteristic had to be carried out on quite a small test set, due to technical limitations in the implementation of the BI. Except for the Subtype U sequence AY046058, all SIV and Subtype U sequences were classified as unknown by USF as well as the BI. Since both tools detect the same sequence as not completely unknown (although different segments were detected as known), this might be a hint that the classification of AY046058 as a pure Subtype U sequence is questionable. To conclude, our analysis does not reveal any significant differences between USF and the BI with respect to their sensitivity.
Versatility of USF & BI
With respect to versatility, the BI seems to be slightly inferior to USF (at least in their current versions). As it is not possible to determine η by a theoretical approach (as described in paragraph 'Results - Calibration - Theoretical determination of η'), both methods require a parameter calibration on training data when applied to a new species, respectively taxon. Regarding breakpoint positions, the BI only provides a graph from which the user would have to deduce the breakpoint positions by visual inspection. Hence, it is not possible to run any automated procedures on the BI if breakpoint positions are required.
In the near future, we plan to incorporate our method in the jpHMM. This would lead to a tool capable not only of assigning the known subtypes of HIV-1 (or subfamilies of other viruses or species) to a query sequence (or parts of it) but also of detecting segments of the genome stemming from a subtype yet unknown. Moreover, we are currently working on the implementation of an information sharing protocol for the jpHMM, which then would attenuate the poor performance of USF when applied to the small size subtypes.
In addition, it has been discussed whether the core gene of some D/E-recombinants of Hepatitis B virus (HBV) might stem from a clade which became rare or extinct . We will apply USF to HBV data in order to investigate this question.
Furthermore, it has been claimed that the HBV genotype G is a recombinant between i) an ancestor comparable in divergence to those between the genotypes A-E, contributing the S gene, and ii) an HBV variant which is much more divergent, contributing the rest of the genome . In the face of this finding, we plan to incorporate more than one unknown subtype in our model so that different degrees of divergence can be modelled.
As mentioned in the section 'Background', it is not possible to find unknown sequence segments by identifying regions of small a posterior probabilities for all of the known subfamilies when applying the jpHMM for example. That is easily exemplified as follows: Let us assume there were only two subtypes A and B and we examined a sequence stemming from an unknown subtype which is genetically considerably closer to Subtype A than to Subtype B. Then this sequence would achieve very large a posteriori probabilities for Subtype A and very small ones for Subtype B. Thus, it would falsely be classified as known.
USF is implemented in C++ and the source code is freely available (see additional file 1).
We would like to thank Thomas Leitner for the encouragement to develop USF and Heinrich Hering for proofreading. This work was supported by the Deutsche Forschungsgemeinschaft (STA 1009/5-1).
- Korber B, Gaschen B, Yusim K, Thakallapally R, Kesmir C, Detours V: Evolutionary and immunological implications of contemporary HIV-1 variation. Br Med Bull 2001, 58: 19–42. 10.1093/bmb/58.1.19View ArticlePubMedGoogle Scholar
- Leitner T: The molecular epidemiology of human viruses. Springer Berlin; 2002.View ArticleGoogle Scholar
- Hraber P, Fischer W, Bruno W, Leitner T, Kuiken C: Comparative analysis of hepatitis C virus phylogenies from coding and non-coding regions: the 5' untranslated region (UTR) fails to classify subtypes. Virology Journal 2006, 3: 103. 10.1186/1743-422X-3-103PubMed CentralView ArticlePubMedGoogle Scholar
- Rhodes T, Wargo H, Hu WS: High Rates of Human Immunodeficiency Virus Type 1 Recombination: Near-Random Segregation of Markers One Kilobase Apart in One Round of Viral Replication. J Virol 2003, 77(20):11193–11200. 10.1128/JVI.77.20.11193-11200.2003PubMed CentralView ArticlePubMedGoogle Scholar
- Hahn BH, Shaw GM, De KM, Sharp PM: AIDS as a Zoonosis: Scientific and Public Health Implications. Science 2000, 287(5453):607–614. 10.1126/science.287.5453.607View ArticlePubMedGoogle Scholar
- Robertson DL, Anderson JP, Bradac JA, Carr JK, Foley B, Funkhouser RK, Gao F, Hahn BH, Kalish ML, Kuiken C, Learn GH, Leitner T, McCutchan F, Osmanov S, Peeters M, Pieniazek D, Salminen M, Sharp PM, Wolinsky S, Korber B: HIV-1 nomenclature proposal. Science 2000, 288: 55–57. 10.1126/science.288.5463.55dView ArticlePubMedGoogle Scholar
- Hoelscher M, Dowling WE, Sanders-Buell E, Carr JK, Harris ME, Thomschke A, Robb ML, Birx DL, McCutchan FE: Detection of HIV-1 subtypes, recombinants, and dual infections in East Africa by a multi-region hybridization assay. AIDS 2002, 16: 2055–2064. 10.1097/00002030-200210180-00011View ArticlePubMedGoogle Scholar
- LANL HIV Databases: CRFs2011. [Http://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html]
- de Oliveira T, Deforche K, Cassol S, Salminen M, Paraskevis D, Seebregts C, Snoeck J, van Rensburg EJ, Wensing AMJ, van de Vijver DA, Boucher CA, Camacho R, Vandamme AM: An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics 2005, 21(19):3797–3800. 10.1093/bioinformatics/bti607View ArticlePubMedGoogle Scholar
- Recombinant Identification Program Web Interface[http://www.hiv.lanl.gov/content/sequence/RIP/RIP.html]
- Zhang M, Schultz AK, Calef C, Kuiken C, Leitner T, Korber B, Morgenstern B, Stanke M: jpHMM at GOBICS: a web server to detect genomic recombinations in HIV-1. Nucleic Acids Res 2006, 34(S2):W463–465. 10.1093/nar/gkl255PubMed CentralView ArticlePubMedGoogle Scholar
- Schultz AK, Zhang M, Bulla I, Leitner T, Korber B, Morgenstern B, Stanke M: jpHMM: Improving the reliability of recombination prediction in HIV-1. Nucl Acids Res 2009, (37 Web Server):W647–651. 10.1093/nar/gkp371Google Scholar
- Maydt J, Lengauer T: Recco: recombination analysis using cost optimization. Bioinformatics 2006, 22(9):1064–1071. 10.1093/bioinformatics/btl057View ArticlePubMedGoogle Scholar
- Pandit A, Sinha S: Using genomic signatures for HIV-1 sub-typing. BMC Bioinformatics 2010, 11(Suppl 1):S26. 10.1186/1471-2105-11-S1-S26PubMed CentralView ArticlePubMedGoogle Scholar
- Wilbe K, Salminen M, Laukkanen T, McCutchan F, Ray SC, Albert J, Leitner T: Characterization of novel recombinant HIV-1 genomes using the branching index. Virology 2003, 316: 116–125. 10.1016/j.virol.2003.08.004View ArticlePubMedGoogle Scholar
- Hraber P, Kuiken C, Waugh M, Geer S, Bruno WJ, Leitner T: Classification of hepatitis C virus and human immunodeficiency virus-1 sequences with the branching index. J Gen Virol 2008, 89(9):2098–2107. 10.1099/vir.0.83657-0PubMed CentralView ArticlePubMedGoogle Scholar
- Schultz AK, Zhang M, Leitner T, Kuiken C, Korber B, Morgenstern B, Stanke M: A jumping profile Hidden Markov Model and applications to recombination sites in HIV and HCV genomes. BMC Bioinformatics 2006, 7: 265. 10.1186/1471-2105-7-265PubMed CentralView ArticlePubMedGoogle Scholar
- Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucl Acids Res 2003, 31(13):3497–3500. 10.1093/nar/gkg500PubMed CentralView ArticlePubMedGoogle Scholar
- Krogh A, Brown M, Mian IS, Sjölander K, Haussler D: Hidden Markov Models in Computational Biology: Applications to Protein Modeling. Journal of Molecular Biology 1994, 235(5):1501–1531. 10.1006/jmbi.1994.1104View ArticlePubMedGoogle Scholar
- Eddy S: Profile hidden Markov models. Bioinformatics 1998, 14(9):755–763. 10.1093/bioinformatics/14.9.755View ArticlePubMedGoogle Scholar
- Viterbi A: Error bounds for convolutional codes and an asymptotically optimum decoding algorithm. Information Theory, IEEE Transactions 1967, 13(2):260–269. 10.1109/TIT.1967.1054010View ArticleGoogle Scholar
- Sjölander K, Karplus K, Brown M, Hughey R, Krogh A, Mian I, Haussler D: Dirichlet mixtures: a method for improved detection of weak but significant protein sequence homology. Comput Appl Biosci 1996, 12(4):327–345.PubMedGoogle Scholar
- Brown DP, Krishnamurthy N, Sjölander K: Automated Protein Subfamily Identification and Classification. PLoS Comput Biol 2007, 3(8):e160. 10.1371/journal.pcbi.0030160PubMed CentralView ArticlePubMedGoogle Scholar
- Akaike H: A new look at the statistical model identification. Automatic Control, IEEE Transactions 1974, 19(6):716–723. 10.1109/TAC.1974.1100705View ArticleGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree: Computing Large Minimum Evolution Trees with Profiles instead of a Distance Matrix. Mol Biol Evol 2009, 26(7):1641–1650. 10.1093/molbev/msp077PubMed CentralView ArticlePubMedGoogle Scholar
- Korber B, Foley B, Kuiken C, Pillai S, Sodroski J: Numbering Positions in HIV Relative to HXB2CG. Human Retroviruses and AIDS 1998, Los Alamos, NM: Theoretical Biology and Biophysics Group, Los Alamos National Laboratory 1998, 102–111.Google Scholar
- Branching Index Web Interface[http://www.hiv.lanl.gov/content/sequence/phyloplace/PhyloPlace.html]
- Simmonds P, Midgley S: Recombination in the Genesis and Evolution of Hepatitis B Virus Genotypes. J Virol 2005, 79(24):15467–15476. 10.1128/JVI.79.24.15467-15476.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Kato H, Orito E, Gish RG, Sugauchi F, Suzuki S, Ueda R, Miyakawa Y, Mizokami M: Characteristics of Hepatitis B Virus Isolates of Genotype G and Their Phylogenetic Differences from the Other Six Genotypes (A through F). J Virol 2002, 76(12):6131–6137. 10.1128/JVI.76.12.6131-6137.2002PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.