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.

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.

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.

A 175 million year history of T cell regulatory molecules reveals widespread selection, with adaptive evolution of disease alleles.

T cell activation plays a central role in immune response and in the maintenance of self-tolerance. We analyzed the evolutionary history of T cell regulatory molecules. Nine genes involved in triggering T cell activation or in regulating the ensuing response evolved adaptively in mammals. Several positively selected sites overlap with positions interacting with the binding partner or with cellular components. Population genetic analysis in humans revealed a complex scenario of local (FASLG, CD40LG, HAVCR2) and worldwide (FAS, ICOSLG) adaptation and H. sapiens-to-Neandertal gene flow (gene transfer between populations). Disease variants in these genes are preferential targets of pathogen-driven selection, and a Crohn's disease risk polymorphism targeted by bacterial-driven selection modulates the expression of ICOSLG in response to a bacterial superantigen. Therefore, we used evolutionary information to generate experimentally testable hypotheses concerning the function of specific genetic variants and indicate that adaptation to infection underlies the maintenance of autoimmune risk alleles.

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.

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.

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.

Key role of SMN/SYNCRIP and RNA-Motif 7 in spinal muscular atrophy: RNA-Seq and motif analysis of human motor neurons.

Spinal muscular atrophy is a motor neuron disorder caused by mutations in SMN1. The reasons for the selective vulnerability of motor neurons linked to SMN (encoded by SMN1) reduction remain unclear. Therefore, we performed deep RNA sequencing on human spinal muscular atrophy motor neurons to detect specific altered gene splicing/expression and to identify the presence of a common sequence motif in these genes. Many deregulated genes, such as the neurexin and synaptotagmin families, are implicated in critical motor neuron functions. Motif-enrichment analyses of differentially expressed/spliced genes, including neurexin2 (NRXN2), revealed a common motif, motif 7, which is a target of SYNCRIP. Interestingly, SYNCRIP interacts only with full-length SMN, binding and modulating several motor neuron transcripts, including SMN itself. SYNCRIP overexpression rescued spinal muscular atrophy motor neurons, due to the subsequent increase in SMN and their downstream target NRXN2 through a positive loop mechanism and ameliorated SMN-loss-related pathological phenotypes in Caenorhabditis elegans and mouse models. SMN/SYNCRIP complex through motif 7 may account for selective motor neuron degeneration and represent a potential therapeutic target.

Evolutionary rates of mammalian telomere-stability genes correlate with karyotype features and female germline expression.

Telomeres protect the ends of eukaryotic chromosomes and are essential for cell viability. In mammals, telomere dynamics vary with life history traits (e.g. body mass and longevity), suggesting differential selection depending on physiological characteristics. Telomeres, in analogy to centromeric regions, also represent candidate meiotic drivers and subtelomeric DNA evolves rapidly. We analyzed the evolutionary history of mammalian genes implicated in telomere homeostasis (TEL genes). We detected widespread positive selection and we tested two alternative hypotheses: (i) fast evolution is driven by changes in life history traits; (ii) a conflict with selfish DNA elements at the female meiosis represents the underlying selective pressure. By accounting for the phylogenetic relationships among mammalian species, we show that life history traits do not contribute to shape diversity of TEL genes. Conversely, the evolutionary rate of TEL genes correlates with expression levels during meiosis and episodes of positive selection across mammalian species are associated with karyotype features (number of chromosome arms). We thus propose a telomere drive hypothesis, whereby (sub)telomeres and telomere-binding proteins are engaged in an intra-genomic conflict similar to the one described for centromeres.

Positive Selection Drives Evolution at the Host-Filovirus Interaction Surface.

Filovirus infection is mediated by engagement of the surface-exposed glycoprotein (GP) by its cellular receptor, NPC1 (Niemann-Pick C1). Two loops in the C domain of NPC1 (NPC1-C) bind filovirus GP. Herein, we show that filovirus GP and NPC1-C evolve under mutual selective pressure. Analysis of a large mammalian phylogeny indicated that strong functional/structural constraints limit the NPC1 sequence space available for adaptive change and most sites at the contact interface with GP are under negative selection. These constraints notwithstanding, we detected positive selection at NPC1-C in all mammalian orders, from Primates to Xenarthra. Different codons evolved adaptively in distinct mammals, and most selected sites are located within the two NPC1-C loops that engage GP, or at their anchor points. In Homininae, NPC1-C was a preferential selection target, and the T419I variant possibly represents a human-specific adaptation to filovirus infection. On the other side of the arms-race, GP evolved adaptively during filovirus speciation. One of the selected sites (S142Q) establishes several atom-to-atom contacts with NPC1-C. Additional selected sites are located within epitopes recognized by neutralizing antibodies, including the 14G7 epitope, where sites selected during the recent EBOV epidemic also map. Finally, pairs of co-evolving sites in Marburgviruses and Ebolaviruses were found to involve antigenic determinants. These findings suggest that the host humoral immune response was a major selective pressure during filovirus speciation. The S142Q variant may contribute to determine Ebolavirus host range in the wild. If this were the case, EBOV/BDBV (S142) and SUDV (Q142) may not share the same reservoir(s).

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.

Evolutionary analysis of Old World arenaviruses reveals a major adaptive contribution of the viral polymerase.

