Computational and Functional Analysis of the Virus-Receptor Interface Reveals Host Range Trade-Offs in New World Arenaviruses.

UNLABELLED: Animal viruses frequently cause zoonotic disease in humans. As these viruses are highly diverse, evaluating the threat that they pose remains a major challenge, and efficient approaches are needed to rapidly predict virus-host compatibility. Here, we develop a combined computational and experimental approach to assess the compatibility of New World arenaviruses, endemic in rodents, with the host TfR1 entry receptors of different potential new host species. Using signatures of positive selection, we identify a small motif on rodent TfR1 that conveys species specificity to the entry of viruses into cells. However, we show that mutations in this region affect the entry of each arenavirus differently. For example, a human single nucleotide polymorphism (SNP) in this region, L212V, makes human TfR1 a weaker receptor for one arenavirus, Machupo virus, but a stronger receptor for two other arenaviruses, Junin and Sabia viruses. Collectively, these findings set the stage for potential evolutionary trade-offs, where natural selection for resistance to one virus may make humans or rodents susceptible to other arenavirus species. Given the complexity of this host-virus interplay, we propose a computational method to predict these interactions, based on homology modeling and computational docking of the virus-receptor protein-protein interaction. We demonstrate the utility of this model for Machupo virus, for which a suitable cocrystal structural template exists. Our model effectively predicts whether the TfR1 receptors of different species will be functional receptors for Machupo virus entry. Approaches such at this could provide a first step toward computationally predicting the "host jumping" potential of a virus into a new host species. IMPORTANCE: We demonstrate how evolutionary trade-offs may exist in the dynamic evolutionary interplay between viruses and their hosts, where natural selection for resistance to one virus could make humans or rodents susceptible to other virus species. We present an algorithm that predicts which species have cell surface receptors that make them susceptible to Machupo virus, based on computational docking of protein structures. Few molecular models exist for predicting the risk of spillover of a particular animal virus into humans or new animal populations. Our results suggest that a combination of evolutionary analysis, structural modeling, and experimental verification may provide an efficient approach for screening and assessing the potential spillover risks of viruses circulating in animal populations.

Widespread sequence variations in VAMP1 across vertebrates suggest a potential selective pressure from botulinum neurotoxins.

Botulinum neurotoxins (BoNT/A-G), the most potent toxins known, act by cleaving three SNARE proteins required for synaptic vesicle exocytosis. Previous studies on BoNTs have generally utilized the major SNARE homologues expressed in brain (VAMP2, syntaxin 1, and SNAP-25). However, BoNTs target peripheral motor neurons and cause death by paralyzing respiratory muscles such as the diaphragm. Here we report that VAMP1, but not VAMP2, is the SNARE homologue predominantly expressed in adult rodent diaphragm motor nerve terminals and in differentiated human motor neurons. In contrast to the highly conserved VAMP2, BoNT-resistant variations in VAMP1 are widespread across vertebrates. In particular, we identified a polymorphism at position 48 of VAMP1 in rats, which renders VAMP1 either resistant (I48) or sensitive (M48) to BoNT/D. Taking advantage of this finding, we showed that rat diaphragms with I48 in VAMP1 are insensitive to BoNT/D compared to rat diaphragms with M48 in VAMP1. This unique intra-species comparison establishes VAMP1 as a physiological toxin target in diaphragm motor nerve terminals, and demonstrates that the resistance of VAMP1 to BoNTs can underlie the insensitivity of a species to members of BoNTs. Consistently, human VAMP1 contains I48, which may explain why humans are insensitive to BoNT/D. Finally, we report that residue 48 of VAMP1 varies frequently between M and I across seventeen closely related primate species, suggesting a potential selective pressure from members of BoNTs for resistance in vertebrates.

Dual host-virus arms races shape an essential housekeeping protein.

