Simplexviruses Successfully Adapt to Their Host by Fine-Tuning Immune Responses.

Primate herpes simplex viruses are species-specific and relatively harmless to their natural hosts. However, cross-species transmission is often associated with severe disease, as exemplified by the virulence of macacine herpesvirus 1 (B virus) in humans. We performed a genome-wide scan for signals of adaptation of simplexviruses to their hominin hosts. Among core genes, we found evidence of episodic positive selection in three glycoproteins, with several selected sites located in antigenic determinants. Positively selected noncore genes were found to be involved in different immune-escape mechanisms. The herpes simplex virus (HSV)-1/HSV-2 encoded product (ICP47) of one of these genes is known to down-modulate major histocompatibility complex class I expression. This feature is not shared with B virus, which instead up-regulates Human Leukocyte Antigen (HLA)-G, an immunomodulatory molecule. By in vitro expression of different ICP47 mutants, we functionally characterized the selection signals. Results indicated that the selected sites do not represent the sole determinants of binding to the transporter associated with antigen processing (TAP). Conversely, the amino acid status at these sites was sufficient to determine HLA-G up-regulation. In fact, both HSV-1 and HSV-2 ICP47 induced HLA-G when mutated to recapitulate residues in B virus, whereas the mutated version of B virus ICP47 failed to determine HLA-G expression. These differences might contribute to the severity of B virus infection in humans. Importantly, they indicate that the evolution of ICP47 in HSV-1/HSV-2 led to the loss of an immunosuppressive effect. Thus, related simplexviruses finely tune the balance between immunosuppressive and immunostimulatory pathways to promote successful co-existence with their primate hosts.

Whole-genome sequencing of hMPXV1 in five Italian cases confirms the occurrence of the predominant epidemic lineage.

The ongoing outbreak of monkeypox virus (hMPXV1) is the largest recorded in historically nonendemic countries. Genomic surveillance has emerged as a pivotal tool to track the spread and monitor the evolution of viral pathogens. Therefore, to assess the genetic diversity of circulating hMPXV1 in northern Italy in June to July 2022, we sequenced and analyzed five complete genomes of viruses sampled from patients presenting with a typical course of hMPXV1 infection. Phylogenetic analysis confirmed that all five genomes belong to the predominant epidemic lineage (B.1). Inspection of genetic changes and comparison with the reference sequence showed the presence of 12 nucleotide substitutions. Seven are nonsynonymous mutations leading to amino acid changes in six proteins belonging to different functional classes. Moreover, 11 of these 12 nucleotide mutations involve GA>AA or TC>TT replacements, suggesting that host APOBEC3 enzymes are responsible for the generation of substitutions in circulating viruses. Finally, metagenomic analysis evidenced bacterial superinfection (Streptococcus pyogenes) in one patient. Through this study, we contributed to expand the number of complete genomes of viruses circulating in Italy and characterize them as belonging to the predominant outbreak lineage.

Past and ongoing adaptation of human cytomegalovirus to its host.

Cytomegaloviruses (order Herpesvirales) display remarkable species-specificity as a result of long-term co-evolution with their mammalian hosts. Human cytomegalovirus (HCMV) is exquisitely adapted to our species and displays high genetic diversity. We leveraged information on inter-species divergence of primate-infecting cytomegaloviruses and intra-species diversity of clinical isolates to provide a genome-wide picture of HCMV adaptation across different time-frames. During adaptation to the human host, core viral genes were commonly targeted by positive selection. Functional characterization of adaptive mutations in the primase gene (UL70) indicated that selection favored amino acid replacements that decrease viral replication in human fibroblasts, suggesting evolution towards viral temperance. HCMV intra-species diversity was largely governed by immune system-driven selective pressure, with several adaptive variants located in antigenic domains. A significant excess of positively selected sites was also detected in the signal peptides (SPs) of viral proteins, indicating that, although they are removed from mature proteins, SPs can contribute to viral adaptation. Functional characterization of one of these SPs indicated that adaptive variants modulate the timing of cleavage by the signal peptidase and the dynamics of glycoprotein intracellular trafficking. We thus used evolutionary information to generate experimentally-testable hypotheses on the functional effect of HCMV genetic diversity and we define modulators of viral phenotypes.

Homology-based classification of accessory proteins in coronavirus genomes uncovers extremely dynamic evolution of gene content.

Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super-groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike-like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion-associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross-communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.

Antigenic variation of SARS-CoV-2 in response to immune pressure.

