Positive selection in dNTPase SAMHD1 throughout mammalian evolution

Positive selection analyses identify a single WWE domain residue that shapes ZAP into a super restriction factor

The host interferon pathway upregulates intrinsic restriction factors in response to viral infection. Many of them block a diverse range of viruses, suggesting that their antiviral functions might have been shaped by multiple viral families during evolution. Virus-host conflicts have led to the rapid adaptation of viral and host proteins at their interaction hotspots. Hence, we can use evolutionary genetic analyses to elucidate antiviral mechanisms and domain functions of restriction factors. Zinc finger antiviral protein (ZAP) is a restriction factor against RNA viruses such as alphaviruses, in addition to other RNA, retro-, and DNA viruses, yet its precise antiviral mechanism is not fully characterized. Previously, an analysis of 13 primate ZAP identified 3 positively selected residues in the poly(ADP-ribose) polymerase-like domain. However, selective pressure from ancient alphaviruses and others likely drove ZAP adaptation in a wider representation of mammals. We performed positive selection analyses in 261 mammalian ZAP using more robust methods with complementary strengths and identified 7 positively selected sites in all domains of the protein. We generated ZAP inducible cell lines in which the positively selected residues of ZAP are mutated and tested their effects on alphavirus replication and known ZAP activities. Interestingly, the mutant in the second WWE domain of ZAP (N658A) is dramatically better than wild-type ZAP at blocking replication of Sindbis virus and other ZAP-sensitive alphaviruses due to enhanced viral translation inhibition. The N658A mutant inhabits the space surrounding the previously reported poly(ADP-ribose) (PAR) binding pocket, but surprisingly has reduced binding to PAR. In summary, the second WWE domain is critical for engineering a super restrictor ZAP and fluctuations in PAR binding modulate ZAP antiviral activity. Our study has the potential to unravel the role of ADP-ribosylation in the host innate immune defense and viral evolutionary strategies that antagonize this post-translational modification.

Author summary

Host proteins and viral proteins that encounter one another are locked in a perpetual genetic arms race. In this evolutionary race, a mutation that confers a survival advantage will become more frequent in the population. By looking at the sequences of genes that are known to have antiviral roles in mammals, we can identify the exact sites where a host and viral protein have interacted and gain insight into how an antiviral protein works. Here, we identified these sites in zinc finger antiviral protein (ZAP), a host protein that blocks many different viruses. We found that changing one of the sites from the original amino acid to another dramatically improves ZAP's antiviral activity against Sindbis virus, an alphavirus, due to improved inhibition of viral translation. Our mutation is also better at inhibiting other members in the Alphavirus genus. We observed that our mutant ZAP has reduced ability to bind poly(ADP-ribose), a post-translational modification that is targeted by alphaviruses for productive infection. Our findings help us better understand how viruses have shaped the evolution of broad-spectrum host antiviral proteins, with great implications for the engineering of super restriction factors.

Gene editing of SAMHD1 in macrophage-like cells reveals complex relationships between SAMHD1 phospho-regulation, HIV-1 restriction, and cellular dNTP levels

IMPORTANCE

We introduce BLaER1 cells as an alternative myeloid cell model in combination with CRISPR/Cas9-mediated gene editing to study the influence of sterile α motif and HD domain-containing protein 1 (SAMHD1) T592 phosphorylation on anti-viral restriction and the control of cellular dNTP levels in an endogenous, physiologically relevant context. A proper understanding of the mechanism of the anti-viral function of SAMHD1 will provide attractive strategies aiming at selectively manipulating SAMHD1 without affecting other cellular functions. Even more, our toolkit may inspire further genetic analysis and investigation of restriction factors inhibiting retroviruses and their cellular function and regulation, leading to a deeper understanding of intrinsic anti-viral immunity.

Gene editing of SAMHD1 in macrophage-like cells reveals complex relationships between SAMHD1 phospho-regulation, HIV-1 restriction and cellular dNTP levels

Sterile α motif (SAM) and HD domain-containing protein 1 (SAMHD1) is a dNTP triphosphate triphosphohydrolase (dNTPase) and a potent restriction factor for immunodeficiency virus 1 (HIV-1), active in myeloid and resting CD4+ T cells. The anti-viral activity of SAMHD1 is regulated by dephosphorylation of the residue T592. However, the impact of T592 phosphorylation on dNTPase activity is still under debate. Whether additional cellular functions of SAMHD1 impact anti-viral restriction is not completely understood. We report BLaER1 cells as a novel human macrophage HIV-1 infection model combined with CRISPR/Cas9 knock-in (KI) introducing specific mutations into the SAMHD1 locus to study mutations in a physiological context. Transdifferentiated BLaER1 cells harbor active dephosphorylated SAMHD1 that blocks HIV-1 reporter virus infection. As expected, homozygous T592E mutation, but not T592A, relieved a block to HIV-1 reverse transcription. Co-delivery of VLP-Vpx to SAMHD1 T592E KI mutant cells did not further enhance HIV-1 infection indicating the absence of an additional SAMHD1-mediated antiviral activity independent of T592 de-phosphorylation. T592E KI cells retained dNTP levels similar to WT cells indicating uncoupling of anti-viral and dNTPase activity of SAMHD1. The integrity of the catalytic site in SAMHD1 was critical for anti-viral activity, yet poor correlation of HIV-1 restriction and global cellular dNTP levels was observed in cells harboring catalytic core mutations. Together, we emphasize the complexity of the relationship between HIV-1 restriction, SAMHD1 enzymatic function and T592 phospho-regulation and provide novel tools for investigation in an endogenous and physiological context.

