Neutral quasispecies evolution and the maximal entropy random walk

Even if they have no impact on phenotype, neutral mutations are not equivalent in the eyes of evolution: A robust neutral variant—one which remains functional after further mutations—is more likely to spread in a large, diverse population than a fragile one. Quasispecies theory shows that the equilibrium frequency of a genotype is proportional to its eigenvector centrality in the neutral network. This paper explores the link between the selection for mutational robustness and the navigability of neutral networks. I show that sequences of neutral mutations follow a “maximal entropy random walk,” a canonical Markov chain on graphs with nonlocal, nondiffusive dynamics. I revisit M. Smith’s word-game model of evolution in this light, finding that the likelihood of certain sequences of substitutions can decrease with the population size. These counterintuitive results underscore the fertility of the interface between evolutionary dynamics, information theory, and physics.

Compartmentalized Replication of SARS-Cov-2 in Upper vs. Lower Respiratory Tract Assessed by Whole Genome Quasispecies Analysis

We report whole-genome and intra-host variability of SARS-Cov-2 assessed by next generation sequencing (NGS) in upper (URT) and lower respiratory tract (LRT) from COVID-19 patients. The aim was to identify possible tissue-specific patterns and signatures of variant selection for each respiratory compartment. Six patients, admitted to the Intensive Care Unit, were included in the study. Thirteen URT and LRT were analyzed by NGS amplicon-based approach on Ion Torrent Platform. Bioinformatic analysis was performed using both realized in-house and supplied by ThermoFisher programs. Phylogenesis showed clade V clustering of the first patients diagnosed in Italy, and clade G for later strains. The presence of quasispecies was observed, with variants uniformly distributed along the genome and frequency of minority variants spanning from 1% to ~30%. For each patient, the patterns of variants in URT and LRT were profoundly different, indicating compartmentalized virus replication. No clear variant signature and no significant difference in nucleotide diversity between LRT and URT were observed. SARS-CoV-2 presents genetic heterogeneity and quasispecies compartmentalization in URT and LRT. Intra-patient diversity was low. The pattern of minority variants was highly heterogeneous and no specific district signature could be identified, nevertheless, analysis of samples, longitudinally collected in patients, supported quasispecies evolution.

Evolution of Avian coronavirus (AvCoV) in BHK-21 and VERO cells

Avian coronavirus (AvCoV) infects a range of tissues in chickens and several other avian species. Although the virus can be isolated in chicken embryos, only a few strains of the 6 genotypes/33 lineages can grow in cell lines, with the Beaudette strain (GI-1 lineage) being the most used for in vitro studies. Considering the differences between cell lines and chicken embryos as habitats for AvCoV, this study aimed to assess the diversity of the genes coding for the nonstructural protein 3 (nsp3) and spike envelope protein (S) after serial passages in BHK-21 and Vero cells. After 14 passages of an embryo-adapted Beaudette strain, the virus loads fluctuated in both cell lines, with the highest loads being 8.72 log genome copies/µL for Vero and 6.36 log genome copies/µL for BHK-21 cells. No polymorphisms were found for nsp3; regarding S, not only aa substitutions (Vero: 8th passage A150S, and 14th S150A; BHK-21: 4th S53F, 8th F53Y, and 8th S95R), but also minor variants could be detected on chromatograms with fluctuating intensities. As the regions of these aa substitutions are within the receptor-binding domain of S, it can be speculated that differences in cell receptors between Vero and BHK-21 cells and the speed of cell death led to the selection of different dominant strains, while the stability of nsp3 supports its function as a protease involved in AvCoV replication. In conclusion, AvCoV quasispecies evolution is influenced by the biological model under consideration, and a gradual transition is seen for minor and major variants.

Comparison of Serum Hepatitis B Virus RNA Levels and Quasispecies Evolution Patterns between Entecavir and Pegylated-Interferon Mono-treatment in Chronic Hepatitis B Patients

ABSTRACT Hepatitis B virus (HBV) RNA may independently predict virological and serological response. This study aimed to compare dynamic changes in serum HBV RNA levels and HBV quasispecies evolution patterns between entecavir and pegylated-interferon mono-treatment in chronic hepatitis B patients and to determine the clinical significance during treatment. TaqMan real-time PCR was used for quantitative analysis. HBV RNA levels were retrospectively determined in serial serum samples from 178 chronic hepatitis B patients who received either entecavir or pegylated-interferon treatment. Both serum HBV DNA and RNA quasispecies were analyzed via next-generation sequencing. Receiver operating characteristics (ROC) analysis was performed to evaluate the prediction value of individual biomarkers for hepatitis B e antigen (HBeAg) seroconversion. Patients who received pegylated-interferon treatment showed stronger declines in HBV RNA levels than did those who received entecavir treatment. Serum HBV RNA levels were lower in patients with subsequent HBeAg seroconversion. At baseline, the level of HBV RNA was better than other indicators in predicting HBeAg seroconversion. Moreover, the predictive value of serum HBV RNA levels was better in the entecavir group. Baseline HBV RNA exhibited a significantly higher genetic diversity than HBV DNA and had a significant decline after 4 weeks of entecavir treatment. Higher baseline genetic diversity may result in a better outcome in pegylated-interferon-treated patients. Serum HBV RNA levels showed different decline kinetics, and HBV RNA quasispecies showed different evolution patterns in entecavir and pegylated-interferon mono-treatment. Taken together, serum HBV RNA may serve as a promising biomarker of HBeAg seroconversion in patients during antiviral treatment.

Nonlocal Reaction-Diffusion Model of Viral Evolution: Emergence of Virus Strains

The work is devoted to the investigation of virus quasispecies evolution and diversification due to mutations, competition for host cells, and cross-reactive immune responses. The model consists of a nonlocal reaction-diffusion equation for the virus density depending on the genotype considered as a continuous variable and on time. This equation contains two integral terms corresponding to the nonlocal effects of virus interaction with host cells and with immune cells. In the model, a virus strain is represented by a localized solution concentrated around some given genotype. Emergence of new strains corresponds to a periodic wave propagating in the space of genotypes. The conditions of appearance of such waves and their dynamics are described.

Characterization of EIAV env Quasispecies during Long-Term Passage In Vitro: Gradual Loss of Pathogenicity

As the only widely used live lentiviral vaccine, the equine infectious anima virus (EIAV) attenuated vaccine was developed by in vitro passaging of a virulent strain for 121 generations. In our previous study, we observed that the attenuated vaccine was gradually selected under increased environmental pressure at the population level (termed a quasispecies). To further elucidate the potential correlation between viral quasispecies evolution and pathogenesis, a systematic study was performed by sequencing env using several methods. Some key mutations were identified within Env, and we observed that increased percentages of these mutations were accompanied by an increased passage number and attenuated virulence. Phylogenetic analysis revealed that env mutations related to the loss of virulence might have occurred evolutionarily. Among these mutations, deletion of amino acid 236 in the V4 region of Env resulted in the loss of one N-glycosylation site that was crucial for virulence. Notably, the 236-deleted sequence represented a “vaccine-specific” mutation that was also found in wild EIAVLN40 strains based on single genome amplification (SGA) analysis. Therefore, our results suggest that the EIAV attenuated vaccine may originate from a branch of quasispecies of EIAVLN40. Generally, the presented results may increase our understanding of the attenuation mechanism of the EIAV vaccine and provide more information about the evolution of other lentiviruses.