Analysis of the bat viruses most closely related to SARS-CoV-2 indicated that the virus probably required limited adaptation to spread in humans. Nonetheless, since its introduction in human populations, SARS-CoV-2 must have been subject to the selective pressure imposed by the human immune system. We exploited the availability of a large number of high-quality SARS-CoV-2 genomes, as well as of validated epitope predictions, to show that B cell epitopes in the spike glycoprotein (S) and in the nucleocapsid protein (N) have higher diversity than nonepitope positions. Similar results were obtained for other human coronaviruses and for sarbecoviruses sampled in bats. Conversely, in the SARS-CoV-2 population, epitopes for CD4+ and CD8+ T cells were not more variable than nonepitope positions. A significant reduction in epitope variability was instead observed for some of the most immunogenic proteins (S, N, ORF8 and ORF3a). Analysis over longer evolutionary time frames indicated that this effect is not due to differential constraints. These data indicate that SARS-CoV-2 evolves to elude the host humoral immune response, whereas recognition by T cells is not actively avoided by the virus. However, we also found a trend of lower diversity of T cell epitopes for common cold coronaviruses, indicating that epitope conservation per se is not directly linked to disease severity. We suggest that conservation serves to maintain epitopes that elicit tolerizing T cell responses or induce T cells with regulatory activity.

Role of NEU3 Overexpression in the Prediction of Efficacy of EGFR-Targeted Therapies in Colon Cancer Cell Lines.

The epidermal growth factor receptor (EGFR), through the MAP kinase and PI3K-Akt-mTOR axis, plays a pivotal role in colorectal cancer (CRC) pathogenesis. The membrane-associated NEU3 sialidase interacts with and desialylates EGFR by promoting its dimerization and downstream effectors' activation. Among the targeted therapies against EGFR, the monoclonal antibody cetuximab is active only in a subgroup of patients not carrying mutations in the MAP kinase pathway. In order to better understand the EGFR-NEU3 interplay and the mechanisms of pharmacological resistance, we investigated the role of NEU3 deregulation in cetuximab-treated CRC cell lines transiently transfected with NEU3 using Western blot analysis. Our results indicate that NEU3 overexpression can enhance EGFR activation only if EGFR is overexpressed, indicating the existence of a threshold for NEU3-mediated EGFR activation. This enhancement mainly leads to the constitutive activation of the MAP kinase pathway. Consequently, we suggest that the evaluation of NEU3 expression cannot entirely substitute the evaluation of EGFR because EGFR-negative cases cannot be stimulated by NEU3. Furthermore, NEU3-mediated hyperactivation of EGFR is counterbalanced by the administration of cetuximab, hypothesizing that a combined treatment of NEU3- and EGFR-targeted therapies may represent a valid option for CRC patients, which must be investigated in the future.

Extensive Positive Selection Drives the Evolution of Nonstructural Proteins in Lineage C Betacoronaviruses.

UNLABELLED: Middle East respiratory syndrome-related coronavirus (MERS-CoV) spreads to humans via zoonotic transmission from camels. MERS-CoV belongs to lineage C of betacoronaviruses (betaCoVs), which also includes viruses isolated from bats and hedgehogs. A large portion of the betaCoV genome consists of two open reading frames (ORF1a and ORF1b) that are translated into polyproteins. These are cleaved by viral proteases to generate 16 nonstructural proteins (nsp1 to nsp16) which compose the viral replication-transcription complex. We investigated the evolution of ORF1a and ORF1b in lineage C betaCoVs. Results indicated widespread positive selection, acting mostly on ORF1a. The proportion of positively selected sites in ORF1a was much higher than that previously reported for the surface-exposed spike protein. Selected sites were unevenly distributed, with nsp3 representing the preferential target. Several pairs of coevolving sites were also detected, possibly indicating epistatic interactions; most of these were located in nsp3. Adaptive evolution at nsp3 is ongoing in MERS-CoV strains, and two selected sites (G720 and R911) were detected in the protease domain. While position 720 is variable in camel-derived viruses, suggesting that the selective event does not represent a specific adaptation to humans, the R911C substitution was observed only in human-derived MERS-CoV isolates, including the viral strain responsible for the recent South Korean outbreak. It will be extremely important to assess whether these changes affect host range or other viral phenotypes. More generally, data herein indicate that CoV nsp3 represents a major selection target and that nsp3 sequencing should be envisaged in monitoring programs and field surveys. IMPORTANCE: Both severe acute respiratory syndrome coronavirus (SARS-CoV) and MERS-CoV originated in bats and spread to humans via an intermediate host. This clearly highlights the potential for coronavirus host shifting and the relevance of understanding the molecular events underlying the adaptation to new host species. We investigated the evolution of ORF1a and ORF1b in lineage C betaCoVs and in 87 sequenced MERS-CoV isolates. Results indicated widespread positive selection, stronger in ORF1a than in ORF1b. Several selected sites were found to be located in functionally relevant protein regions, and some of them corresponded to functional mutations in other coronaviruses. The proportion of selected sites we identified in ORF1a is much higher than that for the surface-exposed spike protein. This observation suggests that adaptive evolution in ORF1a might contribute to host shifts or immune evasion. Data herein also indicate that genetic diversity at nonstructural proteins should be taken into account when antiviral compounds are developed.