Identifying novel host-based diagnostic biomarker panels for COVID-19: a whole-blood/nasopharyngeal transcriptome meta-analysis

BACKGROUND

Regardless of improvements in controlling the COVID-19 pandemic, the lack of comprehensive insight into SARS-COV-2 pathogenesis is still a sophisticated challenge. In order to deal with this challenge, we utilized advanced bioinformatics and machine learning algorithms to reveal more characteristics of SARS-COV-2 pathogenesis and introduce novel host response-based diagnostic biomarker panels.

METHODS

In the present study, eight published RNA-Seq datasets related to whole-blood (WB) and nasopharyngeal (NP) swab samples of patients with COVID-19, other viral and non-viral acute respiratory illnesses (ARIs), and healthy controls (HCs) were integrated. To define COVID-19 meta-signatures, Gene Ontology and pathway enrichment analyses were applied to compare COVID-19 with other similar diseases. Additionally, CIBERSORTx was executed in WB samples to detect the immune cell landscape. Furthermore, the optimum WB- and NP-based diagnostic biomarkers were identified via all the combinations of 3 to 9 selected features and the 2-phases machine learning (ML) method which implemented k-fold cross validation and independent test set validation.

RESULTS

The host gene meta-signatures obtained for SARS-COV-2 infection were different in the WB and NP samples. The gene ontology and enrichment results of the WB dataset represented the enhancement in inflammatory host response, cell cycle, and interferon signature in COVID-19 patients. Furthermore, NP samples of COVID-19 in comparison with HC and non-viral ARIs showed the significant upregulation of genes associated with cytokine production and defense response to the virus. In contrast, these pathways in COVID-19 compared to other viral ARIs were strikingly attenuated. Notably, immune cell proportions of WB samples altered in COVID-19 versus HC. Moreover, the optimum WB- and NP-based diagnostic panels after two phases of ML-based validation included 6 and 8 markers with an accuracy of 97% and 88%, respectively.

CONCLUSIONS

Based on the distinct gene expression profiles of WB and NP, our results indicated that SARS-COV-2 function is body-site-specific, although according to the common signature in WB and NP COVID-19 samples versus controls, this virus also induces a global and systematic host response to some extent. We also introduced and validated WB- and NP-based diagnostic biomarkers using ML methods which can be applied as a complementary tool to diagnose the COVID-19 infection from non-COVID cases.

Recurrent evolution of an inhibitor of ESCRT-dependent virus budding and LINE-1 retrotransposition in primates

Most antiviral proteins recognize specific features of viruses. In contrast, the recently described antiviral factor retroCHMP3 interferes with the "host endosomal complexes required for transport" (ESCRT) pathway to inhibit the budding of enveloped viruses. RetroCHMP3 arose independently on multiple occasions via duplication and truncation of the gene encoding the ESCRT-III factor CHMP3. However, since the ESCRT pathway is essential for cellular membrane fission reactions, ESCRT inhibition is potentially cytotoxic. This raises fundamental questions about how hosts can repurpose core cellular functions into antiviral functions without incurring a fitness cost due to excess cellular toxicity. We reveal the evolutionary process of detoxification for retroCHMP3 in New World monkeys using a combination of ancestral reconstructions, cytotoxicity, and virus release assays. A duplicated, full-length copy of retroCHMP3 in the ancestors of New World monkeys provides modest inhibition of virus budding while exhibiting subtle cytotoxicity. Ancient retroCHMP3 then accumulated mutations that reduced cytotoxicity but preserved virus inhibition before a truncating stop codon arose in the more recent ancestors of squirrel monkeys, resulting in potent inhibition. In species where full-length copies of retroCHMP3 still exist, their artificial truncation generated potent virus-budding inhibitors with little cytotoxicity, revealing the potential for future antiviral defenses in modern species. In addition, we discovered that retroCHMP3 restricts LINE-1 retrotransposition, revealing how different challenges to genome integrity might explain multiple independent origins of retroCHMP3 in different species to converge on new immune functions.

