Sequencing an Ashkenazi reference panel supports population-targeted personal genomics and illuminates Jewish and European origins.

The Ashkenazi Jewish (AJ) population is a genetic isolate close to European and Middle Eastern groups, with genetic diversity patterns conducive to disease mapping. Here we report high-depth sequencing of 128 complete genomes of AJ controls. Compared with European samples, our AJ panel has 47% more novel variants per genome and is eightfold more effective at filtering benign variants out of AJ clinical genomes. Our panel improves imputation accuracy for AJ SNP arrays by 28%, and covers at least one haplotype in ≈ 67% of any AJ genome with long, identical-by-descent segments. Reconstruction of recent AJ history from such segments confirms a recent bottleneck of merely ≈ 350 individuals. Modelling of ancient histories for AJ and European populations using their joint allele frequency spectrum determines AJ to be an even admixture of European and likely Middle Eastern origins. We date the split between the two ancestral populations to ≈ 12-25 Kyr, suggesting a predominantly Near Eastern source for the repopulation of Europe after the Last Glacial Maximum.

Efficacy of re-treatment with TMC435 as combination therapy in hepatitis C virus-infected patients following TMC435 monotherapy.

In the TMC435-C101 study, 6 patients infected with hepatitis C virus genotype 1 were treated with the protease inhibitor TMC435 (200 mg once daily) as monotherapy for 5 days. Approximately 1.5 years later, 5 of these patients were re-treated with TMC435 (200 mg once daily) plus pegylated interferon alfa-2a and ribavirin (PegIFNα-2a and RBV) for 4 weeks, followed by PegIFNα-2a and RBV until week 48 (in the Optimal Protease inhibitor Enhancement of Response to therApy [OPERA-1] study). TMC435-resistant variants, which emerged in all 5 patients during the TMC435-C101 study, were no longer detected at the beginning of the OPERA-1 study based on virus population sequencing. During the OPERA-1 study, 3 patients had a sustained virologic response; deep sequencing indicated low-level persistence of resistant variants in the remaining 2 patients, which might have affected their response to re-treatment.

Study of genotypic and phenotypic HIV-1 dynamics of integrase mutations during raltegravir treatment: a refined analysis by ultra-deep 454 pyrosequencing.

BACKGROUND: The dynamics of raltegravir-resistant variants and their impact on virologic response in 23 HIV-1-infected patients, who started a salvage raltegravir-containing regimen, were investigated. METHODS: Integrase population sequencing and Ultra-Deep-454 Pyrosequencing (UDPS) were performed on plasma samples at baseline and at raltegravir failure. All integrase mutations detected at a frequency ≥1% were considered to be reliable for the UDPS analyses. Phylogenetic and phenotypic resistance analyses were also performed. RESULTS: At baseline, primary resistance mutations were not detected by both population and UDPS genotypic assays; few secondary mutations (T97A-V151I-G163R) were rarely detected and did not show any statistically association either with virologic response at 24-weeks or with the development of resistant variants at failure. At UDPS, not all resistant variants appearing early during treatment evolved as major populations during failure; only specific resistance pathways (Y143R-Q148H/R-N155H) associated with an increased rate of fitness and phenotypic resistance were selected. CONCLUSIONS: Resistance to raltegravir in integrase strand transfer inhibitor-naive patients remains today a rare event, which might be changed by future extensive use of such drugs. In our study, pathways of resistance at failure were not predicted by baseline mutations, suggesting that evolution plus stochastic selection plays a major role in the appearance of integrase-resistance mutations, whereas fitness and resistance are dominant factors acting for the late selection of resistant quasispecies.

TMC310911, a novel human immunodeficiency virus type 1 protease inhibitor, shows in vitro an improved resistance profile and higher genetic barrier to resistance compared with current protease inhibitors.