A common genetic variant in FOXP2 is associated with language-based learning (dis)abilities: Evidence from two Italian independent samples.

Language-based Learning Disabilities (LLDs) encompass a group of complex, comorbid, and developmentally associated deficits in communication. Language impairment and developmental dyslexia (DD) represent the most recognized forms of LLDs. Substantial genetic correlations exist between language and reading (dis)abilities. Common variants in the FOXP2 gene were consistently associated with language- and reading-related neuropsychological and neuroanatomical phenotypes. We tested the effect of a FOXP2 common variant, that is, rs6980093 (A/G), on quantitative measures of language and reading in two independent Italian samples: a population-based cohort of 699 subjects (3-11 years old) and a sample of 572 children with DD (6-18 years old). rs6980093 modulates expressive language in the general population sample, with an effect on fluency scores. In the DD sample, the variant showed an association with the accuracy in the single word reading task. rs6980093 shows distinct genetic models of association in the two cohorts, with a dominant effect of the G allele in the general population sample and heterozygote advantage in the DD cohort. We provide preliminary evidence that rs6980093 associates with language and reading (dis)abilities in two independent Italian cohorts. rs6980093 is an intronic SNP, suggesting that it (or a linked variant) modulates phenotypic association via regulation of FOXP2 expression. Because FOXP2 brain expression is finely regulated, both temporally and spatially, it is possible that the two alleles at rs6980093 differentially modulate expression levels in a developmental stage- or brain area-specific manner. This might help explaining the heterozygote advantage effect and the different genetic models in the two cohorts.

The mammalian complement system as an epitome of host-pathogen genetic conflicts.

The complement system is an innate immunity effector mechanism; its action is antagonized by a wide array of pathogens and complement evasion determines the virulence of several infections. We investigated the evolutionary history of the complement system and of bacterial-encoded complement-interacting proteins. Complement components targeted by several pathogens evolved under strong selective pressure in primates, with selection acting on residues at the contact interface with microbial/viral proteins. Positively selected sites in CFH and C4BPA account for the human specificity of gonococcal infection. Bacterial interactors, evolved adaptively as well, with selected sites located at interaction surfaces with primate complement proteins. These results epitomize the expectation under a genetic conflict scenario whereby the host's and the pathogen's genes evolve within binding avoidance-binding seeking dynamics. In silico mutagenesis and protein-protein docking analyses supported this by showing that positively selected sites, both in the host's and in the pathogen's interacting partner, modulate binding.

NEU3 activity enhances EGFR activation without affecting EGFR expression and acts on its sialylation levels.

Several studies performed over the last decade have focused on the role of sialylation in the progression of cancer and, in particular, on the association between deregulation of sialidases and tumorigenic transformation. The plasma membrane-associated sialidase NEU3 is often deregulated in colorectal cancer (CRC), and it was shown that this enzyme co-immunoprecipitates in HeLa cells with epidermal growth factor receptor (EGFR), the molecular target of most recent monoclonal antibody-based therapies against CRC. To investigate the role of NEU3 sialidase on EGFR deregulation in CRC, we first collected data on NEU3 gene expression levels from a library of commercial colon cell lines, demonstrating that NEU3 transcription is upregulated in these cell lines. We also found EGFR to be hyperphosphorylated in all cell lines, with the exception of SW620 cells and the CCD841 normal intestinal cell line. By comparing the effects induced by overexpression of either the wild-type or the inactive mutant form of NEU3 on EGFR, we demonstrated that the active form of NEU3 enhanced receptor activation without affecting EGFR mRNA or protein expression. Moreover, through western blots and mass spectrometry analysis, we found that EGFR immunoprecipitated from cells overexpressing active NEU3, unlike the receptor from mock cells and cells overexpressing inactive NEU3, is desialylated. On the whole, our data demonstrate that, besides the already reported indirect EGFR activation through GM3, sialidase NEU3 could also play a role on EGFR activation through its desialylation.

