Search

Our aim was to develop an accurate, highly sensitive method for HBV genotype determination and detection of genotype mixtures. We examined the preS and 5′ end of the HBV X gene (5X) regions of the HBV genome using next-generation sequencing (NGS). The 1852 haplotypes obtained were subjected to genotyping via the Distance-Based discrimination method (DB Rule) using two sets of 95 reference sequences of genotypes A–H. In clinical samples from 125 patients, the main genotypes were A, D, F and H in Caucasian, B and C in Asian and A and E in Sub-Saharan patients. Genotype mixtures were identified in 28 (22.40%) cases, and potential intergenotypic recombination was observed in 29 (23.20%) cases. Furthermore, we evaluated sequence conservation among haplotypes classified into genotypes A, C, D, and E by computing the information content. The preS haplotypes exhibited limited shared conserved regions, whereas the 5X haplotypes revealed two groups of conserved regions across the genotypes assessed. In conclusion, we developed an NGS-based HBV genotyping method utilizing the DB Rule for genotype classification. We identified two regions conserved across different genotypes at 5X, offering promising targets for RNA interference-based antiviral therapies. [ABSTRACT FROM AUTHOR]

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 FROM AUTHOR]

In quasispecies diversity studies, the comparison of two samples of varying sizes is a common necessity. However, the sensitivity of certain diversity indices to sample size variations poses a challenge. To address this issue, rarefaction emerges as a crucial tool, serving to normalize and create fairly comparable samples. This study emphasizes the imperative nature of sample size normalization in quasispecies diversity studies using next-generation sequencing (NGS) data. We present a thorough examination of resampling schemes using various simple hypothetical cases of quasispecies showing different quasispecies structures in the sense of haplotype genomic composition, offering a comprehensive understanding of their implications in general cases. Despite the big numbers implied in this sort of study, often involving coverages exceeding 100,000 reads per sample and amplicon, the rarefaction process for normalization should be performed with repeated resampling without replacement, especially when rare haplotypes constitute a significant fraction of interest. However, it is noteworthy that different diversity indices exhibit distinct sensitivities to sample size. Consequently, some diversity indicators may be compared directly without normalization, or instead may be resampled safely with replacement. [ABSTRACT FROM AUTHOR]

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 FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) are potential diagnostic and prognostic biomarkers. However, discrepancies in miRNA patterns and their validation are still frequent due to differences in sample origin, EV isolation, and miRNA sequencing methods. The aim of the present study is to find a reliable EV isolation method for miRNA sequencing, adequate for clinical application. To this aim, two comparative studies were performed in parallel with the same human plasma sample: (i) isolation and characterization of EVs obtained using three procedures: size exclusion chromatography (SEC), iodixanol gradient (GRAD), and its combination (SEC+GRAD) and (ii) evaluation of the yield of miRNA sequences obtained using NextSeq 500 (Illumina) and three miRNA library preparation protocols: NEBNext, NEXTFlex, and SMARTer smRNA-seq. The conclusion of comparison (i) is that recovery of the largest amount of EVs and reproducibility were attained with SEC, but GRAD and SEC+GRAD yielded purer EV preparations. The conclusion of (ii) is that the NEBNext library showed the highest reproducibility in the number of miRNAs recovered and the highest diversity of miRNAs. These results render the combination of GRAD EV isolation and NEBNext library preparation for miRNA retrieval as adequate for clinical applications using plasma samples. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

