HPV genotyping in clinical samples using long-read single-molecule real-time sequencing.

Human papillomavirus (HPV) genotyping is widely used, particularly in combination with high-risk (HR) HPV tests for cervical cancer screening. We developed a genotyping method using sequences of approximately 800 bp in the E6/E7 region obtained by PacBio single molecule real-time sequencing (SMRT) and evaluated its performance against MY09-11 L1 sequencing and after the APTIMA HPV genotyping assay. The levels of concordance of PacBio E6/E7 SMRT sequencing with MY09-11 L1 sequencing and APTIMA HPV genotyping were 100% and 90.8%, respectively. The sensitivity of PacBio E6/EA7 SMRT was slightly greater than that of L1 sequencing and, as expected, lower than that of HR-HPV tests. In the context of cervical cancer screening, PacBio E6/E7 SMRT is then best used after a positive HPV test. PacBio E6/E7 SMRT genotyping is an attractive alternative for HR and LR-HPV genotyping of clinical samples. PacBio SMRT sequencing provides unbiased genotyping and can detect multiple HPV infections and haplotypes within a genotype.

HIV-1 resistance genotyping by ultra-deep sequencing and 6-month virological response to first-line treatment.

OBJECTIVES: To evaluate the routine use of the Sentosa ultra-deep sequencing (UDS) system for HIV-1 polymerase resistance genotyping in treatment-naïve individuals and to analyse the virological response (VR) to first-line antiretroviral treatment. METHODS: HIV drug resistance was determined on 237 consecutive samples from treatment-naïve individuals using the Sentosa UDS platform with two mutation detection thresholds (3% and 20%). VR was defined as a plasma HIV-1 virus load <50 copies/mL after 6 months of treatment. RESULTS: Resistance to at least one antiretroviral drug with a mutation threshold of 3% was identified in 29% and 16% of samples according to ANRS and Stanford algorithms, respectively. The ANRS algorithm also revealed reduced susceptibility to at least one protease inhibitor (PI) in 14.3% of samples, to one reverse transcriptase inhibitor in 12.7%, and to one integrase inhibitor (INSTI) in 5.1%. For a mutation threshold of 20%, resistance was identified in 24% and 13% of samples according to ANRS and Stanford algorithms, respectively. The 6 months VR was 87% and was similar in the 58% of patients given INSTI-based treatment, in the 16% given PI-based treatment and in the 9% given NNRTI-based treatment. Multivariate analysis indicated that the VR was correlated with the baseline HIV virus load and resistance to at least one PI at both 3% and 20% mutation detection thresholds (ANRS algorithm). CONCLUSIONS: The Vela UDS platform is appropriate for determining antiretroviral resistance in patients on a first-line antiretroviral treatment. Further studies are needed on the use of UDS for therapeutic management.

Whole-genome single molecule real-time sequencing of SARS-CoV-2 Omicron.

New variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome can only be identified using accurate sequencing methods. Single molecule real-time (SMRT) sequencing has been used to characterize Alpha and Delta variants, but not Omicron variants harboring numerous mutations in the SARS-CoV-2 genome. This study assesses the performance of a target capture SMRT sequencing protocol for whole genome sequencing (WGS) of SARS-CoV-2 Omicron variants and compared it to that of an amplicon SMRT sequencing protocol optimized for Omicron variants. The failure rate of the target capture protocol (6%) was lower than that of the amplicon protocol (34%, p < 0.001) on our data set, and the median genome coverage with the target capture protocol (98.6% [interquartile range (IQR): 86-99.4]) was greater than that with the amplicon protocol (76.6% [IQR: 66-89.6], [p < 0.001]). The percentages of samples with >95% whole genome coverage were 64% with the target capture protocol and 19% with the amplicon protocol (p < 0.05). The clades of 96 samples determined with both protocols were 93% concordant and the lineages of 59 samples were 100% concordant. Thus, target capture SMRT sequencing appears to be an efficient method for WGS, genotyping and detecting mutations of SARS-CoV-2 Omicron variants.

THETA: a new genotypic approach for predicting HIV-1 CRF02-AG coreceptor usage.

