The runaway evolution of SARS-CoV-2 leading to the highly evolved Delta strain

In new epidemics after the host shift, the pathogens may experience accelerated evolution driven by novel selective pressures. When the accelerated evolution enters a positive feedback loop with the expanding epidemics, the pathogen's runaway evolution may be triggered. To test this possibility in COVID-19, we analyze the extensive databases and identify 5 major waves of strains, one replacing the previous one in 2020-2021. The mutations differ entirely between waves and the number of mutations continues to increase, from 3-4 to 21-31. The latest wave is the Delta strain which accrues 31 new mutations to become highly prevalent. Interestingly, these new mutations in Delta strain emerge in multiple stages with each stage driven by 6-12 coding mutations that form a fitness group. In short, the evolution of SARS-CoV-2 from the oldest to the youngest wave, and from the earlier to the later stages of the Delta wave, is a process of acceleration with more and more mutations. The global increase in the viral population size (M(t), at time t) and the mutation accumulation (R(t)) may have indeed triggered the runaway evolution in late 2020, leading to the highly evolved Alpha and then Delta strain. To suppress the pandemic, it is crucial to break the positive feedback loop between M(t) and R(t), neither of which has yet to be effectively dampened by late 2021. New waves beyond Delta, hence, should not be surprising.

Viral population genomics reveals host and infectivity impact on SARS-CoV-2 adaptive landscape

Public health surveillance, drug treatment development, and optimization of immunological interventions all depend on understanding pathogen adaptation, which differ for specific pathogens. SARS-CoV-2 is an exceptionally successful human pathogen, yet complete understanding of the forces driving its evolution is lacking. Here, we leveraged almost four million SARS-CoV-2 sequences originating mostly from non-vaccinated naive patients to investigate the impact of functional constraints and natural immune pressures on the sequence diversity of the SARS-CoV-2 genome. Overall, we showed that the SARS-CoV-2 genome is under strong and intensifying levels of purifying selection with a minority of sites under diversifying pressure. With a particular focus on the spike protein, we showed that sites under selection were critical for protein stability and virus fitness related to increased infectivity and/or reduced neutralization by convalescent sera. We investigated the genetic diversity of SARS-CoV-2 B and T cell epitopes and determined that the currently known T cell epitope sequences were highly conserved. Outside of the spike protein, we observed that mutations under selection in variants of concern can be associated to beneficial outcomes for the virus. Altogether, the results yielded a comprehensive map of all sites under selection across the entirety of SARS-CoV-2 genome, highlighting targets for future studies to better understand the virus spread, evolution and success.

Spatiotemporal Associations and Molecular Evolution of Highly Pathogenic Avian Influenza A H7N9 Virus in China from 2017 to 2021

Highly pathogenic (HP) H7N9 avian influenza virus (AIV) emerged in China in 2016. HP H7N9 AIV caused at least 33 human infections and has been circulating in poultry farms continuously since wave 5. The genetic divergence, geographic patterns, and hemagglutinin adaptive and parallel molecular evolution of HP H7N9 AIV in China since 2017 are still unclear. Here, 10 new strains of HP H7N9 AIVs from October 2019 to April 2021 were sequenced. We found that HP H7N9 was primarily circulating in Northern China, particularly in the provinces surrounding the Bohai Sea (Liaoning, Hebei, and Shandong) since wave 6. Of note, HP H7N9 AIV phylogenies exhibit a geographical structure compatible with high levels of local transmission after unidirectional rapid geographical expansion towards the north of China in 2017. In addition, we showed that two major subclades were continually expanding with the viral population size undergoing a sharp increase after 2018 with an obvious seasonal tendency. Notably, the hemagglutinin gene showed signs of parallel evolution and positive selection. Our research sheds light on the current epidemiology, evolution, and diversity of HP H7N9 AIV that can help prevent and control the spreading of HP H7N9 AIV.

Porcine Reproductive and Respiratory Syndrome (PRRS) Epidemiology in an Integrated Pig Company of Northern Italy: A Multilevel Threat Requiring Multilevel Interventions

