Characterization of virus‒host recombinant variants of the hepatitis E virus.

Hepatitis E virus is a single-strand, positive-sense RNA virus that can lead to chronic infection in immunocompromised patients. Virus-host recombinant variants (VHRVs) have been described in such patients. These variants integrate part of human genes into the polyproline-rich region that could introduce new post-translational modifications (PTMs), such as ubiquitination. The aim of this study was to characterize the replication capacity of different VHRVs, namely, RNF19A, ZNF787, KIF1B, EEF1A1, RNA18, RPS17, and RPL6. We used a plasmid encoding the Kernow strain, in which the fragment encoding the S17 insertion was deleted (Kernow p6 delS17) or replaced by fragments encoding the different insertions. The HEV RNA concentrations in the supernatants and the HepG2/C3A cell lysates were determined via RT-qPCR. The capsid protein ORF2 was immunostained. The effect of ribavirin was also assessed. The HEV RNA concentrations in the supernatants and the cell lysates were higher for the variants harboring the RNF19A, ZNF787, KIF1B, RPS17, and EEF1A1 insertions than for the Kernow p6 del S17, while it was not with RNA18 or RPL6 fragments. The number of ORF2 foci was higher for RNF19A, ZNF787, KIF1B, and RPS17 than for Kernow p6 del S17. VHRVs with replicative advantages were less sensitive to the antiviral effect of ribavirin. No difference in PTMs was found between VHRVs with a replicative advantage and those without. In conclusion, our study showed that insertions did not systematically confer a replicative advantage in vitro. Further studies are needed to determine the mechanisms underlying the differences in replicative capacity. IMPORTANCE: Hepatitis E virus (HEV) is a major cause of viral hepatitis. HEV can lead to chronic infection in immunocompromised patients. Ribavirin treatment is currently used to treat such chronic infections. Recently, seven virus-host recombinant viruses were characterized in immunocompromised patients. These viruses have incorporated a portion of a human gene fragment into their genome. We studied the consequences of these insertions on the replication capacity. We found that these inserted fragments could enhance virus replication for five of the seven recombinant variants. We also showed that the recombinant variants with replicative advantages were less sensitive to ribavirin in vitro. Finally, we found that the mechanisms leading to such a replicative advantage do not seem to rely on the post-translational modifications introduced by the human gene fragment that could have modified the function of the viral protein. The mechanisms involved in improving the replication of such recombinant viruses remain to be explored.

Comparison of short-read and long-read next-generation sequencing technologies for determining HIV-1 drug resistance.

Accurate HIV-1 genome sequencing is necessary to identify drug resistance mutations (DRMs) in people with HIV-1 (PWH). Next-generation-sequencing (NGS) allows the detection of minor variants and is now available in many laboratories. Our study aimed to compare two NGS approaches, a "short read" sequencing protocol using DeepChek® Whole Genome HIV-1 Assay on Illumina, and a "long read" sequencing protocol of HIV-1 pol and env single-molecule real-time sequencing (SMRT) on Pacific Biosciences (PacBio). We analyzed 16 plasma samples and 13 cellular samples from PWH. HIV-1 whole genome was amplified into five amplicons using DeepChek® Whole Genome HIV-1 Assay and sequenced on an iSeq. 100. In parallel, HIV-1 pol and env genes were separately amplified and sequenced using PacBio SMRT system with the circular consensus sequencing mode on a Sequel IIe. Concordance rates for determining DRMs with both approaches varied depending on the HIV-1 region, with higher concordance in the integrase region compared to the reverse transcriptase and protease regions. DeepChek® Whole Genome HIV-1 Assay exhibited better sensitivity in HIV-1 RNA sequencing of plasmas with lower viral loads. In cell HIV-1 DNA sequencing, the DeepChek® Whole Genome HIV-1 Assay performed better in pol and env sequencing but detected more APOBEC-induced DRMs, which can represent defective proviruses. Our findings indicate that both DeepChek® Whole Genome HIV-1 Assay and PacBio SMRT sequencing exhibit good performance for subtype determination, detection, and quantification of DRMs of the HIV-1 genome. However, some discrepancies were found in cellular samples, highlighting the challenges of interpreting HIV-1 DNA DRMs.

