"Outlaw" mutations in quasispecies of SARS-CoV-2 inhibit replication.

The evolution of SARS-CoV-2, the agent of COVID-19, has been remarkable for its high mutation potential, leading to the appearance of variants. Some mutations have never appeared in the published genomes, which represent consensus, or bona fide genomes. Here we tested the hypothesis that mutations that did not appear in consensus genomes were, in fact, as frequent as the mutations that appeared during the various epidemic episodes, but were not expressed because lethal. To identify these mutations, we analysed the genomes of 90 nasopharyngeal samples and the quasispecies determined by next-generation sequencing. Mutations observed in the quasispecies and not in the consensus genomes were considered to be lethal, what we called "outlaw" mutations. Among these mutations, we analysed the 21 most frequent. Eight of these "outlaws" were in the RNA polymerase and we were able to use a structural biology model and molecular dynamics simulations to demonstrate the functional incapacity of these mutated RNA polymerases. Three other mutations affected the spike, a major protein involved in the pathogenesis of COVID-19. Overall, by analysing the SARS-CoV-2 quasispecies obtained during sequencing, this method made it possible to identify "outlaws," showing areas that could potentially become the target of treatments.

Virus Quasispecies Rarefaction: Subsampling with or without Replacement?

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.

Characterization of BCP/PreC/C region quasispecies in treatment-naive patients with different phases of HBV infection using next-generation sequencing.

BACKGROUND: To analysis of quasispecies (QS) changes and high-frequency mutations in the BCP/PreC/C region of patients at different phases of hepatitis B virus (HBV) infection and provides novel biomarkers for the diagnosis of chronic hepatitis B (CHB) patients. METHODS: With the application of next-generation sequencing technology, we were able to sequence the HBV BCP/PreC/C regions in 40 patients, each at different phases of the HBV infection. The heterogeneity of QS and the frequency of mutations were calculated using MEGA 7 software. RESULTS: Our results show that the complexity and diversity of the BCP/PreC/C QS in HBeAg-positive CHB patients are significantly higher than those in HBeAg-positive chronic infection patients, while HBeAg-negative chronic infection patients had significantly higher QS complexity and diversity than HBeAg-negative CHB patients. In addition, HBeAg-negative patients showed reduced complexity but increased diversity compared with HBeAg-positive patients. Receiver operating characteristic curves showed that G1764A, C2102T, dN and complexity of QS could be used as potential biomarkers for diagnosing HBeAg-positive CHB, while the A2189C, dS and complexity of QS could be used as potential biomarkers for diagnosing HBeAg-negative chronic hepatitis. Finally, our study also found that G1896A and A2159G may be hotspot mutations affecting HBeAg seroconversion. CONCLUSION: Our research elucidates the evolution of HBV by analyzing QS heterogeneity and mutation patterns, offering novel serum biomarkers for enhancing clinical diagnosis and disease prognosis. This comprehensive approach sheds light on the intricate dynamics of HBV progression and paves the way for more precise medical interventions.

RNA three-dimensional structure drives the sequence organization of potato spindle tuber viroid quasispecies.

RNA viruses and viroids exist and evolve as quasispecies due to error-prone replication. Quasispecies consist of a few dominant master sequences alongside numerous variants that contribute to genetic diversity. Upon environmental changes, certain variants within quasispecies have the potential to become the dominant sequences, leading to the emergence of novel infectious strains. However, the emergence of new infectious variants remains unpredictable. Using mutant pools prepared by saturation mutagenesis of selected stem and loop regions, our study of potato spindle tuber viroid (PSTVd) demonstrates that mutants forming local three-dimensional (3D) structures similar to the wild type (WT) are more likely to accumulate in PSTVd quasispecies. The selection mechanisms underlying this biased accumulation are likely associated with cell-to-cell movement and long-distance trafficking. Moreover, certain trafficking-defective PSTVd mutants can be spread by functional sister genomes in the quasispecies. Our study reveals that the RNA 3D structure of stems and loops constrains the evolution of viroid quasispecies. Mutants with a structure similar to WT have a higher likelihood of being maintained within the quasispecies and can potentially give rise to novel infectious variants. These findings emphasize the potential of targeting RNA 3D structure as a more robust approach to defend against viroid infections.

Immune System Deficiencies Do Not Alter SARS-CoV-2 Evolutionary Rate but Favour the Emergence of Mutations by Extending Viral Persistence.

