Search

The repeated failure to treat patients chronically infected with hepatitis E (HEV) and C (HCV) viruses, despite the absence of resistance-associated substitutions (RAS), particularly in response to prolonged treatments with the mutagenic agents of HEV, suggests that quasispecies structure may play a crucial role beyond single point mutations. Quasispecies structured in a flat-like manner (referred to as flat-like) are considered to possess high average fitness, occupy a significant fraction of the functional genetic space of the virus, and exhibit a high capacity to evade specific or mutagenic treatments. In this paper, we studied HEV and HCV samples using high-depth next-generation sequencing (NGS), with indices scoring the different properties describing flat-like quasispecies. The significance of these indices was demonstrated by comparing the values obtained from these samples with those from acute infections caused by respiratory viruses (betacoronaviruses, enterovirus, respiratory syncytial viruses, and metapneumovirus). Our results revealed that flat-like quasispecies in HEV and HCV chronic infections without RAS are characterized by numerous low-frequency haplotypes with no dominant one. Surprisingly, these low-frequency haplotypes (at the nucleotide level) exhibited a high level of synonymity, resulting in much lower diversity at the phenotypic level. Currently, clinical approaches for managing flat-like quasispecies are lacking. Here, we propose methods to identifying flat-like quasispecies, which represents an essential initial step towards exploring alternative treatment protocols for viruses resistant to conventional therapies.

Phage therapy is still in its infancy, but it is increasingly promising as a future alternative for treating antibiotic-resistant bacteria. To investigate the effect of phages on Mycobacterium abscessus complex (MABC), we isolated 113 environmental phages, grown them to high titres, and assayed them on MABC clinical strains through the spot test. Of all the phages, only 16 showed killing activity. Their activity was so temperate to MABC that they could not generate any plaque-forming units (PFUs). The Appelmans method of directed evolution was carried out to evolve these 16 phages into more lytic ones. After only 11 of 30 rounds of evolution, every single clinical strain in our collection, including those that were unsusceptible up to this point, could be lysed by at least one phage. The evolved phages were able to form PFUs on the clinical strains tested. Still, they are temperate at best and require further training. The genomes of one random parental phage and three random evolved phages from Round 13 were sequenced, revealing a diversity of clusters and genes of a variety of evolutionary origins, mostly of unknown function. These complete annotated genomes will be key for future molecular characterisations.

Here, we report the in-host hepatitis E virus (HEV) quasispecies evolution in a chronically infected patient who was treated with three different regimens of ribavirin (RBV) for nearly 6 years. Sequential plasma samples were collected at different time points and subjected to RNA extraction and deep sequencing using the MiSeq Illumina platforms. Specifically, we RT-PCR amplified a single amplicon from the core region located in the open-reading frame 2 (ORF2). At the nucleotide level (genotype), our analysis showed an increase in the number of rare haplotypes and a drastic reduction in the frequency of the master (most represented) sequence during the period when the virus was found to be insensitive to RBV treatment. Contrarily, at the amino acid level (phenotype), our study revealed conservation of the amino acids, which is represented by a high prevalence of the master sequence. Our findings suggest that using mutagenic antivirals concomitant with high viral loads can lead to the selection and proliferation of a rich set of synonymous haplotypes that express the same phenotype. This can also lead to the selection and proliferation of conservative substitutions that express fitness-enhanced phenotypes. These results have important clinical implications, as they suggest that using mutagenic agents as a monotherapy treatment regimen in the absence of sufficiently effective viral inhibitors can result in diversification and proliferation of a highly diverse quasispecies resistant to further treatment. Therefore, such approaches should be avoided whenever possible.

Shortly after the beginning of the SARS-CoV-2 pandemic, many countries implemented sewage sentinel systems to monitor the circulation of the virus in the population. A fundamental part of these surveillance programs is the variant tracking through sequencing approaches to monitor and identify new variants or mutations that may be of importance. Two of the main sequencing platforms are Illumina and Oxford Nanopore Technologies. Here, we compare the performance of MiSeq (Illumina) and MinION (Oxford Nanopore Technologies), as well as two different data processing pipelines, to determine the effect they may have on the results. MiSeq showed higher sequencing coverage, lower error rate, and better capacity to detect and accurately estimate variant abundances than MinION R9.4.1 flow cell data. The use of different variant callers (LoFreq and iVar) and approaches to calculate the variant proportions had a remarkable impact on the results generated from wastewater samples. Freyja, coupled with iVar, may be more sensitive and accurate than LoFreq, especially with MinION data, but it comes at the cost of having a higher error rate. The analysis of MinION R10.4.1 flow cell data using Freyja combined with iVar narrows the gap with MiSeq performance in terms of read quality, accuracy, sensitivity, and number of detected mutations. Although MiSeq should still be considered as the standard method for SARS-CoV-2 variant tracking, MinION's versatility and rapid turnaround time may represent a clear advantage during the ongoing pandemic.