Influence of Nasopharyngeal Viral Load on the Spread of the Omicron BA.2 Variant.

We used variant typing polymerase chain reaction to describe the evolution of severe acute respiratory syndrome coronavirus 2 Omicron sublineages between December 2021 and mid-March 2022. The selective advantage of the BA.2 variant over BA.1 is not due to greater nasopharyngeal viral loads.

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.

Whole-genome sequencing of SARS-CoV-2: Comparison of target capture and amplicon single molecule real-time sequencing protocols.

Fast, accurate sequencing methods are needed to identify new variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome. Single-molecule real-time (SMRT) Pacific Biosciences (PacBio) provides long, highly accurate sequences by circular consensus reads. This study compares the performance of a target capture SMRT PacBio protocol for whole-genome sequencing (WGS) of SARS-CoV-2 to that of an amplicon PacBio SMRT sequencing protocol. The median genome coverage was higher (p < 0.05) with the target capture protocol (99.3% [interquartile range, IQR: 96.3-99.5]) than with the amplicon protocol (99.3% [IQR: 69.9-99.3]). The clades of 65 samples determined with both protocols were 100% concordant. After adjusting for Ct values, S gene coverage was higher with the target capture protocol than with the amplicon protocol. After stratification on Ct values, higher S gene coverage with the target capture protocol was observed only for samples with Ct > 17 (p < 0.01). PacBio SMRT sequencing protocols appear to be suitable for WGS, genotyping, and detecting mutations of SARS-CoV-2.

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.

Evolution of HIV-1 quasispecies and coreceptor use in cell reservoirs of patients on suppressive antiretroviral therapy.

OBJECTIVES: To track changes in the V3 env region of HIV-1 quasispecies and determine virus coreceptor use in cell reservoirs of patients on long-term suppressive antiretroviral therapy (ART). PATIENTS AND METHODS: Ten patients whose plasma viraemia had been suppressed for a median of 5.5 years were followed for 5 years. The V3 env regions of viruses in peripheral blood mononuclear cells were analysed by ultra-deep sequencing (UDS). HIV-1 tropism was predicted using the geno2pheno 5.75 algorithm and a phenotypic assay. RESULTS: The UDS and phenotypic assay data were concordant for predicting HIV-1 tropism. CXCR4-using viruses detected by UDS accounted for 14.7%-76.5% of the virus populations in samples from five patients at enrolment. Five patients harboured pure R5 virus populations and no X4 viruses emerged during the 5 years. The selection pressures estimated by the dN/dS ratio were acting on the V3 region to produce diversification of the quasispecies in CXCR4-infected patients and purification of the quasispecies in R5-infected patients on effective ART. CONCLUSIONS: UDS showed that the virus quasispecies in cell reservoirs of patients on long-term suppressive ART continued to evolve. CXCR4-using variants became more diversified. Analysis of the selection pressures on the virus quasispecies could provide a clearer picture of virus persistence in patients on effective ART.

HIV-1 genotypic resistance testing using single molecule real-time sequencing.

BACKGROUND: HIV-1 resistance testing is recommended in clinical management and next-generation sequencing (NGS) methods are now available in many virology laboratories. OBJECTIVES: To evaluate the diagnostic performance of Long-Read Single Molecule Real-time (SMRT) sequencing (Sequel, PacBio) for HIV-1 polymerase genotyping. STUDY DESIGN: 111 prospective clinical samples (83 plasma and 28 leukocyte-enriched blood fraction) were analyzed for routine HIV-1 resistance genotyping using Sanger sequencing, Vela NGS, and SMRT sequencing. We developed a SMRT sequencing protocol and a bio-informatics pipeline to infer antiretroviral resistance on both haplotype and variant calling approaches. RESULTS: The polymerase was successfully sequenced by the three platforms in 98 % of plasma RNA samples for viral loads above 4 log copies/mL. The success rate decreased to 83 % using Sanger or Vela sequencing and to 67 % using SMRT sequencing for viral loads of 3 to 4 log copies/mL. Sensitivities of 50 %, 54 % and 61 % were obtained using SMRT, Vela, and Sanger sequencing, respectively, in cellular DNA from patients with prolonged undetectable plasma HIV-1 RNA. Ninety-eight percent of resistance-associated mutations (RAMs) identified with Sanger sequencing were detected using SMRT sequencing. Furthermore, 91 % of RAMs (> 5 % threshold) identified with Vela NGS were detected using SMRT sequencing. RAM quantification using Vela and SMRT sequencing was well correlated (Spearman correlation ρ = 0.82; P < 0.0001). CONCLUSIONS: SMRT sequencing of the full-length HIV-1 polymerase appeared performant for characterizing HIV-1 genotypic resistance on both RNA and DNA clinical samples. Long-read sequencing is a new tool for mutation haplotyping and resistance analysis.