TMC310911 is a novel human immunodeficiency virus type 1 (HIV-1) protease inhibitor (PI) structurally closely related to darunavir (DRV) but with improved virological characteristics. TMC310911 has potent activity against wild-type (WT) HIV-1 (median 50% effective concentration [EC(50)], 14 nM) and a wide spectrum of recombinant HIV-1 clinical isolates, including multiple-PI-resistant strains with decreased susceptibility to currently approved PIs (fold change [FC] in EC(50), >10). For a panel of 2,011 recombinant clinical isolates with decreased susceptibility to at least one of the currently approved PIs, the FC in TMC310911 EC(50) was ≤ 4 for 82% of isolates and ≤ 10 for 96% of isolates. The FC in TMC310911 EC(50) was ≤ 4 and ≤ 10 for 72% and 94% of isolates with decreased susceptibility to DRV, respectively. In vitro resistance selection (IVRS) experiments with WT virus and TMC310911 selected for mutations R41G or R41E, but selection of resistant virus required a longer time than IVRS performed with WT virus and DRV. IVRS performed with r13025, a multiple-PI-resistant recombinant clinical isolate, and TMC310911 selected for mutations L10F, I47V, and L90M (FC in TMC310911 EC(50) = 16). IVRS performed with r13025 in the presence of DRV required less time and resulted in more PI resistance-associated mutations (V32I, I50V, G73S, L76V, and V82I; FC in DRV EC(50) = 258). The activity against a comprehensive panel of PI-resistant mutants and the limited in vitro selection of resistant viruses under drug pressure suggest that TMC310911 represents a potential drug candidate for the management of HIV-1 infection for a broad range of patients, including those with multiple PI resistance.

Minor variant detection in amplicons using 454 massive parallel pyrosequencing: experiences and considerations for successful applications.

Ultra-deep sequencing (UDS) of amplicons is a major application for next-generation sequencing technologies, even more so for the 454 Genome Sequencer FLX. Especially for this application, errors that might be introduced during any of the sample processing or data analysis steps should be avoided or at least recognized, as they might lead to aberrant sequence variant calling. Since 454 pyrosequencing relies on PCR-driven target amplification, it is key to differentiate errors introduced during the amplification step from genuine minority variants. Thereto, optimal primer design is imperative because primer selection, primer dimer formation, and nonspecific binding may all affect the quality and outcome of amplicon-based deep sequencing. Also, other intrinsic PCR characteristics including amplification drift and the formation of secondary structures may influence sequencing data quality. We illustrate these phenomena using real life case studies and propose experimental and analytical evidence-based solutions for effective practice. Furthermore, because accuracy of the DNA polymerase is vital for reliable UDS results, a comparative analysis of error profiles from seven different DNA polymerases was performed and experimentally assessed in parallel by 454 sequencing. Finally, intra and interrun variability evaluation of the 454 sequencing protocol revealed highly reproducible results in amplicon-based UDS.

HIV-1 nucleotide mixture detection in the virco(®)TYPE HIV-1 genotyping assay: a comparison between Sanger sequencing and 454 pyrosequencing.