Epidemics and pandemics have occurred since the beginning of time, resulting in millions of deaths. Many such disease outbreaks are caused by viruses. Some viruses, particularly RNA viruses, are characterized by their high genetic variability, and this can affect certain phenotypic features: tropism, antigenicity, and susceptibility to antiviral drugs, vaccines, and the host immune response. The best strategy to face the emergence of new infectious genomes is prompt identification. However, currently available diagnostic tests are often limited for detecting new agents. High-throughput next-generation sequencing technologies based on metagenomics may be the solution to detect new infectious genomes and properly diagnose certain diseases. Metagenomic techniques enable the identification and characterization of disease-causing agents, but they require a large amount of genetic material and involve complex bioinformatic analyses. A wide variety of analytical tools can be used in the quality control and pre-processing of metagenomic data, filtering of untargeted sequences, assembly and quality control of reads, and taxonomic profiling of sequences to identify new viruses and ones that have been sequenced and uploaded to dedicated databases. Although there have been huge advances in the field of metagenomics, there is still a lack of consensus about which of the various approaches should be used for specific data analysis tasks. In this review, we provide some background on the study of viral infections, describe the contribution of metagenomics to this field, and place special emphasis on the bioinformatic tools (with their capabilities and limitations) available for use in metagenomic analyses of viral pathogens.

The SARS-CoV-2 Omicron variant emerged showing higher transmissibility and possibly higher resistance to current COVID-19 vaccines than other variants dominating the global pandemic. In March 2020 we performed a study in clinical samples, where we found that a portion of genomes in the SARS-CoV-2 viral population accumulated deletions immediately before the S1/S2 cleavage site (furin-like cleavage site, PRRAR/S) of the spike gene, generating a frameshift and appearance of a premature stop codon. The main aim of this study was to determine the frequency of defective deletions in prevalent variants from the first to sixth pandemic waves in our setting and discuss whether the differences observed might support epidemiological proposals. The complete SARS-CoV-2 spike gene was deeply studied by next-generation sequencing using the MiSeq platform. More than 90 million reads were obtained from respiratory swab specimens of 78 COVID-19 patients with mild infection caused by the predominant variants circulating in the Barcelona city area during the six pandemic waves: B.1.5, B.1.1, B.1.177, Alpha, Beta, Delta, and Omicron. The frequency of defective genomes found in variants dominating the first and second waves was similar to that seen in Omicron, but differed from the frequencies seen in the Alpha, Beta and Delta variants. The changing pattern of mutations seen in the various SARS-CoV-2 variants driving the pandemic waves over time can affect viral transmission and immune escape. Here we discuss the putative biological effects of defective deletions naturally occurring before the S1/S2 cleavage site during adaption of the virus to human infection., (© 2022. The Author(s).)

The changes occurring in viral quasispecies populations during infection have been monitored using diversity indices, nucleotide diversity, and several other indices to summarize the quasispecies structure in a single value. In this study, we present a method to partition quasispecies haplotypes into four fractions according to their fitness: the master haplotype, rare haplotypes at two levels (those present at <0.1%, and those at 0.1−1%), and a fourth fraction that we term emerging haplotypes, present at frequencies >1%, but less than that of the master haplotype. We propose that by determining the changes occurring in the volume of the four quasispecies fitness fractions together with those of the Hill number profile we will be able to visualize and analyze the molecular changes in the composition of a quasispecies with time. To develop this concept, we used three data sets: a technical clone of the complete SARS-CoV-2 spike gene, a subset of data previously used in a study of rare haplotypes, and data from a clinical follow-up study of a patient chronically infected with HEV and treated with ribavirin. The viral response to ribavirin mutagenic treatment was selection of a rich set of synonymous haplotypes. The mutation spectrum was very complex at the nucleotide level, but at the protein (phenotypic/functional) level the pattern differed, showing a highly prevalent master phenotype. We discuss the putative implications of this observation in relation to mutagenic antiviral treatment.

Background: Hepatitis C virus (HCV) is a major cause of chronic liver infection with 71 million people infected worldwide. Pakistan has the second highest prevalence of HCV infection and more than half (52%) of Pakistani living in Spain reside in Barcelona. The aim of this study was to analyse the seroprevalence and viraemic rate and determine the genotypes and subtypes of HCV among Pakistanis living in the southern metropolitan area of Barcelona., Methods: We included all Pakistani patients seeking primary healthcare in the southern metropolitan area of Barcelona from August 2011 to July 2014. Serum samples were screened for HCV antibodies. HCV viral load was determined by reverse transcription polymerase chain reaction and genotypes and subtypes were performed using Versant HCV Genotype and/or deep-sequencing. Screening for hepatitis B virus (HBV) was also carried out., Results: Among 5877 Pakistani patients, 565 (9.61%) were screened for anti-HCV antibodies, with 68 (12.04%) being positive. The viral load was determined in 65, with 31 presenting active infection and the viraemic rate was 47.69% (95% confidence interval 36.02-59.62). HCV genotyping and subtyping were performed in 24 individuals. Most infections corresponded to HCV genotype 3 (91.67%), and high resolution HCV subtyping was performed in 18 samples, 16 of which presented subtype 3a. One subject presented HBV coinfection with undetectable HBV DNA. During the study period, we identified a possible case of HCV vertical transmission followed by spontaneous viraemia clearance in a chronically infected mother with a C/T IL28B genetic polymorphism., Conclusion: These results suggest that general HCV screening protocols in patients from high prevalence countries, such as Pakistan, would be helpful to identify and treat active HCV infections. This could avoid further transmission and contribute to building targeted health policies for micro-elimination of HCV infection in specific communities., Competing Interests: The authors report no conflicts of interest in this work., (© 2022 Dopico et al.)

