Search

ABSTRACTThe COVID-19 pandemic demonstrated the poor ability of body temperature to reliably identify SARS-CoV-2-infected individuals, an observation that has been made before in the context of other infectious diseases. While acute infection does not always cause fever, it does reliably drive host transcriptional responses as the body responds at the site of infection. These transcriptional changes can occur both in cells that are directly harboring replicating pathogens and in cells elsewhere that receive a molecular signal that infection is occurring. Here, we identify a core set of approximately 70 human genes that are together upregulated in cultured human cells infected by a broad array of viral, bacterial, and fungal pathogens. We have named these “core response” genes. In theory, transcripts from these genes could serve as biomarkers of infection in the human body, in a way that is agnostic to the specific pathogen causing infection. As such, we perform human studies to show that these infection-induced human transcripts can be measured in the saliva of people harboring different types of infections. The number of these transcripts in saliva can correctly classify infection status (whether a person harbors an infection) 91% of the time. Furthermore, in the case of SARS-CoV-2 specifically, the number of core response transcripts in saliva correctly identifies infectious individuals even when enrollees, themselves, are asymptomatic and do not know they are infected.IMPORTANCEThere are a variety of clinical and laboratory criteria available to clinicians in controlled healthcare settings to help them identify whether an infectious disease is present. However, in situations such as a new epidemic caused by an unknown infectious agent, in health screening contexts performed within communities and outside of healthcare facilities or in battlefield or potential biowarfare situations, this gets more difficult. Pathogen-agnostic methods for rapid screening and triage of large numbers of people for infection status are needed, in particular methods that might work on an easily accessible biospecimen like saliva. Here, we identify a small, core set of approximately 70 human genes whose transcripts serve as saliva-based biomarkers of infection in the human body, in a way that is agnostic to the specific pathogen causing infection.

Abstract Background Over the course of the COVID-19 pandemic, colleges and universities have focused on creating policies, such as mask mandates, to minimize COVID-19 transmission both on their campuses and in the surrounding community. Adherence to and opinions about these policies remain largely unknown. Methods The Centers for Disease Control and Prevention (CDC) developed a cross-sectional study, the Mask Adherence and Surveillance at Colleges and Universities Project (MASCUP!), to objectively and inconspicuously measure rates of mask use at institutes of higher education via direct observation. From February 15 through April 11, 2021 the University of Colorado Boulder (CU, n = 2,808 observations) and Colorado State University Fort Collins (CSU, n = 3,225 observations) participated in MASCUP! along with 52 other institutes of higher education (n = 100,353 observations) spanning 21 states and the District of Columbia. Mask use was mandatory at both Colorado universities and student surveys were administered to assess student beliefs and attitudes. Results We found that 91.7%, 93.4%, and 90.8% of persons observed at indoor locations on campus wore a mask correctly at University of Colorado, Colorado State University, and across the 52 other schools, respectively. Student responses to questions about masking were in line with these observed rates of mask use where 92.9% of respondents at CU and 89.8% at CSU believe that wearing masks can protect the health of others. Both Colorado universities saw their largest surges in COVID-19 cases in the fall of 2020, with markedly lower case counts during the mask observation window in the spring of 2021. Conclusion High levels of mask use at Colorado’s two largest campuses aligned with rates observed at other institutes across the country. These high rates of use, coupled with positive student attitudes about mask use, demonstrate that masks were widely accepted and may have contributed to reduced COVID-19 case counts. This study supports an emerging body of literature substantiating masks as an effective, low-cost measure to reduce disease transmission and establishes masking (with proper education and promotion) as a viable tactic to reduce respiratory disease transmission on college campuses.