HIV-1 Protease (PR) and Reverse Transcriptase (RT) genotyping is well established for the management of antiretroviral (ARV) drug therapy, as it is able to detect gene mutations encoding resistance to ARV compounds or drug classes, that are associated with reduced drug susceptibility (i.e. phenotype). A correct phenotypic interpretation from the derived PR-RT genotype (i.e. virtual phenotype), requires a well characterized geno-phenotype correlative database and appropriate statistical predictive models. The applicability of the virtual phenotype for the patient, will, however, not only depend on the accuracy of the statistical models and the database they rely on, but also depend largely on the sequence information that is provided. Since HIV-1 evolves as a complex of closely related but non-identical viral genomes (i.e. quasispecies) it is crucial that the sequencing method used, is able to characterize most of the genetic mixtures that make up the different quasispecies within a single patient. US regulatory agencies require that developers of HIV-1 genotyping assays, determine and report the HIV-1 mixture detection level of their assay. Hence, the mixture scoring sensitivity of the population-based Sanger sequencing method, along with the defined mixture scoring rules, used to drive the virco(®)TYPE HIV-1 virtual phenotype, was investigated by comparing it to the 454 pyrosequencing technique, which is able to generate the complete viral population sequence. To this end the PR-RT coding sequence of 20 clinical isolates was determined by both sequencing methodologies. The genotyping assay which feeds the virco(®)TYPE HIV-1 virtual phenotype was able to call automatically 97.5% (i.e. 268 mixtures) and 95.3% (i.e. 326 mixtures) of the mixtures that were present between 25 and 75% and between 20 and 80% in the viral population, as detected by 454. From the not called mixtures, all but one did present a mixture sequence in the Sanger DNA chromatograms, however, with a peak surface area for the second peak that was below the threshold setting for automatic mixture calling in the basecaller software (i.e. 25%). Viral loads ranged from 470 to 629,000 copies/mL and exerted no effect on the mixture calling relationship between both sequencing methodologies (R(2)=0.92). In some occasions (i.e. 55 mixtures) the genotyping assay would detect automatically mixtures that were present below 20% in the viral population, when measured by 454. Hence, the mixture scoring sensitivity of the automated high throughput virco(®)TYPE HIV-1 genotyping assay is currently set at 97.5% and 95.3%, for mixtures present at 25 and 20% in the viral population and may identify occasionally mutations that are present at lower frequencies. These findings were not influenced by the viral load of the examined samples.

HIV-1 dual/mixed tropic isolates show different genetic and phenotypic characteristics and response to maraviroc in vitro.

Dual/mixed-tropic HIV-1 strains are predominant in a significative proportion of patients, though few information is available regarding the genetic characteristics, quasispecies composition, and susceptibility against CCR5-antagonists of the primary-isolates. For this reason, we investigated in deep details, both phenotypically and genotypically, the characteristics of 54 HIV-1 primary-isolates obtained from HIV-infected patients. Tropism was assessed by multiple-cycles phenotypic-assay on U87MG-CD4(+)-CCR5(+)-/CXCR4(+)-expressing cells. In vitro selection in PBMCs of X4-tropic viral strains following maraviroc-treatment was also performed. Phenotypic-assay reported pure R5-tropic viruses in 31 (57.4%) isolates, dual/mixed-tropic viruses in 22 (40.7%), and pure X4-tropic virus in only 1 (1.8%). Among dual/mixed-tropic isolates, 12 showed a remarkably higher replication-efficacy in CCR5-expressing cells (R5(+)/X4), and 2 in CXCR4-expressing cells (R5/X4(+)). Genotypic-tropism testing showed a correlation between PSSM-scores, geno2pheno false-positive-rate, and V3-net-charge with both CCR5-usage and syncytium-inducing ability. Moreover, specific gp120- and gp41-mutations were significantly associated with tropism and/or syncytium-inducing ability. Ultra-deep V3-pyrosequencing showed the presence of a swarm of genetically distinct species with a preference for CCR5-coreceptor not only in all pure R5-isolates, but also in 6/7 R5(+)/X4-tropic isolates. In both pure-X4 and R5/X4(+)-isolates, we observed extensive prevalence of X4-using species. In vitro selection-experiments with CCR5-inhibitor maraviroc (up to 2 months) showed no-emergence of X4-tropic variants for all R5- and R5(+)/X4-isolates tested (while X4-virus remained fully-resistant). In conclusion, our study shows that dual/mixed-tropic viruses are constituted by different species, whereby those with characteristics R5(+)/X4 are genotypically and phenotypically similar to the pure-R5 isolates; thus the use of CCR5-antagonists in patients with R5(+)/X4-tropic viruses may be a therapeutic-option that deserves further investigations.

Comparison of phenotypic and genotypic tropism determination in triple-class-experienced HIV patients eligible for maraviroc treatment.