Porcine reproductive and respiratory syndrome (PRRS) is probably the most relevant viral disease affecting pig farming. Despite the remarkable efforts paid in terms of vaccination administration and biosecurity, eradication and long-term control have often been frustrated. Unfortunately, few studies are currently available that objectively link, using a formal statistical approach, viral molecular epidemiology to the risk factors determining the observed scenario. The purpose of the present study is to contribute to filling this knowledge gap taking advantage of the advancements in the field of phylodynamics. Approximately one-thousand ORF7 sequences were obtained from strains collected between 2004 and 2021 from the largest Italian pig company, which implements strict compartmentalization among independent three-sites (i.e., sow herds, nurseries and finishing units) pig flows. The history and dynamics of the viral population and its evolution over time were reconstructed and linked to managerial choices. The viral fluxes within and among independent pig flows were evaluated, and the contribution of other integrated pig companies and rurally risen pigs in mediating such spreading was investigated. Moreover, viral circulation in Northern Italy was reconstructed using a continuous phylogeographic approach, and the impact of several environmental features on PRRSV strain persistence and spreading velocity was assessed. The results demonstrate that PRRSV epidemiology is shaped by a multitude of factors, including pig herd management (e.g., immunization strategy), implementation of strict-independent pig flows, and environmental features (e.g., climate, altitude, pig density, road density, etc.) among the others. Small farms and rurally raised animals also emerged as a potential threat for larger, integrated companies. These pieces of evidence suggest that none of the implemented measures can be considered effective alone, and a multidimensional approach, ranging from individual herd management to collaboration and information sharing among different companies, is mandatory for effective infection control.

Spontaneous Mutation in the Movement Protein of Citrus Leprosis Virus C2, in a Heterologous Virus Infection Context, Increases Cell-to-Cell Transport and Generates Fitness Advantage

Previous results using a movement defective alfalfa mosaic virus (AMV) vector revealed that citrus leprosis virus C (CiLV-C) movement protein (MP) generates a more efficient local movement, but not more systemic transport, than citrus leprosis virus C2 (CiLV-C2) MP, MPs belonging to two important viruses for the citrus industry. Here, competition experiment assays in transgenic tobacco plants (P12) between transcripts of AMV constructs expressing the cilevirus MPs, followed by several biological passages, showed the prevalence of the AMV construct carrying the CiLV-C2 MP. The analysis of AMV RNA 3 progeny recovered from P12 plant at the second viral passage revealed the presence of a mix of progeny encompassing the CiLV-C2 MP wild type (MPWT) and two variants carrying serines instead phenylalanines at positions 72 (MPS72F) or 259 (MPS259F), respectively. We evaluated the effects of each modified residue in virus replication, and cell-to-cell and long-distance movements. Results indicated that phenylalanine at position 259 favors viral cell-to-cell transport with an improvement in viral fitness, but has no effect on viral replication, whereas mutation at position 72 (MPS72F) has a penalty in the viral fitness. Our findings indicate that the prevalence of a viral population may be correlated with its greater efficiency in cell-to-cell and systemic movements.

Encapsidation of Staufen-2 Enhances Infectivity of HIV-1

Staufen, the RNA-binding family of proteins, affects various steps in the Human Immuno-Deficiency Virus (HIV-1) replication cycle. While our previous study established Staufen-2–HIV-1 Rev interaction and its role in augmenting nucleocytoplasmic export of RRE-containing viral RNA, viral incorporation of Staufen-2 and its effect on viral propagation were unknown. Here, we report that Staufen-2 interacts with HIV-1 Gag and is incorporated into virions and that encapsidated Staufen-2 boosted viral infectivity. Further, Staufen-2 gets co-packaged into virions, possibly by interacting with host factors Staufen-1 or antiviral protein APOBEC3G, which resulted in different outcomes on the infectivity of Staufen-2-encapsidated virions. These observations suggest that encapsidated host factors influence viral population dynamics and infectivity. With the explicit identification of the incorporation of Staufen proteins into HIV-1 and other retroviruses, such as Simian Immunodeficiency Virus (SIV), we propose that packaging of RNA binding proteins, such as Staufen, in budding virions of retroviruses is probably a general phenomenon that can drive or impact the viral population dynamics, infectivity, and evolution.

The mucosal barrier and anti-viral immune responses can eliminate portions of the viral population during transmission and early viral growth

Little is known about how specific individual viral lineages replicating systemically during acute Human Immunodeficiency Virus or Simian Immunodeficiency Virus (HIV/SIV) infection persist into chronic infection. In this study, we use molecularly barcoded SIV (SIVmac239M) to track distinct viral lineages for 12 weeks after intravenous (IV) or intrarectal (IR) challenge in macaques. Two Mafa-A1*063+ cynomolgus macaques (Macaca fascicularis, CM) were challenged IV, and two Mamu-A1*001+ rhesus macaques (Macaca mulatta, RM) were challenged IR with 200,000 Infectious Units (IU) of SIVmac239M. We sequenced the molecular barcode of SIVmac239M from all animals over the 12 weeks of the study to characterize the diversity and persistence of virus lineages. During the first three weeks post-infection, we found ~70–560 times more unique viral lineages circulating in the animals challenged IV compared to those challenged IR, which is consistent with the hypothesis that the challenge route is the primary driver restricting the transmission of individual viral lineages. We also characterized the sequences of T cell epitopes targeted during acute SIV infection, and found that the emergence of escape variants in acutely targeted epitopes can occur on multiple virus templates simultaneously, but that elimination of some of these templates is likely a consequence of additional host factors. These data imply that virus lineages present during acute infection can still be eliminated from the systemic virus population even after initial selection.