Diverse selective regimes shape genetic diversity at ADAR genes and at their coding targets.

A-to-I RNA editing operated by ADAR enzymes is extremely common in mammals. Several editing events in coding regions have pivotal physiological roles and affect protein sequence (recoding events) or function. We analyzed the evolutionary history of the 3 ADAR family genes and of their coding targets. Evolutionary analysis indicated that ADAR evolved adaptively in primates, with the strongest selection in the unique N-terminal domain of the interferon-inducible isoform. Positively selected residues in the human lineage were also detected in the ADAR deaminase domain and in the RNA binding domains of ADARB1 and ADARB2. During the recent history of human populations distinct variants in the 3 genes increased in frequency as a result of local selective pressures. Most selected variants are located within regulatory regions and some are in linkage disequilibrium with eQTLs in monocytes. Finally, analysis of conservation scores of coding editing sites indicated that editing events are counter-selected within regions that are poorly tolerant to change. Nevertheless, a minority of recoding events occurs at highly conserved positions and possibly represents the functional fraction. These events are enriched in pathways related to HIV-1 infection and to epidermis/hair development. Thus, both ADAR genes and their targets evolved under variable selective regimes, including purifying and positive selection. Pressures related to immune response likely represented major drivers of evolution for ADAR genes. As for their coding targets, we suggest that most editing events are slightly deleterious, although a minority may be beneficial and contribute to antiviral response and skin homeostasis.

Looking at human cytosolic sialidase NEU2 structural features with an interdisciplinary approach.

Circular dichroism (CD) spectra at variable temperatures have been recorded for human cytosolic sialidase NEU2 in buffered water solutions and in the presence of divalent cations. The results show the prevalence of ß-strands together with a considerable amount of α-helical structure, while in the solid state, from both previous X-ray diffraction analysis and our CD data on film samples, the content of ß-strands is higher. In solution, a significant change in CD spectra occurs with an increase in temperature, related to a decrease in α-helix content and a slight increase in ß-strand content. In the same range of temperatures, the enzymatic activity decreases. Although the overall structure of the protein appears to be particularly stable, molecular dynamics simulations performed at various temperatures evidence local conformational changes possibly relevant for explaining the relative lability of enzymatic activity.

A proline-rich loop mediates specific functions of human sialidase NEU4 in SK-N-BE neuronal differentiation.

Human sialidase NEU4 long (N4L) is a membrane-associated enzyme that has been shown to be localized in the outer mitochondrial membrane. A role in different cellular processes has been suggested for this enzyme, such as apoptosis, neuronal differentiation and tumorigenesis. However, the molecular bases for these roles, not found in any of the other highly similar human sialidases, are not understood. We have found that a proline-rich sequence of 81 amino acids, unique to NEU4 sequence, contains potential Akt and Erk1 kinase motifs. Molecular modeling, based on the experimentally determined three-dimensional structure of cytosolic human NEU2, showed that the proline-rich sequence is accommodated in a loop, thus preserving the typical beta-barrel structure of sialidases. In order to investigate the role of this loop in neuronal differentiation, we obtained SK-N-BE neuroblastoma cells stably overexpressing either human wild-type N4L or a deletion mutant lacking the proline-rich loop. Our results demonstrate that the proline-rich region can also enhance cell proliferation and retinoic acid (RA)-induced neuronal differentiation and it is also involved in NEU4 interaction with Akt, as well as in substrate recognition, modifying directly or through the interaction with other protein(s) the enzyme specificity toward sialylated glycoprotein(s). On the whole, our results suggest that N4L could be a downstream component of the PI3K/Akt signaling pathway required for RA-induced differentiation of neuroblastoma SK-N-BE cells.

Evolution of the orthopoxvirus core genome.