Positive natural selection in primate genes of the type I interferon response

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01783-z.

Selective Disruption of SERINC5 Antagonism by Nef Impairs SIV Replication in Primary CD4+ T Cells

The Nef proteins of HIV-1 and SIV enhance viral infectivity by preventing the incorporation of the multipass transmembrane protein serine incorporator 5 (SERINC5), and to a lesser extent SERINC3, into virions. In addition to counteracting SERINCs, SIV Nef also downmodulates several transmembrane proteins from the surface of virus-infected cells, including simian tetherin, CD4 and MHC class I (MHC I) molecules. From a systematic analysis of alanine substitutions throughout the SIV_mac239 Nef protein, we identified residues that are required to counteract SERINC5. This information was used to engineer an infectious molecular clone of SIV (SIV_mac239nef _AV), which differs by two amino acids in the N-terminal domain of Nef that make the virus sensitive to SERINC5 while retaining other activities of Nef. SIV_mac239nef _AV downmodulates CD3, CD4, MHC I and simian tetherin, but cannot counteract SERINC5. In primary rhesus macaque CD4⁺ T cells, SIV_mac239nef _AV exhibits impaired infectivity and replication compared to wild-type SIV_mac239. These results demonstrate that SERINC5 antagonism can be separated from other Nef functions and reveal the impact of SERINC5 on lentiviral replication.Importance: SERINC5, a multipass transmembrane protein, is incorporated into retroviral particles during assembly. This leads to a reduction of particle infectivity by inhibiting virus fusion with the target cell membrane. The Nef proteins of HIV-1 and SIV enhance viral infectivity by preventing the incorporation of SERINC5 into virions. However, the relevance of this restriction factor in viral replication has not been elucidated. Here we report a systematic mapping of Nef residues required for SERINC5 antagonism. Counter screens for three other functions of Nef helped identify two residues in the N-terminal domain of Nef, which when mutated make Nef selectively susceptible to SERINC5. Since Nef is multi-functional, genetic separation of SERINC5 antagonism from its other functions affords comparison of the replication of isogenic viruses that are or are not sensitive to SERINC5. Such a strategy revealed the impact of SERINC5 on SIV replication in primary rhesus macaque CD4⁺ T-cells.

SAMHD1 … and Viral Ways around It

The SAM and HD domain-containing protein 1 (SAMHD1) is a dNTP triphosphohydrolase that plays a crucial role for a variety of different cellular functions. Besides balancing intracellular dNTP concentrations, facilitating DNA damage repair, and dampening excessive immune responses, SAMHD1 has been shown to act as a major restriction factor against various virus species. In addition to its well-described activity against retroviruses such as HIV-1, SAMHD1 has been identified to reduce the infectivity of different DNA viruses such as the herpesviruses CMV and EBV, the poxvirus VACV, or the hepadnavirus HBV. While some viruses are efficiently restricted by SAMHD1, others have developed evasion mechanisms that antagonize the antiviral activity of SAMHD1. Within this review, we summarize the different cellular functions of SAMHD1 and highlight the countermeasures viruses have evolved to neutralize the restriction factor SAMHD1.

Retroviral Restriction Factors and Their Viral Targets: Restriction Strategies and Evolutionary Adaptations

The evolutionary conflict between retroviruses and their vertebrate hosts over millions of years has led to the emergence of cellular innate immune proteins termed restriction factors as well as their viral antagonists. Evidence accumulated in the last two decades has substantially increased our understanding of the elaborate mechanisms utilized by these restriction factors to inhibit retroviral replication, mechanisms that either directly block viral proteins or interfere with the cellular pathways hijacked by the viruses. Analyses of these complex interactions describe patterns of accelerated evolution for these restriction factors as well as the acquisition and evolution of their virus-encoded antagonists. Evidence is also mounting that many restriction factors identified for their inhibition of specific retroviruses have broader antiviral activity against additional retroviruses as well as against other viruses, and that exposure to these multiple virus challenges has shaped their adaptive evolution. In this review, we provide an overview of the restriction factors that interfere with different steps of the retroviral life cycle, describing their mechanisms of action, adaptive evolution, viral targets and the viral antagonists that evolved to counter these factors.

SAMHD1 as the Potential Link Between SARS-CoV-2 Infection and Neurological Complications

