The Missing Expression Level-Evolutionary Rate Anticorrelation in Viruses Does Not Support Protein Function as a Main Constraint on Sequence Evolution

Population immunity enhances the evolution of SARS-CoV-2 in Beijing revealed by wastewater genomic surveillance

This study aims to elucidate the impact of population immunity on the regional evolution of SARS-CoV-2. A total of 3701 wastewater SARS-CoV-2 concentration values and 168 wastewater whole genomes of SARS-CoV-2 were obtained in Beijing over 11 months following the implementation of the "dynamic zero-COVID" policy adjustments in December 2022. The findings indicate that the number of variant strains identified through wastewater surveillance was 2.46 times greater than that detected by clinical monitoring, with single nucleotide polymorphisms showing an increase of up to 7.14 times. This enhanced surveillance facilitates a more comprehensive analysis of regional virus evolution patterns. Following the adjustment of epidemic measure, Beijing experienced three distinct waves of epidemics, and the dominant variant transitioned directly from BA.5 in the first wave to XBB after six months in the second one. During this period, strong population immunity formed by centralized infection in over 90 % of the population blocked the outbreak of internationally prevalent and concerning variants BQ.1 and CH.1.1, resulting in a 12.5 % faster regional evolution of SARS-CoV-2 strains in Beijing compared to the international context. Subsequently, in August 2023, EG.5 became the dominant variant in the third wave, aligning with international trends. The epidemics in Beijing have caused significant positive selection pressure on SARS-CoV-2 strains, favoring those with enhanced antigenic escape mutations in spike gene. These results underscore that the extensive infection after the adjustment of epidemic prevention policies has accelerated the evolution of SARS-CoV-2 in Beijing and been conducive to antigenic escape evolution, which can effectively inform decision making for epidemic control and preemptive vaccine design.

Research advances in mechanisms of arsenic hyperaccumulation of Pteris vittata: Perspectives from plant physiology, molecular biology, and phylogeny

Pteris vittata, as the firstly discovered arsenic (As) hyperaccumulator, has great application value in As-contaminated soil remediation. Currently, the genes involved in As hyperaccumulation in P. vittata have been mined continuously, while they have not been used in practice to enhance phytoremediation efficiency. Aiming to better assist the practice of phytoremediation, this review collects 130 studies to clarify the progress in research into the As hyperaccumulation process in P. vittata from multiple perspectives. Antioxidant defense, rhizosphere activities, vacuolar sequestration, and As efflux are important physiological activities involved in As hyperaccumulation in P. vittata. Among related 19 genes, PHT, TIP, ACR3, ACR2 and HAC family genes play essential roles in arsenate (AsⅤ) transport, arsenite (AsⅢ) transport, vacuole sequestration of AsⅢ, and the reduction of AsⅤ to AsⅢ, respectively. Gene ontology enrichment analysis indicated it is necessary to further explore genes that can bind to related ions, with transport activity, or with function of transmembrane transport. Phylogeny analysis results implied ACR2, HAC and ACR3 family genes with rapid evolutionary rate may be the decisive factors for P. vittata as an As hyperaccumulator. A deeper understanding of the As hyperaccumulation network and key gene components could provide useful tools for further bio-engineered phytoremediation.

Genetic ancestry and population structure of vaccinia virus

Vaccinia virus (VACV) was used for smallpox eradication, but its ultimate origin remains unknown. The genetic relationships among vaccine stocks are also poorly understood. We analyzed 63 vaccine strains with different origin, as well horsepox virus (HPXV). Results indicated the genetic diversity of VACV is intermediate between variola and cowpox viruses, and that mutation contributed more than recombination to VACV evolution. STRUCTURE identified 9 contributing subpopulations and showed that the lowest drift was experienced by the ancestry components of Tian Tan and HPXV/Mütter/Mulford genomes. Subpopulations that experienced very strong drift include those that contributed the ancestry of MVA and IHD-W, in good agreement with the very long passage history of these vaccines. Another highly drifted population contributed the full ancestry of viruses sampled from human/cattle infections in Brazil and, partially, to IOC clones, strongly suggesting that the recurrent infections in Brazil derive from the spillback of IOC to the feral state.

Ongoing Positive Selection Drives the Evolution of SARS-CoV-2 Genomes