Orthopoxviruses comprise several relevant pathogens, including the causative agent of smallpox and monkeypox virus. Analysis of orthopoxvirus genome evolution mainly focused on gene gains/losses. We instead analyzed core genes, which are conserved in all orthopoxviruses. We show that, despite their strong constraint, some genes involved in viral morphogenesis and transcription/replication were targets of pervasive positive selection, which was relatively uncommon in immunomodulatory genes. However at least three of the positively selected genes, E3L, A24R, and H3L, might have evolved in response to immune selection. Episodic positive selection was particularly common on the internal branches of the orthopox phylogeny and on the monkeypox virus lineage. The latter showed evidence of episodic positive selection at the D14L gene, which encodes a modulator of complement activation (MOPICE). Notably, two genes (B1R and A33R) targeted by episodic selection on more than one branch are involved in forms of intra-genomic conflict. Finally, we found that, in orthopoxvirus proteomes, intrinsically disordered regions (IDRs) tend to be less constrained and are common targets of positive selection. Extension of our analysis to all poxviruses showed no evidence that the IDR fraction differs with host range. Conversely, we found a strong effect of base composition, which was however not sufficient to explain IDR fraction. We thus suggest that, in poxviruses, the IDR fraction is maintained by modulating GC content to accommodate disorder-promoting codons. Overall, our data provide novel insight in orthopoxvirus evolution and provide a list of genes and sites that are expected to modulate viral phenotypes.

SARS-CoV-2 ORF3c impairs mitochondrial respiratory metabolism, oxidative stress, and autophagic flux.

Coronaviruses encode a variable number of accessory proteins that are involved in host-virus interaction, suppression of immune responses, or immune evasion. SARS-CoV-2 encodes at least twelve accessory proteins, whose roles during infection have been studied. Nevertheless, the role of the ORF3c accessory protein, an alternative open reading frame of ORF3a, has remained elusive. Herein, we show that the ORF3c protein has a mitochondrial localization and alters mitochondrial metabolism, inducing a shift from glucose to fatty acids oxidation and enhanced oxidative phosphorylation. These effects result in increased ROS production and block of the autophagic flux. In particular, ORF3c affects lysosomal acidification, blocking the normal autophagic degradation process and leading to autolysosome accumulation. We also observed different effect on autophagy for SARS-CoV-2 and batCoV RaTG13 ORF3c proteins; the 36R and 40K sites are necessary and sufficient to determine these effects.

Molecular insight into substrate recognition by human cytosolic sialidase NEU2.

Sialidases or neuramidases are glycoside hydrolases removing terminal sialic acid residues from sialo-glycoproteins and sialo-glycolipids. Viral neuraminidases (NAs) have been extensively characterized and represent an excellent target for antiviral therapy through the synthesis of a series of competitive inhibitors that block the release of newly formed viral particles from infected cells. The human cytosolic sialidase NEU2 is the only mammalian enzyme structurally characterized and represents a valuable model to study the specificity of novel NA inhibitory drugs. Moreover, the availability of NEU2 3D structure represents a pivotal step toward the characterization of the molecular basis of natural substrates recognition by the enzyme. In this perspective, we have carried out a study of molecular docking of NEU2 active site using natural substrates of increasing complexity. Moreover, selective mutations of the residues putatively involved into substrate(s) interaction/recognition have been performed, and the resulting mutant enzymes have been preliminary tested for their catalytic activity and substrate specificity. We found that Q270 is involved in the binding of the disaccharide α(2,3) sialyl-galactose, whereas K45 and Q112 bind the distal glucose of the trisaccharide α(2,3) sialyl-lactose, corresponding to the oligosaccharide moiety of GM3 ganglioside. In addition, E218, beside D46, is proved to be a key catalytic residue, being, together with Y334, the second member of the nucleophile pair required for the catalysis. Overall, our results point out the existence of a dynamic network of interactions that are possibly involved in the recognition of the glycans bearing sialic acid.

Molecular cloning of soluble trehalase from Chironomus riparius larvae, its heterologous expression in Escherichia coli and bioinformatic analysis.

Trehalase is involved in the control of trehalose concentration, the main blood sugar in insects. Here, we describe the molecular cloning of the cDNA encoding for the soluble form of the trehalase from the midge larvae of Chironomus riparius, a well-known bioindicator of the quality of freshwater environments. Molecular cloning was achieved through multiple alignment of Diptera trehalase sequences, allowing the synthesis of internal homology-based primers; the complete open reading frame(ORF) was subsequently obtained through RACE-PCR(where RACE is rapid amplification of cDNA ends). The cDNA contained the 5' untranslated region (UTR), the 3' UTR including a poly(A) tail and the ORF of 1,725 bp consisting of 574 amino acid residues with a predicted molecular mass of 65,778 Da. Recombinant trehalase was successfully expressed in Escherichia coli as a His-tagged protein and purified on Ni-NTA affinity chromatography. Primary structure analysis showed a series of characteristic features shared by all insect trehalases, while three-dimensional structure prediction yielded the typical glucosidase fold, the two key residues involved in the catalytic mechanism being conserved. Production of recombinant insect trehalases opens the way to structural characterizations of the catalytic site, which might represent, among others, an element for reconsidering the enzyme as a target in pest insects' control.