Transferrin Receptor (TfR1) is the cell-surface receptor that regulates iron uptake into cells, a process that is fundamental to life. However, TfR1 also facilitates the cellular entry of multiple mammalian viruses. We use evolutionary and functional analyses of TfR1 in the rodent clade, where two families of viruses bind this receptor, to mechanistically dissect how essential housekeeping genes like TFR1 successfully balance the opposing selective pressures exerted by host and virus. We find that while the sequence of rodent TfR1 is generally conserved, a small set of TfR1 residue positions has evolved rapidly over the speciation of rodents. Remarkably, all of these residues correspond to the two virus binding surfaces of TfR1. We show that naturally occurring mutations at these positions block virus entry while simultaneously preserving iron-uptake functionalities, both in rodent and human TfR1. Thus, by constantly replacing the amino acids encoded at just a few residue positions, TFR1 divorces adaptation to ever-changing viruses from preservation of key cellular functions. These dynamics have driven genetic divergence at the TFR1 locus that now enforces species-specific barriers to virus transmission, limiting both the cross-species and zoonotic transmission of these viruses.

Rapid evolution of BRCA1 and BRCA2 in humans and other primates.

BACKGROUND: The maintenance of chromosomal integrity is an essential task of every living organism and cellular repair mechanisms exist to guard against insults to DNA. Given the importance of this process, it is expected that DNA repair proteins would be evolutionarily conserved, exhibiting very minimal sequence change over time. However, BRCA1, an essential gene involved in DNA repair, has been reported to be evolving rapidly despite the fact that many protein-altering mutations within this gene convey a significantly elevated risk for breast and ovarian cancers. RESULTS: To obtain a deeper understanding of the evolutionary trajectory of BRCA1, we analyzed complete BRCA1 gene sequences from 23 primate species. We show that specific amino acid sites have experienced repeated selection for amino acid replacement over primate evolution. This selection has been focused specifically on humans and our closest living relatives, chimpanzees (Pan troglodytes) and bonobos (Pan paniscus). After examining BRCA1 polymorphisms in 7 bonobo, 44 chimpanzee, and 44 rhesus macaque (Macaca mulatta) individuals, we find considerable variation within each of these species and evidence for recent selection in chimpanzee populations. Finally, we also sequenced and analyzed BRCA2 from 24 primate species and find that this gene has also evolved under positive selection. CONCLUSIONS: While mutations leading to truncated forms of BRCA1 are clearly linked to cancer phenotypes in humans, there is also an underlying selective pressure in favor of amino acid-altering substitutions in this gene. A hypothesis where viruses are the drivers of this natural selection is discussed.

Evolutionary reconstructions of the transferrin receptor of Caniforms supports canine parvovirus being a re-emerged and not a novel pathogen in dogs.

Parvoviruses exploit transferrin receptor type-1 (TfR) for cellular entry in carnivores, and specific interactions are key to control of host range. We show that several key mutations acquired by TfR during the evolution of Caniforms (dogs and related species) modified the interactions with parvovirus capsids by reducing the level of binding. These data, along with signatures of positive selection in the TFRC gene, are consistent with an evolutionary arms race between the TfR of the Caniform clade and parvoviruses. As well as the modifications of amino acid sequence which modify binding, we found that a glycosylation site mutation in the TfR of dogs which provided resistance to the carnivore parvoviruses which were in circulation prior to about 1975 predates the speciation of coyotes and dogs. Because the closely-related black-backed jackal has a TfR similar to their common ancestor and lacks the glycosylation site, reconstructing this mutation into the jackal TfR shows the potency of that site in blocking binding and infection and explains the resistance of dogs until recent times. This alters our understanding of this well-known example of viral emergence by indicating that canine parvovirus emergence likely resulted from the re-adaptation of a parvovirus to the resistant receptor of a former host.

Species-specific features of DARC, the primate receptor for Plasmodium vivax and Plasmodium knowlesi.

The DARC (Duffy antigen/receptor for chemokines) gene, also called Duffy or FY, encodes a membrane-bound chemokine receptor. Two malaria parasites, Plasmodium vivax and Plasmodium knowlesi, use DARC to trigger internalization into red blood cells. Although much has been reported on the evolution of DARC null alleles, little is known about the evolution of the coding portion of this gene or the role that protein sequence divergence in this receptor may play in disease susceptibility or zoonosis. Here, we show that the Plasmodium interaction domain of DARC is nearly invariant in the human population, suggesting that coding polymorphism there is unlikely to play a role in differential susceptibility to infection. However, an analysis of DARC orthologs from 35 simian primate species reveals high levels of sequence divergence in the Plasmodium interaction domain. Signatures of positive selection in this domain indicate that species-specific mutations in the protein sequence of DARC could serve as barriers to the transmission of Plasmodium between primate species.

