ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
Subject(s)
COVID-19 , Nanopores , Humans , SARS-CoV-2/genetics , Wastewater , COVID-19/epidemiology , High-Throughput Nucleotide Sequencing/methodsABSTRACT
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.
Subject(s)
Antiviral Agents , Hepatitis E virus , Humans , Antiviral Agents/pharmacology , Antiviral Agents/therapeutic use , Hepatitis E virus/genetics , Mutagens , Quasispecies/genetics , Ribavirin/pharmacology , Ribavirin/therapeutic useABSTRACT
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.
Subject(s)
Metagenomics , Viruses , Antiviral Agents , Computational Biology , Pandemics , Viruses/geneticsABSTRACT
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.
Subject(s)
Hepatitis E virus , Hepatitis E , Ribavirin , Humans , Follow-Up Studies , Mutagens , Nucleotides , Quasispecies/genetics , Ribavirin/therapeutic use , SARS-CoV-2/genetics , Hepatitis E/drug therapy , Hepatitis E virus/drug effects , Hepatitis E virus/geneticsABSTRACT
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.
Subject(s)
COVID-19 Vaccines , COVID-19 , Humans , COVID-19/epidemiology , COVID-19/genetics , SARS-CoV-2/genetics , Codon, NonsenseABSTRACT
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.
ABSTRACT
Background: Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the main concern is whether reinfections are possible, and which are the associated risk factors. This study aims to describe the clinical and molecular characteristics of 24 sequence-confirmed reinfection SARS-CoV-2 cases over 1 year in Barcelona (Catalonia, Spain). Methods: Patients with > 45 days between two positive PCR tests regardless of symptoms and negative tests between episodes were initially considered as suspected reinfection cases from November 2020 to May 2021. Whole-genome sequencing (WGS) was performed to confirm genetic differences between consensus sequences and for phylogenetic studies based on PANGOLIN nomenclature. Reinfections were confirmed by the number of mutations, change in lineage, or epidemiological criteria. Results: From 39 reported suspected reinfection cases, complete viral genomes could be sequenced from both episodes of 24 patients, all were confirmed as true reinfections. With a median age of 44 years (interquartile range [IQR] 32-65), 66% were women and 58% were healthcare workers (HCWs). The median days between episodes were 122 (IQR 72-199), occurring one-third within 3 months. Reinfection episodes were frequently asymptomatic and less severe than primary infections. The absence of seroconversion was associated with symptomatic reinfections. Only one case was reinfected with a variant of concern (VOC). Conclusion: Severe acute respiratory syndrome coronavirus 2 reinfections can occur in a shorter time than previously reported and are mainly found in immunocompetent patients. Surveillance through WGS is useful to identify viral mutations associated with immune evasion.
ABSTRACT
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.
Subject(s)
Hepatitis C, Chronic , Viruses , Animals , Antiviral Agents , High-Throughput Nucleotide Sequencing , PandemicsABSTRACT
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.
Subject(s)
COVID-19 , SARS-CoV-2 , COVID-19/epidemiology , Humans , SARS-CoV-2/genetics , Wastewater , Wastewater-Based Epidemiological MonitoringABSTRACT
Herein, we describe the genetic diversity of circulating SARS-CoV-2 viruses by whole-genome sequencing (WGS) in Barcelona city (Catalonia, Spain) throughout the first four pandemic waves. From weeks 11/2020-24/2021, SARS-CoV-2-positive respiratory samples were randomly selected per clinical setting (80% from primary care or 20% from the hospital), age group, and week. WGS was performed following the ARTICv3 protocol on MiSeq or NextSeq2000 Illumina platforms. Nearly complete consensus sequences were used for genetic characterization based on GISAID and PANGOLIN nomenclatures. From 2475 samples, 2166 (87%) were fully sequenced (78% from primary care and 22% from hospital settings). Multiple genetic lineages were co-circulating, but four were predominant at different periods. While B.1.5 (50.68%) and B.1.1 (32.88%) were the major lineages during the first pandemic wave, B.1.177 (66.85%) and B.1.1.7 (83.80%) were predominant during the second, third, and fourth waves, respectively. Almost all (96.4%) were carrying D614G mutation in the S protein, with additional mutations that define lineages or variants. But some mutations of concern, such as E484K from B.1.351 and P.1 lineages are currently under monitoring, together with those observed in the receptor-binding domain or N-terminal domain, such as L452R and T478K from B.1.617.2 lineage. The fact that a predominant lineage was observed in each pandemic wave suggests advantageous properties over other contemporary co-circulating variants. This genetic variability should be monitored, especially when a massive vaccination campaign is ongoing because the potential selection and emergence of novel antigenic SARS-CoV-2 strains related to immunological escapement events.
Subject(s)
COVID-19/epidemiology , Genome, Viral , Mutation , SARS-CoV-2/classification , SARS-CoV-2/genetics , Spike Glycoprotein, Coronavirus/genetics , Adolescent , Adult , Aged , COVID-19/prevention & control , COVID-19/transmission , COVID-19/virology , COVID-19 Vaccines/administration & dosage , Child , Child, Preschool , Computational Biology/methods , Epidemiological Monitoring , Female , Humans , Infant , Infant, Newborn , Male , Middle Aged , Nasopharynx/virology , Physical Distancing , Prevalence , SARS-CoV-2/pathogenicity , Spain/epidemiology , Vaccination/methods , Whole Genome SequencingABSTRACT
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.
Subject(s)
Hepacivirus/classification , Hepatitis C, Chronic/therapy , Liver Cirrhosis/virology , Liver Transplantation/adverse effects , Antiviral Agents/pharmacology , Antiviral Agents/therapeutic use , Disease Progression , Hepacivirus/genetics , Hepacivirus/isolation & purification , Humans , Liver Cirrhosis/prevention & control , Phylogeny , Quasispecies , Time-to-TreatmentABSTRACT
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_52 (year_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.
Subject(s)
COVID-19 , SARS-CoV-2 , Humans , Pandemics , WastewaterABSTRACT
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.
Subject(s)
COVID-19/virology , Genetic Variation , SARS-CoV-2/genetics , Spike Glycoprotein, Coronavirus/genetics , Adult , Amino Acid Sequence , Base Sequence , Female , Host-Pathogen Interactions , Humans , Male , Middle Aged , Protein Conformation , SARS-CoV-2/physiology , Virus ReplicationABSTRACT
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.
Subject(s)
Hepatitis C, Chronic , Sofosbuvir , Aged , Aminoisobutyric Acids , Antiviral Agents/therapeutic use , Carbamates , Cyclopropanes , Drug Resistance, Viral , Drug Therapy, Combination , Genotype , Hepacivirus/genetics , Hepatitis C, Chronic/drug therapy , Heterocyclic Compounds, 4 or More Rings , Humans , Lactams, Macrocyclic , Leucine/analogs & derivatives , Male , Middle Aged , Proline/analogs & derivatives , Quinoxalines , Sofosbuvir/therapeutic use , Sulfonamides , Treatment Failure , Viral Nonstructural Proteins/geneticsABSTRACT
Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed.
Subject(s)
Hepatitis C, Chronic , Hepatitis C , Amino Acid Substitution , Antiviral Agents/pharmacology , Antiviral Agents/therapeutic use , Drug Resistance, Viral/genetics , Genotype , Hepacivirus/genetics , Hepatitis C/drug therapy , Hepatitis C, Chronic/drug therapy , Humans , Treatment Failure , Viral Nonstructural Proteins/geneticsABSTRACT
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.
Subject(s)
Betacoronavirus/genetics , Coronavirus Infections/virology , Genome, Viral/genetics , Pneumonia, Viral/virology , Quasispecies/genetics , Respiratory Tract Infections/virology , Spike Glycoprotein, Coronavirus/genetics , Adult , Aged , COVID-19 , Computational Biology , Female , Gene Deletion , High-Throughput Nucleotide Sequencing , Humans , Male , Middle Aged , Pandemics , RNA Cleavage , SARS-CoV-2 , Sequence Analysis, RNAABSTRACT
A percentage of hepatitis C virus (HCV)-infected patients fail direct acting antiviral (DAA)-based treatment regimens, often because of drug resistance-associated substitutions (RAS). The aim of this study was to characterize the resistance profile of a large cohort of patients failing DAA-based treatments, and investigate the relationship between HCV subtype and failure, as an aid to optimizing management of these patients. A new, standardized HCV-RAS testing protocol based on deep sequencing was designed and applied to 220 previously subtyped samples from patients failing DAA treatment, collected in 39 Spanish hospitals. The majority had received DAA-based interferon (IFN) α-free regimens; 79% had failed sofosbuvir-containing therapy. Genomic regions encoding the nonstructural protein (NS) 3, NS5A, and NS5B (DAA target regions) were analyzed using subtype-specific primers. Viral subtype distribution was as follows: genotype (G) 1, 62.7%; G3a, 21.4%; G4d, 12.3%; G2, 1.8%; and mixed infections 1.8%. Overall, 88.6% of patients carried at least 1 RAS, and 19% carried RAS at frequencies below 20% in the mutant spectrum. There were no differences in RAS selection between treatments with and without ribavirin. Regardless of the treatment received, each HCV subtype showed specific types of RAS. Of note, no RAS were detected in the target proteins of 18.6% of patients failing treatment, and 30.4% of patients had RAS in proteins that were not targets of the inhibitors they received. HCV patients failing DAA therapy showed a high diversity of RAS. Ribavirin use did not influence the type or number of RAS at failure. The subtype-specific pattern of RAS emergence underscores the importance of accurate HCV subtyping. The frequency of "extra-target" RAS suggests the need for RAS screening in all three DAA target regions.
Subject(s)
Antiviral Agents/therapeutic use , Drug Resistance, Multiple, Viral/genetics , Hepacivirus/drug effects , Hepacivirus/genetics , Mutation , Antiviral Agents/pharmacology , Cohort Studies , Drug Therapy, Combination , Genotype , Hepatitis C/drug therapy , High-Throughput Nucleotide Sequencing , Humans , Spain , Treatment FailureABSTRACT
Previous studies documented that long-term hepatitis C virus (HCV) replication in human hepatoma Huh-7.5 cells resulted in viral fitness gain, expansion of the mutant spectrum, and several phenotypic alterations. In the present work, we show that mutational waves (changes in frequency of individual mutations) occurred continuously and became more prominent as the virus gained fitness. They were accompanied by an increasing proportion of heterogeneous genomic sites that affected 1 position in the initial HCV population and 19 and 69 positions at passages 100 and 200, respectively. Analysis of biological clones of HCV showed that these dynamic events affected infectious genomes, since part of the fluctuating mutations became incorporated into viable genomes. While 17 mutations were scored in 3 biological clones isolated from the initial population, the number reached 72 in 3 biological clones from the population at passage 200. Biological clones differed in their responses to antiviral inhibitors, indicating a phenotypic impact of viral dynamics. Thus, HCV adaptation to a specific constant environment (cell culture without external influences) broadens the mutant repertoire and does not focus the population toward a limited number of dominant genomes. A retrospective examination of mutant spectra of foot-and-mouth disease virus passaged in cell cultures suggests a parallel behavior here described for HCV. We propose that virus diversification in a constant environment has its basis in the availability of multiple alternative mutational pathways for fitness gain. This mechanism of broad diversification should also apply to other replicative systems characterized by high mutation rates and large population sizes.IMPORTANCE The study shows that extensive replication of an RNA virus in a constant biological environment does not limit exploration of sequence space and adaptive options. There was no convergence toward a restricted set of adapted genomes. Mutational waves and mutant spectrum broadening affected infectious genomes. Therefore, profound modifications of mutant spectrum composition and consensus sequence diversification are not exclusively dependent on environmental alterations or the intervention of population bottlenecks.