A model-based information sharing protocol for profile Hidden Markov Models used for HIV-1 recombination detection
© Bulla et al.; licensee BioMed Central Ltd. 2014
Received: 19 November 2013
Accepted: 4 June 2014
Published: 19 June 2014
In many applications, a family of nucleotide or protein sequences classified into several subfamilies has to be modeled. Profile Hidden Markov Models (pHMMs) are widely used for this task, modeling each subfamily separately by one pHMM. However, a major drawback of this approach is the difficulty of dealing with subfamilies composed of very few sequences. One of the most crucial bioinformatical tasks affected by the problem of small-size subfamilies is the subtyping of human immunodeficiency virus type 1 (HIV-1) sequences, i.e., HIV-1 subtypes for which only a small number of sequences is known.
To deal with small samples for particular subfamilies of HIV-1, we introduce a novel model-based information sharing protocol. It estimates the emission probabilities of the pHMM modeling a particular subfamily not only based on the nucleotide frequencies of the respective subfamily but also incorporating the nucleotide frequencies of all available subfamilies. To this end, the underlying probabilistic model mimics the pattern of commonality and variation between the subtypes with regards to the biological characteristics of HI viruses. In order to implement the proposed protocol, we make use of an existing HMM architecture and its associated inference engine.
We apply the modified algorithm to classify HIV-1 sequence data in the form of partial HIV-1 sequences and semi-artificial recombinants. Thereby, we demonstrate that the performance of pHMMs can be significantly improved by the proposed technique. Moreover, we show that our algorithm performs significantly better than Simplot and Bootscanning.
Information sharing protocol
Profile Hidden Markov Models (pHMMs) are widely employed to model nucleotide or protein sequence families, in particular for database search [1, 2]. In numerous applications (see e.g. [3–5]), a sequence family classified into subfamilies is given, with each subfamily being modeled separately by one pHMM (i.e. only the information comprehended in the respective subfamily is used). If enough information is available about each subfamily (i.e. it consists of a sufficient number of sequences) this approach is feasible. Otherwise, the model is unable to detect sequences stemming from subfamilies composed of too few sequences .
To overcome this drawback, one can use an information sharing protocol which makes use of the information available about other subfamilies to model subfamilies composed of only a few sequences. Although they are usually named differently, information sharing protocols are a widely used concept in bioinformatics and other fields. For example, deducing pseudocounts  or Dirichlet mixtures  from a complete multiple sequence alignment to model a prior distribution for position-wise nucleotide or protein emission probabilities, constitutes a protocol for sharing information between different columns of the multiple sequence alignment (MSA). In speech recognition, using the same parameters for different HMMs within complex systems of HMMs in order to reduce the complexity of the system is called parameter tying . In fact, if a prior distribution in a Bayesian framework  is deduced from a superordinate entity, one can regard the estimation of probabilities characterizing a subordinate part of this entity by application of the estimated prior distribution as an information sharing protocol.
Sjölander et al.  suggested such a protocol for automated protein subfamily identification and classification. Their approach to model the emission frequencies of a protein subfamily S is based on the idea to add the nucleotide frequencies of other subfamilies T (in weighted form) to the ones of S as done with pseudocounts. Hereby, the added frequencies of T are weighted by the probability of the amino acid frequencies of T under a model whose emission probabilities are constructed based only on the frequencies of S.
HIV-1 is classified into the four main phylogenetic groups M, N, O, and P. All four were introduced into humans by separate zoonotic events, whereby the sources were different simian immunodeficiency viruses from chimpanzees  and gorilla , respectively. The M group is the most intensively studied one as it is responsible for the HIV pandemic. This group is further divided into ten subtypes, whereby Subtype A and F are further subdivided into sub-subtypes . Inter-subtype recombination is very frequent among HIV-1 subtypes : Up to now, 55 circulating recombinant forms (CRFs) have been reported , i.e., recombinant forms exhibiting at least three epidemiological independent sequences.
The subtypes of HIV-1 Group M can be divided into two groups: On the one hand, the ones for which at least 30 full-length sequences are available (subtypes A–G), on the other hand, subtypes H, J, and K, for which up to now only four (subtypes H and J) or two (Subtype K), respectively, complete-length sequences have been sampled. Thus, we refer to the first group – subtypes A–G – as “large-size” subtypes and to the second one – subtypes H–K – as “small-size” subtypes in what follows. Moreover, 11 full-length and about 1000 partial sequences have been labeled as “unknown”. That is, these sequences do not form an independent subtype or CRF. One can expect that by further sampling activities new small-size subtypes will be found, by finding sequences which are closely related to sequences of an unknown subtype already known, as well as by discovering sets of epidemiological independent sequences which are unrelated to the sequences now available .
Up to now, roughly 50 tools have been developed for the purpose of recognition of recombinants and breakpoint detection in HIV-1 sequences. Examples are Bootscanning , Simplot , Recco , the REGA HIV-1 Subtyping Tool , and MaxChi2 . In 2006, Schultz et al. developed the jumping profile Hidden Markov Model (jpHMM), an algorithm for subtyping, recombination analysis, and breakpoint detection, which we applied to HIV-1 and hepatitis C virus sequences . The jpHMM is composed of one pHMM for each viral subtype. In addition to the usual transitions within these pHMMs, it allows for so-called jumps between different pHMMs at nearby positions. Thus, it is possible to jump between states associated with different subtypes, depending on the local similarity of the query sequence to the subtypes. Although performing well for most input sequences, it did not correctly predict the HIV-1 subtypes of smaller size H, J, and K. Further analysis revealed that this deficit was due to the lack of an appropriate inter-subtype information sharing protocol, i.e., no information from the large-size subtypes was used to model the small-size subtypes.
The aim of this paper is to develop an information sharing protocol suitable for the application to HIV-1 Group M sequences. The introduced protocol is incorporated into the HIV recombination detection tool jumping profile Hidden Markov Model [4, 20, 21], short jpHMM. We investigate the influence of this modification on the performance of jpHMM. This application is of particular interest as reliable and accurate classification of HIV-1 (and other viral sequence data) facilitates to understand the influence of genetic diversity on host immune response, is crucial for epidemiological studies, and provides therapeutic decision support [22–24].
Small-size HIV subtypes
The effects of a small subtype size
Since a vast majority of columns in the HIV-1 MSA are strongly conserved over subtypes, it is very probable that in fact the most realistic way to model the emission probabilities in this example consists of using a joint emission probability for all subtypes. Hence, the modeled emission probabilities of Subtype K should be much closer to the ones of Subtype A and C. They should still be smaller (but to a lesser extent) since one has to account for the following (rather improbable) scenario: The nucleotide frequencies of Subtype K only coincidentally suggests that all subtypes should be modeled jointly, but the real emission probabilities are, in fact, considerably different.
A workaround for this kind of situation consists of using the same emission probabilities for all subtypes in case a column is completely conserved (notice that in our example this does not apply to the column under consideration). But given that there are about 2000 non-recombinant full-length HIV-1 Group M genomes reported by now, almost all columns in an alignment composed of all these sequences are not 100% conserved.
In a preparatory step, Sjölander’s approach performs a sequence weighting based on the total numbers of independent counts. Unfortunately, this procedure results in sequence weights that are far too low if applied to HIV nucleotide sequence data instead of protein amino acid data. Trying to circumvent this problem by leaving this step out causes the influence of other subtypes to become too strong. Since Sjölander’s approach is, in addition, an ad hoc one and aims at protein sequence data which are considerably more variable than HIV-1 sequences, we developed an alternative algorithm. This new algorithm is suitable for HIV-1 nucleotide sequence data and is derived from a probabilistic model of the interdependence of different HIV-1 subtypes.
Our model of the emission probabilities is based on the observation that for almost all sites in an HIV-1 MSA 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 them, the phenomenon that some but not all of the subtypes exhibit equal emission probabilities could be explained biologically as follows: Assume a site allows for more than one nucleotide to be present (i.e. at least two alleles are observed). Then, there are only very few characteristics of the virus which determine which degree of fitness the virus has for the different nucleotides possible at the respective site. Moreover, these characteristics take only very few values. As the characteristics and their values at a particular site are small in number, the number of different nucleotide distributions observed at the respective site is also small.
To clarify this reasoning, let us assume that for a site i in the MSA 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. For example, it could be the case that some part of the virus can assume two distinct shapes, influencing which nucleotides can be present at site isuch that the virus is able to survive. Then, the following might hold for the nucleotide distributions of the subtypes:
If the characteristics takes the value 0, it might be that the virus can only survive if adenine is present at site i. This leads to a nucleotide distribution where adenine has a probability very near to one.
In case it takes the value 1 the virus might also be able to survive if cytosine is present, but with a significant disadvantage with respect to its fitness. This might lead to a nucleotide distribution where adenine has a probability of, say, about 90% and the one of cytosine is about 10%.
In this case we would observe two different nucleotide distributions at site i: One part of the subtypes will show one nucleotide distribution and the other subtypes will show another. In the following, we will call the different nucleotide distributions (resp. emission probabilities) at the site “source”. That is, in the example just given there are two sources.
In order to estimate the emission probabilities of the subtypes, we iterate through all possible source combinations for each position. For each source combination we estimate the emission probabilities of each source based on the nucleotide frequencies of the subtypes assigned to it using a Dirichlet priora. We obtain the emission probabilities of a subtype by averaging over all source combinations, taking for each source combination the emission probabilities of the source the subtype is assigned to. Hereby, we weight each source combination by its likelihood. Further details are provided in the Subsection ‘Methods - Details on the information sharing protocol’.
In HIV, the distribution of gaps in the MSA provides little help in distinguishing different subtypes. Most parts of the genome are highly conserved, possessing nearly no insertions and deletions. In contrast, small parts are highly variable (possessing lots of gaps), but they vary so strongly that they provide little information for distinguishing different subtypes. For this reason, we estimate the transition probabilities (using a Dirichlet distribution as prior) for all subtypes jointly from the complete MSA. The same approach was applied by  for subfamilies belonging to the same protein family.
In order to determine the recombination pattern of a query sequence, the so-called Viterbi path  is computed, i.e., the most probable path of hidden states through the model generating this sequence. In this way, each position of the query sequence is assigned to exactly one parental subtype since i) each state of the model is assigned to exactly one pHMM (ignoring the special states at the begin and end of the model) and ii) each position of the query sequence is generated by one state of the model. The recombination breakpoints correspond to the positions of the jumps between different subtypes.
Details on the information sharing protocol
Let 1,…,F be the subtype (or, more generally, subfamily) indices. The individual sources in a source combination are indexed by 1,…,R. The space of all source combinations is denoted by S, the source of subtype i by s i . For each source j of a source combination , we denote the subtypes assigned to source j by . That is, if F=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 indicate that source combination by an additional index, in case several source combinations are considered). The number of nucleotides, denoted by N, is equal to 4.
Prior probability of number of sources
We denote the probability of a given number of sources j by P(R = j) = ρ j . These probabilities are estimated as described in : For each column in the alignment of all full-length HIV-1 Group M sequences, we decide based on the AIC (Akaike Information Criterion) by how many sources the nucleotide frequencies of that column are modeled in the most likely way. Then, the prior probability of a given number of sources R is estimated by the relative frequency of for how many alignment columns R was found to be the most likely number of sources. Since the frequencies for R>6 are zero or near to zero, we set the maximum number of sources R m a x to 6. The estimated values are (ρ j )j=1,…, 6=(0.37,0.30,0.20,0.097,0.029,0.0047).
Restricting the search space
Estimation of emission probabilities
In case we only have one subtype, we use the following Bayesian approach to estimate the emission probabilities of the subtype at a fixed site given the nucleotide frequencies at that site: We assume that the a priori distribution on the emission probabilities 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.
That is, we estimate the emission probability of Subtype 1 for each fixed source combination by . Then we sum these estimates over all possible source combinations, weighting the summands by the probability of the respective source combination given the observed nucleotide frequencies.
- 1.In  (last equation on p. 6), we showed that(2)
- 2.Applying Bayes’ rule, we obtain(3)
The sequences in T r e c are semi-artificial recombinants of real-world subtype sequences. The purpose of T r e c is to test the sensitivity of our algorithm with respect to the detection of segments originating from small-size subtypes. The sequences in T r e c are composed of three segments, each of length 250 bps. The first and last segment stem from a regular-size subtype (A-D, F, and G), the middle part stems from one of the small-size subtypes (H, J, and K). The choice to group the subtypes this way into normal- and small-size ones, is canonical since the subtypes H, J, and K contain only 2 resp. 3 sequences, whereas the other subtypes contain 30+ sequences. For each pair composed of a normal- and a small-size subtype, T r e c contains 10 semi-artificial recombinants composed of sequences from the respective pair of subtypes. That is, we have 180 sequences in total. The first position of the 750 bps long genome part from which a semi-artificial recombinant is constructed is chosen randomly between position 1007 and 8501 with respect to the HIV-1 reference sequence HXB2 (see ).
The set T p u r e is composed of parts of sequences from regular-size pure subtypes, i.e. partial subtype sequences. The purpose of T p u r e is to test the specificity of our algorithm with respect to the detection of segments originating from small-size and regular-size subtypes. Each sequence in T p u r e has a length of 2000 bps. For each subtype, T p u r e contains 20 partial sequences, i.e. 180 sequences in total. The starting position of the sequence part taken is chosen randomly between position 1007 and 8501 with respect to HXB2.
The motivation for using the two measures s p o s and s s e g m consists in each concentrating on a different type of classification error: the primary purpose of s p o s is to identify smaller shifts of the breakpoints, under the assumption that the predicted recombination patterns are correct. On the contrary, s s e g m does not aim for these errors, but is designed for detecting different recombination patterns. Furthermore, these scores are suitable for being used in an automatic evaluation. In contrast to that, using the distances between the correct and the predicted breakpoints (which is a more canonical approach) would require to deal with differences in the recombination patterns (like the recombination patterns AGA and AGAGA) manually. We use the arithmetic mean of s p o s and s s e g m as an overall score for a single test sequence, denoting it s t o t .
Notice that applying the introduced measures on the two test sets does not exactly measure what we stated as purpose of the test sets (sensitivity resp. specificity to detect segments originating from small-size subtypes): For T r e c we also consider misclassification in the first and last segment and for T p u r e we do not distinguish between wrongful inserts of segments of normal- and small-size subtypes. One can rather view T r e c and T p u r e as boundaries between which all reasonable application scenarios lie: No standard user will want to emphasize sensitivity more than measured by one of the introduced scores on T r e c and, analogously, no one specificity more than measured on T p u r e .
Versions of jpHMM
The evaluated models
jpHMM p r o b
jpHMM m l
jpHMM s c a l
jpHMM s e m i
jpHMM o r i g
Comparison of different versions of jpHMM
On T p u r e , jpHMM s e m i classifies 82.2% of the sequences correctly, jpHMM s c a l 85.6%, jpHMM m l 97.2%, and jpHMM p r o b 98.9%. Applying the sign test yields that jpHMM p r o b is better than jpHMM s e m i with p-value p=9·10−10, than jpHMM s c a l with p=6·10−8, and than jpHMM m l with p=0.125.
On T r e c , jpHMM s e m i achieves an average s t o t of 0.74, jpHMM m l of 0.83, and jpHMM p r o b of 0.90. Based on the Wilcoxon signed rank test jpHMM p r o b performs better than jpHMM s e m i with p<2·10−16 and outperforms jpHMM m l with p=3·10−6.
Comparison of jpHMM with Simplot and Bootscanning
Simplot employs a distance-based method, whereas Bootscanning is based on a phylogenetic approach. Both can be applied to both single sequences and groups of sequences, such as HIV-1 subtypes. Simplot as well as Bootscanning require a certain number of sequences (respectively, groups) as input (three for Simplot, four for Bootscanning, including the query sequence).
As neither Simplot nor Bootscanning allow for automatic execution, we restrict our comparison to 84 sequences of T p u r e and 36 sequences of T r e c . Since it turned out during our evaluation that correctly classifying segments from large subtypes is not a challenge for any algorithm, we compared the methods i) for each large-size subtype on four randomly selected sequences from T p u r e and ii) on all sequences of T p u r e from small-size subtypes. For the comparison on semi-artificial recombinants, we randomly picked two sequences for each combination of a large- and a small-size subtype from T r e c .
All four methods classify the 24 pure segments from large-size subtypes correctly. Out of the 60 segments from the small-size subtypes, jpHMM p r o b misclassifies 2, jpHMM m l 5, Simplot 11, and Bootscanning 6. Applying the sign test yields that jpHMM p r o b performs better than Simplot with p=0.006 and better than Bootscanning with p=0.11.
On the 36 semi-artificial recombinants taken from T r e c , jpHMM p r o b achieves an average s p o s of 0.926, jpHMM m l of 0.863, Simplot of 0.833, and Bootscanning of 0.839. Hereby, jpHMM p r o b achieves a higher score than Simplot for 28 sequences with Simplot better performing for the remaining 8 sequences. Comparing jpHMM p r o b with Bootscannnig, jpHMM p r o b performs better in 32 out of 36 cases, with Bootscanning being better for the rest. Applying the Wilcoxon signed-rank test yields that jpHMM p r o b performs better than Simplot with p=2·10−4 and better than Bootscanning with p=2·10−6.
Influence of the degree of conservation on the performance of jpHMM
The degree of conservation has an influence on the performance of jpHMM p r o b . The score achieved by jpHMM p r o b for sequences with identity lower than 70% is significantly lower than the one for sequences with identity higher than 70%. Above the threshold of 70%, the degree of conservation does not seem to have a significant influence on the prediction performance.
Application to CRF04
Notice that subtypes F and K as well as subtypes A, G, and J, respectively, are closely related . Hence, one can expect different algorithms to classify a particular segment as belonging to two different subtypes of one of these two phylogenetic groups. Furthermore, the original classification assigns parts of the sequence to Subtype U (“unknown”). Since the input alignment of jpHMM does not contain sequences of an unknown subtype, jpHMM cannot designate any part of a query sequence as Subtype U.
We also explored how the results of jpHMM p r o b change with different parameter values. Increasing p j u m p to 10−5 does not change the prediction at all, while for p j u m p =10−4 the number of segments is increased by four. Decreasing p j u m p to 10−7 or 10−8 yields the same result. The fragment from Subtype G, predicted in the vif gene, vanishes and instead Subtype A is predicted in this location.
On an Intel Core 2 Duo E8400 with 3 GHz and 4 GB of RAM, the application of jpHMM p r o b to a sequence from T r e c (length 750 bps) takes about 16 minutes. The lion’s share of the computational time is due to the calculation of the emission probabilities. It could be decreased considerably by implementing further heuristics to neglect source combinations with very low likelihood.
Discussion and conclusions
In this paper, we presented a novel information sharing protocol, realized in the form of a probabilistic model, for estimating the emission probabilities of pHMMs. We incorporated this protocol into the recombination detection program jpHMM and applied it to semi-artificial and real-world HIV-1 Group M sequences. As we were able to confirm an improvement of the prediction accuracy of jpHMM by our protocol, we demonstrated the value of the protocol for a crucial biological and medical application employing pHMMs.
The work of biological and medical researchers can be supported in various ways by improving the precision of recombination and breakpoint detection for HIV-1. For example, discovering recombination-free genome regions allows for the identification of barriers to genetic crossing between different viral strains. This kind of barriers then may help identifying vulnerable aspects of HIV biology in form of combinations of mutations strongly decreasing the survivability of an HI virus. This can facilitate finding promising targets for antiviral strategies (Talk of M. Negroni on HIV Dynamics & Evolution 2011).
Furthermore, recombination analysis and breakpoint detection do not only play a role for numerous viruses, e.g., flaviviruses, coronaviruses, alphaviruses, avian oncoviruses, rotaviruses, and influenza viruses [23, 33–39]. It is also needed for bacteria [40, 41] and within the scope of the detection of 16S rRNA chimera, which evolve during PCR [42–44]. Moreover, our approach could also be employed to improve software for remote homology detection of proteins. A heuristic probabilistic approach addressing the problem of small protein families has been introduced in .
We have applied i) the novel version of jpHMM incorporating our information sharing protocol, ii) former versions of jpHMM, and iii) the established recombination breakpoint detection tools Simplot and Bootscanning to i) sequences designated (in the LANL HIV database) as originating from a pure subtype, ii) semi-artificial recombinants of two real-world subtype sequences, iii) a full-length sequence of CRF04. In total, we compared the novel version of jpHMM with three formerly introduced ones, differing with regards to how they estimate the emission probabilities: i) a machine learning approach using linear discriminant learning, ii) a heuristic approach of scaling the pseudocounts, iii) the original approach without adaptions for small-size subtypes.
The results were compared using appropriate performance measures. On all test sequence sets, the novel version of jpHMM either significantly outperformed the other algorithms or achieved better results, but without reaching statistical significance. The latter is mostly a result of the fact that for most algorithms the task of correctly classifying segments of large-size subtypes is not a challenge, thus yielding very similar results for the corresponding type of test.
Compared with the machine-learning approach (jpHMM m l ) introduced in , the method presented in this paper (jpHMM p r o b ) seems to be considerably better performing on very small subtypes: On the one hand, jpHMM m l achieves an average score (s t o t ) of 0.87 on the semi-artificial recombinants constructed from Subtype H or J, respectively, whereas the average score on recombinants constructed from Subtype K is 0.76. That is, there is a difference in regard to the average score of 0.11 between recombinants built from Subtype H and J, respectively, and Subtype K. On the other hand, for jpHMM p r o b the average score is 0.92 for recombinants built from Subtype H and J, respectively, and 0.87 for Subtype K. That is, a difference of 0.05 can be observed. Thus, our model-driven information sharing approach seems in particular more suitable to model very small subtypes than a machine learning approach based on regularized linear discriminant learning and a transformation of the resulting weights into valid probabilities via a softmax function.
Astonishingly, for the tests on semi-artificial recombinants, the phylogenetic distance between the respective small- and large-size subtypes does not seem to have a notable influence on the performance of jpHMM p r o b . The phylogenetically nearest large-size subtype of Subtype H is Subtype C, of Subtype J is Subtype A, and of Subtype K is Subtype F. Nevertheless, only semi-artificial recombinants composed of Subtype J and A achieve considerably poorer results than recombinants composed of Subtype J and another large-size subtype. For Subtypes H and K, recombinants built from Subtypes C and F, respectively, do not yield a notably low performance. Since Subtypes K and F are phylogenetically closer than H and C (J and A, respectively), we can conclude that the low performance for recombinants built from Subtype J and A is probably simply coincidental and that the influence of the phylogenetic distance is at most weak for recombinants between HIV-1 Group M subtypes.
The relatively weak performance of jpHMM p r o b on partial sequences stemming from weakly conserved genome parts can easily be explained by the fact that jpHMM is not well suited for this kind of sequences: pHMMs assume that the group of sequences to be modeled can be well aligned. If not, the resulting pHMM consists mostly of insert states which do not constitute a sophisticated model for unalignable sequences.
One might wonder whether our test setting, with very short segments, is appropriate if the normal application requires considerably lower jumping probabilities. That is, one could suppose that our test setting does not prove that the algorithm performs well in correctly classifying a very long fragment (i.e. several thousand bps) stemming entirely from one subtype. But if the algorithm tended to assign such a fragment to more than one subtype, it would also misclassify shorter subfragments. Hence, by testing it on partial length sequences from one subtype we have in fact also covered these long fragments from one subtype.
Application to CRF04
In order to verify that jpHMM p r o b is capable of handling real world data, we applied it to a CRF04 sequence. This application shows – to some extent – why we stuck to semi-artificial recombinants instead of evaluating jpHMM p r o b on CRFs. For CRFs, it is unclear which classification is the “correct” one. This results from the scientific community not coming to any consensus on which subtype classifications are the correct ones. Beyond that, it is even discussed whether the entire classification system exhibits flaws: For example, it has been debated whether Subtype G is in fact a recombinant and CRF02 a pure subtype [45–48]. Moreover, the total number of CRFs is still quite small and a large fraction of them possesses only a few breakpoints, i.e., the number of potential test sequences is limited. Furthermore, it is a widespread approach to evaluate the performance of recombination detection tools using (semi-)artificial recombinants [21, 46, 49, 50]. Therefore, we chose to test jpHMM p r o b on semi-artificial recombinants, despite one could perceive testing on CRFs as more “realistic”.
We recently have developed a probabilistic algorithm called the “Unknown Subtype Finder (USF)”  for the purpose of automatically detecting parts of an input sequence which stem from a subtype yet unknown. For USF, the emission probabilities of an unknown subtype are also calculated employing the concept of sources which is the base of the information sharing protocol presented in this paper. Consequently, we plan the integration of USF and jpHMM into one tool. This tool then could be used to assign the known subtypes of HIV-1 to a (full- or part-length) query sequence as well as to detect segments of the genome originating from a subtype yet unknown.
Beyond this, one might think of several extensions to the probabilistic model presented in this paper:
The probability of whether the emission probabilities of a given set of subtypes are modeled jointly might depend on the subtypes in the set. Like this, it could be accounted for the fact that phylogenetic close subtypes should tend more to be modeled jointly than phylogenetic distant ones. As the phylogeny of the subtypes, which subtypes tend to be jointly modeled might change along the genome.
The information sharing protocol could be based on a hierarchical classification system (as it is with protein families, superfamilies, folds, and classes ) instead of a flat classification system (for example, HÌV-1 Group M subtypes). Similarly, how strongly the emission probabilities of an entity from the classification system (e.g. a protein family) is influenced by another entity could depend on how close the two entities are in the classification system. Analogously, one could also allow for unknown entities differing with regards to how distant they are to the known entities, e.g., unknown protein families and unknown protein superfamilies.
The probability of jumping between subtypes could be modeled as varying along the genome, motivated by the fact that the distribution of recombination breakpoints is not uniform among the genome .
The C++ source code of jpHMM and Matlab code for evaluating the different versions of jpHMM is available at http://jphmm.gobics.de.
a The explanation here emphasizes understandability. Due to computational efficiency, the actual computations are carried out in a quite different way.
b The first ten Bell numbers B(1) through B(10) are 1, 2, 5, 15, 52, 203, 877, 4140, 21147, 115975.
This work was supported by the Deutsche Forschungsgemeinschaft (STA 1009/5-1 and BU 2685/4-1) and the National Institutes of Health (AI087520). We thank Thomas Leitner for his comments on the draft.
- Krogh A, Brown M, Mian IS, Sjölander K, Haussler D: Hidden Markov Models in computational biology: applications to protein modeling. J Mol Biol. 1994, 235 (5): 1501-1531.View ArticlePubMedGoogle Scholar
- Eddy S: Profile Hidden Markov Models. Bioinformatics. 1998, 14 (9): 755-763.View ArticlePubMedGoogle Scholar
- Brown DP, Krishnamurthy N, Sjölander K: Automated protein subfamily identification and classification. PLoS Comput Biol. 2007, 3 (8): e160-View ArticlePubMed CentralPubMedGoogle 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-View ArticlePubMed CentralPubMedGoogle Scholar
- Laht S, Koua D, Kaplinski L, Lisacek F, Stöcklin R, Remm M: Identification and classification of conopeptides using profile Hidden Markov Models. Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics. 2012, 1824 (3): 488-492. [http://www.sciencedirect.com/science/article/pii/S1570963911003219],View ArticleGoogle Scholar
- Durbin R, Eddy SR, Krogh A: Biological Sequence Analysis: Probabilistic Models of Proteins and Nucleic Acids. 1998, UK: Cambridge University PressView 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
- Young S: The general use of tying in phoneme-based HMM speech recognisers. Comput. Appl. Biosci, Volume 3. 1992, 569-572.. [http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=225844],Google Scholar
- Bolstad WM: Introduction to Bayesian Statistics. 2007, New York: John Wiley & SonsView ArticleGoogle Scholar
- Hahn BH, Shaw GM, De Cock KM, Sharp PM: AIDS as a Zoonosis: scientific and public health implications. Science. 2000, 287 (5453): 607-614.View ArticlePubMedGoogle Scholar
- Plantier JC, Leoz M, Dickerson JE, De Oliveira F, Cordonnier F, Lemee V, Damond F, Robertson DL, Simon F: A new human immunodeficiency virus derived from gorillas. Nat Med. 2009, 15 (8): 871-872.View 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.View 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.View ArticlePubMedGoogle Scholar
- LANL HIV Databases: CRFs. [http://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html],
- Salminen MO, Carr JK, Burke DS, McCutchan FE: Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses. 1995, 11: 1423-1425.View ArticlePubMedGoogle Scholar
- Robertson DL, Sharp PM, McCutchan FE, Hahn BH: Recombination in HIV-1. Nature. 1995, 374: 124-126.View ArticlePubMedGoogle Scholar
- Maydt J, Lengauer T: Recco: recombination analysis using cost optimization. Bioinformatics. 2006, 22 (9): 1064-1071.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Smith JM: Analyzing the mosaic structure of genes. J Mol Evol. 1992, 34: 126-129.PubMedGoogle Scholar
- 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-W465.View ArticlePubMed CentralPubMedGoogle 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 issue): W647-W651.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Leitner T: The Molecular Epidemiology of Human Viruses. 2002, Berlin: SpringerView 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. Virol J. 2006, 3: 103-View ArticlePubMed CentralPubMedGoogle Scholar
- Seifert M, Gohr A, Strickert M, Grosse I: Parsimonious higher-order Hidden Markov Models for improved array-CGH analysis with applications toArabidopsis thaliana. PLoS Comput Biol. 2012, 8: e1.002286-View ArticleGoogle Scholar
- Viterbi A: Error bounds for convolutional codes and an asymptotically optimum decoding algorithm. IEEE Trans Inform Theor. 1967, 13 (2): 260-269.View ArticleGoogle Scholar
- Unterthiner T, Schultz AK, Bulla J, Morgenstern B, Stanke M, Bulla I: Detection of viral sequence fragments of HIV-1 subfamilies yet unknown. BMC Bioinformatics. 2011, 12: 93-View ArticlePubMed CentralPubMedGoogle Scholar
- LANL HIV Databases: Numbering Positions in HIV Relative to HXB2CG. Http://www.hiv.lanl.gov/content/sequence/HIV/REVIEWS/HXB2.html,
- Bulla I, Schultz AK, Meinicke P: Improving Hidden Markov Models for classification of human immunodeficiency virus-1 subtypes through linear classifier learning. Stat Appl Genet Mol. 2012,, 11. http://www.ncbi.nlm.nih.gov/pubmed/22499688,Google Scholar
- Gao F, Robertson DL, Carruthers CD, Li Y, Bailes E, Kostrikis LG, Salminen MO, Bibollet-Ruche F, Peeters M, Ho DD, Shaw GM, Sharp PM, Hahn BH: An isolate of human immunodeficiency virus type 1 originally classified as subtype I represents a complex mosaic comprising three different group M subtypes (A, G, and I). J Virol. 1998, 72 (12): 10234-10241.PubMed CentralPubMedGoogle Scholar
- Paraskevis D, Magiorkinis M, Vandamme AM, Kostrikis LG, Hatzakis A: Re-analysis of human immunodeficiency virus type 1 isolates from Cyprus and Greece, initially designated ‘subtype I’, reveals a unique complex A/G/H/K/? Mosaic pattern. J Gen Virol. 2001, 82 (3): 575-580.View ArticlePubMedGoogle Scholar
- McCutchan FE: Understanding the genetic diversity of HIV-1. AIDS. 2000, 14 (S3): S31-S44.PubMedGoogle Scholar
- Worobey M, Rambaut A, Holmes EC: Widespread intra-serotype recombination in natural populations of dengue virus. Proc Natl Acad Sci U S A. 1999, 96 (13): 7352-7357.View ArticlePubMed CentralPubMedGoogle Scholar
- Lai MM: RNA recombination in animal and plant viruses. Microbiol Mol Biol Rev. 1992, 56: 61-79.Google Scholar
- Hahn CS, Lustig S, Strauss EG, Strauss JH: Western equine encephalitis virus is a recombinant virus. Proc Natl Acad Sci U S A. 1988, 85 (16): 5997-6001.View ArticlePubMed CentralPubMedGoogle Scholar
- Shaikh R, Linial M, Coffin J, Eisenman R: Recombinant avian oncoviruses I. Alterations in the precursor to the internal structural proteins. Virology. 1978, 87 (2): 326-338.View ArticlePubMedGoogle Scholar
- Bergmann M, Garcia-Sastre A, Palese P: Transfection-mediated recombination of influenza A virus. J Virol. 1992, 66 (12): 7576-7580.PubMed CentralPubMedGoogle Scholar
- Orlich M, Gottwald H, Rott R: Nonhomologous recombination between the hemagglutinin gene and the nucleoprotein gene of an influenza virus. Virology. 1994, 204: 462-465.View ArticlePubMedGoogle Scholar
- Jarvis TC, Kirkegaard K: Poliovirus RNA recombination: mechanistic studies in the absence of selection. Virology. 1992, 11 (8): 3135-3145.Google Scholar
- Feil EJ, Spratt BG: Recombination and the population structures of bacterial pathogens. Annu Rev Microbiol. 2001, 55: 561-590.View ArticlePubMedGoogle Scholar
- Goss EM, Kreitman M, Bergelson J: Genetic diversity, recombination and cryptic clades in pseudomonas viridiflava infecting natural populations of Arabidopsis thaliana. Genetics. 2005, 169: 21-35.View ArticlePubMed CentralPubMedGoogle Scholar
- Hugenholtz P, Huber T: Chimeric 16S rDNA sequences of diverse origin are accumulating in the public databases. Int J Syst Evol Microbiol. 2003, 53: 289-293.View ArticlePubMedGoogle Scholar
- Ashelford KE, Chuzhanova NA, Fry JC, Jones AJ, Weightman AJ: At least 1 in 20 16S rRNA sequence records currently held in public repositories is estimated to contain substantial anomalies. Appl Environ Microbiol. 2005, 71 (12): 7724-7736.View ArticlePubMed CentralPubMedGoogle Scholar
- Bolling S: Chimärenerkennung bei 16S RNA. Diploma Thesis, University of Greifswald;. 2012,Google Scholar
- Abecasis AB, Lemey P, Vidal N, de Oliveira T, Peeters M, Camacho R, Shapiro B, Rambaut A, Vandamme AM: Recombination confounds the early evolutionary history of human immunodeficiency virus type 1: subtype G is a circulating recombinant form. J Virol. 2007, 81 (16): 8543-8551.View ArticlePubMed CentralPubMedGoogle Scholar
- Lemey P, Lott M, Martin D, Moulton V: Identifying recombinants in human and primate immunodeficiency virus sequence alignments using quartet scanning. BMC Bioinformatics. 2009, 10: 126-View ArticlePubMed CentralPubMedGoogle Scholar
- Bulla I, Schultz AK, Schreiber F, Zhang M, Leitner T, Korber B, Morgenstern B, Stanke M: HIV classification using the coalescent theory. Bioinformatics. 2010, 26 (11): 1409-1415.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang M, Foley B, Schultz AK, Macke J, Bulla I, Stanke M, Morgenstern B, Korber B, Leitner T: The role of recombination in the emergence of a complex and dynamic HIV epidemic. Retrovirology. 2010, 7: 25-View ArticlePubMed CentralPubMedGoogle Scholar
- Bruen TC, Philippe H, Bryant D: A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006, 172 (4): 2665-2681. [http://www.genetics.org/content/172/4/2665.abstract],View ArticlePubMed CentralPubMedGoogle Scholar
- Kosakovsky Pond SL, Posada D, Stawiski E, Chappey C, Poon AF, Hughes G, Fearnhill E, Gravenor MB, Leigh Brown AJ, Frost SD: An evolutionary model-based algorithm for accurate phylogenetic breakpoint mapping and subtype prediction in HIV-1. PLoS Comput Biol. 2009, 5 (11): e1000581-View ArticlePubMed CentralPubMedGoogle Scholar
- Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol. 1995, 247 (4): 536-540.PubMedGoogle Scholar
- Fan J, Negroni M, Robertson DL: The distribution of HIV-1 recombination breakpoints. Infect Genet Evol. 2007, 7 (6): 717-723.View 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 credited.