Summary: We present a protocol to detect cells that have been infected by RNA viruses. The method, RNA fluorescence in situ hybridization flow cytometry (RNA FISH-Flow), uses 48 fluorescently labeled DNA probes that hybridize in tandem to viral RNA. RNA FISH-Flow probes can be synthesized to match any RNA virus genome, in either sense or anti-sense, enabling detection of genomes or replication intermediates within cells. Flow cytometry enables high-throughput analysis of infection dynamics within a population at the single cell level.For complete details on the use and execution of this protocol, please refer to Warren et al. (2022).1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Abstract Background The Type I interferon response is an important first-line defense against viruses. In turn, viruses antagonize (i.e., degrade, mis-localize, etc.) many proteins in interferon pathways. Thus, hosts and viruses are locked in an evolutionary arms race for dominance of the Type I interferon pathway. As a result, many genes in interferon pathways have experienced positive natural selection in favor of new allelic forms that can better recognize viruses or escape viral antagonists. Here, we performed a holistic analysis of selective pressures acting on genes in the Type I interferon family. We initially hypothesized that the genes responsible for inducing the production of interferon would be antagonized more heavily by viruses than genes that are turned on as a result of interferon. Our logic was that viruses would have greater effect if they worked upstream of the production of interferon molecules because, once interferon is produced, hundreds of interferon-stimulated proteins would activate and the virus would need to counteract them one-by-one. Results We curated multiple sequence alignments of primate orthologs for 131 genes active in interferon production and signaling (herein, “induction” genes), 100 interferon-stimulated genes, and 100 randomly chosen genes. We analyzed each multiple sequence alignment for the signatures of recurrent positive selection. Counter to our hypothesis, we found the interferon-stimulated genes, and not interferon induction genes, are evolving significantly more rapidly than a random set of genes. Interferon induction genes evolve in a way that is indistinguishable from a matched set of random genes (22% and 18% of genes bear signatures of positive selection, respectively). In contrast, interferon-stimulated genes evolve differently, with 33% of genes evolving under positive selection and containing a significantly higher fraction of codons that have experienced selection for recurrent replacement of the encoded amino acid. Conclusion Viruses may antagonize individual products of the interferon response more often than trying to neutralize the system altogether.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in over 6 million deaths worldwide. The high variability in COVID-19 symptoms remains one of the most interesting mysteries of the pandemic. Genetic and environmental factors are likely to be key determinants of COVID-19 symptomatology. Here, we explored ACE2 as a genetic determinant for SARS-CoV-2 infection and COVID-19 symptomatology. Each human genome encodes two alleles of ACE2, which encodes the cell entry receptor for SARS-CoV-2. Here, we determined whether naturally occurring human ACE2 (hACE2) polymorphisms in the human population affect SARS-CoV-2 infection and the severity of COVID-19 symptoms. ACE2 variants S19P, I21V, E23K, K26R, K31R, N33I, H34R, E35K, and T92I showed increased virus infection compared to wild-type ACE2; thus, these variants could increase the risk for COVID-19. In contrast, variants D38V, Y83H, I468V, and N638S showed reduced infection, indicating a potential protective effect. hACE2 variants K26R and T92I increased infection by three-fold without changing the levels of ACE2 on the surface of the cells, suggesting that these variants may increase the risk of severe COVID-19. On the contrary, hACE2 variants D38V and Y83H decreased SARS-CoV-2 infection by four- and ten-fold, respectively, without changing surface expression, suggesting that these variants may protect against severe COVID-19. Remarkably, all protective hACE2 Polymorphisms were found almost exclusively in Asian populations, which may provide a partial explanation for the low COVID-19 mortality rates in Asian countries. Thus, hACE2 polymorphisms may modulate susceptibility to SARS-CoV-2 in the host and partially account for the differences in severity of COVID-19 among different ethnic groups.

Summary: Tripartite motif-containing protein 5α (TRIM5α) is a cellular antiviral restriction factor that prevents early events in retrovirus replication. The activity of TRIM5α is thought to be limited to retroviruses as a result of highly specific interactions with capsid lattices. In contrast to this current understanding, we show that both human and rhesus macaque TRIM5α suppress replication of specific flaviviruses. Multiple viruses in the tick-borne encephalitis complex are sensitive to TRIM5α-dependent restriction, but mosquito-borne flaviviruses, including yellow fever, dengue, and Zika viruses, are resistant. TRIM5α suppresses replication by binding to the viral protease NS2B/3 to promote its K48-linked ubiquitination and proteasomal degradation. Importantly, TRIM5α contributes to the antiviral function of IFN-I against sensitive flaviviruses in human cells. Thus, TRIM5α possesses remarkable plasticity in the recognition of diverse virus families, with the potential to influence human susceptibility to emerging flaviviruses of global concern. : The antiviral activity of TRIM5α is thought to be limited to retroviruses as a result of highly specific interactions with capsid lattices. Here, Chiramel et al. demonstrate that TRIM5α restricts replication of specific flaviviruses by binding and degrading the viral protease. Keywords: TRIM5α, flavivirus, retrovirus, interferon, restriction factor, retrovirus, interferon stimulated genes, tick-borne encephalitis virus

