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.

OAS (2′–5′-oligoadenylate synthases) proteins and cyclic GMP–AMP synthase (cGAS, gene symbol: MB21D1) patrol the cytoplasm for the presence of foreign nucleic acids. Upon binding to double-stranded RNA or double-stranded DNA, OAS proteins and cGAS produce nucleotide second messengers to activate RNase L and STING (stimulator of interferon genes, gene symbol: TMEM173), respectively; this leads to the initiation of antiviral responses. We analyzed the evolutionary history of the MB21D1–TMEM173 and OAS–RNASEL axes in primates and bats and found evidence of widespread positive selection in both orders. In TMEM173, residue 230, a major determinant of response to natural ligands and to mimetic drugs (e.g., DMXAA), was positively selected in Primates and Chiroptera. In both orders, selection also targeted an α-helix/loop element in RNase L that modulates the enzyme preference for single-stranded RNA versus stem loops. Analysis of positively selected sites in OAS1, OAS2, and MB21D1 revealed parallel evolution, with the corresponding residues being selected in different genes. As this cannot result from gene conversion, these data suggest that selective pressure acting on OAS and MB21D1 genes is related to nucleic acid recognition and to the specific mechanism of enzyme activation, which requires a conformational change. Finally, a population genetics-phylogenetics analysis in humans, chimpanzees, and gorillas detected several positively selected sites in most genes. Data herein shed light into species-specific differences in infection susceptibility and in response to synthetic compounds, with relevance for the design of synthetic compounds as vaccine adjuvants.

Genome-wide association studies (GWASs) have been widely applied to identify genetic factors that affect complex diseases or traits. Presently, the GWAS Catalog includes >2800 human studies. Of these, only a minority have investigated the susceptibility to infectious diseases or the response to therapies for the treatment or prevention of infections. Despite their limited application in the field, GWASs have provided valuable insights by pinpointing associations to both innate and adaptive immune response loci, as well as novel unexpected risk factors for infection susceptibility. Herein, we discuss some issues and caveats of GWASs for infectious diseases, we review the most recent findings ensuing from these studies, and we provide a brief summary of selected GWASs for infections in non-human mammals. We conclude that, although the general trend in the field of complex traits is to shift from GWAS to next-generation sequencing, important knowledge on infectious disease-related traits can be still gained by GWASs, especially for those conditions that have never been investigated using this approach. We suggest that future studies will benefit from the leveraging of information from the host's and pathogen's genomes, as well as from the exploration of models that incorporate heterogeneity across populations and phenotypes. Interactions within HLA genes or among HLA variants and polymorphisms located outside the major histocompatibility complex may also play an important role in shaping the susceptibility and response to invading pathogens.

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.

Plasmodium falciparum, the causative agent of the deadliest form of malaria, is a member of the Laverania subgenus, which includes ape-infecting parasites. P. falciparum is thought to have originated in gorillas, although infection is now restricted to humans. Laverania parasites display remarkable host-specificity, which is partially mediated by the interaction between parasite ligands and host receptors. We analyse the evolution of BSG (basigin) and GYPA (glycophorin A) in primates/hominins, as well as of their Plasmodium-encoded ligands, PfRH5 and PfEBA175. We show that, in primates, positive selection targeted two sites in BSG (F27 and H102), both involved in PfRH5 binding. A population genetics-phylogenetics approach detected the strongest selection for the gorilla lineage: one of the positively selected sites (K191) is a major determinant of PfRH5 binding affinity. Analysis of RH5 genes indicated episodic selection on the P. falciparum branch; the positively selected W447 site is known to stabilize the interaction with human basigin. Conversely, we detect no selection in the receptor-binding region of EBA175 in the P. falciparum lineage. Its host receptor, GYPA, shows evidence of positive selection in all hominid lineages; selected codons include glycosylation sites that modulate PfEBA175 binding affinity. Data herein provide an evolutionary explanation for species-specific binding of the PfRH5-BSG ligand-receptor pair and support the hypothesis that positive selection at these genes drove the host shift leading to the emergence of P. falciparum as a human pathogen.