Evaluation of Automatic Analysis of Ultradeep Pyrosequencing Raw Data to Determine Percentages of HIV Resistance Mutations in Patients Followed-Up in Hospital

Phylogenetic analysis of HIV-1 archived DNA in blood and gut-associated lymphoid tissue in two patients under antiretroviral therapy

One of the approaches to cure human immunodeficiency virus (HIV) is the use of therapeutic vaccination. We have launched the Provir/Latitude 45 study to identify conserved CTL epitopes in archived HIV-1 DNA according to the HLA class I alleles in aviremic patients under antiretroviral therapy (ART). A HIV-1 polypeptidic therapeutic vaccine based on viral sequence data obtained from circulating blood was proposed; here, our aim was to compare the proviral DNA in blood and gut-associated lymphoid tissue (GALT). Peripheral blood mononuclear cells and gut biopsies were obtained from two HIV-1 infected patients under successful antiretroviral therapy. Total DNA was extracted including the proviral DNA. The HIV-1 reverse transcriptase was sequenced in both compartments using next generation sequencing followed by single genome sequencing; phylogenetic trees were established and compared. The proviral sequences of both compartments intra-patient exhibited a very low genetic divergence while it was possible to differentiate the sequences inter-patients; single genome sequencing analysis of two couples of samples confirmed that there was no compartmentalization of the sequences intra-patient. We conclude that, considering these two cases, the proviral DNA sequences in blood and GALT are similar and that the epitope analysis of HIV-1 provirus in blood should be considered as relevant to that observed in the GALT, a hard-to-reach major compartment, and can therefore be used for therapeutic vaccine approaches.

Provir/Latitude 45 study: A step towards a multi-epitopic CTL vaccine designed on archived HIV-1 DNA and according to dominant HLA I alleles

One of the approaches by which the scientific community is seeking to cure HIV is the use of therapeutic vaccination. Previous studies have highlighted the importance of the virus-specific CD8+ T cell cytotoxic responses for the immune control of HIV and have oriented research on vaccine constructs based on CTL epitopes from circulating HIV-1 strains. The clinical trials with therapeutic vaccines to date have had limited success likely due to (i) a discrepancy between archived CTL epitopes in the viral reservoir and those in circulating viruses before antiretroviral therapy (ART) initiation and (ii) the lack of strong affinity between the selected CTL epitopes and the HLA grooves for presentation to CD8+ cells. To overcome these limitations, we launched the Provir/Latitude 45 study to identify conserved CTL epitopes in archived HIV-1 DNA according to the HLA class I alleles of aviremic patients, most of whom are under ART. The near full-length genomes or Gag, Pol and Nef regions of proviral DNA were sequenced by Sanger and/or Next Generation Sequencing (NGS). The HLA-A and B alleles were defined by NGS or molecular analysis. The TuTuGenetics software, which moves a sliding window of 8 to 10 amino acids through the amino acid alignment, was combined with the Immune Epitope Data Base (IEDB) to automatically calculate the theoretical binding affinity of identified epitopes to the HLA alleles for each individual. We identified 15 conserved epitopes in Pol (11), Gag (3), and Nef (1) according to their potential presentation by the dominant HLA-A and B alleles and now propose to use the corresponding conserved peptides in a multi-epitopic vaccine (HLA-fitted VAC, HFVAC).

Evaluation of different analysis pipelines for the detection of HIV-1 minority resistant variants

OBJECTIVE

Reliable detection of HIV minority resistant variants (MRVs) requires bioinformatics analysis with specific algorithms to obtain good quality alignments. The aim of this study was to analyze ultra-deep sequencing (UDS) data using different analysis pipelines.

METHODS

HIV-1 protease, reverse transcriptase (RT) and integrase sequences from antiretroviral-naïve patients were obtained using GS-Junior® (Roche) and MiSeq® (Illumina) platforms. MRVs were defined as variants harbouring resistance-mutation present at a frequency of 1%-20%. Reads were analyzed using different alignment algorithms: Amplicon Variant Analyzer®, Geneious® compared to SmartGene® NGS HIV-1 module.

RESULTS

101 protease and 51 RT MRVs identified in 139 protease and 124 RT sequences generated with a GS-Junior® platform were analyzed using AVA® and SmartGene® software. The correlation coefficients for the MRVs were R2 = 0.974 for protease and R2 = 0.972 for RT. Discordances (n = 13 in protease and n = 15 in RT) mainly concerned low-level MRVs (i.e., with frequencies of 1%-2%, n = 18/28) and they were located in homopolymeric regions (n = 10/15). Geneious® and SmartGene® software were used to analyze 143 protease, 45 RT and 26 integrase MRVs identified in 172 protease, 69 RT, and 72 integrase sequences generated with a MiSeq® platform. The correlation coefficients for the MRVs were R2 = 0.987 for protease, R2 = 0.995 for RT and R2 = 0.993 for integrase. Discordances (n = 9 in protease, n = 3 in RT, and n = 3 in integrase) mainly concerned low-level MRVs (n = 13/15).

CONCLUSION

We found an excellent correlation between the various UDS analysis pipelines that we tested. However, our results indicate that specific attention should be paid to low-level MRVs, for which the use of two different analysis pipelines and visual inspection of sequences alignments might be beneficial. Thus, our results argue for use of a 2% threshold for MRV detection, rather than the 1% threshold, to minimize misalignments and time-consuming sight reading steps essential to ensure accurate results for MRV frequencies below 2%.

Sanger and Next-Generation Sequencing data for characterization of CTL epitopes in archived HIV-1 proviral DNA

One of the strategies for curing viral HIV-1 is a therapeutic vaccine involving the stimulation of cytotoxic CD8-positive T cells (CTL) that are Human Leucocyte Antigen (HLA)-restricted. The lack of efficiency of previous vaccination strategies may have been due to the immunogenic peptides used, which could be different from a patient's virus epitopes and lead to a poor CTL response. To counteract this lack of specificity, conserved epitopes must be targeted. One alternative is to gather as many data as possible from a large number of patients on their HIV-1 proviral archived epitope variants, taking into account their genetic background to select the best presented CTL epitopes. In order to process big data generated by Next-Generation Sequencing (NGS) of the DNA of HIV-infected patients, we have developed a software package called TutuGenetics. This tool combines an alignment derived either from Sanger or NGS files, HLA typing, target gene and a CTL epitope list as input files. It allows automatic translation after correction of the alignment obtained between the HxB2 reference and the reads, followed by automatic calculation of the MHC IC50 value for each epitope variant and the HLA allele of the patient by using NetMHCpan 3.0, resulting in a csv file as output result. We validated this new tool by comparing Sanger and NGS (454, Roche) sequences obtained from the proviral DNA of patients at success of ART included in the Provir Latitude 45 study and showed a 90% correlation between the quantitative results of NGS and Sanger. This automated analysis combined with complementary samples should yield more data regarding the archived CTL epitopes according to the patients' HLA alleles and will be useful for screening epitopes that in theory are presented efficiently to the HLA groove, thus constituting promising immunogenic peptides for a therapeutic vaccine.