RNA viruses are a major source of emerging and re-emerging infectious diseases around the world. We developed a method to identify RNA viruses that is based on the fact that RNA viruses produce double-stranded RNA (dsRNA) while replicating. Purifying and sequencing dsRNA from the total RNA isolated from infected tissue allowed us to recover dsRNA virus sequences and replicated sequences from single-stranded RNA (ssRNA) viruses. We refer to this approach as dsRNA-Seq. By assembling dsRNA sequences into contigs we identified full length or partial RNA viral genomes of varying genome types infecting mammalian culture samples, identified a known viral disease agent in laboratory infected mice, and successfully detected naturally occurring RNA viral infections in reptiles. Here, we show that dsRNA-Seq is a preferable method for identifying viruses in organisms that don’t have sequenced genomes and/or commercially available rRNA depletion reagents. In addition, a significant advantage of this method is the ability to identify replicated viral sequences of ssRNA viruses, which is useful for distinguishing infectious viral agents from potential noninfectious viral particles or contaminants.

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 (435VIIC438) 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.

In many biological applications, we would like to be able to computationally predict mutational effects on affinity in protein–protein interactions. However, many commonly used methods to predict these effects perform poorly in important test cases. In particular, the effects of multiple mutations, non alanine substitutions, and flexible loops are difficult to predict with available tools and protocols. We present here an existing method applied in a novel way to a new test case; we interrogate affinity differences resulting from mutations in a host–virus protein–protein interface. We use steered molecular dynamics (SMD) to computationally pull the machupo virus (MACV) spike glycoprotein (GP1) away from the human transferrin receptor (hTfR1). We then approximate affinity using the maximum applied force of separation and the area under the force-versus-distance curve. We find, even without the rigor and planning required for free energy calculations, that these quantities can provide novel biophysical insight into the GP1/hTfR1 interaction. First, with no prior knowledge of the system we can differentiate among wild type and mutant complexes. Moreover, we show that this simple SMD scheme correlates well with relative free energy differences computed via free energy perturbation. Second, although the static co-crystal structure shows two large hydrogen-bonding networks in the GP1/hTfR1 interface, our simulations indicate that one of them may not be important for tight binding. Third, one viral site known to be critical for infection may mark an important evolutionary suppressor site for infection-resistant hTfR1 mutants. Finally, our approach provides a framework to compare the effects of multiple mutations, individually and jointly, on protein–protein interactions.

Here, we develop a simple molecular test for SARS-CoV-2 in saliva based on reverse transcription loop-mediated isothermal amplification. The test has two steps: (1) heat saliva with a stabilization solution and (2) detect virus by incubating with a primer/enzyme mix. After incubation, saliva samples containing the SARS-CoV-2 genome turn bright yellow. Because this test is pH dependent, it can react falsely to some naturally acidic saliva samples. We report unique saliva stabilization protocols that rendered 295 healthy saliva samples compatible with the test, producing zero false positives. We also evaluated the test on 278 saliva samples from individuals who were infected with SARS-CoV-2 but had no symptoms at the time of saliva collection, and from 54 matched pairs of saliva and anterior nasal samples from infected individuals. The Saliva TwoStep test described herein identified infections with 94% sensitivity and >99% specificity in individuals with sub-clinical (asymptomatic or pre-symptomatic) infections.

Human dengue viruses emerged from primate reservoirs, yet paradoxically dengue does not reach high titers in primate models. This presents a unique opportunity to examine the genetics of spillover versus reservoir hosts. The dengue virus 2 (DENV2) - encoded protease cleaves human STING, reducing type I interferon production and boosting viral titers in humans. We find that both human and sylvatic (reservoir) dengue viruses universally cleave human STING, but not the STING of primates implicated as reservoir species. The special ability of dengue to cleave STING is thus specific to humans and a few closely related ape species. Conversion of residues 78/79 to the human-encoded ‘RG’ renders all primate (and mouse) STINGs sensitive to viral cleavage. Dengue viruses may have evolved to increase viral titers in the dense and vast human population, while maintaining decreased titers and pathogenicity in the more rare animals that serve as their sustaining reservoir in nature.

Viruses are responsible for many of the adaptive mutations in the human genome.

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.

The H3N2 canine influenza virus – which originally came from birds – is evolving to become more transmissible between dogs.