Phylodynamic and Recombination Analyses of Avian Infectious Bronchitis GI-23 Reveal a Widespread Recombinant Cluster and New Among-Countries Linkages

Infectious bronchitis virus GI-23 lineage, although described approximately two decades ago in the Middle East, has recently drawn remarkable attention and is considered an “emerging” lineage due to its current spread to several other regions, including Europe. Despite the relevance, no comprehensive studies are available investigating its epidemiologic and evolutionary pattern. The present phylodynamic study was designed to fill this gap, benefitting from a collection of freely available GI-23 sequences and ad-hoc generated European ones. After a relatively ancient origin in the Middle East, likely in the first half of the previous century, GI-23 circulated largely undetected or underdiagnosed for a long time in this region, likely causing little damage, potentially because of low virulence coupled with limited development of avian industry in the considered years and regions and insufficient diagnostic activity. The following development of the poultry industry and spread to other countries led to a progressive but slow increase of viral population size between the late ‘90s and 2010. An increase in viral virulence could also be hypothesized. Of note, a big recombinant cluster, likely originating in the Middle East but spreading thereafter, especially to Europe through Turkey, demonstrated a much-marked increase in viral population size compared to previously circulating variants. The extensive available GI-23 sequence datasets allowed to demonstrate several potential epidemiological links among African, Asian, and European countries, not described for other IBV lineages. However, differently from previously investigated IBV lineages, its spread appears to primarily involve neighbouring countries and those with strong economic and political relationships. It could thus be speculated that frequent effective contacts among locations are necessary for efficient strain transmission. Some countries appear to play a major role as a “bridge” among less related locations, being Turkey the most relevant example. The role of vaccination in controlling the viral population was also tentatively evaluated. However, despite some evidence suggesting such an effect, the bias in sequence and data availability and the variability in the applied vaccination protocols prevent robust conclusions and warrant further investigations.

Influence of GlycoGag on the Incorporation of Host Membrane Proteins Into the Envelope of the Moloney Murine Leukemia Virus

Analysis of viral particle heterogeneity produced from infected cells has been limited by the inefficiency of traditional analytical methods to characterize large populations of viruses at an individual particle level. Flow virometry (FVM) is an emerging technique based on flow cytometry principles that enables a high throughput, multiparametric, and phenotypic characterization of viruses at a single particle resolution. Here, we performed FVM to analyze surface markers found on Murine Leukemia Virus (MLV) and glycosylated Gag-deficient (glycoGag) MLV. The glycoGag viral accessory protein has several roles in the MLV viral infection cycle including directing retroviral assembly and particle release at lipid rafts. Based on previous studies, we hypothesize that glycoGag modulates host protein incorporation into the viral envelope during viral assembly and budding. Here, by using FVM, we reveal that glycoGag is associated with an increased incorporation of the host-derived tetraspanins CD81 and CD63 along with the lipid raft marker and immune antigen Thy1.2 during the assembly and release of viral particles from infected NIH 3T3, EL4, and primary CD4+ T cells. Moreover, we show differences in the uptake of host proteins by viruses that are released from the two cell lines and primary T lymphocytes. Additionally, at the individual viral particle level, we observed a degree of expression heterogeneity of host-derived antigens within the viral population. Finally, certain cellular antigens can show either enrichment or exclusion from the viral envelope depending on whether glycoGag is expressed by the virus. This suggests that glycoGag is involved in a mechanism of selective host protein incorporation into the viral envelope.

Replication and transmission of an influenza A(H3N2) virus harboring the polymerase acidic I38T substitution, in guinea pigs.

The polymerase acidic (PA) I38T substitution is a dominant marker of resistance to baloxavir. We evaluated the impact of I38T on the fitness of a contemporary influenza A(H3N2) virus. Influenza A/Switzerland/9715293/2013 (H3N2) wild-type (WT) virus and its I38T mutant were rescued by reverse genetics. Replication kinetics were compared using ST6GalI-MDCK and A549 cells and infectivity/contact transmissibility were evaluated in guinea pigs. Nasal wash (NW) viral titres were determined by TCID50 ml−1 in ST6GalI-MDCK cells. Competition experiments were performed and the evolution of viral population was assessed by droplet digital RT-PCR. I38T did not alter in vitro replication. I38T induced comparable titres vs the WT in guinea pigs NWs and the two viruses transmitted equally by direct contact. However, a 50 %:50 % mixture inoculum evolved to mean WT/I38T ratios of 71 %:29 % and 66.4 %:33.6 % on days 4 and 6 p.i., respectively. Contemporary influenza A(H3N2)-I38T PA variants may conserve a significant level of viral fitness.