Wastewater based epidemiology (WBE) offers an overview of the SARS-CoV-2 variants circulating among the population thereby serving as a proper surveillance method. The variant of concern (VOC) Alpha was first identified in September 2020 in the United Kingdom, and rapidly became dominant across Europe. Our objective was to elucidate the Alpha VOC outcompetition rate and identify mutations in the spike glycoprotein (S) gene, indicative of the circulation of the Alpha VOC and/or other variants in the population through wastewater analysis. In the period covered by this study (November 2020-April 2021), forteen wastewater treatment plants (WWTPs) were weekly sampled. The total number of SARS-CoV-2 genome copies per L (GC/L) was determined with a Real-Time qPCR, targeting the N gene. Surveillance of the Alpha VOC circulation was ascertained using a duplex RT-qPCR, targeting and discriminating the S gene. Our results showed that in a period of 6 weeks the Alpha VOC was present in all the studied WWTPs, and became dominant in 11 weeks on average. The outcompetition rates of the Alpha VOC were estimated, and their relationship with different parameters statistically analyzed. The rapid spread of the Alpha VOC was influenced by its initial input and by the previous circulation of SARS-COV-2 in the population. This latter point could be explained by its higher transmissibility, particularly advantadgeous when a certain degree of herd immunity exists. Moreover, the presence of signature mutations of SARS-COV-2 variants were established by deep-sequencing of the complete S gene. The circulation of the Alpha VOC in the area under study was confirmed, and additionally two combinations of mutations in the S glycoprotein (T73A and D253N, and S477N and A522S) that could affect antibody binding were identified., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

A common trait among RNA viruses is their high capability to acquire genetic variability due to viral and host mechanisms. Next-generation sequencing (NGS) analysis enables the deep study of the viral quasispecies in samples from infected individuals. In this study, the viral quasispecies complexity and single nucleotide polymorphisms of the SARS-CoV-2 spike gene of coronavirus disease 2019 (COVID-19) patients with mild or severe disease were investigated using next-generation sequencing (Illumina platform). SARS-CoV-2 spike variability was higher in patients with long-lasting infection. Most substitutions found were present at frequencies lower than 1%, and had an A → G or T → C pattern, consistent with variants caused by adenosine deaminase acting on RNA-1 (ADAR1). ADAR1 affected a small fraction of replicating genomes, but produced multiple, mainly non-synonymous mutations. ADAR1 editing during replication rather than the RNA-dependent RNA polymerase (nsp12) was the predominant mechanism generating SARS-CoV-2 genetic variability. However, the mutations produced are not fixed in the infected human population, suggesting that ADAR1 may have an antiviral role, whereas nsp12-induced mutations occurring in patients with high viremia and persistent infection are the main source of new SARS-CoV-2 variants.

Cirrhosis derived from chronic hepatitis C virus (HCV) infection is still a common indication for liver transplantation (LT). Reinfection of the engrafted liver is universal in patients with detectable viral RNA at the time of transplant and causes fast progression of cirrhosis (within 5 years) in around one-third of these patients. To prevent damage to the liver graft, effective direct-acting antiviral (DAA) therapy is required as soon as possible. However, because of post-LT clinical instability, it is difficult to determine the optimal time to start DAAs with a low risk of complications. Evaluate changes in quasispecies complexity following LT and seek a predictive index of fast liver damage progression to determine the timing of DAA initiation. HCV genomes isolated from pre-LT and 15-day post-LT serum samples of ten patients, who underwent orthotopic LT, were quantified and sequenced using a next-generation sequencing platform. Sequence alignments, phylogenetic trees, quasispecies complexity measures, biostatistics analyses, adjusted R2 values, and analysis of variance (ANOVA) were carried out. Three different patterns of reinfection were observed (viral bottlenecking, conserved pre-LT population, and mixed populations), suggesting that bottlenecking or homogenization of the viral population is not a generalized effect after liver graft reinfection. None of the quasispecies complexity measures predicted the future degree of liver damage. Higher and more uniform viral load (VL) values were observed in all pre-LT samples, but values were more dispersed in post-LT samples. However, VL increased significantly from the pre-LT to 15-day post-LT samples in patients with advanced fibrosis at 1-year post-LT, suggesting that a VL increase on day 15 may be a predictor of fast liver fibrosis progression. HCV kinetics after LT differ between patients and are not fibrosis-dependent. Higher VL at day 15 post-LT versus pre-LT samples may predict fast liver fibrosis progression.