The antiviral endoribonuclease, RNase L, is activated by the mammalian innate immune response to destroy host and viral RNA to ultimately reduce viral gene expression. Herein, we show that RNase L and RNase L-mediated mRNA decay are primarily localized to the cytoplasm. Consequently, RNA-binding proteins (RBPs) translocate from the cytoplasm to the nucleus upon RNase L activation due to the presence of intact nuclear RNA. The re-localization of RBPs to the nucleus coincides with global alterations to RNA processing in the nucleus. While affecting many host mRNAs, these alterations are pronounced in mRNAs encoding type I and type III interferons and correlate with their retention in the nucleus and reduction in interferon protein production. Similar RNA processing defects also occur during infection with either dengue virus or SARS-CoV-2 when RNase L is activated. These findings reveal that the distribution of RBPs between the nucleus and cytosol is dictated by the availability of RNA in each compartment. Thus, viral infections that trigger RNase L-mediated cytoplasmic RNA in the cytoplasm also alter RNA processing in the nucleus, resulting in an ingenious multi-step immune block to protein biogenesis.

AbstractThe gene encoding the ubiquitous DNA repair protein, Ku70p, has undergone extensive copy number expansion during primate evolution. Gene duplications of KU70KU70NUKUNUKUKU70NUKUNUKU2NUKUs

Individuals chronically infected with HIV-1 harbor complex viral populations within their bloodstreams. Recently, it has come to light that when these people infect others, the new infection is typically established by only one or a small number of virions from within this complex viral swarm. An important goal is to characterize the biological properties of HIV-1 virions that seed and exist early in new human infections because these are potentially the only viruses against which a prophylactic HIV-1 vaccine would need to elicit protection. This includes understanding how the Envelope (Env) protein of these virions interacts with the T-cell receptor CD4, which supports attachment and entry of HIV-1 into target cells. We examined early HIV-1 isolates for their ability to infect cells via the CD4 receptor of 15 different primate species. Primates were the original source of HIV-1 and now serve as valuable animal models for studying HIV-1. We find that most primary isolates of HIV-1 from the blood, including early isolates, are highly selective and enter cells through some primate CD4 receptor orthologs but not others. This phenotype is remarkably consistent, regardless of route of transmission, viral subtype, or time of isolation post infection. We show that the weak CD4 binding affinity of blood-derived HIV-1 isolates is what makes them sensitive to the small sequence differences in CD4 from one primate species to the next. To substantiate this, we engineered an early HIV-1 Env to have high, medium, or low binding affinity to CD4, and we show that it loses the ability to enter cells via the CD4 receptor of many primate species as the binding affinity gets weaker. Based on the phenotype of selective use of primate CD4, we find that weak CD4 binding appears to be a nearly universal property of HIV-1 circulating in the bloodstream. Therefore, weak binding to CD4 must be a selected and important property in the biology of HIV-1 in the body. We identify six primate species that encode CD4 receptors that fully support the entry of early HIV-1 isolates despite their low binding affinity for CD4. These findings will help inform long-standing efforts to model HIV-1 transmission and early disease in primates.

Viruses of wild and domestic animals can infect humans in a process called zoonosis, and these events can give rise to explosive epidemics such as those caused by the HIV and Ebola viruses. While humans are constantly exposed to animal viruses, those that can successfully infect and transmit between humans are exceedingly rare. The key event in zoonosis is when an animal virus begins to replicate (one virion making many) in the first human subject. Only at this point will the animal virus first experience the selective environment of the human body, rendering possible viral adaptation and refinement for humans. In addition, appreciable viral titers in this first human may enable infection of a second, thus initiating selection for viral variants with increased capacity for spread. We assert that host genetics plays a critical role in defining which animal viruses in nature will achieve this key event of replication in a first human host. This is because animal viruses that pose the greatest risk to humans will have few (or no) genetic barriers to replicating themselves in human cells, thus requiring minimal mutations to make this jump. Only experimental virology provides a path to identifying animal viruses with the potential to replicate themselves in humans because this information will not be evident from viral sequencing data alone.

Yeasts serve as hosts to several types of genetic parasites. Few studies have addressed the evolutionary trajectory of yeast genes that control the stable co-existence of these parasites with their host cell. In Saccharomyces yeasts, the retrovirus-like Ty retrotransposons must access the nucleus. We show that several genes encoding components of the yeast nuclear pore complex have experienced natural selection for substitutions that change the encoded protein sequence. By replacing these S. cerevisiae genes with orthologs from other Saccharomyces species, we discovered that natural sequence changes have affected the mobility of Ty retrotransposons. Specifically, changing the genetic sequence of NUP84 or NUP82 to match that of other Saccharomyces species alters the mobility of S. cerevisiae Ty1 and Ty3. Importantly, all tested housekeeping functions of NUP84 and NUP82 remained equivalent across species. Signatures of natural selection, resulting in altered interactions with viruses and parasitic genetic elements, are common in host defense proteins. Yet, few instances have been documented in essential housekeeping proteins. The nuclear pore complex is the gatekeeper of the nucleus. This study shows how the evolution of this large, ubiquitous eukaryotic complex can alter the replication of a molecular parasite, but concurrently maintain essential host functionalities regarding nucleocytoplasmic trafficking.