What takes decades, centuries or millennia to happen with a natural ecosystem, it takes only days, weeks or months with a replicating viral quasispecies in a host, especially when under treatment. Some methods to quantify the evolution of a quasispecies are introduced and discussed, along with simple simulated examples to help in the interpretation and understanding of the results. The proposed methods treat the molecules in a quasispecies as individuals of competing species in an ecosystem, where the haplotypes are the competing species, and the ecosystem is the quasispecies in a host, and the evolution of the system is quantified by monitoring changes in haplotype frequencies. The correlation between the proposed indices is also discussed, and the R code used to generate the simulations, the data and the plots is provided. The virtues of the proposed indices are finally shown on a clinical case. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: In this unicentric prospective study, we analyzed the most prevalent mutations in HCC (TERT promoter, TP53, CTNNB1, AXIN1 and ARID1A) in plasma cfDNA by next-generation sequencing, aiming to elucidate their value as prognostic noninvasive biomarkers. Total cfDNA (cut-off value 2 ng/µL), number of mutated genes and number of detectable mutations on cfDNA were significantly related to mortality. Number of mutated genes and number of detected mutations in cfDNA and the ratio between number of mutations and total amount of cfDNA were also significantly associated with recurrence. Detection of more than four mutations in cfDNA correlated with a higher risk of death. Dynamic changes in cfDNA mutations were detected prior to radiological detection of HCC recurrence. We believe that these results support the proof of principle and launching of validation studies to confirm that total cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patients. Background: Cell-free DNA (cfDNA) concentrations have been described to be inversely correlated with prognosis in cancer. Mutations in HCC-associated driver genes in cfDNA have been reported, but their relation with patient's outcome has not been described. Our aim was to elucidate whether mutations found in cfDNA could be representative from those present in HCC tissue, providing the rationale to use the cfDNA to monitor HCC. Methods: Tumoral tissue, paired nontumor adjacent tissue and blood samples were collected from 30 HCC patients undergoing curative therapies. Deep sequencing targeting HCC driver genes was performed. Results: Patients with more than 2 ng/µL of cfDNA at diagnosis had higher mortality (mean OS 24.6 vs. 31.87 months, p = 0.01) (AUC = 0.782). Subjects who died during follow-up, had a significantly higher number of mutated genes (p = 0.015) and number of mutations (p = 0.015) on cfDNA. Number of mutated genes (p = 0.001), detected mutations (p = 0.001) in cfDNA and ratio (number of mutations/cfDNA) (p = 0.003) were significantly associated with recurrence. However, patients with a ratio (number of mutations/cfDNA) above 6 (long-rank p = 0.0003) presented a higher risk of recurrence than those with a ratio under 6. Detection of more than four mutations in cfDNA correlated with higher risk of death (long-rank p = 0.042). Conclusions: In summary, cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patients [ABSTRACT FROM AUTHOR]

Deletions in the 3′ end region of the hepatitis B virus (HBV) X open reading frame (HBX) may affect the core promoter (Cp) and have been frequently associated with hepatocellular carcinoma (HCC). The aim of this study was to investigate the presence of variants with deletions and/or insertions (Indels) in this region in the quasispecies of 50 chronic hepatitis B (CHB) patients without HCC. We identified 103 different Indels in 47 (94%) patients, in a median of 3.4% of their reads (IQR, 1.3–8.4%), and 25% (IQR, 13.1–40.7%) of unique sequences identified in each quasispecies (haplotypes). Of those Indels, 101 (98.1%) caused 44 different altered stop codons, the most commonly observed were at positions 128, 129, 135, and 362 (putative position). Moreover, 39 (37.9%) Indels altered the TATA-like box (TA) sequences of Cp; the most commonly observed caused TA2 + TA3 fusion, creating a new putative canonical TATA box. Four (8%) patients developed negative clinical outcomes after a median follow-up of 9.4 (8.7–12) years. In conclusion, we observed variants with Indels in the HBX 3′ end in the vast majority of our CHB patients, some of them encoding alternative versions of HBx with potential functional roles, and/or alterations in the regulation of transcription. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Selection of inhibitor-resistant viral mutants is universal for viruses that display quasi-species dynamics, and hepatitis C virus (HCV) is no exception. Here we review recent results on drug resistance in HCV, with emphasis on resistance to the newly-developed, directly-acting antiviral agents, as they are increasingly employed in the clinic. We put the experimental observations in the context of quasi-species dynamics, in particular what the genetic and phenotypic barriers to resistance mean in terms of exploration of sequence space while HCV replicates in the liver of infected patients or in cell culture. Strategies to diminish the probability of viral breakthrough during treatment are briefly outlined. [ABSTRACT FROM AUTHOR]