Intact proviruses are enriched in the colon and associated with PD-1+TIGIT- mucosal CD4+ T cells of people with HIV-1 on antiretroviral therapy.

BACKGROUND: The persistence of intact replication-competent HIV-1 proviruses is responsible for the virological rebound off treatment. The gut could be a major reservoir of HIV-1 due to the high number of infected target cells. METHODS: We collected blood samples and intestinal biopsies (duodenum, ileum, colon) from 42 people with HIV-1 receiving effective antiretroviral therapy. We used the Intact Proviral DNA Assay to estimate the frequency of intact HIV-1 proviruses in the blood and in the intestinal mucosa of these individuals. We analyzed the genetic complexity of the HIV-1 reservoir by performing single-molecule next-generation sequencing of HIV-1 env DNA. The activation/exhaustion profile of mucosal T lymphocytes was assessed by flow cytometry. FINDINGS: Intact proviruses are particularly enriched in the colon. Residual HIV-1 transcription in the gut is associated with persistent mucosal and systemic immune activation. The HIV-1 intestinal reservoir appears to be shaped by the proliferation of provirus-hosting cells. The genetic complexity of the viral reservoir in the colon is positively associated with TIGIT expression but negatively with PD-1, and inversely related to its intact content. The size of the intact reservoir in the colon is associated with PD-1+TIGIT- mucosal CD4+ T cells, particularly in CD27+ memory cells, whose proliferation and survival could contribute to the enrichment of the viral reservoir by intact proviruses. INTERPRETATION: Enrichment in intact proviruses makes the gut a key compartment for HIV-1 persistence on antiretroviral therapy. FUNDING: This project was supported by grants from the ANRS-MIE (ANRS EP61 GALT), Sidaction, and the Institut Universitaire de France.

Impact of Casirivimab-Imdevimab on Severe Acute Respiratory Syndrome Coronavirus 2 Delta Variant Nasopharyngeal Virus Load and Spike Quasispecies.

Background: The increasing use of monoclonal antibodies (mAbs) to treat coronavirus disease 2019 raises questions about their impact on the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mAb-resistant variants. We assessed the impact of Casirivimab-Imdevimab on SARS-CoV-2 mutations associated with reduced mAb activity in treated patients. Methods: We measured the nasopharyngeal (NP) viral load and sequenced the haplotypes of spike gene of 50 patients infected with the SARS-CoV-2 delta variant and treated with Casirivimab-Imdevimab using single-molecule real-time sequencing. Results: The NP SARS-CoV-2 viral load of patients treated with Casirivimab-Imdevimab decreased from 8.13 (interquartile range [IQR], 7.06-8.59) log10 copies/mL pretreatment to 3.67 (IQR, 3.07-5.15) log10 copies/mL 7 days later (P < .001). Of the 36 patients for whom follow-up timepoints Spike sequencing were available, none of the Spike mutations that reduced mAb activity were detected. Conclusions: Casirivimab-Imdevimab is an effective treatment for patients infected with the SARS-CoV-2 delta variant. Despite selective pressure on SARS-CoV-2 Spike quasispecies, we detected no key mutations that reduced mAb activity in our patients.

Influence of treatment with neutralizing monoclonal antibodies on the SARS-CoV-2 nasopharyngeal load and quasispecies.