[This corrects the article DOI: 10.1371/journal.ppat.1006906.].

HIV-1 arose as the result of spillover of simian immunodeficiency viruses (SIVs) from great apes in Africa, namely from chimpanzees and gorillas. Chimpanzees and gorillas were, themselves, infected with SIV after virus spillover from African monkeys. During spillover events, SIV is thought to require adaptation to the new host species. The host barriers that drive viral adaptation have predominantly been attributed to restriction factors, rather than cofactors (host proteins exploited to promote viral replication). Here, we consider the role of one cofactor, RanBP2, in providing a barrier that drove viral genome evolution during SIV spillover events. RanBP2 (also known as Nup358) is a component of the nuclear pore complex known to facilitate nuclear entry of HIV-1. Our data suggest that transmission of SIV from monkeys to chimpanzees, and then from chimpanzees to gorillas, both coincided with changes in the viral capsid that allowed interaction with RanBP2 of the new host species. However, human RanBP2 subsequently provided no barrier to the zoonotic transmission of SIV from chimpanzees or gorillas, indicating that chimpanzee- and gorilla-adapted SIVs are pre-adapted to humans in this regard. Our observations are in agreement with RanBP2 driving virus evolution during cross-species transmissions of SIV, particularly in the transmissions to and between great ape species.

Schlafen11 (encoded by the SLFN11 gene) has been shown to inhibit the accumulation of HIV-1 proteins. We show that the SLFN11 gene is under positive selection in simian primates and is species-specific in its activity against HIV-1. The activity of human Schlafen11 is relatively weak compared to that of some other primate versions of this protein, with the versions encoded by chimpanzee, orangutan, gibbon, and marmoset being particularly potent inhibitors of HIV-1 protein production. Interestingly, we find that Schlafen11 is functional in the absence of infection and reduces protein production from certain non-viral (GFP) and even host (Vinculin and GAPDH) transcripts. This suggests that Schlafen11 may just generally block protein production from non-codon optimized transcripts. Because Schlafen11 is an interferon-stimulated gene with a broad ability to inhibit protein production from many host and viral transcripts, its role may be to create a general antiviral state in the cell. Interestingly, the strong inhibitors such as marmoset Schlafen11 consistently block protein production better than weak primate Schlafen11 proteins, regardless of the virus or host target being analyzed. Further, we show that the residues to which species-specific differences in Schlafen11 potency map are distinct from residues that have been targeted by positive selection. We speculate that the positive selection of SLFN11 could have been driven by a number of different factors, including interaction with one or more viral antagonists that have yet to be identified.

In eukaryotes, the degradation of cellular mRNAs is accomplished by Xrn1 and the cytoplasmic exosome. Because viral RNAs often lack canonical caps or poly-A tails, they can also be vulnerable to degradation by these host exonucleases. Yeast lack sophisticated mechanisms of innate and adaptive immunity, but do use RNA degradation as an antiviral defense mechanism. One model is that the RNA of yeast viruses is subject to degradation simply as a side effect of the intrinsic exonuclease activity of proteins involved in RNA metabolism. Contrary to this model, we find a highly refined, species-specific relationship between Xrn1p and the "L-A" totiviruses of different Saccharomyces yeast species. We show that the gene XRN1 has evolved rapidly under positive natural selection in Saccharomyces yeast, resulting in high levels of Xrn1p protein sequence divergence from one yeast species to the next. We also show that these sequence differences translate to differential interactions with the L-A virus, where Xrn1p from S. cerevisiae is most efficient at controlling the L-A virus that chronically infects S. cerevisiae, and Xrn1p from S. kudriavzevii is most efficient at controlling the L-A-like virus that we have discovered within S. kudriavzevii. All Xrn1p orthologs are equivalent in their interaction with another virus-like parasite, the Ty1 retrotransposon. Thus, the activity of Xrn1p against totiviruses is not simply an incidental consequence of the enzymatic activity of Xrn1p, but rather Xrn1p co-evolves with totiviruses to maintain its potent antiviral activity and limit viral propagation in Saccharomyces yeasts. Consistent with this, we demonstrated that Xrn1p physically interacts with the Gag protein encoded by the L-A virus, suggesting a host-virus interaction that is more complicated than just Xrn1p-mediated nucleolytic digestion of viral RNAs.

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.

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.

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.