The hepatitis delta virus (HDV) genome has an autocatalytic region called the ribozyme, which is essential for viral replication. The aim of this study was to use next-generation sequencing (NGS) to analyze the ribozyme quasispecies (QS) in order to study its evolution and identify highly conserved regions potentially suitable for a gene-silencing strategy. HDV RNA was extracted from 2 longitudinal samples of chronic HDV patients and the ribozyme (nucleotide, nt 688–771) was analyzed using NGS. QS conservation, variability and genetic distance were analyzed. Mutations were identified by aligning sequences with their specific genotype consensus. The main relevant mutations were tested in vitro. The ribozyme was conserved overall, with a hyper-conserved region between nt 715–745. No difference in QS was observed over time. The most variable region was between nt 739–769. Thirteen mutations were observed, with three showing a higher frequency: T23C, T69C and C64 deletion. This last strongly reduced HDV replication by more than 1 log in vitro. HDV Ribozyme QS was generally highly conserved and was maintained during follow-up. The most conserved portion may be a valuable target for a gene-silencing strategy. The presence of the C64 deletion may strongly impair viral replication, as it is a potential mechanism of viral persistence. [ABSTRACT FROM AUTHOR]

Universal history is characterized by continuous evolution, in which civilizations are born and die. This evolution is associated with multiple factors, among which the role of microorganisms is often overlooked. Viruses and bacteria have written or decisively contributed to terrible episodes of history, such as the Black Death in 14th century Europe, the annihilation of pre-Columbian American civilizations, and pandemics such as the 1918 Spanish flu or the current COVID-19 pandemic caused by the coronavirus SARS-CoV-2. Nevertheless, it is clear that we could not live in a world without these tiny beings. Endogenous retroviruses have been key to our evolution and for the regulation of gene expression, and the gut microbiota helps us digest compounds that we could not otherwise process. In addition, we have used microorganisms to preserve or prepare food for millennia and more recently to obtain drugs such as antibiotics or to develop recombinant DNA technologies. Due to the enormous importance of microorganisms for our survival, they have significantly influenced the population genetics of different human groups. This paper will review the role of microorganisms as "villains" who have been responsible for tremendous mortality throughout history but also as "friends" who help us survive and evolve. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed. [ABSTRACT FROM AUTHOR]

Diagnosis and clinical management of people infected with hepatitis C virus (HCV) relies on results from a combination of serological and virological tests. The aim of this study was to compare the performance of dried plasma spots (DPS), prepared using the cobas® Plasma Separation Card (PSC), to plasma and serum from venipuncture, for HCV diagnosis. We carried out a prospective study using DPS and paired plasma or serum samples. Serum and DPS samples were analyzed by immunoassay using Elecsys® Anti-HCV II (Roche). Plasma and DPS samples were analyzed using the cobas® HCV viral load and cobas® HCV genotyping tests (Roche). All DPS samples that had high anti-HCV antibody titers in serum were also antibody-positive, as were five of eight samples with moderate titers. Eight samples with low titers in serum were negative with DPS. Among 80 samples with plasma HCV viral loads between 61.5 and 2.2 × 108 IU/mL, 74 were RNA-positive in DPS. The mean viral load difference between plasma and DPS was 2.65 log10 IU/mL. The performance of DPS for detection of serological and virological markers of hepatitis C virus infection was comparable to that of the conventional specimen types. However, the limits of detection were higher for DPS. [ABSTRACT FROM AUTHOR]

The influence of quasispecies dynamics on long-term virus diversification in nature is a largely unexplored question. Specifically, whether intra-host nucleotide and amino acid variation in quasispecies fit the variation observed in consensus sequences or data bank alignments is unknown. Genome conservation and dynamics simulations are used for the computational design of universal vaccines, therapeutic antibodies and pan-genomic antiviral agents. The expectation is that selection of escape mutants will be limited when mutations at conserved residues are required. This strategy assumes long-term (epidemiologically relevant) conservation but, critically, does not consider short-term (quasispecies-dictated) residue conservation. We calculated mutant frequencies of individual loci from mutant spectra of hepatitis C virus (HCV) populations passaged in cell culture and from infected patients. Nucleotide or amino acid conservation in consensus sequences of the same populations, or in the Los Alamos HCV data bank did not match residue conservation in mutant spectra. The results relativize the concept of sequence conservation in viral genetics and suggest that residue invariance in data banks is an insufficient basis for the design of universal viral ligands for clinical purposes. Our calculations suggest relaxed mutational restrictions during quasispecies dynamics, which may contribute to higher calculated short-term than long-term viral evolutionary rates. [ABSTRACT FROM AUTHOR]