Development of a novel human immunodeficiency virus type 1 subtyping tool, subtype analyzer (STAR): Analysis of subtype distribution in London

We have developed a high throughput computational tool for assigning subtype to HIV-1, based solely on protease and reverse transcriptase (PR-RT) amino acid sequence, generated routinely for clinical assessment of genotypic drug resistance. Subtype-specific profiles were created by generation of position-specific scoring matrices (PSSMs) from multiple amino acids alignments of HIV-1 sequence data from GenBank, phylogenetically divided into subtypes A, AG, B, C, D, F/K, G, H, and J and the separate groups N and O. Query sequences of unknown subtype are aligned with these profiles and a score is derived by comparing each amino acid position in the unknown sequence to the normalized frequency distribution of amino acids at the corresponding positions in the subtype alignments. The highest score is used to assign subtype to the query sequence. Leave one out cross-validation analysis showed the Subtype Analyzer (STAR) was 99% accurate in subtype assignation. STAR can be updated with additional subtype-specific sequence data from sequence databases. STAR was used to classify HIV-1 PR-RT sequences from 843 HIV-1 clinical isolates submitted for drug resistance profiling in London. Within this dataset 26.9% of sequences were classified by STAR as non-B subtypes.