MOTIVATION: The circulating recombinant form of HIV-1 CRF02-AG is the most frequent non-B subtype in Europe. Anti-HIV therapy and pathophysiological studies on the impact of HIV-1 tropism require genotypic determination of HIV-1 tropism for non-B subtypes. But genotypic approaches based on analysis of the V3 envelope region perform poorly when used to determine the tropism of CRF02-AG. We, therefore, designed an algorithm based on information from the gp120 and gp41 ectodomain that better predicts the tropism of HIV-1 subtype CRF02-AG. RESULTS: We used a bio-statistical method to identify the genotypic determinants of CRF02-AG coreceptor use. Toulouse HIV Extended Tropism Algorithm (THETA), based on a Least Absolute Shrinkage and Selection Operator method, uses HIV envelope sequence from phenotypically characterized clones. Prediction of R5X4/X4 viruses was 86% sensitive and that of R5 viruses was 89% specific with our model. The overall accuracy of THETA was 88%, making it sufficiently reliable for predicting the tropism of subtype CRF02-AG sequences. AVAILABILITY AND IMPLEMENTATION: Binaries are freely available for download at https://github.com/viro-tls/THETA. It was implemented in Matlab and supported on MS Windows platform. The sequence data used in this work are available from GenBank under the accession numbers MK618182-MK618417.

Impact of the mutational load on the virological response to a first-line rilpivirine-based regimen.

OBJECTIVES: To determine how the load of rilpivirine-resistant variants (mutational load) influences the virological response (VR) of HIV-1-infected patients to a rilpivirine-based first-line regimen. PATIENTS AND METHODS: Four hundred and eighty-nine patients infected with HIV-1 whose reverse transcriptase gene had been successfully resistance genotyped using next-generation sequencing were given a first-line regimen containing rilpivirine. Variables associated with the VR at 12 months were identified using a logistic model. The results were used to build a multivariate model for each mutational load threshold and the R2 variations were analysed to identify the mutational load threshold that best predicted the VR. RESULTS: The mutational load at baseline was the only variable linked to the VR at 12 months (P  < 0.01). The VR at 12 months decreased from 96.9% to 83.4% when the mutational load was >1700 copies/mL and to 50% when the mutational load was > 9000 copies/mL. The threshold of 9000 copies/mL was associated with the VR at 12 months with an OR of 36.7 (95% CI 4.7-285.1). The threshold of 1700 copies/mL was associated with the VR at 12 months with an OR of 7.2 (95% CI 1.4-36.8). CONCLUSIONS: There is quantifiable evidence that determining a mutational load threshold can be used to identify those patients on a first-line regimen containing rilpivirine who are at risk of virological failure. The clinical management of HIV-infected patients can be improved by evaluating the frequency of mutant variants at a threshold of < 20% together with the plasma HIV-1 viral load at the time of resistance genotyping.

Diversity of hepatitis E virus genotype 3.

Hepatitis E virus genotype 3 (HEV-3) can lead to chronic infection in immunocompromised patients, and ribavirin is the treatment of choice. Recently, mutations in the polymerase gene have been associated with ribavirin failure but their frequency before treatment according to HEV-3 subtypes has not been studied on a large data set. We used single-molecule real-time sequencing technology to sequence 115 new complete genomes of HEV-3 infecting French patients. We analyzed phylogenetic relationships, the length of the polyproline region, and mutations in the HEV polymerase gene. Eighty-five (74%) were in the clade HEV-3efg, 28 (24%) in HEV-3chi clade, and 2 (2%) in HEV-3ra clade. Using automated partitioning of maximum likelihood phylogenetic trees, complete genomes were classified into subtypes. Polyproline region length differs within HEV-3 clades (from 189 to 315 nt). Investigating mutations in the polymerase gene, distinct polymorphisms between HEV-3 subtypes were found (G1634R in 95% of HEV-3e, G1634K in 56% of HEV-3ra, and V1479I in all HEV-3efg, clade HEV-3ra, and HEV-3k strains). Subtype-specific polymorphisms in the HEV-3 polymerase have been identified. Our study provides new complete genome sequences of HEV-3 that could be useful for comparing strains circulating in humans and the animal reservoir.