BACKGROUND: Determination of HIV-1 tropism is a pre-requisite to the use of CCR5 antagonists. This study evaluated the potential of population genotypic tropism tests (GTTs) in clinical practice, and the correlation with phenotypic tropism tests (PTTs) in patients accessing routine HIV care. METHODS: Forty-nine consecutive plasma samples for which an original Trofile(TM) assay was performed were obtained from triple-class-experienced patients in need of a therapy change. Viral tropism was defined as the consensus of three or more tropism calls obtained from the combination of two independent population PTT assays (Trofile Biosciences, San Francisco, CA, USA, and Virco, Beerse, Belgium), population GTTs and GTTs based on ultra-deep sequencing. If no consensus was reached, a clonal PTT was performed in order to finalize the tropism call. This two-step approach allowed the definition of a reference tropism call. RESULTS: According to the reference tropism result, 35/49 samples were CCR5 tropic (R5) (patients eligible for maraviroc treatment) and 14/49 were assigned as non-R5 tropic. The non-R5 samples [patients not eligible for maraviroc treatment according to the FDA/European Medicines Agency (EMEA) label] group included both the CXCR4 (X4) samples and the dual and mixed CCR5/CXCR4 (R5/X4) samples. Compared with Trofile(TM) population PTTs, population GTTs showed a higher sensitivity (97%) and a higher negative predictive value (91%), but almost equal specificity and an equal positive predictive value. CONCLUSIONS: In line with recent reports from clinical trial data, our data support the use of population genotypic tropism testing as a tool for tropism determination before the start of maraviroc.

HIV-1 V3 envelope deep sequencing for clinical plasma specimens failing in phenotypic tropism assays.

BACKGROUND: HIV-1 infected patients for whom standard gp160 phenotypic tropism testing failed are currently excluded from co-receptor antagonist treatment. To provide patients with maximal treatment options, massively parallel sequencing of the envelope V3 domain, in combination with tropism prediction tools, was evaluated as an alternative tropism determination strategy. Plasma samples from twelve HIV-1 infected individuals with failing phenotyping results were available. The samples were submitted to massive parallel sequencing and to confirmatory recombinant phenotyping using a fraction of the gp120 domain. RESULTS: A cut-off for sequence reads interpretation of 5 to10 times the sequencing error rate (0.2%) was implemented. On average, each sample contained 7 different V3 haplotypes. V3 haplotypes were submitted to tropism prediction algorithms, and 4/14 samples returned with presence of a dual/mixed (D/M) tropic virus, respectively at 3%, 10%, 11%, and 95% of the viral quasispecies. V3 tropism prediction was confirmed by gp120 phenotyping, except for two out of 4 D/M predicted viruses (with 3 and 95%) which were phenotypically R5-tropic. In the first case, the result was discordant due to the limit of detection for the phenotyping technology, while in the latter case the prediction algorithms were not computing the viral tropism correctly. CONCLUSIONS: Although only demonstrated on a limited set of samples, the potential of the combined use of "deep sequencing + prediction algorithms" in cases where routine gp160 phenotype testing cannot be employed was illustrated. While good concordance was observed between gp120 phenotyping and prediction of R5-tropic virus, the results suggest that accurate prediction of X4-tropic virus would require further algorithm development.

Tracking the evolution of multiple in vitro hepatitis C virus replicon variants under protease inhibitor selection pressure by 454 deep sequencing.

Resistance to hepatitis C virus (HCV) inhibitors targeting viral enzymes has been observed in in vitro replicon studies and during clinical trials. The factors determining the emergence of resistance and the changes in the viral quasispecies population under selective pressure are not fully understood. To assess the dynamics of variants emerging in vitro under various selective pressures with TMC380765, a potent macrocyclic HCV NS3/4A protease inhibitor, HCV genotype 1b replicon-containing cells were cultured in the presence of a low, high, or stepwise-increasing TMC380765 concentration(s). HCV replicon RNA from representative samples thus obtained was analyzed using (i) population, (ii) clonal, and (iii) 454 deep sequencing technologies. Depending on the concentration of TMC380765, distinct mutational patterns emerged. In particular, culturing with low concentrations resulted in the selection of low-level resistance mutations (F43S and A156G), whereas high concentrations resulted in the selection of high-level resistance mutations (A156V, D168V, and D168A). Clonal and 454 deep sequencing analysis of the replicon RNA allowed the identification of low-frequency preexisting mutations possibly contributing to the mutational pattern that emerged. Stepwise-increasing TMC380765 concentrations resulted in the emergence and disappearance of multiple replicon variants in response to the changing selection pressure. Moreover, two different codons for the wild-type amino acids were observed at certain NS3 positions within one population of replicons, which may contribute to the emerging mutational patterns. Deep sequencing technologies enabled the study of minority variants present in the HCV quasispecies population present at baseline and during antiviral drug pressure, giving new insights into the dynamics of resistance acquisition by HCV.

