Viral quasispecies refers to the complex and dynamic mutant distributions (also termed mutant spectra, clouds, or swarms) that arise as a result of high error rates during RNA genome replication. The mutant spectrum of individual RNA virus populations is modified by continuous generation of variant genomes, competition and interactions among them, environmental influences, bottleneck events, and bloc transmission of viral particles. Quasispecies dynamics provides a new perspective on how viruses adapt, evolve, and cause disease, and sheds light on strategies to combat them. Molecular flexibility, together with ample opportunity of mutant cloud traffic in our global world, are key ingredients of viral disease emergences, as exemplified by the recent COVID-19 pandemic. In the present article, we present a brief overview of the molecular basis of mutant swarm formation and dynamics, and how the latter relates to viral disease and epidemic spread. We outline future challenges derived of the highly diverse cellular world in which viruses are necessarily installed.

Defective genomes are part of SARS-CoV-2 quasispecies. High-resolution, ultra-deep sequencing of bulk RNA from viral populations does not distinguish RNA mutations, insertions, and deletions in viable genomes from those in defective genomes. To quantify SARS-CoV-2 infectious variant progeny, virus from four individual plaques (biological clones) of a preparation of isolate USA-WA1/2020, formed on Vero E6 cell monolayers, was subjected to further biological cloning to yield 9 second-generation and 15 third-generation sub-clones. Consensus genomic sequences of the biological clones and sub-clones included an average of 2.8 variations per viable genome, relative to the consensus sequence of the parental USA-WA1/2020 virus. This value is 6.5-fold lower than the estimates for biological clones of other RNA viruses such as bacteriophage Qβ, foot-and-mouth disease virus, or hepatitis C virus in cell culture. The mutant spectrum complexity of the nsp12 (polymerase)- and spike (S)-coding region was unique in the progeny of each of 10 third-generation sub-clones; they shared 2.4% of the total of 164 different mutations and deletions scored in the 3,719 genomic residues that were screened. The presence of minority out-of-frame deletions revealed the ease of defective genome production from an individual infectious genome. Several low-frequency point mutations and deletions were clade-discordant in that they were not typical of USA-WA1/2020 but served to define the consensus sequences of future SARS-CoV-2 clades. Implications for SARS-CoV-2 adaptability and COVID-19 control of the viable genome heterogeneity and the generation of complex mutant spectra from individual genomes are discussed. IMPORTANCE: Sequencing of biological clones is a means to identify mutations, insertions, and deletions located in viable genomes. This distinction is particularly important for viral populations, such as those of SARS-CoV-2, that contain large proportions of defective genomes. By sequencing biological clones and sub-clones, we quantified the heterogeneity of the viable complement of USA-WA1/2020 to be lower than exhibited by other RNA viruses. This difference may be due to a reduced mutation rate or to limited tolerance of the large coronavirus genome to incorporate mutations and deletions and remain functional or a combination of both influences. The presence of clade-discordant residues in the progeny of individual biological sub-clones suggests limitations in the occupation of sequence space by SARS-CoV-2. However, the complex and unique mutant spectra that are rapidly generated from individual genomes suggest an aptness to confront selective constraints.

Understanding how viral mutant spectra organize and explore genotype space is essential for elucidating the mechanisms that drive molecular evolution. Here, we use deep-sequencing data of an amplicon in the A2 protein of the RNA bacteriophage Q β to reconstruct genotype networks comprising tens of thousands of haplotypes. The study of populations evolved under different temperature regimes reveals robust and reproducible patterns that arise from the interplay between fundamental geometrical motifs of sequence spaces and population dynamics. Mutant swarms exhibit a self-similar, hierarchical organization in which sequences cluster around highly connected, abundant cores that continuously regenerate diversity during evolution. The immediate neighborhood of these cores is rapidly rebuilt and extensively sampled, while a few mutations away sampling becomes dynamical and sparse. This population structure emerges from a dynamic, out-of-equilibrium balance between replication and mutational exploration and suggests that Q β populations do not rely primarily on neutral networks to navigate genotype space or to generate diversity. Combining genotype networks from populations adapted to different temperatures reveals early evolutionary divergence, with partially overlapping yet distinct populations that remain connected through short mutational paths. Even at the time scale of these experiments, evolutionary trajectories remain multiple, preventing the backward reconstruction of unique trajectories once mutations have been fixed. Together, this analysis provides a detailed view of the fine-scale processes shaping the evolution of heterogeneous viral populations and establishes genotype networks as a powerful framework for visualizing and interpreting the organization and diversification of viral quasispecies.

Understanding how ecological factors shape viral evolution is essential for predicting adaptation in RNA viruses. In this study, we investigated the evolutionary dynamics of bacteriophage Qβ under varying host densities, focusing on two nonsynonymous mutations—A1930G and C2011A—located in the A1 protein. Using experimental evolution, phenotypic assays, and competition experiments, we found that C2011A is consistently selected at low bacterial densities, enhancing viral entry but reducing burst size. In contrast, A1930G is fixed at high densities, despite similar phenotypic effects, suggesting its advantage arises from interactions with additional mutations. Clonal analysis revealed that compensatory or beneficial mutations modulate the fitness of A1930G, enabling its fixation. The absence of both mutations in the same genome points to negative epistasis, confirmed by the poor performance of the double mutant generated by site-directed mutagenesis. Sequencing of intermediate transfers showed early emergence of A1930G, but its fixation was prevented by clonal interference with C2011A. These findings highlight how host availability, fitness trade-offs, epistasis, and competition among variants shape the adaptive landscape of RNA viruses.

RNA viruses form genetically diverse populations structured as mutant spectra, or quasispecies, whose internal organization influences their evolutionary and adaptive dynamics. While genetic diversity has been extensively characterized, the structural organization of viral populations in sequence space remains less explored. Here, we compare genotype network architectures in two RNA viruses with markedly different evolutionary contexts: bacteriophage $Q\beta$ evolving in controlled laboratory conditions and severe acute respiratory syndrome coronavirus 2 evolving within infected human hosts. Using deep sequencing data, we reconstruct the genotype network of mutationally coupled variants within viral populations and analyse their topological properties. Despite large differences in genome size, mutation rate, and ecological setting, both viruses exhibit a common organization: a highly abundant central haplotype surrounded by layers of variants of diminishing abundance as Hamming distance to the central haplotype increases. All reconstructed networks share qualitative and quantitative topological features, displaying a hierarchical structure. The robust organization of both populations under multiple conditions suggests that RNA viruses may share a common genotype network architecture governed by fundamental properties of sequence space and the generic mechanisms of replication and mutation. Genotype networks provide a unifying framework to describe viral population structure beyond conventional diversity measures and, by revealing how local constraints shape mutational search, offers insights into the predictability of viral evolution.