Evolution of Virus-like Features and Intrinsically Disordered Regions in Retrotransposon-derived Mammalian Genes.

Several mammalian genes have originated from the domestication of retrotransposons, selfish mobile elements related to retroviruses. Some of the proteins encoded by these genes have maintained virus-like features; including self-processing, capsid structure formation, and the generation of different isoforms through -1 programmed ribosomal frameshifting. Using quantitative approaches in molecular evolution and biophysical analyses, we studied 28 retrotransposon-derived genes, with a focus on the evolution of virus-like features. By analyzing the rate of synonymous substitutions, we show that the -1 programmed ribosomal frameshifting mechanism in three of these genes (PEG10, PNMA3, and PNMA5) is conserved across mammals and originates alternative proteins. These genes were targets of positive selection in primates, and one of the positively selected sites affects a B-cell epitope on the spike domain of the PNMA5 capsid, a finding reminiscent of observations in infectious viruses. More generally, we found that retrotransposon-derived proteins vary in their intrinsically disordered region content and this is directly associated with their evolutionary rates. Most positively selected sites in these proteins are located in intrinsically disordered regions and some of them impact protein posttranslational modifications, such as autocleavage and phosphorylation. Detailed analyses of the biophysical properties of intrinsically disordered regions showed that positive selection preferentially targeted regions with lower conformational entropy. Furthermore, positive selection introduces variation in binary sequence patterns across orthologues, as well as in chain compaction. Our results shed light on the evolutionary trajectories of a unique class of mammalian genes and suggest a novel approach to study how intrinsically disordered region biophysical characteristics are affected by evolution.

Depletion of CpG dinucleotides in bacterial genomes may represent an adaptation to high temperatures.

Dinucleotide biases have been widely investigated in the genomes of eukaryotes and viruses, but not in bacteria. We assembled a dataset of bacterial genomes (>15 000), which are representative of the genetic diversity in the kingdom Eubacteria, and we analyzed dinucleotide biases in relation to different traits. We found that TpA dinucleotides are the most depleted and that CpG dinucleotides show the widest dispersion. The abundances of both dinucleotides vary with genomic G + C content and show a very strong phylogenetic signal. After accounting for G + C content and phylogenetic inertia, we analyzed different bacterial lifestyle traits. We found that temperature preferences associate with the abundance of CpG dinucleotides, with thermophiles/hyperthemophiles being particularly depleted. Conversely, the TpA dinucleotide displays a bias that only depends on genomic G + C composition. Using predictions of intrinsic cyclizability we also show that CpG depletion may associate with higher DNA bendability in both thermophiles/hyperthermophiles and mesophiles, and that the former are predicted to have significantly more flexible genomes than the latter. We suggest that higher bendability is advantageous at high temperatures because it facilitates DNA positive supercoiling and that, through modulation of DNA mechanical properties, local or global CpG depletion controls genome organization, most likely not only in bacteria.

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts.

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Dinucleotide biases in RNA viruses that infect vertebrates or invertebrates.

IMPORTANCE: Akin to a molecular signature, dinucleotide composition can be exploited by the zinc-finger antiviral protein (ZAP) to restrict CpG-rich (and UpA-rich) RNA viruses. ZAP evolved in tetrapods, and it is not encoded by invertebrates and fish. Because a systematic analysis is missing, we analyzed the genomes of RNA viruses that infect vertebrates or invertebrates. We show that vertebrate single-stranded (ss) RNA(+) viruses and, to a lesser extent, double-stranded RNA viruses tend to have stronger CpG bias than invertebrate viruses. Conversely, ssRNA(-) viruses have similar dinucleotide composition whether they infect vertebrates or invertebrates. Analysis of ssRNA(+) viruses that infect mammals, reptiles, and fish indicated that ZAP is unlikely to be a major driver of CpG depletion. We also show that, compared to other coronaviruses, the genome of SARS-CoV-2 is not homogeneously CpG-depleted. Our study provides new insights into virus evolution and strategies for recoding RNA virus genomes.