Resistance profile of darunavir: combined 24-week results from the POWER trials.

The resistance profile of darunavir (TMC114) in treatment-experienced patients was explored using pooled week 24 data from POWER 1, 2, and 3 at the recommended dose of darunavir with low-dose ritonavir (darunavir/r, 600/100 mg bid, N = 458). Baseline darunavir fold change in EC(50) was a strong predictor of virological response at week 24. Preliminary phenotypic clinical cut-offs of 10 and 40 were established. Virological response to darunavir/r was maintained in the presence at baseline of a high number of IAS-USA PI resistance-associated mutations (IAS-USA PI RAMS); a diminished response occurred with >or=14. Eleven protease mutations associated with diminished darunavir/r virological response were identified (V11I, V32I, L33F, I47V, I50V, I54L/M, G73S, L76V, I84V, and L89V). These darunavir resistance-associated mutations (DRV RAMS) occurred in the presence of a high number of IAS-USA PI RAMS. Virological response was diminished with three or more DRV RAMS in the background of a high number of IAS-USA PI RAMS. Incremental numbers of DRV RAMS were more predictive of outcome than were IAS-USA PI RAMS. Mutations developing during darunavir/r virological failure (V32I, L33F, I47V, I54L, and L89V) were also featured in the DRV RAMS list. Site-directed mutants carrying these five mutations, or any one of these mutations either alone or together with one or two IAS-USA PI RAMS, showed no reduced darunavir susceptibility, suggesting that a high number of additional background mutations is required for darunavir resistance. In this population of treatment-experienced patients, darunavir/r demonstrated significantly greater efficacy than investigator-selected control PIs of trials POWER 1 and 2, regardless of baseline viral genotype or phenotype, while exhibiting a high genetic barrier to the development of resistance.

Indolopyridones inhibit human immunodeficiency virus reverse transcriptase with a novel mechanism of action.

We have discovered a novel class of human immunodeficiency virus (HIV) reverse transcriptase (RT) inhibitors that block the polymerization reaction in a mode distinct from those of the nucleoside or nucleotide RT inhibitors (NRTIs) and nonnucleoside RT inhibitors (NNRTIs). For this class of indolopyridone compounds, steady-state kinetics revealed competitive inhibition with respect to the nucleotide substrate. Despite substantial structural differences with classical chain terminators or natural nucleotides, these data suggest that the nucleotide binding site of HIV RT may accommodate this novel class of RT inhibitors. To test this hypothesis, we have studied the mechanism of action of the prototype compound indolopyridone-1 (INDOPY-1) using a variety of complementary biochemical tools. Time course experiments with heteropolymeric templates showed "hot spots" for inhibition following the incorporation of pyrimidines (T>C). Moreover, binding studies and site-specific footprinting experiments revealed that INDOPY-1 traps the complex in the posttranslocational state, preventing binding and incorporation of the next complementary nucleotide. The novel mode of action translates into a unique resistance profile. While INDOPY-1 susceptibility is unaffected by mutations associated with NNRTI or multidrug NRTI resistance, mutations M184V and Y115F are associated with decreased susceptibility, and mutation K65R confers hypersusceptibility to INDOPY-1. This resistance profile provides additional evidence for active site binding. In conclusion, this class of indolopyridones can occupy the nucleotide binding site of HIV RT by forming a stable ternary complex whose stability is mainly dependent on the nature of the primer 3' end.