Evidence for ACE2-utilizing coronaviruses (CoVs) related to severe acute respiratory syndrome CoV in bats.

In 2002, severe acute respiratory syndrome (SARS)-coronavirus (CoV) appeared as a novel human virus with high similarity to bat coronaviruses. However, while SARS-CoV uses the human angiotensin-converting enzyme 2 (ACE2) receptor for cellular entry, no coronavirus isolated from bats appears to use ACE2. Here we show that signatures of recurrent positive selection in the bat ACE2 gene map almost perfectly to known SARS-CoV interaction surfaces. Our data indicate that ACE2 utilization preceded the emergence of SARS-CoV-like viruses from bats.

Ancient and recent adaptive evolution of primate non-homologous end joining genes.

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes.

Identification and dissection of a complex DNA repair sensitivity phenotype in Baker's yeast.

Complex traits typically involve the contribution of multiple gene variants. In this study, we took advantage of a high-density genotyping analysis of the BY (S288c) and RM strains of Saccharomyces cerevisiae and of 123 derived spore progeny to identify the genetic loci that underlie a complex DNA repair sensitivity phenotype. This was accomplished by screening hybrid yeast progeny for sensitivity to a variety of DNA damaging agents. Both the BY and RM strains are resistant to the ultraviolet light-mimetic agent 4-nitroquinoline 1-oxide (4-NQO); however, hybrid progeny from a BYxRM cross displayed varying sensitivities to the drug. We mapped a major quantitative trait locus (QTL), RAD5, and identified the exact polymorphism within this locus responsible for 4-NQO sensitivity. By using a backcrossing strategy along with array-assisted bulk segregant analysis, we identified one other locus, MKT1, and a QTL on Chromosome VII that also link to the hybrid 4-NQO-sensitive phenotype but confer more minor effects. This work suggests an additive model for sensitivity to 4-NQO and provides a strategy for mapping both major and minor QTL that confer background-specific phenotypes. It also provides tools for understanding the effect of genetic background on sensitivity to genotoxic agents.

Incompatibilities involving yeast mismatch repair genes: a role for genetic modifiers and implications for disease penetrance and variation in genomic mutation rates.

Genetic background effects underlie the penetrance of most genetically determined phenotypes, including human diseases. To explore how such effects can modify a mutant phenotype in a genetically tractable system, we examined an incompatibility involving the MLH1 and PMS1 mismatch repair genes using a large population sample of geographically and ecologically diverse Saccharomyces cerevisiae strains. The mismatch repair incompatibility segregates into naturally occurring yeast strains, with no strain bearing the deleterious combination. In assays measuring the mutator phenotype conferred by different combinations of MLH1 and PMS1 from these strains, we observed a mutator phenotype only in combinations predicted to be incompatible. Surprisingly, intragenic modifiers could be mapped that specifically altered the strength of the incompatibility over a 20-fold range. Together, these observations provide a powerful model in which to understand the basis of disease penetrance and how such genetic variation, created through mating, could result in new mutations that could be the raw material of adaptive evolution in yeast populations.

A putative SUMO interacting motif in the B30.2/SPRY domain of rhesus macaque TRIM5α important for NF-κB/AP-1 signaling and HIV-1 restriction.

