Detection of novel HIV-1 drug resistance mutations by support vector analysis of deep sequence data and experimental validation

The global HIV-1 pandemic comprises many genetically divergent subtypes. Most of our understanding of drug resistance in HIV-1 derives from subtype B, which predominates in North America and western Europe. However, about 90% of the pandemic represents non-subtype B infections. Here, we use deep sequencing to analyze HIV-1 from infected individuals in Uganda who were either treatment-n aive or who experienced virologic failure on ART without the expected patterns of drug resistance. Our objective was to detect potentially novel associations between mutations in HIV-1 integrase and treatment outcomes in Uganda, where most infections are subtypes A or D. We retrieved a total of 380 archived plasma samples from patients at the Joint Clinical Research Centre (Kampala), of which 328 were integrase inhibitor-n aive and 52 were raltegravir (RAL)-based treatment failures.  Next, we developed a bioinformatic pipeline for alignment and variant calling of the deep sequence data obtained from these samples from a MiSeq platform (Illumina). To detect associations between within-patient polymorphisms and treatment outcomes, we used a support vector machine (SVM) for feature selection with multiple imputation to account for partial reads and low quality base calls. Candidate point mutations of interest were experimentally introduced into the HIV-1 subtype B NL4-3 backbone to determine susceptibility to RAL in U87.CD4.CXCR4 cells.  Finally, we carried out replication capacity experiments with wild-type and mutant viruses in TZM-bl cells in the presence and absence of RAL. Our analyses not only identified the known major mutation N155H and accessory mutations G163R and V151I, but also novel mutations I203M and I208L as most highly associated with RAL failure.  The I203M and I208L mutations resulted in significantly decreased susceptibility to RAL (44.0-fold and 54.9-fold, respectively) compared to wild-type virus (EC50=0.32 nM), and may represent novel pathways of HIV-1 resistance to modern treatments.