Selective events at individual sites underlie the evolution of monkeypox virus clades.

In endemic regions (West Africa and the Congo Basin), the genetic diversity of monkeypox virus (MPXV) is geographically structured into two major clades (Clades I and II) that differ in virulence and host associations. Clade IIb is closely related to the B.1 lineage, which is dominating a worldwide outbreak initiated in 2022. Lineage B.1 has however accumulated mutations of unknown significance that most likely result from apolipoprotein B mRNA editing catalytic polypeptide-like 3 (APOBEC3) editing. We applied a population genetics-phylogenetics approach to investigate the evolution of MPXV during historical viral spread in Africa and to infer the distribution of fitness effects. We observed a high preponderance of codons evolving under strong purifying selection, particularly in viral genes involved in morphogenesis and replication or transcription. However, signals of positive selection were also detected and were enriched in genes involved in immunomodulation and/or virulence. In particular, several genes showing evidence of positive selection were found to hijack different steps of the cellular pathway that senses cytosolic DNA. Also, a few selected sites in genes that are not directly involved in immunomodulation are suggestive of antibody escape or other immune-mediated pressures. Because orthopoxvirus host range is primarily determined by the interaction with the host immune system, we suggest that the positive selection signals represent signatures of host adaptation and contribute to the different virulence of Clade I and II MPXVs. We also used the calculated selection coefficients to infer the effects of mutations that define the predominant human MPXV1 (hMPXV1) lineage B.1, as well as the changes that have been accumulating during the worldwide outbreak. Results indicated that a proportion of deleterious mutations were purged from the predominant outbreak lineage, whose spread was not driven by the presence of beneficial changes. Polymorphic mutations with a predicted beneficial effect on fitness are few and have a low frequency. It remains to be determined whether they have any significance for ongoing virus evolution.

An APOBEC3 Mutational Signature in the Genomes of Human-Infecting Orthopoxviruses.

The ongoing worldwide monkeypox outbreak is caused by viral lineages (globally referred to as hMPXV1) that are related to but distinct from clade IIb MPXV viruses transmitted within Nigeria. Analysis of the genetic differences has indicated that APOBEC-mediated editing might be responsible for the unexpectedly high number of mutations observed in hMPXV1 genomes. Here, using 1,624 publicly available hMPXV1 sequences, we analyzed the mutations that accrued between 2017 and the emergence of the current predominant variant (B.1), as well as those that that have been accumulating during the 2022 outbreak. We confirmed an overwhelming prevalence of C-to-T and G-to-A mutations, with a sequence context (5'-TC-3') consistent with the preferences of several human APOBEC3 enzymes. We also found that mutations preferentially occur in highly expressed viral genes, although no transcriptional asymmetry was observed. A comparison of the mutation spectrum and context was also performed against the human-specific variola virus (VARV) and the zoonotic cowpox virus (CPXV), as well as fowlpox virus (FWPV). The results indicated that in VARV genomes, C-to-T and G-to-A changes were more common than the opposite substitutions, although the effect was less marked than for hMPXV1. Conversely, no preference toward C-to-T and G-to-A changes was observed in CPXV and FWPV. Consistently, the sequence context of C-to-T changes confirmed a preference for a T in the -1 position for VARV, but not for CPXV or FWPV. Overall, our results strongly support the view that, irrespective of the transmission route, orthopoxviruses infecting humans are edited by the host APOBEC3 enzymes. IMPORTANCE Analysis of the viral lineages responsible for the 2022 monkeypox outbreak suggested that APOBEC enzymes are driving hMPXV1 evolution. Using 1,624 public sequences, we analyzed the mutations that accumulated between 2017 and the emergence of the predominant variant and those that characterize the last outbreak. We found that the mutation spectrum of hMPXV1 has been dominated by TC-to-TT and GA-to-AA changes, consistent with the editing activity of human APOBEC3 proteins. We also found that mutations preferentially affect highly expressed viral genes, possibly because transcription exposes single-stranded DNA (ssDNA), a target of APOBEC3 editing. Notably, analysis of the human-specific variola virus (VARV) and the zoonotic cowpox virus (CPXV) indicated that in VARV genomes, TC-to-TT and GA-to-AA changes are likewise extremely frequent. Conversely, no preference toward TC-to-TT and GA-to-AA changes is observed in CPXV. These results suggest that APOBEC3 proteins have an impact on the evolution of different human-infecting orthopoxviruses.

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.

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.