SARS-CoV-2 is a new RNA virus affecting humans and spreads extensively through world populations since its first outbreak in December, 2019. Whether the transmissibility and pathogenicity of SARS-CoV-2 in humans after zoonotic transfer are actively evolving, and driven by adaptation to the new host and environments is still under debate. Understanding the evolutionary mechanism underlying epidemiological and pathological characteristics of COVID-19 is essential for predicting the epidemic trend, and providing guidance for disease control and treatments. Interrogating novel strategies for identifying natural selection using within-species polymorphisms and 3,674,076 SARS-CoV-2 genome sequences of 169 countries as of December 30, 2021, we demonstrate with population genetic evidence that during the course of SARS-CoV-2 pandemic in humans, 1) SARS-CoV-2 genomes are overall conserved under purifying selection, especially for the 14 genes related to viral RNA replication, transcription, and assembly; 2) ongoing positive selection is actively driving the evolution of 6 genes (e.g., S, ORF3a, and N) that play critical roles in molecular processes involving pathogen-host interactions, including viral invasion into and egress from host cells, and viral inhibition and evasion of host immune response, possibly leading to high transmissibility and mild symptom in SARS-CoV-2 evolution. According to an established haplotype phylogenetic relationship of 138 viral clusters, a spatial and temporal landscape of 556 critical mutations is constructed based on their divergence among viral haplotype clusters or repeatedly increase in frequency within at least 2 clusters, of which multiple mutations potentially conferring alterations in viral transmissibility, pathogenicity, and virulence of SARS-CoV-2 are highlighted, warranting attentions.

Evidence for a mouse origin of the SARS-CoV-2 Omicron variant

The rapid accumulation of mutations in the SARS-CoV-2 Omicron variant that enabled its outbreak raises questions as to whether its proximal origin occurred in humans or another mammalian host. Here, we identified 45 point mutations that Omicron acquired since divergence from the B.1.1 lineage. We found that the Omicron spike protein sequence was subjected to stronger positive selection than that of any reported SARS-CoV-2 variants known to evolve persistently in human hosts, suggesting a possibility of host-jumping. The molecular spectrum of mutations (i.e., the relative frequency of the 12 types of base substitutions) acquired by the progenitor of Omicron was significantly different from the spectrum for viruses that evolved in human patients but resembled the spectra associated with virus evolution in a mouse cellular environment. Furthermore, mutations in the Omicron spike protein significantly overlapped with SARS-CoV-2 mutations known to promote adaptation to mouse hosts, particularly through enhanced spike protein binding affinity for the mouse cell entry receptor. Collectively, our results suggest that the progenitor of Omicron jumped from humans to mice, rapidly accumulated mutations conducive to infecting that host, then jumped back into humans, indicating an inter-species evolutionary trajectory for the Omicron outbreak.

Host-specific asymmetric accumulation of mutation types reveals that the origin of SARS-CoV-2 is consistent with a natural process

The capacity of RNA viruses to adapt to new hosts and rapidly escape the host immune system is largely attributable to de novo genetic diversity that emerges through mutations in RNA. Although the molecular spectrum of de novo mutations-the relative rates at which various base substitutions occur-are widely recognized as informative toward understanding the evolution of a viral genome, little attention has been paid to the possibility of using molecular spectra to infer the host origins of a virus. Here, we characterize the molecular spectrum of de novo mutations for SARS-CoV-2 from transcriptomic data obtained from virus-infected cell lines, enabled by the use of sporadic junctions formed during discontinuous transcription as molecular barcodes. We find that de novo mutations are generated in a replication-independent manner, typically on the genomic strand, and highly dependent on mutagenic mechanisms specific to the host cellular environment. De novo mutations will then strongly influence the types of base substitutions accumulated during SARS-CoV-2 evolution, in an asymmetric manner favoring specific mutation types. Consequently, similarities between the mutation spectra of SARS-CoV-2 and the bat coronavirus RaTG13, which have accumulated since their divergence strongly suggest that SARS-CoV-2 evolved in a host cellular environment highly similar to that of bats before its zoonotic transfer into humans. Collectively, our findings provide data-driven support for the natural origin of SARS-CoV-2.

G-quadruplexes in genomes of viruses infecting eukaryotes or prokaryotes are under different selection pressures from hosts

G-quadruplexes in viral genomes can be applied as the targets of antiviral therapies, which has attracted wide interest. However, it is still not clear whether the pervasive number of such elements in the viral world is the result of natural selection for functionality. In this study, we identified putative quadruplex-forming sequences (PQSs) across the known viral genomes and analyzed the abundance, structural stability, and conservation of viral PQSs. A Viral Putative G-quadruplex Database (ViPGD,http://jsjds.hzau.edu.cn/MBPC/ViPGD/index.php/home/index) was constructed to collect the details of each viral PQS, which provides guidance for selecting the desirable PQS. The PQS with two putative G-tetrads (G₂-PQS) was significantly enriched in both eukaryotic viruses and prokaryotic viruses, while the PQSs with three putative G-tetrads (G₃-PQS) were only enriched in eukaryotic viruses and depleted in prokaryotic viruses. The structural stability of PQSs in prokaryotic viruses was significantly lower than that in eukaryotic viruses. Conservation analysis showed that the G₂-PQS, instead of G₃-PQS, was highly conserved within the genus. This suggested that the G₂-quadruplex might play an important role in viral biology, and the difference in the occurrence of G-quadruplex between eukaryotic viruses and prokaryotic viruses may result from the different selection pressures from hosts.