During the COVID-19 pandemic, immunosuppressed patients showed prolonged SARS-CoV-2 infections, with several studies reporting the accumulation of mutations in the viral genome. The weakened immune system present in these individuals, along with the effect of antiviral therapies, are thought to create a favourable environment for intra-host viral evolution and have been linked to the emergence of new viral variants which strongly challenged containment measures and some therapeutic treatments. To assess whether impaired immunity could lead to the increased instability of viral genomes, longitudinal nasopharyngeal swabs were collected from eight immunocompromised patients and fourteen non-immunocompromised subjects, all undergoing SARS-CoV-2 infection. Intra-host viral evolution was compared between the two groups through deep sequencing, exploiting a probe-based enrichment method to minimise the possibility of artefactual mutations commonly generated in amplicon-based methods, which heavily rely on PCR amplification. Although, as expected, immunocompromised patients experienced significantly longer infections, the acquisition of novel intra-host viral mutations was similar between the two groups. Moreover, a thorough analysis of viral quasispecies showed that the variability of viral populations in the two groups is comparable not only at the consensus level, but also when considering low-frequency mutations. This study suggests that a compromised immune system alone does not affect SARS-CoV-2 within-host genomic variability.

Quasispecies productivity.

The quasispecies theory is a helpful concept in the explanation of RNA virus evolution and behaviour, with a relevant impact on methods used to fight viral diseases. It has undergone some adaptations to integrate new empirical data, especially the non-deterministic nature of mutagenesis, and the variety of behavioural motifs in cooperation, competition, communication, innovation, integration, and exaptation. Also, the consortial structure of quasispecies with complementary roles of memory genomes of minority populations better fits the empirical data than did the original concept of a master sequence and its mutant spectra. The high productivity of quasispecies variants generates unique sequences that never existed before and will never exist again. In the present essay, we underline that such sequences represent really new ontological entities, not just error copies of previous ones. Their primary unique property, the incredible variant production, is suggested here as quasispecies productivity, which replaces the error-replication narrative to better fit into a new relationship between mankind and living nature in the twenty-first century.

The dynamic variation position and predominant quasispecies of hepatitis B virus: Novel predictors of early hepatocarcinoma.

To find the predictors of early HCC based on the dynamic changes of HBV quasispecies, this study utilizing the second-generation sequencing (NGS) and high-order multiplex droplet digital PCR (ddPCR) technology to examine the HBV quasispecies in serum of total 247 subjects recruited from high-incidence area of HCC. In the discovery stage, 15 non-synonymous Single Nucleotide Polymorphisms (SNPs) with higher variant proportion in HCC case group were founded (all P<0.05). Furthermore, the variant proportions in some of these SNPs were observed changing regularly within 5 years before the onset of HCC, and 5 of them located in HBX, 2 in HBS and 2 in HBC. The HBV predominant quasispecies and their consensus sequences were identified by genetic evolution analysis, in which the high HBS and HBC quasispecies heterogeneity were found associated with the forming of multifarious quasispecies clones, and the HBX gene had the highest proportion of predominant quasispecies (46.7 % in HBX vs 12.7 % and 13.8 % in HBS and HBC respectively) with the key variations (G1512A, A1630G, T1753C/G/A, A1762T and G1764A) determined. In the validation stage, we confirmed that the combined double mutations of G1512A+A1630G, A1762T+G1764A, and the combined triple mutations of T1753C/G/A + A1762T+G1764A, all expressed higher in early HCC cases when comparing with control group (all P<0.05). We also demonstrated the advantages of ddPCR using in multi-variations detection in large-sample for early HCC surveillance and screening. So we think that the dynamic of key HBV variation positions and their different combinations determined by quasispecies anlysis in this study can act as the novel predictors of early hepatocarcinoma and suitable to popularize and apply in HCC screening.

Cell-cell communication and initial population composition shape the structure of potato spindle tuber viroid quasispecies.