Since its first identification in the United Kingdom in late 2020, the highly transmissible B.1.1.7 variant of SARS-CoV-2 has become dominant in several countries raising great concern. We developed a duplex real-time RT-qPCR assay to detect, discriminate, and quantitate SARS-CoV-2 variants containing one of its mutation signatures, the ΔHV69/70 deletion, and used it to trace the community circulation of the B.1.1.7 variant in Spain through the Spanish National SARS-CoV-2 Wastewater Surveillance System (VATar COVID-19). The B.1.1.7 variant was detected earlier than clinical epidemiological reporting by the local authorities, first in the southern city of Málaga (Andalucía) in week 20&#95;52 (year&#95;week), and multiple introductions during Christmas holidays were inferred in different parts of the country. Wastewater-based B.1.1.7 tracking showed a good correlation with clinical data and provided information at the local level. Data from wastewater treatment plants, which reached B.1.1.7 prevalences higher than 90% for ≥2 consecutive weeks showed that 8.1 ± 2.0 weeks were required for B.1.1.7 to become dominant. The study highlights the applicability of RT-qPCR-based strategies to track specific mutations of variants of concern as soon as they are identified by clinical sequencing and their integration into existing wastewater surveillance programs, as a cost-effective approach to complement clinical testing during the COVID-19 pandemic.

Direct-acting antivirals (DAAs) resolve chronic HCV infection in >95% of patients, but a small percentage do not respond to DAA-based therapy. These may be difficult to treat because of resistance-associated substitutions (RAS) emerging after treatment failure. Triple therapy with sofosbuvir (SOF)/velpatasvir (VEL)/voxilaprevir (VOX) is the recommended retreatment after DAA-based failure. However, in rare cases, failure to triple therapy occurs, and there is little information characterizing the viruses that relapse. To determine the RAS profile after failing SOF/VEL/VOX, and seek suitable alternatives for retreatment, samples from 5 patients were analysed using MiSeq Illumina deep sequencing before and after triple therapy. All patients were men, aged 59-78 years, 2 HCV genotype (G) 1b and 3 G3a. The most prevalent NS3 substitutions after SOF/VEL/VOX failure were Y56F and A166T. Four patients had the NS5A RAS, Y93H, after triple failure, and Y93H was observed in both G1b patients before retreatment and after SOF/ledipasvir failure. In 2 G3a patients, Y93H appeared at triple failure, and on the other G3a, A30K persisted in 100% of viral genomes. Finally, G1b patients showed C316N in NS5B, associated with SOF failure, but G3a patients had no known NS5B substitutions. HCV RAS analysis identified the following substitutions present at higher rates after triple failure: Y56F in NS3 (G1b), A166T in NS3 (G3a), A30K or Y93H in NS5A, and C316N in NS5B (G1b). A RAS-based salvage treatment (SOF + glecaprevir/pibrentasvir + RBV) was successfully used in one G3a patient., (© 2021 John Wiley & Sons Ltd.)

The SARS-CoV-2 spike (S) protein, the viral mediator for binding and entry into the host cell, has sparked great interest as a target for vaccine development and treatments with neutralizing antibodies. Initial data suggest that the virus has low mutation rates, but its large genome could facilitate recombination, insertions, and deletions, as has been described in other coronaviruses. Here, we deep-sequenced the complete SARS-CoV-2 S gene from 18 patients (10 with mild and 8 with severe COVID-19), and found that the virus accumulates deletions upstream and very close to the S1/S2 cleavage site (PRRAR/S), generating a frameshift with appearance of a stop codon. These deletions were found in a small percentage of the viral quasispecies (2.2%) in samples from all the mild and only half the severe COVID-19 patients. Our results suggest that the virus may generate free S1 protein released to the circulation. We suggest that natural selection has favoured a "Don't burn down the house" strategy, in which free S1 protein may compete with viral particles for the ACE2 receptor, thus reducing the severity of the infection and tissue damage without losing transmission capability.