Dating the Emergence of Human Endemic Coronaviruses.

Four endemic coronaviruses infect humans and cause mild symptoms. Because previous analyses were based on a limited number of sequences and did not control for effects that affect molecular dating, we re-assessed the timing of endemic coronavirus emergence. After controlling for recombination, selective pressure, and molecular clock model, we obtained similar tMRCA (time to the most recent common ancestor) estimates for the four coronaviruses, ranging from 72 (HCoV-229E) to 54 (HCoV-NL63) years ago. The split times of HCoV-229E and HCoV-OC43 from camel alphacoronavirus and bovine coronavirus were dated ~268 and ~99 years ago. The split times of HCoV-HKU1 and HCoV-NL63 could not be calculated, as their zoonoticic sources are unknown. To compare the timing of coronavirus emergence to that of another respiratory virus, we recorded the occurrence of influenza pandemics since 1500. Although there is no clear relationship between pandemic occurrence and human population size, the frequency of influenza pandemics seems to intensify starting around 1700, which corresponds with the initial phase of exponential increase of human population and to the emergence of HCoV-229E. The frequency of flu pandemics in the 19th century also suggests that the concurrence of HCoV-OC43 emergence and the Russian flu pandemic may be due to chance.

Adaptation of the endemic coronaviruses HCoV-OC43 and HCoV-229E to the human host.

Four coronaviruses (HCoV-OC43, HCoV-HKU1, HCoV-NL63, and HCoV-229E) are endemic in human populations. All these viruses are seasonal and generate short-term immunity. Like the highly pathogenic coronaviruses, the endemic coronaviruses have zoonotic origins. Thus, understanding the evolutionary dynamics of these human viruses might provide insight into the future trajectories of SARS-CoV-2 evolution. Because the zoonotic sources of HCoV-OC43 and HCoV-229E are known, we applied a population genetics-phylogenetic approach to investigate which selective events accompanied the divergence of these viruses from the animal ones. Results indicated that positive selection drove the evolution of some accessory proteins, as well as of the membrane proteins. However, the spike proteins of both viruses and the hemagglutinin-esterase (HE) of HCoV-OC43 represented the major selection targets. Specifically, for both viruses, most positively selected sites map to the receptor-binding domains (RBDs) and are polymorphic. Molecular dating for the HCoV-229E spike protein indicated that RBD Classes I, II, III, and IV emerged 3-9 years apart. However, since the appearance of Class V (with much higher binding affinity), around 25 years ago, limited genetic diversity accumulated in the RBD. These different time intervals are not fully consistent with the hypothesis that HCoV-229E spike evolution was driven by antigenic drift. An alternative, not mutually exclusive possibility is that strains with higher affinity for the cellular receptor have out-competed strains with lower affinity. The evolution of the HCoV-OC43 spike protein was also suggested to undergo antigenic drift. However, we also found abundant signals of positive selection in HE. Whereas such signals might result from antigenic drift, as well, previous data showing co-evolution of the spike protein with HE suggest that optimization for human cell infection also drove the evolution of this virus. These data provide insight into the possible trajectories of SARS-CoV-2 evolution, especially in case the virus should become endemic.

Evolution and Genetic Diversity of Primate Cytomegaloviruses.

Cytomegaloviruses (CMVs) infect many mammals, including humans and non-human primates (NHPs). Human cytomegalovirus (HCMV) is an important opportunistic pathogen among immunocompromised patients and represents the most common infectious cause of birth defects. HCMV possesses a large genome and very high genetic diversity. NHP-infecting CMVs share with HCMV a similar genomic organization and coding content, as well as the course of viral infection. Recent technological advances have allowed the sequencing of several HCMV strains from clinical samples and provided insight into the diversity of NHP-infecting CMVs. The emerging picture indicates that, with the exclusion of core genes (genes that have orthologs in all herpesviruses), CMV genomes are relatively plastic and diverse in terms of gene content, both at the inter- and at the intra-species level. Such variability most likely underlies the strict species-specificity of these viruses, as well as their ability to persist lifelong and with relatively little damage to their hosts. However, core genes, despite their strong conservation, also represented a target of adaptive evolution and subtle changes in their coding sequence contributed to CMV adaptation to different hosts. Indubitably, important knowledge gaps remain, the most relevant of which concerns the role of viral genetics in HCMV-associated human disease.