Serotonin transporter gene methylation and emotional regulation in preschool children born preterm: A longitudinal evaluation of the role of negative emotionality in infancy.

The aim of the study was to assess the contribution of negative emotionality at 3 months (T1) and serotonin transporter gene (SLC6A4) DNA methylation at 4.5 years of age (T2) to emotion regulation in pre-schoolers born very preterm and full-term. Forty one children (n = 21 born very preterm, n = 20 born full-term) participated in the study. Fretful behavior was assessed at T1 in response to the Face-to-FaceStill-Face (FFSF) paradigm. At T2, SLC6A4 DNA methylation was analyzed and emotion regulation was assessed using an observational procedure (i.e., the Pre-schooler Regulation of Emotional Stress, PRES). The very preterm group displayed higher emotion dysregulation during the PRES Reactivity phase than the full-term group. Higher levels of fretful behavior at 3 months were associated with greater emotional distress only for very preterm children with higher methylation at T2. No significant associations emerged in the full-term group. Despite current findings cannot be generalized owing to the relatively small sample size, this work provides preliminary longitudinal evidence about the link between negative emotionality during infancy, stress-linked epigenetic status at 4.5 years and emotion dysregulation in preschoolers born preterm.

El propósito del estudio fue evaluar la contribución de la emocionalidad negativa a los 3 meses (T1) y la metilación del ADN en el gen transportador de la serotonina (SLC6A4) a los 4 años y medio de edad (T2) a la regulación de la emoción en prescolares nacidos muy antes de la gestación completa o de gestación completa. Cuarenta y un niños (n = 21 nacidos muy antes de la gestación completa, n = 20 nacidos de gestación completa) participaron en el estudio. El comportamiento irritable se evaluó a T1 como respuesta al Cara-a-Cara del paradigma de la Cara Inmóvil (FFSF). A T2, se analizó la metilación de ADN SLC6A4 y se evaluó la regulación de la emoción usando un procedimiento de observación (v.g. La Regulación del Estrés Emocional del Prescolar, PRES). El grupo nacido muy antes de la gestación completa mostró una más alta desregulación durante la fase de Reactividad PRES que el grupo nacido de gestación completa. Los niveles más altos de comportamiento irritable a los 3 meses se asociaron con una mayor angustia emocional solamente para los niños nacidos muy antes de la gestación completa con más alta metilación al T2. Ninguna asociación significativa surgió del grupo nacido de gestación completa. A pesar de que los actuales resultados no se pueden generalizar debido al tamaño relativamente pequeño del grupo muestra, este trabajo ofrece aporta evidencia longitudinal preliminar acerca de la conexión entre la emocionalidad negativa durante la infancia, el estado epigenético relacionado con el estrés a los 4 años y medio y la desregulación de la emoción en prescolares nacidos antes de la completa gestación.