RNA viruses and viroids replicate with high mutation rates, forming quasispecies, population of variants centered around dominant sequences. The mechanisms governing quasispecies remain unclear. Plasmodesmata regulate viroid movement and were hypothesized to impact viroid quasispecies. Here, we sequenced the progeny of potato spindle tuber viroid intermediate (PSTVd-I) strain from mature guard cells lacking plasmodesmal connections and from in vitro-cultivated mesophyll cell protoplasts from systemic leaves of early-infected tomato (Solanum lycopersicum) plants. Remarkably, more variants accumulated in guard cells compared to whole leaves. Similarly, after extended cell culture, we observed more variants in cultivated mesophyll protoplasts. Coinfection and single-cell sequencing experiments demonstrated that the same plant cell can be infected multiple times by the same or different PSTVd sequences. To study the impact of initial population composition on PSTVd-I quasispecies, we conducted coinfections with PSTVd-I and variants. Two inoculum ratios (10:1 or 1:10) established quasispecies with or without PSTVd-I as the master sequence. In the absence of the master sequence, the percentage of novel variants initially increased. Moreover, a 1:1 PSTVd-I/variant RNA ratio resulted in PSTVd-I dominating (>50%), while the variants reached 20%. After PSTVd-I-only infection, the variants reached around 10%, while after variant-only infection, the variants were significantly more than 10%. These results emphasize the role of cell-to-cell communication and initial population composition in shaping PSTVd quasispecies.

Modelling the structure and evolution of cultural information as quasispecies.

We present a novel mode of cultural evolution whereby some forms of transmission may be modelled as quasispecies. The model incorporates the effect of high rates of error in certain forms of communication; while also building on the structural similarities between biological molecules and written language. Firstly, both written language and key biological molecules, such as RNA and proteins, are modular. Within these molecules, structural domains may be recombined, while retaining their function. Likewise, sentences are structured as combinations of clauses, in which each clause contains a domain of information. The clausal structure permits the recombination of information to adopt different meanings, while allowing each unit to retain its identity. Secondly, by virtue of intrinsically-high error rates, we show that some, but not all, aspects of communicated culture information exists as rapidly evolving clouds within the population. These clouds of cultural information behave as quasispecies, which we model with varying mutation rates and suitable selection coefficients. We then integrate these ideas with the application of Shannon Diversity Index to produce a more holistic view of culture that is centred on the evolution of its information. Re-imagining culture, as evolving clouds of information, unifies the mode in which information is stored culturally and biologically, and opens up new avenues of comparative analysis.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies.

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

In-Host HEV Quasispecies Evolution Shows the Limits of Mutagenic Antiviral Treatments.

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.

From viral democratic genomes to viral wild bunch of quasispecies.

The tremendous majority of RNA genomes from pathogenic viruses analyzed and deposited in databases are consensus or "democratic" genomes. They represent the genomes most frequently found in the clinical samples of patients but do not account for the huge genetic diversity of coexisting genomes, which is better described as quasispecies. A viral quasispecies is defined as the dynamic distribution of nonidentical but closely related mutants, variants, recombinant, or reassortant viral genomes. Viral quasispecies have collective behavior and dynamics and are the subject of internal interactions that comprise interference, complementation, or cooperation. In the setting of SARS-CoV-2 infection, intrahost SARS-CoV-2 genetic diversity was recently notably reported for immunocompromised, chronically infected patients, for patients treated with monoclonal antibodies targeting the viral spike protein, and for different body compartments of a single patient. A question that deserves attention is whether such diversity is generated postinfection from a clonal genome in response to selection pressure or is already present at the time of infection as a quasispecies. In the present review, we summarize the data supporting that hosts are infected by a "wild bunch" of viruses rather than by multiple virions sharing the same genome. Each virion in the "wild bunch" may have different virulence and tissue tropisms. As the number of viruses replicated during host infections is huge, a viral quasispecies at any time of infection is wide and is also influenced by host-specific selection pressure after infection, which accounts for the difficulty in deciphering and predicting the appearance of more fit variants and the evolution of epidemics of novel RNA viruses.

The influence of immune cells on the existence of virus quasi-species.

This article investigate a nonlocal reaction-diffusion system of equations modeling virus distribution with respect to their genotypes in the interaction with the immune response. This study demonstrates the existence of pulse solutions corresponding to virus quasi-species. The proof is based on the Leray-Schauder method, which relies on the topological degree for elliptic operators in unbounded domains and a priori estimates of solutions. Furthermore, linear stability analysis of a spatially homogeneous stationary solution identifies the critical conditions for the emergence of spatial and spatiotemporal structures. Finally, numerical simulations are used to illustrate nonlinear dynamics and pattern formation in the nonlocal model.