OBJECTIVES: To evaluate the impact of neutralizing monoclonal antibody (mAb) treatment and to determine whether the selective pressure of mAbs could facilitate the proliferation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with spike protein mutations that might attenuate mAb effectiveness. PATIENTS AND METHODS: We evaluated the impact of mAbs on the nasopharyngeal (NP) viral load and virus quasispecies of mAb-treated patients using single-molecule real-time sequencing. The mAbs used were: Bamlanivimab alone (four patients), Bamlanivimab/Etesevimab (23 patients) and Casirivimab/Imdevimab (five patients). RESULTS: The NP SARS-CoV-2 viral load of mAb-treated patients decreased from 8.2 log10 copies/mL before administration to 4.3 log10 copies/mL 7 days after administration. Five immunocompromised patients given Bamlanivimab/Etesevimab were found to have mAb activity-reducing spike mutations. Two patients harboured SARS-CoV-2 variants with a Q493R spike mutation 7 days after administration, as did a third patient 14 days after administration. The fourth patient harboured a variant with a Q493K spike mutation 7 days post-treatment, and the fifth patient had a variant with a E484K spike mutation on day 21. The emergence of the spike mutation was accompanied by stabilization or rebound of the NP viral load in three of five patients. CONCLUSION: Two-mAb therapy can drive the selection of resistant SARS-CoV-2 variants in immunocompromised patients. Patients given mAbs should be closely monitored and measures to limit virus spread should be reinforced.

Prediction of SARS-CoV-2 Variant Lineages Using the S1-Encoding Region Sequence Obtained by PacBio Single-Molecule Real-Time Sequencing.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is the causal agent of the COVID-19 pandemic that emerged in late 2019. The outbreak of variants with mutations in the region encoding the spike protein S1 sub-unit that can make them more resistant to neutralizing or monoclonal antibodies is the main point of the current monitoring. This study examines the feasibility of predicting the variant lineage and monitoring the appearance of reported mutations by sequencing only the region encoding the S1 domain by Pacific Bioscience Single Molecule Real-Time sequencing (PacBio SMRT). Using the PacBio SMRT system, we successfully sequenced 186 of the 200 samples previously sequenced with the Illumina COVIDSeq (whole genome) system. PacBio SMRT detected mutations in the S1 domain that were missed by the COVIDseq system in 27/186 samples (14.5%), due to amplification failure. These missing positions included mutations that are decisive for lineage assignation, such as G142D (n = 11), N501Y (n = 6), or E484K (n = 2). The lineage of 172/186 (92.5%) samples was accurately determined by analyzing the region encoding the S1 domain with a pipeline that uses key positions in S1. Thus, the PacBio SMRT protocol is appropriate for determining virus lineages and detecting key mutations.

Influence of SARS-CoV-2 Variant B.1.1.7, Vaccination, and Public Health Measures on the Spread of SARS-CoV-2.

The spread of SARS-CoV-2 and the resulting disease COVID-19 has killed over 2.6 million people as of 18 March 2021. We have used a modified susceptible, infected, recovered (SIR) epidemiological model to predict how the spread of the virus in regions of France will vary depending on the proportions of variants and on the public health strategies adopted, including anti-COVID-19 vaccination. The proportion of SARS-CoV-2 variant B.1.1.7, which was not detected in early January, increased to become 60% of the forms of SARS-CoV-2 circulating in the Toulouse urban area at the beginning of February 2021, but there was no increase in positive nucleic acid tests. Our prediction model indicates that maintaining public health measures and accelerating vaccination are efficient strategies for the sustained control of SARS-CoV-2.

Classification of the Zoonotic Hepatitis E Virus Genotype 3 Into Distinct Subgenotypes.

Hepatitis E virus (HEV) genotype 3 is the most common genotype linked to HEV infections in Europe and America. Three major clades (HEV-3.1, HEV-3.2, and HEV-3.3) have been identified but the overlaps between intra-subtype and inter-subtype p-distances make subtype classification inconsistent. Reference sequences have been proposed to facilitate communication between researchers and new putative subtypes have been identified recently. We have used the full or near full-length HEV-3 genome sequences available in the Genbank database (April 2020; n = 503) and distance analyses of clades HEV-3.1 and HEV-3.2 to determine a p-distance cut-off (0.093 nt substitutions/site) in order to define subtypes. This could help to harmonize HEV-3 genotyping, facilitate molecular epidemiology studies and investigations of the biological and clinical differences between HEV-3 subtypes.