TRIM5α from the rhesus macaque (TRIM5αRh) is a restriction factor that shows strong activity against HIV-1. TRIM5αRh binds specifically to HIV-1 capsid (CA) through its B30.2/PRYSPRY domain shortly after entry of the virus into the cytoplasm. Recently, three putative SUMO interacting motifs (SIMs) have been identified in the PRYSPRY domain of human and macaque TRIM5α. However, structural modeling of this domain suggested that two of them were buried in the hydrophobic core of the protein, implying that interaction with SUMO was implausible, while the third one was not relevant to restriction. In light of these results, we re-analyzed the TRIM5αRh PRYSPRY sequence and identified an additional putative SIM ((435)VIIC(438)) which we named SIM4. This motif is exposed at the surface of the PRYSPRY domain, allowing potential interactions with SUMO or SUMOylated proteins. Introducing a double mutation in SIM4 (V435K, I436K) did not alter stability, unlike mutations in SIM1. SIM4-mutated TRIM5αRh failed to bind HIV-1CA and lost the ability to restrict this virus. Accordingly, SIM4 undergoes significant variation among primates and substituting this motif with naturally occurring SIM4 variants affected HIV-1 restriction by TRIM5αRh, suggesting a direct role in capsid recognition. Interestingly, SIM4-mutated TRIM5αRh also failed to activate NF-κB and AP-1-mediated transcription. Although there is no direct evidence that SIM4 is involved in direct interaction with SUMO or a SUMOylated protein, mutating this motif strongly reduced co-localization of TRIM5αRh with SUMO-1 and with PML, a SUMOylated nuclear protein. In conclusion, this new putative SIM is crucial for both direct interaction with incoming capsids and for NF-κB/AP-1 signaling. We speculate that the latter function is mediated by interactions of SIM4 with a SUMOylated protein involved in the NF-κB/AP-1 signaling pathways.

Filovirus receptor NPC1 contributes to species-specific patterns of ebolavirus susceptibility in bats.

Biological factors that influence the host range and spillover of Ebola virus (EBOV) and other filoviruses remain enigmatic. While filoviruses infect diverse mammalian cell lines, we report that cells from African straw-colored fruit bats (Eidolon helvum) are refractory to EBOV infection. This could be explained by a single amino acid change in the filovirus receptor, NPC1, which greatly reduces the affinity of EBOV-NPC1 interaction. We found signatures of positive selection in bat NPC1 concentrated at the virus-receptor interface, with the strongest signal at the same residue that controls EBOV infection in Eidolon helvum cells. Our work identifies NPC1 as a genetic determinant of filovirus susceptibility in bats, and suggests that some NPC1 variations reflect host adaptations to reduce filovirus replication and virulence. A single viral mutation afforded escape from receptor control, revealing a pathway for compensatory viral evolution and a potential avenue for expansion of filovirus host range in nature.

Detection of heterozygous mutations in the genome of mismatch repair defective diploid yeast using a Bayesian approach.

DNA replication errors that escape polymerase proofreading and mismatch repair (MMR) can lead to base substitution and frameshift mutations. Such mutations can disrupt gene function, reduce fitness, and promote diseases such as cancer and are also the raw material of molecular evolution. To analyze with limited bias genomic features associated with DNA polymerase errors, we performed a genome-wide analysis of mutations that accumulate in MMR-deficient diploid lines of Saccharomyces cerevisiae. These lines were derived from a common ancestor and were grown for 160 generations, with bottlenecks reducing the population to one cell every 20 generations. We sequenced to between 8- and 20-fold coverage one wild-type and three mutator lines using Illumina Solexa 36-bp reads. Using an experimentally aware Bayesian genotype caller developed to pool experimental data across sequencing runs for all strains, we detected 28 heterozygous single-nucleotide polymorphisms (SNPs) and 48 single-nt insertion/deletions (indels) from the data set. This method was evaluated on simulated data sets and found to have a very low false-positive rate (∼6 × 10(-5)) and a false-negative rate of 0.08 within the unique mapping regions of the genome that contained at least sevenfold coverage. The heterozygous mutations identified by the Bayesian genotype caller were confirmed by Sanger sequencing. All of the mutations were unique to a given line, except for a single-nt deletion mutation which occurred independently in two lines. All 48 indels, composed of 46 deletions and two insertions, occurred in homopolymer (HP) tracts [i.e., 47 poly(A) or (T) tracts, 1 poly(G) or (C) tract] between 5 and 13 bp long. Our findings are of interest because HP tracts are present at high levels in the yeast genome (>77,400 for 5- to 20-nt HP tracts), and frameshift mutations in these regions are likely to disrupt gene function. In addition, they demonstrate that the mutation pattern seen previously in mismatch repair defective strains using a limited number of reporters holds true for the entire genome.