Genetic diversity of hepatitis B virus quasispecies in different biological compartments reveals distinct genotypes.

The selection pressure imposed by the host immune system impacts hepatitis B virus (HBV) quasispecies variability. This study evaluates HBV genetic diversity in different biological fluids. Twenty paired serum, oral fluid, and DBS samples from chronic HBV carriers were analyzed using both Sanger and next generation sequencing (NGS). The mean HBV viral load in serum was 5.19 ± 4.3 log IU/mL (median 5.29, IQR 3.01-7.93). Genotype distribution was: HBV/A1 55% (11/20), A2 15% (3/20), D3 10% (2/20), F2 15% (3/20), and F4 5% (1/20). Genotype agreement between serum and oral fluid was 100% (genetic distances 0.0-0.006), while that between serum and DBS was 80% (genetic distances 0.0-0.115). Two individuals presented discordant genotypes in serum and DBS. Minor population analysis revealed a mixed population. All samples displayed mutations in polymerase and/or surface genes. Major population analysis of the polymerase pointed to positions H122 and M129 as the most polymorphic (≥ 75% variability), followed by V163 (55%) and I253 (50%). Neither Sanger nor NGS detected any antiviral primary resistance mutations in the major populations. Minor population analysis, however, demonstrated the rtM204I resistance mutation in all individuals, ranging from 2.8 to 7.5% in serum, 2.5 to 6.3% in oral fluid, and 3.6 to 7.2% in DBS. This study demonstrated that different fluids can be used to assess HBV diversity, nonetheless, genotypic differences according to biological compartments can be observed.

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies.

BACKGROUND: HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering. METHODS: A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance. RESULTS: In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/µL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases. CONCLUSIONS: Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

Ongoing HIV-1 evolution and reservoir reseeding in two elite controllers with genetically diverse peripheral proviral quasispecies.

BACKGROUND: Elite controllers (EC) are human immunodeficiency virus (HIV)-positive individuals who can maintain low viral loads for extended periods without antiretroviral therapy due to multifactorial and individual characteristics. Most have a small HIV-1 reservoir composed of identical proviral sequences maintained by clonal expansion of infected CD4+ T cells. However, some have a more diverse peripheral blood mononuclear cell (PBMC)-associated HIV-1 reservoir with unique sequences. OBJECTIVES: To understand the turnover dynamics of the PBMC-associated viral quasispecies in ECs with relatively diverse circulating proviral reservoirs. METHODS: We performed single genome amplification of the env gene at three time points during six years in two EC with high intra-host HIV DNA diversity. FINDINGS: Both EC displayed quite diverse PBMCs-associated viral quasispecies (mean env diversity = 1.9-4.1%) across all time-points comprising both identical proviruses that are probably clonally expanded and unique proviruses with evidence of ongoing evolution. HIV-1 env glycosylation pattern suggests that ancestral and evolving proviruses may display different phenotypes of resistance to broadly neutralising antibodies consistent with persistent immune pressure. Evolving viruses may progressively replace the ancestral ones or may remain as minor variants in the circulating proviral population. MAIN CONCLUSIONS: These findings support that the high intra-host HIV-1 diversity of some EC resulted from long-term persistence of archival proviruses combined with the continuous reservoir's reseeding and low, but measurable, viral evolution despite undetectable viremia.

Traceability of SARS-CoV-2 transmission through quasispecies analysis.

During COVID-19 pandemic, consensus genomic sequences were used for rapidly monitor the spread of the virus worldwide. However, less attention was paid to intrahost genetic diversity. In fact, in the infected host, SARS-CoV-2 consists in an ensemble of replicating and closely related viral variants so-called quasispecies. Here we show that intrahost single nucleotide variants (iSNVs) represent a target for contact tracing analysis. Our data indicate that in the acute phase of infection, in highly likely transmission links, the number of viral particles transmitted from one host to another (bottleneck size) is large enough to propagate iSNVs among individuals. Furthermore, we demonstrate that, during SARS-CoV-2 outbreaks when the consensus sequences are identical, it is possible to reconstruct the transmission chains by genomic investigations of iSNVs. Specifically, we found that it is possible to identify transmission chains by limiting the analysis of iSNVs to only three well-conserved genes, namely nsp2, ORF3, and ORF7.