Tripartite Motif (TRIM) ubiquitin ligases act in the innate immune response against viruses. One of the best characterized members of this family, TRIM5α, serves as a potent retroviral restriction factor with activity against HIV. Here, we characterize what are likely to be the youngest TRIM genes in the human genome. For instance, we have identified 11 TRIM genes that are specific to humans and African apes (chimpanzees, bonobos, and gorillas) and another 7 that are human-specific. Many of these young genes have never been described, and their identification brings the total number of known human TRIM genes to approximately 100. These genes were acquired through segmental duplications, most of which originated from a single locus on chromosome 11. Another polymorphic duplication of this locus has resulted in these genes being copy number variable within the human population, with a Han Chinese woman identified as having 12 additional copies of these TRIM genes compared to other individuals screened in this study. Recently, this locus was annotated as one of 34 "hotspot" regions that are also copy number variable in the genomes of chimpanzees and rhesus macaques. Most of the young TRIM genes originating from this locus are expressed, spliced, and contain signatures of positive natural selection in regions known to determine virus recognition in TRIM5α. However, we find that they do not restrict the same retroviruses as TRIM5α, consistent with the high degree of divergence observed in the regions that control target specificity. We propose that this recombinationally volatile locus serves as a reservoir from which new TRIM genes arise through segmental duplication, allowing primates to continually acquire new antiviral genes that can be selected to target new and evolving pathogens.

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.

TRIM5alpha provides a cytoplasmic block to retroviral infection, and orthologs encoded by some primates are active against HIV. Here, we present an evolutionary comparison of the TRIM5 gene to its closest human paralogs: TRIM22, TRIM34, and TRIM6. We show that TRIM5 and TRIM22 have a dynamic history of gene expansion and loss during the evolution of mammals. The cow genome contains an expanded cluster of TRIM5 genes and no TRIM22 gene, while the dog genome encodes TRIM22 but has lost TRIM5. In contrast, TRIM6 and TRIM34 have been strictly preserved as single gene orthologs in human, dog, and cow. A more focused analysis of primates reveals that, while TRIM6 and TRIM34 have evolved under purifying selection, TRIM22 has evolved under positive selection as was previously observed for TRIM5. Based on TRIM22 sequences obtained from 27 primate genomes, we find that the positive selection of TRIM22 has occurred episodically for approximately 23 million years, perhaps reflecting the changing pathogenic landscape. However, we find that the evolutionary episodes of positive selection that have acted on TRIM5 and TRIM22 are mutually exclusive, with generally only one of these genes being positively selected in any given primate lineage. We interpret this to mean that the positive selection of one gene has constrained the adaptive flexibility of its neighbor, probably due to genetic linkage. Finally, we find a striking congruence in the positions of amino acid residues found to be under positive selection in both TRIM5alpha and TRIM22, which in both proteins fall predominantly in the beta2-beta3 surface loop of the B30.2 domain. Astonishingly, this same loop is under positive selection in the multiple cow TRIM5 genes as well, indicating that this small structural loop may be a viral recognition motif spanning a hundred million years of mammalian evolution.

Host genomes have adopted several strategies to curb the proliferation of transposable elements and viruses. A recently discovered novel primate defense against retroviral infection involves a single-stranded DNA-editing enzyme, APOBEC3G, that causes hypermutation of HIV. The HIV-encoded virion infectivity factor (Vif) protein targets APOBEC3G for destruction, setting up a genetic conflict between the APOBEC3G and Vif genes. This kind of conflict leads to rapid fixation of mutations that alter amino acids at the protein-protein interface, referred to as positive selection. We show that the APOBEC3G gene has been subject to strong positive selection throughout the history of primate evolution. Unexpectedly, this selection appears more ancient than, and is likely only partially caused by, modern lentiviruses. Furthermore, five additional APOBEC genes in the human genome appear to be engaged in similar genetic conflicts, displaying some of the highest signals for positive selection in the human genome. Despite being only recently discovered, editing of RNA and DNA may thus represent an ancient form of host defense in primate genomes.