Purpose:  To study resistance-associated substitutions using next-generation sequencing in Egyptian hepatitis C virus-infected patients failing direct-acting antiviral treatment., Methods:  The current study describes three cases of treatment failure in patients referred to Zagazig Viral Hepatitis Treatment Center (ZVHTC), Sharkia Governorate, Egypt. RAS were identified and characterized using deep sequencing. The first patient had breakthrough while receiving a daclatasvir (DCV)+sofosbuvir (SOF) regimen, patient 2 relapsed after treatment with DCV+SOF+ribavirin (RBV), and patient 3 relapsed after DCV+SOF therapy. A serum sample was collected from each patient at failure and sent to Vall d'Hebron Research Institute at Hospital Universitari Vall d'Hebron in Barcelona (Spain) for deep-sequencing study to identify and characterize the RAS present in the samples., Results:  The following were identified: L28M, L30S and L28M+L30S in patient 1, L30R in patient 2, and R155C, D168E, L28M, L30H, L30S, L28M+L30H, and L28M+L30S in patient 3., Conclusion: To the best of our knowledge, this is the first report from Egypt of patients failing DAA-based therapy, describing the associated RAS. This information will be of help to understand the natural history of HCV in Egyptian patients and guide the proper choice of retreatment protocols., Competing Interests: The funder, Roche Diagnostics S.L., provided support in the form of salaries for one author (Josep Gregori), but the company did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The authors report no other conflicts of interest in this work., (© 2019 Amer et al.)

Hepatitis C virus (HCV) is a highly variable infectious agent, classified into 8 genotypes and 86 subtypes. Our laboratory has implemented an in-house developed high-resolution HCV subtyping method based on next-generation sequencing (NGS) for error-free classification of the virus using phylogenetic analysis and analysis of genetic distances in sequences from patient samples compared to reference sequences. During routine diagnostic, a sample from an Equatorial Guinea patient could not be classified into any of the existing subtypes. The whole genome was analyzed to confirm that the new isolate could be classified as a new HCV subtype. In addition, naturally occurring resistance-associated substitutions (RAS) were analyzed by NGS. Whole-genome analysis based on p-distances suggests that the sample belongs to a new HCV genotype 1 subtype. Several RAS in the NS3 (S122T, D168E and I170V) and NS5A protein (Q(1b)24K, R(1b)30Q and Y93L+Y93F) were found, which could limit the use of some inhibitors for treating this subtype. RAS studies of new subtypes are of great interest for tailoring treatment, as no data on treatment efficacy are reported. In our case, the patient has not yet been treated, and the RAS report will be used to design the most effective treatment., Competing Interests: The funder “Roche Diagnostics S.L.” provided support in the form of salaries for one author (Josep Gregori), but the company did not have any additional role in the study design, data collection and analysis, decision to publish or preparation of the manuscript. The other authors report no conflicts of interest in this work.

Background: A huge outbreak in the men-having-sex-with-men (MSM) has hit Europe during the years 2016-2018. Outbreak control has been hampered by vaccine shortages in many countries, and to minimize their impact, reduction of antigen doses has been implemented. However, these measures may have consequences on the evolution of hepatitis A virus (HAV), leading to the emergence of antigenic variants. Cases in vaccinated MSM patients have been detected in Barcelona, opening the possibility to study HAV evolution under immune pressure., Methods: We performed deep-sequencing analysis of ten overlapping fragments covering the complete capsid coding region of HAV. A total of 14578255 reads were obtained and used for the analysis of virus evolution in vaccinated versus non-vaccinated patients. We estimated maximum and minimum mutation frequencies, and Shannon entropy in the quasispecies of each patient. Non-synonymous (NSyn) mutations affecting residues exposed in the capsid surface were located, with respect to epitopes, using the recently described crystal structure of HAV, as an indication of its potential role in escaping to the effect of vaccines., Findings: HAV evolution at the quasispecies level, in non-vaccinated and vaccinated patients, revealed higher diversity in epitope-coding regions of the vaccinated group. Although amino acid replacements occurring in and around the epitopes were observed in both groups, their abundance was significantly higher in the quasispecies of vaccinated patients, indicating ongoing processes of fixation., Interpretation: Our data suggest positive selection of antigenic variants in some vaccinated patients, raising concerns for new vaccination polices directed to the MSM group., (Copyright © 2018. Published by Elsevier B.V.)

[This corrects the article DOI: 10.1371/journal.pone.0201850.].

Background: Controversy is ongoing about whether a minority mutant present at frequencies below 15% may be clinically relevant and should be considered to guide treatment., Methods: Resistance-associated substitution (RAS) studies were performed in patients before and at failure of antiviral treatments using Next-generation hepatitis C virus (HCV) sequencing (NGS)., Results: We have found two patients with genotype 1a infection having RAS in 3.5%-7.1% of the viral population at baseline that were selected during ledipasvir + sofosbuvir treatment. Co-selection of RAS located in a region not directly affected by the antiviral treatment also occurred. This observation calls into question, the recommendations to guide RAS-based direct-acting antiviral (DAA) treatment only when RAS are present in >15% of the sequences generated., Conclusion: Our results suggests that RAS study should include all three HCV DAA target proteins and minority mutants should be considered as clinically relevant., Competing Interests: Disclosure The authors report no conflicts of interest in this work.