Le but de cette étude était d'évaluer la contribution de l'émotivité négative à 3 mois (T1) et du gène vecteur de la sérotonine (SLC6A4) méthylation de l'ADN à l'âge de 4,5 ans (T2) à la régulation de l'émotion chez les enfants d'âge préscolaire nés très prématurés et à plein terme. Quarante et un enfant (n = 21 nés très prématurés, n = 20 nés à plein terme) ont participé à l'étude. Le comportement agité a été évalué au T1 en réponse au paradigme face-à-face visage inexpressif (abrégé FFSF en anglais). Au T2, la méthylation de l'ADN SLC6A4 a été analysée et la régulation de l'émotion a été évaluée en utilisant un protocole d'observation (à savoir, la Régulation du Stress Emotionnel de l'Enfant d'Age Préscolaire, abrégé en anglais PRES). Le groupe très prématuré a fait état d'une dysrégulation de l'émotion plus élevée durant la phase de Réactivité PRES que le groupe né à plein terme. Des niveaux plus élevés de comportement agité à 3 mois étaient liés à une détresse émotionnelle plus grande uniquement pour les enfants très prématurés avec une méthylation plus élevée au T2. Aucune association importante n'a émergé dans le groupe à plein terme. En dépit du fait que les résultats actuels ne peuvent pas être généralisés à cause de la taille relativement petite de l'échantillon, ce travail offre des preuves longitudinales préliminaires sur le lien entre l'émotivité négative durant la petite enfant, le statut épigénétique lié au stress à 4,5 ans et la dysrégulation de l'émotion chez les enfants d'âge préscolaires nés avant terme.

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.

Kinetochore proteins and microtubule-destabilizing factors are fast evolving in eutherian mammals.

Centromeres have central functions in chromosome segregation, but centromeric DNA and centromere-binding proteins evolve rapidly in most eukaryotes. The selective pressure(s) underlying the fast evolution of centromere-binding proteins are presently unknown. An attractive possibility is that selfish centromeres promote their preferential inclusion in the oocyte and centromeric proteins evolve to suppress meiotic drive (centromere drive hypothesis). We analysed the selective patterns of mammalian genes that encode kinetochore proteins and microtubule (MT)-destabilizing factors. We show that several of these proteins evolve at the same rate or faster than proteins with a role in centromere specification. Elements of the kinetochore that bind MTs or that bridge the interaction between MTs and the centromere represented the major targets of positive selection. These data are in line with the possibility that the genetic conflict fuelled by meiotic drive extends beyond genes involved in centromere specification. However, we cannot exclude that different selective pressures underlie the rapid evolution of MT-destabilizing factors and kinetochore components. Whatever the nature of such pressures, they must have been constant during the evolution of eutherian mammals, as we found a surprisingly good correlation in dN/dS (ratio of the rate of nonsynonymous and synonymous substitutions) across orders/clades. Finally, when phylogenetic relationships were accounted for, we found little evidence that the evolutionary rates of these genes change with testes size, a proxy for sperm competition. Our data indicate that, in analogy to centromeric proteins, kinetochore components are fast evolving in mammals. This observation may imply that centromere drive plays out at multiple levels or that these proteins adapt to lineage-specific centromeric features.

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.

Hemodynamic and behavioral peculiarities in response to emotional stimuli in children with attention deficit hyperactivity disorder: An fNIRS study.

BACKGROUND: Children with attention deficit hyperactivity disorder (ADHD) exhibit behavioral inhibition deficits, which often lead to emotional dysregulation (ED) affecting individual ability to control emotions and behavioral responses. In ADHD, ED is associated with poor outcomes and comorbidities, with both externalizing and internalizing disorders. This work aims to evaluate sensitivity to emotional stimuli in children with ADHD using functional Near Infrared Spectroscopy (fNIRS). METHODS: During frontal fNIRS recording, 20 children with ADHD and 25 typically developing (TD) peers performed a visual continuous performance task with stimuli of different emotional content (i.e., positive, negative, neutral, and control stimuli without emotional content). This is a cognitive task designed to evaluate the ability to recognize emotional stimuli and to deal with emotional interference. RESULTS: The ADHD sample showed more variability in response time to stimuli and more false alarms compared to TD group. fNIRS data showed between-group differences in right prefrontal and frontal cortices, with wider hemoglobin concentration changes in the TD group, during positive, negative, and neutral conditions. LIMITATIONS: Owing to the limited possibility of near infrared light to penetrate tissue, fNIRS can only measure cortical activations, while it would be of interest to identify the subcortical areas linked to emotional processing, too. CONCLUSIONS: Findings suggest the presence of emotional processing deficits in children with ADHD, as suggested by poor performances on the e-CPT task, and of peculiar sensitivity to emotional stimuli, linked to atypical hemodynamics of right prefrontal and frontal areas.