The recent pandemic of coronavirus infectious illness 2019 (COVID19) triggered by SARS-CoV-2 has rapidly spread around the globe, generating in severe events an acute, highly lethal pneumonia and death. In the past two hitherto similar CoVs, the severe acute respiratory syndrome CoV (SARS-CoV-1) and Middle East respiratory syndrome CoV (MERS-CoV) also gained universal attention as they produced clinical symptoms similar to those of SARS-CoV-2 utilizing angiotensin-converting enzyme 2 (ACE2) receptor and dipeptidyl peptidase 4 (DPP4) to go into the cells. COVID-19 may also present with overtly neurological symptoms. The proper understanding of the expression and dissemination of ACE2 in central and peripheral nerve systems is crucial to understand better the neurological morbidity caused by COVID-19. Using the STRING bioinformatic tool and references through text mining tools associated to Coronaviruses, we identified SAMHD1 as the probable link to neurological symptoms. Paralleled to the response to influenza A virus and, specifically, respiratory syncytial virus, SARS-CoV-2 evokes a response that needs robust induction of a subclass of cytokines, including the Type I and, obviously, Type III interferons as well as a few chemokines. We correlate ACE2 to the pathogenesis and neurologic complications of COVID-19 and found that SAMHD1 links to NF-κB pathway. No correlation was found with other molecules associated with Coronavirus infection, including ADAR, BST2, IRF3, IFITM3, ISG15, MX1, MX2, RNASEL, RSAD2, and VPRBP. We suggest that SAMHD1 is the molecule that may be behind the mechanisms of the neurological complications associated with COVID-19.

The C-terminal domain of feline and bovine SAMHD1 proteins has a crucial role in lentiviral restriction

SAM and HD domain-containing protein 1 (SAMHD1) is a host factor that restricts reverse transcription of lentiviruses such as HIV in myeloid cells and resting T cells through its dNTP triphosphohydrolase (dNTPase) activity. Lentiviruses counteract this restriction by expressing the accessory protein Vpx or Vpr, which targets SAMHD1 for proteasomal degradation. SAMHD1 is conserved among mammals, and the feline and bovine SAMHD1 proteins (fSAM and bSAM) restrict lentiviruses by reducing cellular dNTP concentrations. However, the functional regions of fSAM and bSAM that are required for their biological functions are not well-characterized. Here, to establish alternative models to investigate SAMHD1 in vivo, we studied the restriction profile of fSAM and bSAM against different primate lentiviruses. We found that both fSAM and bSAM strongly restrict primate lentiviruses and that Vpx induces the proteasomal degradation of both fSAM and bSAM. Further investigation identified one and five amino acid sites in the C-terminal domain (CTD) of fSAM and bSAM, respectively, that are required for Vpx-mediated degradation. We also found that the CTD of bSAM is directly involved in mediating bSAM's antiviral activity by regulating dNTPase activity, whereas the CTD of fSAM is not. Our results suggest that the CTDs of fSAM and bSAM have important roles in their antiviral functions. These findings advance our understanding of the mechanism of fSAM- and bSAM-mediated viral restriction and might inform strategies for improving HIV animal models.

Integrated structural and evolutionary analysis reveals common mechanisms underlying adaptive evolution in mammals

Understanding the molecular basis of adaptation to the environment is a central question in evolutionary biology, yet linking detected signatures of positive selection to molecular mechanisms remains challenging. Here we demonstrate that combining sequence-based phylogenetic methods with structural information assists in making such mechanistic interpretations on a genomic scale. Our integrative analysis shows that positively selected sites tend to colocalize on protein structures and that positively selected clusters are found in functionally important regions of proteins, indicating that positive selection can contravene the well-known principle of evolutionary conservation of functionally important regions. This unexpected finding, along with our discovery that positive selection acts on structural clusters, opens previously unexplored strategies for the development of better models of protein evolution. Remarkably, proteins where we detect the strongest evidence of clustering belong to just two functional groups: Components of immune response and metabolic enzymes. This gives a coherent picture of pathogens and xenobiotics as important drivers of adaptive evolution of mammals.

ChromaClade: combined visualisation of phylogenetic and sequence data

BACKGROUND

Studying site-specific amino acid frequencies by eye can reveal biologically significant variability and lineage-specific adaptation. This so-called 'sequence gazing' often informs bioinformatics and experimental research. But it is important to also account for the underlying phylogeny, since similarities may be due to common descent rather than selection pressure, and because it is important to distinguish between founder effects and convergent evolution. We set out to combine phylogenetic and sequence data to produce evolutionarily insightful visualisations.

RESULTS

We present ChromaClade, a convenient tool with a graphical user-interface that works in concert with popular tree viewers to produce colour-annotated phylogenies highlighting residues found in each taxon and at each site in a sequence alignment. Colouring branches according to residues found at descendent tips also quickly identifies lineage-specific residues and those internal branches where key substitutions have occurred. We demonstrate applications of ChromaClade to human immunodeficiency virus and influenza A virus datasets, illustrating cases of conservative, adaptive and convergent evolution.

CONCLUSIONS

We find this to be a powerful approach for visualising site-wise residue distributions and detecting evolutionary patterns, especially in large datasets. ChromaClade is available for Windows, macOS and Unix or Linux; program executables and source code are available at github.com/chrismonit/chroma_clade .