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

Abstract 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-naïve 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-naïve 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 (EC 50 =0.32 nM), and may represent novel pathways of HIV-1 resistance to modern treatments. Author summary There are many different types of HIV-1 around the world. Most of the research on how HIV-1 can become resistant to drug treatment has focused on the type (B) that is the most common in high-income countries. However, about 90% of infections around the world are caused by a type other than B. We used next-generation sequencing to analyze samples of HIV-1 from patients in Uganda (mostly infected by types A and D) for whom drug treatment failed to work, and whose infections did not fit the classic pattern of adaptation based on B. Next, we used machine learning to detect mutations in these virus populations that could explain the treatment outcomes. Finally, we experimentally added two candidate mutations identified by our analysis to a laboratory strain of HIV-1 and confirmed that they conferred drug resistance to the virus. Our study reveals new pathways that other types of HIV-1 may use to evolve resistance to drugs that make up the current recommended treatment for newly diagnosed individuals.