Evolutionary analysis of exogenous and integrated HHV-6A/HHV-6B populations.

Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6A/HHV-6B, iciHHV-6A/iciHHV-6B). We reconstructed the population structures of HHV-6A and HHV-6B, showing that HHV-6A diverged less than HHV-6B genomes from the projected common ancestral population. Thus, HHV-6B genomes experienced stronger drift, as also supported by calculation of nucleotide diversity and Tajima's D. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, ancestry proportions were similar in exogenous HHV-6B viruses and iciHHV-6B. In line with previous indications, this suggests the distinct exogenous viral populations that originated iciHHV-6B in subjects with European and Asian ancestry are still causing infections in the corresponding geographic areas. Notably, for both iciHHV-6A and iciHHV-6B, we found that European and American sequences tend to have high proportions of ancestry from viral populations that experienced considerable drift, suggesting that they underwent one or more bottlenecks followed by population expansion. Finally, analysis of HHV-6B exogenous viruses sampled in Japan indicated that proportions of ancestry components of most of these viruses are different from the majority of those sampled in the USA. More generally, we show that, in both viral species, both integrated and exogenous viral genomes have different ancestry components, partially depending on geographic location. It would be extremely important to determine whether such differences account for the diversity of HHV-6A/HHV-6B-associated clinical symptoms and epidemiology. Also, the sequencing of additional exogenous and integrated viral genomes will be instrumental to confirm and expand our conclusions, which are based on a relatively small number of genomes, sequenced with variable quality, and with unequal sampling in terms of geographic origin.

Intrinsically disordered regions are abundant in simplexvirus proteomes and display signatures of positive selection.

Whereas the majority of herpesviruses co-speciated with their mammalian hosts, human herpes simplex virus 2 (HSV-2, genus Simplexvirus) most likely originated from the cross-species transmission of chimpanzee herpesvirus 1 to an ancestor of modern humans. We exploited the peculiar evolutionary history of HSV-2 to investigate the selective events that drove herpesvirus adaptation to a new host. We show that HSV-2 intrinsically disordered regions (IDRs)-that is, protein domains that do not adopt compact three-dimensional structures-are strongly enriched in positive selection signals. Analysis of viral proteomes indicated that a significantly higher portion of simplexvirus proteins is disordered compared with the proteins of other human herpesviruses. IDR abundance in simplexvirus proteomes was not a consequence of the base composition of their genomes (high G + C content). Conversely, protein function determines the IDR fraction, which is significantly higher in viral proteins that interact with human factors. We also found that the average extent of disorder in herpesvirus proteins tends to parallel that of their human interactors. These data suggest that viruses that interact with fast-evolving, disordered human proteins, in turn, evolve disordered viral interactors poised for innovation. We propose that the high IDR fraction present in simplexvirus proteomes contributes to their wider host range compared with other herpesviruses.

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.

Recent Out-of-Africa Migration of Human Herpes Simplex Viruses.