The Old World (OW) arenavirus complex includes several species of rodent-borne viruses, some of which (i.e., Lassa virus, LASV and Lymphocytic choriomeningitis virus, LCMV) cause human diseases. Most LCMV and LASV infections are caused by rodent-to-human transmissions. Thus, viral evolution is largely determined by events that occur in the wildlife reservoirs. We used a set of human- and rodent-derived viral sequences to investigate the evolutionary history underlying OW arenavirus speciation, as well as the more recent selective events that accompanied LASV spread in West Africa. We show that the viral RNA polymerase (L protein) was a major positive selection target in OW arenaviruses and during LASV out-of-Nigeria migration. No evidence of selection was observed for the glycoprotein, whereas positive selection acted on the nucleoprotein (NP) during LCMV speciation. Positively selected sites in L and NP are surrounded by highly conserved residues, and the bulk of the viral genome evolves under purifying selection. Several positively selected sites are likely to modulate viral replication/transcription. In both L and NP, structural features (solvent exposed surface area) are important determinants of site-wise evolutionary rate variation. By incorporating several rodent-derived sequences, we also performed an analysis of OW arenavirus codon adaptation to the human host. Results do not support a previously hypothesized role of codon adaptation in disease severity for non-Nigerian strains. In conclusion, L and NP represent the major selection targets and possible determinants of disease presentation; these results suggest that field surveys and experimental studies should primarily focus on these proteins.

An evolutionary analysis of antigen processing and presentation across different timescales reveals pervasive selection.

The antigenic repertoire presented by MHC molecules is generated by the antigen processing and presentation (APP) pathway. We analyzed the evolutionary history of 45 genes involved in APP at the inter- and intra-species level. Results showed that 11 genes evolved adaptively in mammals. Several positively selected sites involve positions of fundamental importance to the protein function (e.g. the TAP1 peptide-binding domains, the sugar binding interface of langerin, and the CD1D trafficking signal region). In CYBB, all selected sites cluster in two loops protruding into the endosomal lumen; analysis of missense mutations responsible for chronic granulomatous disease (CGD) showed the action of different selective forces on the very same gene region, as most CGD substitutions involve aminoacid positions that are conserved in all mammals. As for ERAP2, different computational methods indicated that positive selection has driven the recurrent appearance of protein-destabilizing variants during mammalian evolution. Application of a population-genetics phylogenetics approach showed that purifying selection represented a major force acting on some APP components (e.g. immunoproteasome subunits and chaperones) and allowed identification of positive selection events in the human lineage. We also investigated the evolutionary history of APP genes in human populations by developing a new approach that uses several different tests to identify the selection target, and that integrates low-coverage whole-genome sequencing data with Sanger sequencing. This analysis revealed that 9 APP genes underwent local adaptation in human populations. Most positive selection targets are located within noncoding regions with regulatory function in myeloid cells or act as expression quantitative trait loci. Conversely, balancing selection targeted nonsynonymous variants in TAP1 and CD207 (langerin). Finally, we suggest that selected variants in PSMB10 and CD207 contribute to human phenotypes. Thus, we used evolutionary information to generate experimentally-testable hypotheses and to provide a list of sites to prioritize in follow-up analyses.

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.

Serotonin Transporter Gene (SLC6A4) Methylation Associates With Neonatal Intensive Care Unit Stay and 3-Month-Old Temperament in Preterm Infants.

Preterm birth and Neonatal Intensive Care Unit (NICU) stay are early adverse stressful experiences, which may result in an altered temperamental profile. The serotonin transporter gene (SLC6A4), which has been linked to infant temperament, is susceptible to epigenetic regulation associated with early stressful experience. This study examined a moderation model in which the exposure to NICU-related stress and SLC6A4 methylation moderated infant temperament at 3 months of age. SLC6A4 methylation at 20 CpG sites was quantified in preterm infants (N = 48) and full-term infants (N = 30) from Italian middle-class families. Results suggested that in preterm infants NICU-related stress might be associated with alterations of serotonergic tone as a consequence of SLC6A4 methylation, which in turn, might associate with temperamental difficulties assessed at 3 months of age.

Evolutionary analysis identifies an MX2 haplotype associated with natural resistance to HIV-1 infection.

The protein product of the myxovirus resistance 2 (MX2) gene restricts HIV-1 and simian retroviruses. We demonstrate that MX2 evolved adaptively in mammals with distinct sites representing selection targets in distinct branches; selection mainly involved residues in loop 4, previously shown to carry antiviral determinants. Modeling data indicated that positively selected sites form a continuous surface on loop 4, which folds into two antiparallel α-helices protruding from the stalk domain. A population genetics-phylogenetics approach indicated that the coding region of MX2 mainly evolved under negative selection in the human lineage. Nonetheless, population genetic analyses demonstrated that natural selection operated on MX2 during the recent history of human populations: distinct selective events drove the frequency increase of two haplotypes in the populations of Asian and European ancestry. The Asian haplotype carries a susceptibility allele for melanoma; the European haplotype is tagged by rs2074560, an intronic variant. Analyses performed on three independent European cohorts of HIV-1-exposed seronegative individuals with different geographic origin and distinct exposure route showed that the ancestral (G) allele of rs2074560 protects from HIV-1 infection with a recessive effect (combined P = 1.55 × 10(-4)). The same allele is associated with lower in vitro HIV-1 replication and increases MX2 expression levels in response to IFN-α. Data herein exploit evolutionary information to identify a novel host determinant of HIV-1 infection susceptibility.

Positive selection underlies the species-specific binding of Plasmodium falciparum RH5 to human basigin.

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.

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.