HIV-1 genotypic resistance testing using the Vela automated next-generation sequencing platform.

Objectives: To evaluate the diagnostic performance of the Vela next-generation sequencing (NGS) system in conjunction with the Sentosa SQ HIV Genotyping Assay for genotyping HIV-1. Methods: Plasma RNA was extracted and templates prepared with the Sentosa SX instrument before sequencing the HIV-1 polymerase on the Sentosa SQ301 Sequencer (PGM IonTorrent). The Vela NGS System was compared with direct sequencing and the 454 GS-FLX (Roche) and MiSeq (Illumina) systems for genotypic resistance testing on clinical samples. Results: The Vela NGS system detected majority resistance mutations in subtype B and CRF02-AG samples at 500 copies/mL and minority variants with a sensitivity of 5% at 100 000 copies/mL. The Vela NGS system and direct sequencing identified resistance mutations with 97% concordance in 46 clinical samples. Vela identified 1/20 of the 1%-5% mutations identified by 454, 5/12 of the 5%-20% mutations and 60/61 of the >20% mutations. Vela identified 3/14 of the 1%-5% mutations identified by MiSeq, 0/2 of the 5%-20% mutations and 47/47 of the >20% mutations. The resistance mutation quantifications by Vela and 454 were concordant (bias: 2.31%), as were those by Vela and MiSeq (bias: 1.06%). Conclusions: The Vela NGS system provides automated nucleic acid extraction, PCR reagent distribution, library preparation and bioinformatics analysis. The analytical performance was very good when compared with direct sequencing, but was less sensitive than two other NGS platforms for detecting minority variants.

Performance comparison of deep sequencing platforms for detecting HIV-1 variants in the pol gene.

The present study compares the performances of an in-house sequencing protocol developed on MiSeq, the Sanger method, and the 454 GS-FLX for detecting and quantifying drug-resistant mutations (DRMs) in the human immunodeficiency virus polymerase gene (reverse transcriptase [RT] and protease [PR]). MiSeq sequencing identified all the resistance mutations detected by bulk sequencing (n = 84). Both the MiSeq and 454 GS-FLX platforms identified 67 DRMs in the RT and PR regions, but a further 25 DRMs were identified by only one or other of them. Pearson's analysis showed good concordance between the percentage of drug-resistant variants determined by MiSeq and 454 GS-FLX sequencing (ρ = .77, P < .0001). The MiSeq platform is as accurate as the 454 GS-FLX Roche system for determining RT and PR DRMs and could be used for monitoring human immunodeficiency virus type 1 drug resistance.

Mutation in the Hepatitis E Virus Polymerase and Outcome of Ribavirin Therapy.

Hepatitis E virus (HEV) can lead to chronic infection in solid-organ transplant patients. Ribavirin is efficient for treatment of chronically infected patients. Recently, the1634R mutation in the HEV polymerase has been associated with treatment failure. However, it is unclear if this mutation can be used as a prognostic marker of treatment outcome. We studied the prevalence of the 1634R mutation in the HEV polymerase of patients starting ribavirin therapy, the influence of the 1634R variants on the viral response, the frequency of the 1634R mutation in patients whose treatment failed, and its impact on ribavirin retreatment. We analyzed pretreatment samples from 63 solid-organ transplant patients with chronic hepatitis E using deep sequencing; 42 patients had a sustained virologic response (SVR), and 21 were non-SVR patients. We detected the 1634R variant by deep sequencing in 36.5% (23/63) of the patients (proportions, 1.3 to 100%). The 1634R variant was detected in 31.0% (13/42) of baseline plasma samples from patients with SVR and in 47.6% (10/21) in the other patients (P = 0.2). The presence of this mutation did not influence the initial decrease in viral RNA. Lastly, a second prolonged ribavirin treatment led to SVR in 70% of the patients who initially did not have SVR, despite the presence of the 1634R variant. We conclude that the presence of the 1634R variant at ribavirin initiation does not lead to absolute ribavirin resistance. Although its proportion increased in patients whose treatment failed, the presence of the 1634R variant did not compromise the response to a second ribavirin treatment.