Herpes simplex virus types 1 and 2 (HSV-1 and HSV-2) are ubiquitous human pathogens. Both viruses evolved from simplex viruses infecting African primates and they are thus thought to have left Africa during early human migrations. We analyzed the population structure of HSV-1 and HSV-2 circulating strains. Results indicated that HSV-1 populations have limited geographic structure and the most evident clustering by geography is likely due to recent bottlenecks. For HSV-2, the only level of population structure is accounted for by the so-called "worldwide" and "African" lineages. Analysis of ancestry components and nucleotide diversity, however, did not support the view that the worldwide lineage followed early humans during out-of-Africa dispersal. Although phylogeographic analysis confirmed an African origin for both viruses, molecular dating with a method that corrects for the time-dependent rate phenomenon indicated that HSV-1 and HSV-2 migrated from Africa in relatively recent times. In particular, we estimated that the HSV-2 worldwide lineage left the continent in the 18th century, which corresponds to the height of the transatlantic slave trade, possibly explaining the high prevalence of HSV-2 in the Americas (second highest after Africa). The limited geographic clustering of HSV-1 makes it difficult to date its exit from Africa. The split between the basal clade, containing mostly African sequences, and all other strains was dated at ∼5,000 years ago. Our data do not imply that herpes simplex viruses did not infect early humans but show that the worldwide distribution of circulating strains is the result of relatively recent events.

You Will Never Walk Alone: Codispersal of JC Polyomavirus with Human Populations.

JC polyomavirus (JCPyV) is one of the most prevalent human viruses. Findings based on the geographic distribution of viral subtypes suggested that JCPyV codiverged with human populations. This view was however challenged by data reporting a much more recent origin and expansion of JCPyV. We collected information on ∼1,100 worldwide strains and we show that their geographic distribution roughly corresponds to major human migratory routes. Bayesian phylogeographic analysis inferred a Subsaharan origin for JCPyV, although with low posterior probability. High confidence inference at internal nodes provided strong support for a long-standing association between the virus and human populations. In line with these data, pairwise FST values for JCPyV and human mtDNA sampled from the same areas showed a positive and significant correlation. Likewise, very strong relationships were found when node ages in the JCPyV phylogeny were correlated with human population genetic distances (nuclear-marker based FST). Reconciliation analysis detected a significant cophylogenetic signal for the human population and JCPyV trees. Notably, JCPyV also traced some relatively recent migration events such as the expansion of people from the Philippines/Taiwan area into Remote Oceania, the gene flow between North-Eastern Siberian and Ainus, and the Koryak contribution to Circum-Arctic Americans. Finally, different molecular dating approaches dated the origin of JCPyV in a time frame that precedes human out-of-Africa migration. Thus, JCPyV infected early human populations and accompanied our species during worldwide dispersal. JCPyV typing can provide reliable geographic information and the virus most likely adapted to the genetic background of human populations.

Retraction to: A complex evolutionary relationship between HHV-6A and HHV-6B.

[This retracts the article DOI: 10.1093/ve/vez043.].

A complex evolutionary relationship between HHV-6A and HHV-6B.

Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6, iciHHV6). We reconstructed the population structure of HHV-6 and we show that HHV-6A genomes diverged less than HHV-6B genomes from the ancestral common HHV-6A/B population. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, exogenous viral and iciHHV-6B populations were similar in terms of ancestry components, with no evident geographic structuring. Most HHV-6B genomes sampled to date derive from viral populations that experienced considerable drift. However, a population of HHV-6 exogenous viruses, currently classified as HHV-6B and sampled in New York state, formed a separate cluster (NY cluster) and harbored a considerable portion of HHV-6A-like ancestry. Recombination detection methods identified these viruses as interspecies recombinants, but phylogenetic reconstruction indicated that the recombination signals are due to shared ancestry. In analogy to iciHHV-6A, NY cluster viruses have high nucleotide diversity and constant population size. We propose that HHV-6A sequences and the NY cluster population diverged from an ancestral HHV-6A-like population. A relatively recent bottleneck of the NY (or a related) population with subsequent expansion originated most HHV-6B genomes currently sampled. Our findings indicate that the distinction between HHV-6A and -6B is not as clear-cut as previously thought. More generally, epidemiological and clinical surveys would benefit from taking HHV-6 genetic diversity into account.