Noroviruses are the main cause of epidemics of acute gastroenteritis at a global scale. Although chronically infected immunocompromised individuals are regarded as potential reservoirs for the emergence of new viral variants, viral quasispecies distribution and evolution patterns in acute symptomatic and asymptomatic infections has not been extensively studied. Amplicons of 450 nts from the P2 coding capsid domain were studied using next-generation sequencing (454/GS-Junior) platform. Inter-host diversity between symptomatic and asymptomatic acutely infected individuals linked to the same outbreak as well as their viral intra-host diversity over time were characterized. With an average of 2848 reads per sample and a cutoff frequency of 0.1%, minor variant haplotypes were detected in 5 out of 8 specimens. Transmitted variants could not be confirmed in all infected individuals in one outbreak. The observed initial intra-host viral diversity in asymptomatically infected subjects was higher than in symptomatic ones. Viral quasispecies evolution over time within individuals was host-specific, with an average of 2.8 nt changes per day (0.0062 changes per nucleotide per day) in a given symptomatic case. Nucleotide polymorphisms were detected in 28 out of 450 analyzed nucleotide positions, 32.14% of which were synonymous and 67.86% were non-synonymous. Most observed amino acid changes emerged at or near blockade epitopes A, B, D and E. Our results suggest that acutely infected individuals, even in the absence of symptoms, which go underreported and may enhance transmission, may contribute to norovirus genetic variability and evolution., Competing Interests: Competing Interests Statement Dr. Josep Gregori is a member of “Roche Diagnostic, S.L.”. This does not alter our adherence to PLoS ONE policies on sharing data and materials.

Hepatitis C virus (HCV) is a highly divergent virus currently classified into seven major genotypes and 86 subtypes (ICTV, June 2017), which can have differing responses to therapy. Accurate genotyping/subtyping using high-resolution HCV subtyping enables confident subtype identification, identifies mixed infections and allows detection of new subtypes. During routine genotyping/subtyping, one sample from an Equatorial Guinea patient could not be classified into any of the subtypes. The complete genomic sequence was compared to reference sequences by phylogenetic and sliding window analysis. Resistance-associated substitutions (RASs) were assessed by deep sequencing. The unclassified HCV genome did not belong to any of the existing genotype 1 (G1) subtypes. Sliding window analysis along the complete genome ruled out recombination phenomena suggesting that it belongs to a new HCV G1 subtype. Two NS5A RASs (L31V+Y93H) were found to be naturally combined in the genome which could limit treatment possibilities in patients infected with this subtype.

Objectives: This study aimed to characterize the chronically infected general hepatitis C virus (HCV) population in Barcelona using a highly sensitive subtyping method that can identify the 67 recognized HCV subtypes and diagnose mixed infection by various genotypes/subtypes in a single individual. The resulting information has implications for selecting optimal direct-acting antiviral (DAA) treatment for each patient and establishing public healthcare policies in our setting., Methods: Consecutive HCV patients (treatment-naïve or interferon-based failures) attending Vall d'Hebron Hospital outpatient clinics from February 2015 to May 2016 (N=1473) were included in the study. Patient samples were characterized using HCV subtyping by next-generation ultra-deep pyrosequencing., Results: The following genotypes (G) were found: G1 (1126/1473 (76.4%)), G4 (145/1473 (9.8%)), G3 (135/1473 (9.2%)), G2 (51/1473 (3.5%)), and G5 (1/1473 (0.1%)). Twenty-two subtypes were seen: 1b (790/1473 (53.6%)), 1a (332/1473 (22.5%)), 3a (133/1473 (9.0%)), 4d (105/1473 (7.1%)), 4a (29/1473 (2.0%)), and 2c (25/1473 (1.7%)), with 16 low-prevalence subtypes accounting for the remaining 3.0% (44/1473). There was a worrisome 1.0% (15/1473) of mixed infections. G2 (51/1473 (3.5%)) showed a high level of heterogeneity. Analyses by age groups showed a predominance of G1b over G1a (428/506 (84.6%) vs. 24/506 (4.7%)) in patients born before 1950 (N=506/1473), and similar percentages of these subtypes in those born between 1951 and 1975 (N=834/1473) (315/834, 37.8% vs. 266/834, 31.9%) and after 1976 (N=133/1473) (47/133, 35.3% vs. 42/133, 31.6%)., Conclusions: Subtype distribution showed a higher level of heterogeneity than was expected, particularly for G2. Prevalence of mixed infections was around 1%. HCV subtype distribution related to patient age group suggested that patients born from 1936 to 1975 in our setting should undergo screening for the infection. Next-generation sequencing enabled better classification of candidates for DAA-based treatment., (Copyright © 2017 European Society of Clinical Microbiology and Infectious Diseases. Published by Elsevier Ltd. All rights reserved.)

The unequivocal identification of hepatitis C virus (HCV) subtypes 1a/1b and genotypes 2 to 6 is required for optimizing the effectiveness of interferon-free, direct-acting antiviral therapies. We compared the performance of a new real-time HCV genotyping assay used on the Cobas 4800 system (C4800) with that of high-resolution HCV subtyping (HRCS). In total, 502 samples were used, including 184 samples from chronic HCV patients (from routine laboratory activity during April 2016), 5 stored samples with double HCV genotype infections for testing the limitations of the method, and 313 samples from a screening protocol implemented in our hospital (from May to August 2016) based on the new method to further determine its genotyping accuracy. A total of 282 samples, including 171 from April 2016 (the 13 remaining had too low of a viral load for HRCS), 5 selected with double infections, and 106 from screening, were analyzed by both methods, and 220 were analyzed only by the C4800. The C4800 correctly subtyped 125 of 126 1a/1b samples, and the 1 remaining sample was reported as genotype 1. The C4800 correctly genotyped 38 of 45 non-1a/1b samples (classified by HRCS), and it reported the remaining 7 samples as indeterminate. One hundred two of 106 non-1a/1b genotype samples that were identified using the C4800 for screening were confirmed by HRCS. In the 4 remaining samples, 3 were correctly reported as genotype 1 (without defining the subtype) and 1 was reported as indeterminate. None of the samples were misgenotyped. Four of 7 samples with double HCV infections were correctly genotyped by the C4800. Excluding the 5 selected double-infected samples, the C4800 showed 95.7% concordant results for genotyping HCVs 2 to 6 and 1a/1b subtyping, and 99.2% concordance for subtyping 1a/1b single infections in clinical samples. To improve laboratory workflow, we propose using the C4800 as a first-line test for HCV genotyping and 1a/1b classification, followed by transferring non-1a/1b samples to a center where HRCS is available, if further characterization is needed., (Copyright © 2017 Nieto-Aponte et al.)

Aim: To develop a fast, low-cost diagnostic strategy to identify single point mutations in highly variable genomes such as hepatitis C virus (HCV)., Methods: In patients with HCV infection, resistance-associated amino acid substitutions within the viral quasispecies prior to therapy can confer decreased susceptibility to direct-acting antiviral agents and lead to treatment failure and virological relapse. One such naturally occurring mutation is the Q80K substitution in the HCV-NS3 protease gene, which confers resistance to PI inhibitors, particularly simeprevir. Low-cost, highly sensitive techniques enabling routine detection of these single point mutations would be useful to identify patients at a risk of treatment failure. LightCycler methods, based on real-time PCR with sequence-specific probe hybridization, have been implemented in most diagnostic laboratories. However, this technique cannot identify single point mutations in highly variable genetic environments, such as the HCV genome. To circumvent this problem, we developed a new method to homogenize all nucleotides present in a region except the point mutation of interest., Results: Using nucleotide-specific probes Q, K, and R substitutions at position 80 were clearly identified at a sensitivity of 10% (mutations present at a frequency of at least 10% were detected). The technique was successfully applied to identify the Q80K substitution in 240 HCV G1 serum samples, with performance comparable to that of direct Sanger sequencing, the current standard procedure for this purpose. The new method was then validated in a Catalonian population of 202 HCV G1-infected individuals. Q80K was detected in 14.6% of G1a patients and 0% of G1b in our setting., Conclusion: A fast, low-cost diagnostic strategy based on real-time PCR and fluorescence resonance energy transfer probe melting curve analysis has been successfully developed to identify single point mutations in highly variable genomes such as hepatitis C virus. This technique can be adapted to detect any single point mutation in highly variable genomes., Competing Interests: Conflict-of-interest statement: There are no conflicts of interest to report.

Hepatitis C virus (HCV) is classified into seven major genotypes and 67 subtypes. Recent studies have shown that in HCV genotype 1-infected patients, response rates to regimens containing direct-acting antivirals (DAAs) are subtype dependent. Currently available genotyping methods have limited subtyping accuracy. We have evaluated the performance of a deep-sequencing-based HCV subtyping assay, developed for the 454/GS-Junior platform, in comparison with those of two commercial assays (Versant HCV genotype 2.0 and Abbott Real-time HCV Genotype II) and using direct NS5B sequencing as a gold standard (direct sequencing), in 114 clinical specimens previously tested by first-generation hybridization assay (82 genotype 1 and 32 with uninterpretable results). Phylogenetic analysis of deep-sequencing reads matched subtype 1 calling by population Sanger sequencing (69% 1b, 31% 1a) in 81 specimens and identified a mixed-subtype infection (1b/3a/1a) in one sample. Similarly, among the 32 previously indeterminate specimens, identical genotype and subtype results were obtained by direct and deep sequencing in all but four samples with dual infection. In contrast, both Versant HCV Genotype 2.0 and Abbott Real-time HCV Genotype II failed subtype 1 calling in 13 (16%) samples each and were unable to identify the HCV genotype and/or subtype in more than half of the non-genotype 1 samples. We concluded that deep sequencing is more efficient for HCV subtyping than currently available methods and allows qualitative identification of mixed infections and may be more helpful with respect to informing treatment strategies with new DAA-containing regimens across all HCV subtypes., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

We have investigated the reliability and reproducibility of HCV viral quasispecies quantification by ultra-deep pyrosequencing (UDPS) methods. Our study has been divided in two parts. First of all, by UDPS sequencing of clone mixes samples we have established the global noise level of UDPS and fine tuned a data treatment workflow previously optimized for HBV sequence analysis. Secondly, we have studied the reproducibility of the methodology by comparing 5 amplicons from two patient samples on three massive sequencing platforms (FLX+, FLX and Junior) after applying the error filters developed from the clonal/control study. After noise filtering the UDPS results, the three replicates showed the same 12 polymorphic sites above 0.7%, with a mean CV of 4.86%. Two polymorphic sites below 0.6% were identified by two replicates and one replicate respectively. A total of 25, 23 and 26 haplotypes were detected by GS-Junior, GS-FLX and GS-FLX+. The observed CVs for the normalized Shannon entropy (Sn), the mutation frequency (Mf), and the nucleotidic diversity (Pi) were 1.46%, 3.96% and 3.78%. The mean absolute difference in the two patients (5 amplicons each), in the GS-FLX and GS-FLX+, were 1.46%, 3.96% and 3.78% for Sn, Mf and Pi. No false polymorphic site was observed above 0.5%. Our results indicate that UDPS is an optimal alternative to molecular cloning for quantitative study of HCV viral quasispecies populations, both in complexity and composition. We propose an UDPS data treatment workflow for amplicons from the RNA viral quasispecies which, at a sequencing depth of at least 10,000 reads per strand, enables to obtain sequences and frequencies of consensus haplotypes above 0.5% abundance with no erroneous mutations, with high confidence, resistant mutants as minor variants at the level of 1%, with high confidence that variants are not missed, and highly confident measures of quasispecies complexity.

Background: Selection of amino acid substitutions associated with resistance to nucleos(t)ide-analog (NA) therapy in the hepatitis B virus (HBV) reverse transcriptase (RT) and their combination in a single viral genome complicates treatment of chronic HBV infection and may affect the overlapping surface coding region. In this study, the variability of an overlapping polymerase-surface region, critical for NA resistance, is investigated before treatment and under antiviral therapy, with assessment of NA-resistant amino acid changes simultaneously occurring in the same genome (linkage analysis) and their influence on the surface coding region., Methodology/principal Findings: Serum samples obtained from chronic HBV-infected patients at pre-treatment and during sequential NA treatment with lamivudine, adefovir, and entecavir were analyzed by ultra-deep pyrosequencing (UDPS) using the GS-FLX platform (454 Life Sciences-Roche). The pre-treatment HBV quasispecies was not enriched with NA-resistant substitutions. The frequencies of this type of substitutions at pre-treatment did not predict the frequencies observed during lamivudine treatment. On linkage analysis of the RT region studied, NA-resistant HBV variants (except for rtA181T) were present in combinations of amino acid substitutions that increased in complexity after viral breakthrough to entecavir, at which time the combined variant rtL180M-S202G-M204V-V207I predominated. In the overlapping surface region, NA-resistant substitutions caused selection of stop codons in a significant percentage of sequences both at pre-treatment and during sequential treatment; the rtA181T substitution, related to sW172stop, predominated during treatment with lamivudine and adefovir. A highly conserved RT residue (rtL155), even more conserved than the essential residues in the RT catalytic motif YMDD, was identified in all samples., Conclusions: UDPS methodology enabled quantification of HBV quasispecies variants, even those harboring complex combinations of amino acid changes. The high percentage of potentially defective genomes, especially in the surface region, suggests effective trans-complementation of these variants.

Background: HCV is a major cause of morbidity and mortality in HIV-coinfected patients. Several observational studies have suggested that HCV response to pegylated interferon and ribavirin is lower in HIV-coinfected patients treated with abacavir. It has been postulated that abacavir could compete with ribavirin to be phosphorylated, leading to a reduction in the active form of the drug (triphosphorylated ribavirin). Here, we studied the effect of abacavir, tenofovir or lamivudine addition on the suppressive activity of ribavirin in an HCV RNA replicon system., Methods: We used the human hepatoma HuH-7 cell clone 9B containing the HCV genotype 1b replicon I389/NS3-3'. Cells were treated for 24 h with ribavirin (0, 10 and 50 μM) plus abacavir, tenofovir or lamivudine at doses of 0, 10 and 50 μM and HCV RNA production was quantified by real-time PCR in triplicate assays. Results were expressed as mean±SD of the HCV RNA produced per cell (log(10) IU/cell). Means were compared using the Student's t-test., Results: Ribavirin treatment produced a dose-dependent suppression of HCV RNA production by the replicon system. Combination of ribavirin and interferon resulted in an additive antiviral activity. The addition of abacavir did not modify the suppressive activity of ribavirin on the replicon HCV RNA expression. Similar results were obtained when ribavirin was used in combination with tenofovir or lamivudine., Conclusions: In a subgenomic HCV RNA replicon system, the antiviral effect of ribavirin was not modified by the addition of abacavir, tenofovir or lamivudine., (© 2011 International Medical Press)