Next Article in Journal
A Review of the Bromodomain and Extraterminal Domain Epigenetic Reader Proteins: Function on Virus Infection and Cancer
Next Article in Special Issue
Translation Inhibition Mediated by Interferon-Stimulated Genes during Viral Infections
Previous Article in Journal
Clinical Application of Adenovirus (AdV): A Comprehensive Review
Previous Article in Special Issue
Factor-Dependent Internal Ribosome Entry Site and -1 Programmed Frameshifting Signal in the Bemisia-Associated Dicistrovirus 2
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 16
Issue 7
10.3390/v16071095
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tanapox Virus and Yaba Monkey Tumor Virus K3 Orthologs Inhibit Primate Protein Kinase R in a Species-Specific Fashion

by
Dewi Megawati
1,2,
Jeannine N. Stroup
1,
Chorong Park
1,†,
Taylor Clarkson
1,‡,
Loubna Tazi
1,
Greg Brennan
1 and
Stefan Rothenburg
1,*
1
Department of Medical Microbiology and Immunology, School of Medicine, University of California, Davis, CA 95616, USA
2
Department of Microbiology and Parasitology, Faculty of Medicine and Health Sciences, Warmadewa University, Denpasar 80239, Bali, Indonesia
*
Author to whom correspondence should be addressed.
†
Current Address: Department of Microbiology-Immunology, Feinberg School of Medicine, Northwestern University, Chicago, IL 60208, USA.
‡
Current Address: Colorado Department of Public Health and Environment, Denver, CO 80246, USA.
Viruses 2024, 16(7), 1095; https://doi.org/10.3390/v16071095
Submission received: 7 June 2024 / Revised: 2 July 2024 / Accepted: 3 July 2024 / Published: 8 July 2024
(This article belongs to the Special Issue Viral Strategies and Cellular Countermeasures that Regulate mRNA Access to the Translation Apparatus)
Browse Figures
Review Reports Versions Notes

Abstract

:
Yaba monkey tumor virus (YMTV) and Tanapox virus (TPV) are members of the Yatapoxvirus genus and can infect humans and other primates. Despite the threat posed by yatapoxviruses, the factors determining their host range are poorly understood. In this study, we analyzed the ability of YMTV and TPV orthologs of vaccinia virus K3 (called 012 in YMTV and TPV), which share 75% amino acid identity with one another, to inhibit PKR from 15 different primate species. We first used a luciferase-based reporter, and found that YMTV and TPV K3 orthologs inhibited PKR in a species-specific manner and showed distinct PKR inhibition profiles. TPV 012 inhibited PKR from 11 primates, including humans, substantially better than YMTV 012. In contrast, both K3 orthologs inhibited the other four primate PKRs comparably well. Using YMTV 012 and TPV 012 hybrids, we mapped the region responsible for the differential PKR inhibition to the C- terminus of the K3 orthologs. Next, we generated chimeric vaccinia virus strains to investigate whether TPV K3 and YMTV K3 orthologs could rescue the replication of a vaccinia virus strain that lacks PKR inhibitors K3L and E3L. Virus replication in primate-derived cells generally correlated with the patterns observed in the luciferase-based assay. Together, these observations demonstrate that yatapoxvirus K3 orthologs have distinct PKR inhibition profiles and inhibit PKR in a species-specific manner, which may contribute to the differential susceptibility of primate species to yatapoxvirus infections.

1. Introduction

Yatapoxviruses are a small group of Chordopoxvirinae that are pathogenic to primates, including humans. The Yatapoxvirus genus includes Yaba monkey tumor virus (YMTV) and Tanapox virus (TPV) [1]. Previously, Yaba-like disease virus (YLDV) was considered a third member; however, its genome shares 98.68% sequence identity with TPV, and can be considered a different variant of the same virus [2]. The genomes of YMTV and TPV share 75% sequence identity [3]. TPV was first isolated from human skin biopsies during a TPV outbreak in 1957 in Tana River Valley, Kenya [4,5]. YMTV was first isolated from an outbreak of subcutaneous tumors in monkeys in Yaba, Nigeria, in 1958 [6]. Yatapoxvirus genomes are A+T rich and range in size from 135 to 145 kb [2,3,7,8]. TPV and YMTV exhibit differences in disease presentations. TPV infection is characterized by vesicular skin lesions, whereas YMTV infection produces a localized histiocyte-filled tumor (histiocytomas) [9].
While both TPV and YMTV infect a wide range of primates, they appear to exhibit host specificity. Serological surveys have indicated that TPV is endemic in African and Malaysian monkeys, but not in Indian rhesus macaques or New World monkeys [4,5,10]. Human Tanapox disease is considered endemic to several regions of Africa and causes febrile illness and vesicular skin lesions similar to those produced in non-human primates. Sporadic cases have been identified in 30 locations spanning 6000 km from Sierra Leone to Tanzania, with larger outbreaks occurring from time to time [11]. Similarly, antibodies against YMTV have been found in great apes and Old World monkeys but not in New World monkeys or Indian rhesus macaques [10]. A serological survey of 456 primate sera including 26 chimpanzees, 326 Old World monkeys (African green monkeys, patas monkeys, baboons, colobus monkeys, and rhesus macaques), and 104 New World monkeys (spider monkeys, squirrel monkeys, owl monkeys, marmosets, and capuchin monkeys) indicated that antibodies against YMTV were evident in chimpanzees and Old World monkeys but not in any of the New World monkeys [10].
Following experimental subcutaneous inoculation with YMTV, tumor-like masses were detected in Asiatic rhesus macaques (Macaca mulatta, Macaca irus, and Macaca speciosa) but not in African green monkeys, mangabey monkeys, patas monkeys, mice, rats, rabbits, or guinea pigs [6,9,12], indicating host specificity. While no natural infection of YMTV in humans has been reported, the infection of human volunteers as well as an accidental infection of a laboratory worker have resulted in the development of mild histiocytomas [13]. Despite the threat posed to human health, the factors determining the host range of yatapoxviruses are poorly understood. In addition, vaccination with vaccinia virus does not appear to protect against yatapoxvirus infections [13,14]. Thus, it is important to study factors that determine the host range of yatapoxviruses.
A major determinant for the tropism of many viruses is the attachment and entry via host-specific cell receptors [15]. However, as poxviruses use ubiquitous cell surface receptors for entry, it is believed that the evasion of the innate immune response post entry is a more important determinant of poxvirus tropism [16,17]. Protein Kinase R (PKR) is a prominent host restriction factor against poxvirus infection (reviewed in [18]). PKR exists in a monomeric inactive form and is activated upon binding to double-stranded RNA (dsRNA) produced during infection by most viruses, which leads to PKR dimerization and autophosphorylation. Activated PKR phosphorylates the alpha subunit of eukaryotic initiation factor 2 (eIF2) [19]. Phosphorylated eIF2α has a high binding affinity for the regulatory core of the guanine nucleotide exchange factor eIF2B and prevents eIF2B from catalyzing GDP-GTP exchange to eIF2α [20]. As GTP-bound eIF2 is required for translation initiation, the resulting low availability of GTP-eIF2 leads to the shut-off of eIF2-dependent protein synthesis, including that of viral proteins. To overcome the antiviral effects of PKR activity, most poxviruses encode two PKR inhibitors, named E3 and K3 in vaccinia virus. E3 binds dsRNA and prevents PKR dimerization, whereas K3 inhibits PKR by acting as a pseudosubstrate for the eIF2α binding site [21,22,23,24].
It was previously shown that the inhibition of host PKR by K3 orthologs of several poxvirus families was a key determinant of host specificity [22,25,26,27,28]. Vaccinia virus lacking the E3L gene (VACVΔE3L) but maintaining the K3L gene was unable to replicate in HeLa cells but remained replication-competent in BHK21 cells, suggesting a role for K3 in determining host range [22]. This finding was supported by recent studies on leporipoxvirus, capripoxvirus, and orthopoxvirus K3 orthologs. The myxoma virus (MYXV) K3 ortholog M156 specifically inhibited rabbit PKR but failed to inhibit other PKR species [25]. Similarly, capripoxvirus K3 orthologs strongly inhibited human, goat, and sheep PKR, but only weakly inhibited mouse and cow PKR [26]. Our studies on the inhibition of PKR from a panel of mammalian species by orthopoxvirus K3 orthologs exhibited distinct inhibition profiles. Here, we used a luciferase-based assay and infections of mammalian cell lines to examine the ability of TPV and YMTV K3 orthologs to inhibit PKR derived from 15 primate species. Our results demonstrate that TPV and YMTV K3 orthologs have distinct PKR inhibition profiles and that they inhibit primate PKR in a species-specific manner.

2. Materials and Methods

2.1. Cell Lines

Tert immortalized gibbon fibroblasts (gibbon Tert-GF), HeLa cells (human, ATCC #CCL-2), HeLa PKR-knock-out (PKRko) cells, and BSC-40 cells (ATCC CRL-2761) were kindly provided by Dr. Adam Geballe [29]. RK13 cells (rabbit) expressing E3 and K3 (designated RK13+E3L+K3L) were previously described [30]. The cells were grown in Dulbecco’s Modified Eagle’s Medium (DMEM, Life Technologies, Carlsbad, CA, USA) supplemented with 5% fetal bovine serum (FBS, Lonza Bioscience, Walkersville, MD, USA) or 10% FBS and 100 IU/mL penicillin/streptomycin (Gibco, Life Technologies, Carlsbad, CA, USA). The RK13+E3L+K3L cell culture medium contained 500 μg/mL geneticin and 300 μg/mL zeocin (Life Technologies).

2.2. Plasmids

A library of primate PKR expressing plasmids were kindly provided by Dr. Nels Elde [31]. Primate PKR and viral antagonist genes VACV K3L (YP_232916.1), TPV 012 (EF420156.1), YMTV 012 (AY386371.1) were subcloned into the pSG5 expression vector for luciferase-based reporter assays. pGL3 luciferase reporter vector was purchased from Promega, Madison, WI, USA. Cloning was carried out using Gibson assembly techniques (New England Biolabs, Ipswich, MA, USA) [32]. To monitor VACV K3L, TPV 012, and YMTV 012 gene expressions, these genes were cloned into C-terminus DYK-tagged pSG5 vector for transient transfection assay. To generate hybrid NTPV-CYMTV and NYMTV-CTPV K3 orthologs, we PCR-amplified NYMTV 012 (nucleotide 1–99), CYMTV 012 (nucleotide 100–267), NTPV 012 (nucleotide 1–99), and CTPV 012 (nucleotide 100–267) and cloned the TPV 012 and YMTV 012 hybrid genes into pSG5 vector using Gibson assembly. To generate recombinant VACV expressing TPV 012 and YMTV 012 genes, the yatapoxvirus K3 ortholog genes were cloned into p837-GOI-mCherry-E3L, as previously described [33].

2.3. Luciferase-Based Reporter Assays

Luciferase-based assays were performed as previously described [34]. Briefly, 5 × 104 HeLa PKRko cells were seeded per well in 24-well plates overnight. The HeLa PKRko cells were co-transfected with 200 ng of the indicated PKR expression vector, 200 ng of each viral antagonist expression vector, and 50 ng of pGL3 firefly luciferase expression vector (Promega), following the manufacturer’s protocol. Cells were lysed with mammalian lysis buffer (GE Healthcare, Chicago, IL, USA) at 48 h post transfection. Luciferase activity was measured using a GloMax luminometer (Promega) by adding luciferin (Promega) reagent to the cell lysates, as per the manufacturer’s recommendations. Data are presented as relative luciferase activity in which all data were normalized to pSG5 empty vector. Experiments were conducted in triplicate for each of the three independent experiments.

2.4. Virus and Infections

Vaccinia virus variant vP872 [35] was kindly provided by Dr. Bertram Jacobs. VC-R4+VACV K3L and VC-R4+sheeppox 011 (SPPV 011), as well as the generation of VC-R4 from vP872 variant, were previously described [26,33]. The generation of chimeric vaccinia virus, VC-R4+TPV 012, and VC-R4+YMTV 012, was carried out through the scarless integration of the open reading frames of the yatapoxvirus K3L orthologs into the E3L locus [33]. The chimeric viruses were plaque-purified twice, and the K3L ortholog gene integrations were confirmed by means of Sanger sequencing. The viruses were purified by means of zonal sucrose gradient centrifugation, and the virus titer was determined on confluent 12-well plates of RK13+E3L+K3L cells. Plaque assays were performed with confluent 6-well plates of the indicated cell lines, which were infected with 50 plaque-forming units (pfu) of each indicated virus. One hour post infection, the media were replaced with DMEM containing 5% FBS, 1% carboxymethylcellulose (CMC), and 100 U P/S. After 48 h, cells were stained with 0.1% crystal violet, and the excess staining was washed with water. The plates were imaged using an iBright Imaging System (Invitrogen, Waltham, MA, USA). Virus infection assays were performed in confluent 6-well plates of the indicated cells. The cells were infected with each indicated virus at a MOI of 0.01 to assess replication in multiple-cycle replication assays. Cells and supernatants were collected at 30 h post infection and subjected to three rounds of freezing at −80 °C and thawing at 37 °C. Lysates were sonicated twice for 15 s, at 50% amplitude (Qsonica Q500, Newtown, CT, USA). Virus titers were measured in plaque-forming units per mL (pfu/mL) on RK13+E3L+K3L cells. Virus titer data were analyzed using one-way ANOVA followed by Tukey’s multiple comparisons test (GraphPad Prism version 10.2.3).

2.5. PCR of Viral Genomic DNA

HeLa PKRko cells were seeded in 10 cm dishes to a confluency of 90–100%. The cells were then infected with the indicated viruses at a MOI of 0.1 for 24 h. Viral genomic DNA (gDNA) extraction was performed as previously described [36]. About 100 ng of the isolated viral gDNA was used as a template in a PCR targeting K3L ortholog genes using Phusion High Fidelity DNA polymerase (NEB #M0530L). The forward primer sequence was 5′ GACGAACCACCAGAGGATGATG 3′ and the reverse primer sequence was 5′ AGTACTACAATTTGTATTTTTTAATCTATCTCA 3′. PCR products were gel-purified with Monarch DNA Gel Extraction Kit (NEB #T1020) and Sanger sequencing was performed to confirm the correct integration of the K3L orthologs.

2.6. Immunoblot Analyses

To determine the level of K3 ortholog expression in a transient transfection system, 4 × 105 HeLa PKRko cells were seeded in a 6-well plate. On the next day, the cells were transfected with 1 µg of pSG5 empty vector, DYK-tagged VACV K3L, DYK-tagged TPV 012, and DYK-tagged YMTV 012, using a Genjet to DNA ratio of 2:1. At 48 h post transfection, cells were washed with PBS, lysed with 1% sodium dodecyl sulfate (SDS) in DPBS, and sonicated at 50% amplitude for 5 s twice. About 10 µg of proteins were run on 12% SDS polyacrylamide gels (Biorad, Hercules, CA, USA) and transferred to polyvinyl difluoride (PVDF, GE Healthcare) membranes. The membranes were blocked with SuperBlock blocking buffer (Thermofisher, Waltham, MA, USA) for 1 h and probed with primary antibodies against Flag M2 (Sigma F1804, St. Louis, MO, USA) and beta-actin (Sigma A1978) at a dilution of 1:2000 in SuperBlock blocking buffer (Thermofisher) for 1 h at room temperature. After washing with Tris Buffered Saline 0.1% Tween® 20 (TBST, Fisher Bioreagents), the membranes were probed with secondary antibody horseradish peroxidase-conjugated donkey anti-mouse (Jackson ImmunoResearch Laboratories Inc., 715-035-150, West Grove, PA, USA) at a dilution of 1:10,000 in 1% (w/v) non-fat milk dissolved in TBST for 1 h. The membranes were washed three times for 10 min, and proteins were detected with Amersham™ ECL™ (GE Healthcare). Images were taken using the iBright Imaging System (Invitrogen).
To determine the level of K3L ortholog expression by the chimeric vaccinia virus, 8 × 105 HeLa PKRko cells were seeded in a 6-well plate. On the next day, the cells were infected with VC-R4, VC-R4+VACV K3L, VC-R4+TPV 012, and VC-R4+YMTV 012 at a MOI of 3 to ensure that the vast majority of cells were infected. Protein lysates were collected at 24 h post infection. A total of 12 µg of protein was electrophoresed using 12% SDS polyacrylamide gels and transferred to polyvinyl difluoride (PVDF, GE Healthcare) membranes. The membranes were probed with primary antibodies for beta-actin (Sigma A1978) at a dilution of 1:2000 in Thermofisher Superblock blocking buffer, primary antibodies against TPV 012 (Genscript, Piscataway, NJ, USA) at a dilution of 1:1000 in 5% milk dissolved in TBST, and VACV-K3 at a dilution of 1:500 in 5% milk dissolved in TBST overnight at 4 °C. Anti-TPV 012 (cVQVIRTDKLKGYVDVRHIT) and anti-VACV K3 (cKVIRVDYTKGYIDVNYKRM) were custom-produced using peptide–KLH conjugate in New Zealand rabbit (GenScript). After being washed with TBST three times, the membranes were probed with horseradish peroxidase-conjugated donkey anti-rabbit IgG (H+L) secondary antibody (Invitrogen, A16023) at a dilution of 1:10,000 in 1% milk for anti-TPV 012 and anti-VACV K3, and with horseradish peroxidase-conjugated donkey anti-mouse IgG (H+L) (Jackson ImmunoResearch Laboratories Inc., 715-035-150) at a 1:10,000 dilution in TBST containing 1% (w/v) non-fat milk for beta-actin. The membranes were washed three times for 10 min, and proteins were detected with Amersham™ ECL™ (GE Healthcare). Images were taken using the iBright Imaging System (Invitrogen).

2.7. Phylogenetic Analysis

The amino-acid sequences of PKR from 21 different primate species (Table 1) were used to generate the phylogenetic tree. The multiple sequence alignment was performed using MUSCLE [37], and the resulting alignment was used to construct a phylogenetic tree using the maximum likelihood approach, as implemented in the program PhyML 3.0 [38]. The nodal support of the phylogenetic tree was assessed by means of a bootstrap analysis of 100 replicates. The resulting phylogenetic tree was visualized using FigTree v1.4.4 (available at: http://tree.bio.ed.ac.uk/software/figtree/).

3. Results

3.1. Species-Specific Inhibition of Primate PKRs by Yatapoxvirus K3 Orthologs

The inhibition of PKR by K3 orthologs has been shown to be a key determinant for the host range of poxviruses [24,25,26,27,28,34]. In this study, we examined the role of the K3 orthologs TPV 012 and YMTV 012 in inhibiting PKR derived from a diverse set of primate species. The comparison of the amino acid sequences between TPV and YMTV 012 revealed 22 amino acid differences (75% sequence identity). Yatapoxvirus K3 orthologs and VACV K3 showed about 36% sequence identity (Figure 1A,B). We cloned these three orthologs into the pSG5 expression vector with two C-terminal FLAG tags and transfected them into HeLa PKRko cells. The Western blot analysis of whole cell lysates showed that vaccinia K3 and TPV 012 were expressed at comparable levels, whereas YMTV 012 was expressed at higher level than the former (Figure 1C). We also raised polyclonal antibodies against a conserved epitope of TPV 012 (Figure 1A). Western blot HeLa PKRko cell lysates transiently transfected with untagged K3L orthologs showed a higher expression of YMTV 012 (Figure 1D). A previously generated antibody raised against VACV K3 detected both VACV K3 and YMTV 012 (Figure 1D), which could be explained by shared epitopes.
We next tested the sensitivities of 15 different primate species to VACV K3 and the yatapoxvirus K3 orthologs in an established luciferase-based expression assay in which PKR-deficient HeLa cells are co-transfected with a luciferase reporter plasmid, PKR, and PKR inhibitors [25,34]. In this assay, PKR is activated by dsRNA that is formed by overlapping transcripts generated from the transfected plasmids [40] and suppresses the translation of reporter mRNA, which serves as a proxy to quantify relative translation. Luciferase activity is thus inversely correlated with PKR activity, and the inhibition of PKR leads to elevated luciferase activity. We first co-transfected cells with PKRs from 15 primate species representing apes (human, gorilla, chimpanzee, orangutan, and white-cheeked gibbon), Old World monkeys (colobus monkey, François’ leaf monkey, baboon, sooty mangabey, rhesus macaques, talapoin monkey, African green monkey, and patas monkey) and New World monkeys (tamarin and dusky titi), together with the luciferase reporter plasmid to assess PKR activity. All PKRs showed a strong and comparable inhibition of luciferase activity (Figure 2A). We next co-transfected each PKR with either K3 ortholog and the luciferase reporter plasmid (Figure 2B). Most primate PKRs tested were largely resistant to VACV K3 (less than 2-fold inhibition). Orangutan and white-cheeked gibbon PKR were weakly inhibited (about 2-fold), while rhesus macaque and tamarin PKR were inhibited at intermediate low levels (about 4-fold). In total, 11 of the 15 primate PKRs were inhibited substantially better by TPV 012 than by YMTV 012. PKRs from chimpanzee, white-cheeked gibbon, colobus monkey, and François’ leaf monkey were inhibited comparably well by TPV 012 and YMTV 012. In addition to different inhibition profiles when comparing TPV 012 with YMTV 012, species-specific differences in PKR susceptibility were observed, with orangutan, François leaf monkey, and rhesus macaque PKR being most strongly inhibited, and nine other PKRs being intermediately inhibited by TPV 012. In contrast, nine tested PKRs were only weakly inhibited by YMTV 012.
We also projected the relative sensitivities of the tested PKRs to the K3 orthologs on a phylogenetic tree generated with PKR amino acid sequences from a panel of 21 primate species using rabbit PKR as an outgroup to visualize the relatedness among the primate species PKRs (Figure 3). This tree recapitulates the relatedness of primates PKRs from an earlier study. The closer relatedness of human PKR to gorilla PKR than to chimpanzee PKR, which is different from the generally accepted primate phylogeny, likely is a result of positive selection during PKR evolution [31]. We used a previously described scale to grade relative PKR inhibition [28], with 1 for no or very weak inhibition (<2-fold inhibition), 2 for weak inhibition (2- to 3-fold inhibition), 3 for intermediate low inhibition (3- to 5-fold inhibition), 4 for intermediate high inhibition (5- to 7-fold inhibition), 5 for high inhibition (7- to 10-fold inhibition), and 6 for very high inhibition (>10-fold inhibition).

3.2. Differential PKR Inhibition Was Governed by the C-Terminus of Yatapoxvirus K3 Orthologs

The amino acid sequence alignment of YMTV and TPV K3 orthologs identified 22 amino acid differences between the two orthologs, which are distributed throughout the protein sequence. To investigate which regions are important for the differential PKR inhibition by yatapoxvirus K3 orthologs, we generated constructs encoding hybrid YMTV-TPV K3 orthologs in which the C-terminal regions of the two inhibitors were swapped. We took advantage of the shared region in the middle part of the gene (FTVFLPEFG) to separate the N-and C-termini, which leaves 11 amino acid differences in each terminus (Figure 4A). We analyzed the ability of these hybrid proteins to inhibit human, colobus monkey, and dusky titi PKR using the luciferase assay. The hybrid NTPV-CYMTV exhibited inhibition profiles that were more similar to YMTV 012 for each PKR tested (Figure 4B). Correspondingly, the NYMTV-CTPV hybrid showed a similar inhibition pattern to wild-type TPV 012. These results indicate that the C-terminus of TPV 012 and YMTV 012 is important for the observed differential PKR inhibition, although we cannot exclude a contribution of the N-terminus.

3.3. Chimeric Viruses Expressing TPV 012 or YMTV 012 Displayed Cell-Type-Specific Differences in Plaque Formation and Virus Replication

Next, we investigated whether the ability of TPV 012 and YMTV 012 orthologs to inhibit primate PKRs correlated with the ability of the viral genes to rescue the replication of a VACV strain that lacks PKR inhibitors in primate-derived cell lines. The VACV strain VC-R4 (VACVΔE3LΔK3L) is a highly attenuated vaccinia virus variant that can only replicate in PKR-deficient cell lines or cell lines expressing PKR antagonists [33]. We inserted either TPV 012 or YMTV 012 into VC-R4 at the E3L locus using the scarless integration method as previously described, designating the chimeric viruses VC-R4-TPV 012 or VC-R4-YMTV 012, respectively [33]. To analyze the expression of the K3 orthologs by these chimeric viruses, we infected HeLa PKRko cells and performed Western blots using custom-made antibodies at 24 h post infection. In agreement with the transfection data (Figure 1), YMTV 012 was detected at a higher level compared to TPV 012 (Figure 5A). We next examined the plaque formation of these viruses in comparison to VACV vP872 (which lacks K3L but contains the second VACV PKR inhibitor E3), VC-R4+VACV K3L, and VC-R4+SPPV 011 which expresses the sheeppox virus K3 ortholog [26], in different primate-derived cell lines: HeLa (human), gibbon Tert-GF (gibbon), and BSC-40 and PRO1190 (African green monkey). As control, we used RK13+K3L+E3L cells (rabbit), which express VACV E3 and K3 and in which VC-R4 is permissive. We infected these cells with 50 PFU/well. As expected, all viruses formed plaques in permissive RK13+K3L+E3L cells. VC-R4 only formed plaques in RK13+K3L+E3L cells, and vP872 and VC-R4+SPPV 011 formed plaques in all cells tested (Figure 5B).
VC-R4+VACV K3L developed small plaques only in BSC-40 cells. VC-R4+TPV 012 infection led to the formation of plaques in all tested cells with comparable sizes to plaques formed by both VC-R4+SPPV 011 and vP872. Relative to VC-R4+TPV 012, VC-R4+YMTV 012 formed smaller plaques in HeLa and gibbon Tert-GF cells and failed to cause plaque formation in both African green monkey cells. In general, based on the plaque assay, TPV 012 rescued VC-R4 replication substantially better than YMTV 012 in primate-derived cell lines, despite the higher expression of YMTV 012.
We next measured virus titers to determine the capacity of yatapoxvirus K3 orthologs to restore the replication of VC-R4. We infected RK13+K3L+E3L and primate-derived cell lines HeLa (human), gibbon Tert-GF, BSC-40 (African green monkeys), and PRO1190 (African green monkeys) cells with either VC-R4, vP872, VC-R4+VACV K3L, VC-R4+SPPV 11, VC-R4+TPV 012, or VC-R4+YMTV 012 at a MOI of 0.01. Cell lysates were collected at 30 h post infection, and the virus titer was determined in RK13+K3L+E3L cells. Consistent with the plaque assays and luciferase assays, all viruses replicated to comparable titers in RK13+K3L+E3L cells (Figure 6). VC-R4 did not replicate in any cells tested (virus titer below 103 pfu/mL), except in RK13+K3L+E3L cells. VACV K3L rescued VC-R4 virus replication in all cells tested to variable extents, with virus titers about 10- to 100-fold lower than vP872. In HeLa and gibbon Tert-GF cells, SPPV 011 and TPV 012 both restored the replication of VC-R4 to a level comparable to vP872. In BSC-40 and PRO1190 cells, TPV 012 and SPPV 011 rescued VC-R4 replication comparably. Consistent with the plaque assays, VC-R4+YMTV 012 replicated to slightly lower levels than VC-R4+TPV 012 in HeLa cells, although this difference was not significant. In contrast, VC-R4+YMTV 012 significantly replicated with less efficiency than VC-R4+TPV 012 in gibbon Tert-GF cells (11-fold lower), BSC-40 cells (1394-fold lower), and PRO1190 cells (121-fold lower).

4. Discussion

Vaccinia virus K3 is a poxvirus host range factor that acts by inhibiting the host restriction factor PKR. Homology between the K3 and the S1 domain of eIF2α allows K3 to function as a pseudosubstrate of PKR, competitively inhibiting eIF2α phosphorylation [19]. This interaction between PKR and K3 has been studied using yeast and mammalian cell assay systems [41,42]. Previously, Elde et al. have studied the interaction of 10 primate PKR variants with VACV K3 and showed that PKR from Old World monkeys and New World monkeys were generally susceptible to VACV K3, whereas PKR from hominoids were resistant to VACV K3, indicating a species-specific inhibition of primate PKR by VACV K3 [31]. We recently described that K3 orthologs from myxoma virus, which infects rabbits and hares, and capripoxviruses, which infect sheep, goats, and cattle, also inhibit host PKRs in a species-specific manner [25,26,43]. Our recent studies revealed that K3 inhibition could not be predicted by the phylogenetic relatedness of PKR alone [28]. Therefore, in this study, we experimentally determined whether K3 orthologs of the Yatapoxvirus genus also exhibit host-specific PKR inhibition.
Yatapoxviruses infect primates and can cause zoonotic infections in humans. Here, we analyzed whether TPV and YMTV K3 orthologs inhibit PKR from 15 different primate species in comparison to VACV K3. Our results confirm and extend the results of species-specific VACV K3 inhibition previously reported in both yeast-based [31] and luciferase-based assays [28]. White-cheeked gibbon PKR was largely resistant to VACV K3, whereas VACV K3 strongly inhibited tamarin PKR and rhesus macaque PKR. Interestingly, dusky titi PKR, which was previously shown to be susceptible to VACV K3 in the yeast-based assay, appeared to be largely resistant to VACV K3 in our luciferase-based assays. Although the yeast assay can be used to predict the inhibition of PKR by viral antagonists, yeast expression systems have very different post-translational modifications relative to mammalian cells, which might affect protein function [44]. In general, our data showed that TPV 012 and YMTV 012 exhibited distinct PKR inhibition profiles and that these K3 orthologs inhibited primate PKR in a species-specific manner. All tested primate PKRs were sensitive to TPV 012, although to different degrees. We found a distinct profile of inhibition by YMTV 012 within the Old World monkeys in which talapoin monkey, African green monkey, and patas monkey PKRs were largely resistant to YMTV 012 inhibition. In contrast, rhesus macaque PKR was largely sensitive, supporting that the interaction between K3 orthologs and PKR could only be partially inferred by phylogenetic relatedness. This finding is in line with subcutaneous inoculation studies, showing that YMTV induced tumor formation in rhesus monkeys, whereas no tumors were observed in sooty mangabey, patas monkey, and African green monkeys [6,9,12].
TPV 012 and YMTV 012 differ by 22 amino acids out of 88 residues (Figure 1). K3 orthologs from the orthopoxvirus genus share a higher sequence identity (80.7–100%) among their members [28] than TPV 012 and YMTV 012. We found that residues within the C-terminus region are important for the differential PKR inhibition by these orthologs. There are two critical conserved motifs in the C-terminus of the K3 orthologs: the PKR recognition motif and helix insert region [19]. The PKR recognition motif, which includes residues K74, Y76, and D78, provides important determinants for high-affinity binding to PKR [19]. TPV 012 and YMTV 012 orthologs have identical residues in the PKR recognition motif “KGYVD”, one amino acid difference compared to the PKR recognition motif “KGYID” of VACV K3. The helix insert of eIF2α appears to undergo conformational change upon binding to PKR, resulting in the Ser51 site of eIF2α becoming fully accessible to the phosphoacceptor binding site of PKR [45]. The K3 helix insert region is structurally homologous to the region adjacent to Ser51 in eIF2α. The unique conformation of the helix insert region of K3L functions as an inhibitor of PKR activation by preventing PKR trans-autophosphorylation and acting as a pseudosubstrate of PKR [45]. In addition, we and others have shown that one or more amino acid residues in orthopoxvirus K3 orthologs are responsible for the differential PKR inhibition and that the amino acid variants are predominantly located in the helix insert region [27,28].
In this study, the function of yatapoxvirus K3 proteins was further examined in the context of virus infection. The E3L and K3L double-knockout vaccinia virus strain VC-R4 (VACVΔE3LΔK3L) can only replicate in PKR-deficient cell lines or cell lines that express PKR antagonists. By inserting K3L orthologs from different poxviruses in the E3L locus of VC-R4, we can analyze the chimeric viruses’ ability to replicate in PKR-competent cell lines. The advantage of using the chimeric virus system is that it allows us to assess the ability of K3L orthologs in rescuing the virus in PKR-competent cells. As expected, VC-R4 could not replicate in all tested cells, except RK13+E3L+K3L cells. VC-R4 formed smaller plaques than vP872, which can be explained by incomplete PKR inhibition by the ectopically expressed E3 and K3, and is consistent with the slightly less efficient replication of VC-R4 than vP872 in this cell line [46]. VC-R4 expressing SPPV 011 and TPV 012 replicated as efficiently as vP872 in all tested cell lines, as shown in plaque assays and virus titer. In contrast, VC-R4+VACV K3L and VC-R4+YMTV 012 formed smaller plaques in all tested cell lines, except in RK13+E3L+K3L where these viruses formed comparable plaque size with other viruses. Together these results indicate that VACV K3L and YMTV 012 are generally weaker inhibitors for the PKR orthologs that we tested. This result was recapitulated in the virus replication assay. There was an 11-fold higher virus titer in VC-R4+TPV 012-infected gibbon Tert-GF cells than in VC-R4+YMTV 012-infected cells. The difference in virus titer was more pronounced in PRO1190 and BSC-40 cells, with about 100- to 1000-fold lower virus titer in VC-R4+YMTV 012 infections. Although both cell lines are derived from African green monkeys, BSC-40 are immortalized epithelial cells and PRO1190 are primary fibroblasts. Differential susceptibility to various VACV recombinants by these two cell lines is well established, and may in part be due to three non-synonymous variants in PRO1190 PKR relative to BSC-40 PKR, although other host proteins, e.g., RNase L, have been implicated in some studies, suggesting that this phenotype is likely multifactorial [47,48,49]. While the general patterns observed in the luciferase-based assay were recapitulated in the infection assay, we observed some differences in the magnitude of inhibition.
Overall, the data presented here support previous observations regarding species-specific PKR inhibition by VACV K3 orthologs and extend them to a new genus of poxviruses [27,28,31]. Within primate species PKR, TPV 012 exhibited a broader range of PKR inhibition than did YMTV 012 and VACV K3. VACV K3 appeared to have a narrow range of PKR inhibition within primate species but efficiently inhibited PKR from a wide range mammals PKR [28]. It is worth mentioning that the ability of a virus to replicate and disseminate efficiently in its host is dependent on the entire interactome between the host and the virus [50]. In addition to K3L, yatapoxviruses encode several host range genes, including E3L, C7L, M11L, Serpins, P28-like, and B5R [51]. The role of the E3L gene of YMTV and C7L of YLDV in host range factors has been previously described [52,53]. However, the role of other genes as host range factors needs to be experimentally confirmed.

5. Conclusions

The possible reemergence of smallpox-like diseases through the zoonosis of animal poxviruses is a major public health concern. Members of the Yatapoxvirus genus infect humans and other primates and have caused outbreaks throughout Africa for the past 60 years. In addition, human infections with yatapoxviruses have been reported in animal handlers at primate centers in the United States and among travelers who visited Africa. Despite the threat posed to human health, the factors determining the host range of these viruses remain poorly understood. Here, we describe the inhibitory profiles of TPV and YMTV K3 orthologs against primate PKRs, and show that they display virus and host species-specific PKR inhibition. The TPV and YMTV K3 orthologs exhibited distinct PKR inhibition profiles, with TPV 012 generally inhibiting PKRs more strongly than YMTV 012. The species-specific PKR inhibition by the TPV and YMTV K3 orthologs may contribute to the differential susceptibility of primate species to TPV and YMTV infections.

Author Contributions

Conceptualization, D.M. and S.R.; methodology, D.M., J.N.S., C.P., T.C., L.T., G.B. and S.R.; validation, D.M., J.N.S., C.P., T.C., L.T., G.B. and S.R.; formal analysis, D.M., J.N.S., C.P., T.C., L.T., G.B. and S.R.; resources, S.R.; data curation, D.M. and L.T.; writing—original draft preparation, D.M.; writing—review and editing, D.M., J.N.S., C.P., L.T., G.B. and S.R.; visualization, D.M., L.T. and S.R.; supervision, S.R.; project administration, S.R.; funding acquisition, S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Grant R01 AI114851 (to S.R.) from the National Institute of Allergy and Infectious Diseases.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data necessary for the interpretation of the findings presented in this work are contained within the manuscript figures.

Acknowledgments

We are grateful to Jeffrey Hastings for technical assistance. We thank Nels Elde, Bertram Jacobs, Charles Samuel, and Adam Geballe for generously providing materials and reagents. Parts of this work were included in the PhD dissertation of D.M.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Haller, S.L.; Peng, C.; McFadden, G.; Rothenburg, S. Poxviruses and the evolution of host range and virulence. Infect. Genet. Evol. 2014, 21, 15–40. [Google Scholar] [CrossRef] [PubMed]
  2. Lee, H.J.; Essani, K.; Smith, G.L. The genome sequence of Yaba-like disease virus, a yatapoxvirus. Virology 2001, 281, 170–192. [Google Scholar] [CrossRef]
  3. Brunetti, C.R.; Amano, H.; Ueda, Y.; Qin, J.; Miyamura, T.; Suzuki, T.; Li, X.; Barrett, J.W.; McFadden, G. Complete genomic sequence and comparative analysis of the tumorigenic poxvirus Yaba monkey tumor virus. J. Virol. 2003, 77, 13335–13347. [Google Scholar] [CrossRef]
  4. Downie, A.W.; Taylor-Robinson, C.H.; Caunt, A.E.; Nelson, G.S.; Manson-Bahr, P.E.; Matthews, T.C. Tanapox: A new disease caused by a pox virus. Br. Med. J. 1971, 1, 363–368. [Google Scholar] [CrossRef]
  5. Jezek, Z.; Arita, I.; Szczeniowski, M.; Paluku, K.M.; Ruti, K.; Nakano, J.H. Human tanapox in Zaire: Clinical and epidemiological observations on cases confirmed by laboratory studies. Bull. World Health Organ. 1985, 63, 1027–1035. [Google Scholar]
  6. Bearcroft, W.G.; Jamieson, M.F. An outbreak of subcutaneous tumours in rhesus monkeys. Nature 1958, 182, 195–196. [Google Scholar] [CrossRef]
  7. Knight, J.C.; Novembre, F.J.; Brown, D.R.; Goldsmith, C.S.; Esposito, J.J. Studies on Tanapox virus. Virology 1989, 172, 116–124. [Google Scholar] [CrossRef]
  8. Nazarian, S.H.; Barrett, J.W.; Frace, A.M.; Olsen-Rasmussen, M.; Khristova, M.; Shaban, M.; Neering, S.; Li, Y.; Damon, I.K.; Esposito, J.J.; et al. Comparative genetic analysis of genomic DNA sequences of two human isolates of Tanapox virus. Virus Res. 2007, 129, 11–25. [Google Scholar] [CrossRef] [PubMed]
  9. Sproul, E.E.; Metzgar, R.S.; Grace, J.T., Jr. The pathogenesis of Yaba virus-induced histiocytomas in primates. Cancer Res. 1963, 23, 671–675. [Google Scholar]
  10. Downie, A.W. Serological evidence of infection with Tana and Yaba pox viruses among several species of monkey. J. Hyg. 1974, 72, 245–250. [Google Scholar] [CrossRef]
  11. Monroe, B.P.; Nakazawa, Y.J.; Reynolds, M.G.; Carroll, D.S. Estimating the geographic distribution of human Tanapox and potential reservoirs using ecological niche modeling. Int. J. Health Geogr. 2014, 13, 34. [Google Scholar] [CrossRef] [PubMed]
  12. Ambrus, J.L.; Strandström, H.V. Susceptibility of Old World monkeys to Yaba virus. Nature 1966, 211, 876. [Google Scholar] [CrossRef]
  13. Grace, J.T., Jr.; Mirand, E.A. Human susceptibility to a simian tumor virus. Ann. N. Y. Acad. Sci. 1963, 108, 1123–1128. [Google Scholar] [CrossRef] [PubMed]
  14. Stich, A.; Meyer, H.; Köhler, B.; Fleischer, K. Tanapox: First report in a European traveller and identification by PCR. Trans. R. Soc. Trop. Med. Hyg. 2002, 96, 178–179. [Google Scholar] [CrossRef] [PubMed]
  15. Dimitrov, D.S. Virus entry: Molecular mechanisms and biomedical applications. Nat. Reviews. Microbiol. 2004, 2, 109–122. [Google Scholar] [CrossRef]
  16. McFadden, G. Poxvirus tropism. Nat. Reviews. Microbiol. 2005, 3, 201–213. [Google Scholar] [CrossRef]
  17. Moss, B. Poxvirus cell entry: How many proteins does it take? Viruses 2012, 4, 688–707. [Google Scholar] [CrossRef]
  18. Jacobs, B.L.; Langland, J.O. When two strands are better than one: The mediators and modulators of the cellular responses to double-stranded RNA. Virology 1996, 219, 339–349. [Google Scholar] [CrossRef]
  19. Dar, A.C.; Sicheri, F. X-ray crystal structure and functional analysis of vaccinia virus K3L reveals molecular determinants for PKR subversion and substrate recognition. Mol. Cell 2002, 10, 295–305. [Google Scholar] [CrossRef]
  20. Adomavicius, T.; Guaita, M.; Zhou, Y.; Jennings, M.D.; Latif, Z.; Roseman, A.M.; Pavitt, G.D. The structural basis of translational control by eIF2 phosphorylation. Nat. Commun. 2019, 10, 2136. [Google Scholar] [CrossRef]
  21. Romano, P.R.; Zhang, F.; Tan, S.L.; Garcia-Barrio, M.T.; Katze, M.G.; Dever, T.E.; Hinnebusch, A.G. Inhibition of double-stranded RNA-dependent protein kinase PKR by vaccinia virus E3: Role of complex formation and the E3 N-terminal domain. Mol. Cell. Biol. 1998, 18, 7304–7316. [Google Scholar] [CrossRef] [PubMed]
  22. Langland, J.O.; Jacobs, B.L. The role of the PKR-inhibitory genes, E3L and K3L, in determining vaccinia virus host range. Virology 2002, 299, 133–141. [Google Scholar] [CrossRef] [PubMed]
  23. Langland, J.O.; Jacobs, B.L. Inhibition of PKR by vaccinia virus: Role of the N- and C-terminal domains of E3L. Virology 2004, 324, 419–429. [Google Scholar] [CrossRef] [PubMed]
  24. Haller, S.L.; Park, C.; Bruneau, R.C.; Megawati, D.; Zhang, C.; Vipat, S.; Peng, C.; Senkevich, T.G.; Brennan, G.; Tazi, L.; et al. Molecular basis for the host range function of the poxvirus PKR inhibitor E3. bioRxiv 2024. [Google Scholar] [CrossRef]
  25. Peng, C.; Haller, S.L.; Rahman, M.M.; McFadden, G.; Rothenburg, S. Myxoma virus M156 is a specific inhibitor of rabbit PKR but contains a loss-of-function mutation in Australian virus isolates. Proc. Natl. Acad. Sci. USA 2016, 113, 3855–3860. [Google Scholar] [CrossRef] [PubMed]
  26. Park, C.; Peng, C.; Brennan, G.; Rothenburg, S. Species-specific inhibition of antiviral protein kinase R by capripoxviruses and vaccinia virus. Ann. N.Y. Acad. Sci. 2019, 1438, 18–29. [Google Scholar] [CrossRef] [PubMed]
  27. Cao, J.; Varga, J.; Deschambault, Y. Poxvirus encoded eIF2α homolog, K3 family proteins, is a key determinant of poxvirus host species specificity. Virology 2020, 541, 101–112. [Google Scholar] [CrossRef] [PubMed]
  28. Park, C.; Peng, C.; Rahman, M.J.; Haller, S.L.; Tazi, L.; Brennan, G.; Rothenburg, S. Orthopoxvirus K3 orthologs show virus- and host-specific inhibition of the antiviral protein kinase PKR. PLoS Pathog. 2021, 17, e1009183. [Google Scholar] [CrossRef] [PubMed]
  29. Carpentier, K.S.; Esparo, N.M.; Child, S.J.; Geballe, A.P. A Single Amino Acid Dictates Protein Kinase R Susceptibility to Unrelated Viral Antagonists. PLoS Pathog. 2016, 12, e1005966. [Google Scholar] [CrossRef]
  30. Rahman, M.M.; Liu, J.; Chan, W.M.; Rothenburg, S.; McFadden, G. Myxoma virus protein M029 is a dual function immunomodulator that inhibits PKR and also conscripts RHA/DHX9 to promote expanded host tropism and viral replication. PLoS Pathog. 2013, 9, e1003465. [Google Scholar] [CrossRef]
  31. Elde, N.C.; Child, S.J.; Geballe, A.P.; Malik, H.S. Protein kinase R reveals an evolutionary model for defeating viral mimicry. Nature 2009, 457, 485–489. [Google Scholar] [CrossRef] [PubMed]
  32. Gibson, D.G.; Young, L.; Chuang, R.Y.; Venter, J.C.; Hutchison, C.A., 3rd; Smith, H.O. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 2009, 6, 343–345. [Google Scholar] [CrossRef] [PubMed]
  33. Vipat, S.; Brennan, G.; Park, C.; Haller, S.L.; Rothenburg, S. Rapid, Seamless Generation of Recombinant Poxviruses using Host Range and Visual Selection. J. Vis. Exp. JoVE 2020, 159, e61049. [Google Scholar] [CrossRef]
  34. Rothenburg, S.; Seo, E.J.; Gibbs, J.S.; Dever, T.E.; Dittmar, K. Rapid evolution of protein kinase PKR alters sensitivity to viral inhibitors. Nat. Struct. Mol. Biol. 2009, 16, 63–70. [Google Scholar] [CrossRef] [PubMed]
  35. Beattie, E.; Tartaglia, J.; Paoletti, E. Vaccinia virus-encoded eIF-2 alpha homolog abrogates the antiviral effect of interferon. Virology 1991, 183, 419–422. [Google Scholar] [CrossRef] [PubMed]
  36. Esposito, J.; Condit, R.; Obijeski, J. The preparation of orthopoxvirus DNA. J. Virol. Methods 1981, 2, 175–179. [Google Scholar] [CrossRef]
  37. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  38. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed]
  39. Sievers, F.; Higgins, D.G. Clustal Omega for making accurate alignments of many protein sequences. Protein Sci. 2018, 27, 135–145. [Google Scholar] [CrossRef]
  40. Nejepinska, J.; Malik, R.; Wagner, S.; Svoboda, P. Reporters transiently transfected into mammalian cells are highly sensitive to translational repression induced by dsRNA expression. PLoS ONE 2014, 9, e87517. [Google Scholar] [CrossRef]
  41. Kawagishi-Kobayashi, M.; Silverman, J.B.; Ung, T.L.; Dever, T.E. Regulation of the protein kinase PKR by the vaccinia virus pseudosubstrate inhibitor K3L is dependent on residues conserved between the K3L protein and the PKR substrate eIF2alpha. Mol. Cell. Biol. 1997, 17, 4146–4158. [Google Scholar] [CrossRef] [PubMed]
  42. Davies, M.V.; Chang, H.W.; Jacobs, B.L.; Kaufman, R.J. The E3L and K3L vaccinia virus gene products stimulate translation through inhibition of the double-stranded RNA-dependent protein kinase by different mechanisms. J. Virol. 1993, 67, 1688–1692. [Google Scholar] [CrossRef] [PubMed]
  43. Yu, H.; Peng, C.; Zhang, C.; Stoian, A.M.M.; Tazi, L.; Brennan, G.; Rothenburg, S. Maladaptation after a virus host switch leads to increased activation of the pro-inflammatory NF-κB pathway. Proc. Natl. Acad. Sci. USA 2022, 119, e2115354119. [Google Scholar] [CrossRef] [PubMed]
  44. Vieira Gomes, A.M.; Souza Carmo, T.; Silva Carvalho, L.; Mendonça Bahia, F.; Parachin, N.S. Comparison of Yeasts as Hosts for Recombinant Protein Production. Microorganisms 2018, 6, 38. [Google Scholar] [CrossRef] [PubMed]
  45. Dar, A.C.; Dever, T.E.; Sicheri, F. Higher-order substrate recognition of eIF2alpha by the RNA-dependent protein kinase PKR. Cell 2005, 122, 887–900. [Google Scholar] [CrossRef]
  46. Hand, E.S.; Haller, S.L.; Peng, C.; Rothenburg, S.; Hersperger, A.R. Ectopic expression of vaccinia virus E3 and K3 cannot rescue ectromelia virus replication in rabbit RK13 cells. PLoS ONE 2015, 10, e0119189. [Google Scholar] [CrossRef] [PubMed]
  47. Brennan, G.; Kitzman, J.O.; Rothenburg, S.; Shendure, J.; Geballe, A.P. Adaptive gene amplification as an intermediate step in the expansion of virus host range. PLoS Pathog. 2014, 10, e1004002. [Google Scholar] [CrossRef] [PubMed]
  48. Brennan, G.; Kitzman, J.O.; Shendure, J.; Geballe, A.P. Experimental Evolution Identifies Vaccinia Virus Mutations in A24R and A35R That Antagonize the Protein Kinase R Pathway and Accompany Collapse of an Extragenic Gene Amplification. J. Virol. 2015, 89, 9986–9997. [Google Scholar] [CrossRef] [PubMed]
  49. Banerjee, S.; Smith, C.; Geballe, A.P.; Rothenburg, S.; Kitzman, J.O.; Brennan, G. Gene amplification acts as a molecular foothold to facilitate cross-species adaptation and evasion of multiple antiviral pathways. Virus Evol. 2022, 8, veac105. [Google Scholar] [CrossRef]
  50. Rothenburg, S.; Brennan, G. Species-Specific Host-Virus Interactions: Implications for Viral Host Range and Virulence. Trends Microbiol. 2020, 28, 46–56. [Google Scholar] [CrossRef]
  51. Bratke, K.A.; McLysaght, A.; Rothenburg, S. A survey of host range genes in poxvirus genomes. Infect. Genet. Evol. 2013, 14, 406–425. [Google Scholar] [CrossRef] [PubMed]
  52. Myskiw, C.; Arsenio, J.; Hammett, C.; van Bruggen, R.; Deschambault, Y.; Beausoleil, N.; Babiuk, S.; Cao, J. Comparative analysis of poxvirus orthologues of the vaccinia virus E3 protein: Modulation of protein kinase R activity, cytokine responses, and virus pathogenicity. J. Virol. 2011, 85, 12280–12291. [Google Scholar] [CrossRef] [PubMed]
  53. Meng, X.; Chao, J.; Xiang, Y. Identification from diverse mammalian poxviruses of host-range regulatory genes functioning equivalently to vaccinia virus C7L. Virology 2008, 372, 372–383. [Google Scholar] [CrossRef] [PubMed]
Viruses 16 01095 g001
Figure 1. Amino acid differences among TPV 012, YMTV 012, and VACV K3. (A). Residues differing from TPV 012 are highlighted in red. The regions to the left and right of the dotted line were considered the N-terminal and C-terminal regions of the inhibitors, respectively. Sequences used to generate antibodies against TPV 012 and VACV K3 are boxed. (B). Percent identities between tested K3 orthologs were calculated from the multiple sequence alignment using Clustal Omega [39]. (C,D). HeLa PKRko cells were transfected with 1 µg of either pSG5 empty vector, pSG5 encoding FLAG-tag VACV K3, FLAG-tag TPV 012, or FLAG-tag YMTV 012. Proteins were harvested at 48 h post transfection, and protein expressions were detected with anti-Flag M2, anti-TPV 012, anti-VACV K3, and anti-β actin antibodies as indicated.
Figure 1. Amino acid differences among TPV 012, YMTV 012, and VACV K3. (A). Residues differing from TPV 012 are highlighted in red. The regions to the left and right of the dotted line were considered the N-terminal and C-terminal regions of the inhibitors, respectively. Sequences used to generate antibodies against TPV 012 and VACV K3 are boxed. (B). Percent identities between tested K3 orthologs were calculated from the multiple sequence alignment using Clustal Omega [39]. (C,D). HeLa PKRko cells were transfected with 1 µg of either pSG5 empty vector, pSG5 encoding FLAG-tag VACV K3, FLAG-tag TPV 012, or FLAG-tag YMTV 012. Proteins were harvested at 48 h post transfection, and protein expressions were detected with anti-Flag M2, anti-TPV 012, anti-VACV K3, and anti-β actin antibodies as indicated.
Viruses 16 01095 g001
Viruses 16 01095 g002
Figure 2. Differential sensitivities of tested PKRs to VACV K3, TPV 012, and YMTV 012 orthologs. (A). HeLa PKRko cells were co-transfected with expression vectors encoding firefly luciferase (50 ng) and with either empty vector or PKR (200 ng) from the indicated species. Luciferase activities were measured 48 h after transfection and normalized to empty vector transfected cells to obtain relative luciferase activities. (B). HeLa PKRko cells were transfected with expression vectors encoding firefly luciferase (50 ng), PKR (200 ng), and PKR inhibitor (200 ng) from the indicated species. Luciferase activities were measured 48 h after transfection and normalized to PKR-only transfected cells to obtain relative luciferase activities. Error bars represent the standard deviations from three independent transfections. The luciferase assay results shown are representative of at least three independent experiments.
Figure 2. Differential sensitivities of tested PKRs to VACV K3, TPV 012, and YMTV 012 orthologs. (A). HeLa PKRko cells were co-transfected with expression vectors encoding firefly luciferase (50 ng) and with either empty vector or PKR (200 ng) from the indicated species. Luciferase activities were measured 48 h after transfection and normalized to empty vector transfected cells to obtain relative luciferase activities. (B). HeLa PKRko cells were transfected with expression vectors encoding firefly luciferase (50 ng), PKR (200 ng), and PKR inhibitor (200 ng) from the indicated species. Luciferase activities were measured 48 h after transfection and normalized to PKR-only transfected cells to obtain relative luciferase activities. Error bars represent the standard deviations from three independent transfections. The luciferase assay results shown are representative of at least three independent experiments.
Viruses 16 01095 g002
Viruses 16 01095 g003
Figure 3. Sensitivities of tested PKRs projected on a PKR phylogenetic tree. We generated a phylogenetic tree with PKR protein sequences from a panel of 21 primate species and used rabbit PKR as an outgroup. The tree was built by PhyML with 100 bootstrap replicates and based on the maximum-likelihood method. Relative sensitivities of PKRs to either VACV K3, TPV 012, or YMTV 012 from multiple luciferase-based reporter assays are shown on a scale from 1 to 6. Scale 1 = no or very weak inhibition (≤2-fold increase in luciferase readout); 2 = weak inhibition (2- to 3-fold); 3 = intermediate low inhibition (3- to 5-fold); 4 = intermediate high inhibition (5- to 7-fold); 5 = high inhibition (7- to 10-fold); 6=very high inhibition (≥10-fold).
Figure 3. Sensitivities of tested PKRs projected on a PKR phylogenetic tree. We generated a phylogenetic tree with PKR protein sequences from a panel of 21 primate species and used rabbit PKR as an outgroup. The tree was built by PhyML with 100 bootstrap replicates and based on the maximum-likelihood method. Relative sensitivities of PKRs to either VACV K3, TPV 012, or YMTV 012 from multiple luciferase-based reporter assays are shown on a scale from 1 to 6. Scale 1 = no or very weak inhibition (≤2-fold increase in luciferase readout); 2 = weak inhibition (2- to 3-fold); 3 = intermediate low inhibition (3- to 5-fold); 4 = intermediate high inhibition (5- to 7-fold); 5 = high inhibition (7- to 10-fold); 6=very high inhibition (≥10-fold).
Viruses 16 01095 g003
Viruses 16 01095 g004
Figure 4. Identification of region of TPV 012 and YMTV 012 orthologs that confer differential PKR inhibition. (A). The predicted amino acid sequences of TPV 012 and YMTV 012 hybrid proteins. Constructs were designed such that the C-terminal regions of the two inhibitors were swapped to generate hybrid proteins. Different residues among K3 orthologs are highlighted in red. (B). Human HeLa-PKRko cells were transfected with expression vector encoding firefly luciferase, indicated PKRs, and either hybrid NTPV-CYMTV or NYMTV-CTPV inhibitors. The inhibition of indicated PKRs by TPV 012 and YMTV 012 serves as control. Luciferase activity was measured 48 h after transfection. Luciferase light units were normalized to PKR-only transfected cells to obtain relative luciferase activities. Error bars represent triplicate transfections and data are representative of three independent experiments. The data were analyzed using the unpaired t-test (GraphPad). p-values are indicated as follows: ns (non-significant) p > 0.05, * p < 0.05, ** p < 0.01.
Figure 4. Identification of region of TPV 012 and YMTV 012 orthologs that confer differential PKR inhibition. (A). The predicted amino acid sequences of TPV 012 and YMTV 012 hybrid proteins. Constructs were designed such that the C-terminal regions of the two inhibitors were swapped to generate hybrid proteins. Different residues among K3 orthologs are highlighted in red. (B). Human HeLa-PKRko cells were transfected with expression vector encoding firefly luciferase, indicated PKRs, and either hybrid NTPV-CYMTV or NYMTV-CTPV inhibitors. The inhibition of indicated PKRs by TPV 012 and YMTV 012 serves as control. Luciferase activity was measured 48 h after transfection. Luciferase light units were normalized to PKR-only transfected cells to obtain relative luciferase activities. Error bars represent triplicate transfections and data are representative of three independent experiments. The data were analyzed using the unpaired t-test (GraphPad). p-values are indicated as follows: ns (non-significant) p > 0.05, * p < 0.05, ** p < 0.01.
Viruses 16 01095 g004
Viruses 16 01095 g005
Figure 5. Generation of chimeric VC-R4 expressing yatapoxvirus K3 orthologs. (A). HeLa PKRko were infected with VC-R4, VC-R4+VACV K3, VC-R4+SPPV 011, VC-R4+TPV 012, or VC-R4+YMTV 012. Protein lysates were collected at 24 h post infection. The expression of K3 orthologs was determined by immunoblotting with the indicated antibodies. (B). Plaque formation of VC-R4 expressing K3 orthologs. A variety of primate-derived cells and RK13+E3L+K3L were infected with 50 pfu/mL of the indicated viruses, and plaque formation was visualized with crystal violet at 48 h after infection.
Figure 5. Generation of chimeric VC-R4 expressing yatapoxvirus K3 orthologs. (A). HeLa PKRko were infected with VC-R4, VC-R4+VACV K3, VC-R4+SPPV 011, VC-R4+TPV 012, or VC-R4+YMTV 012. Protein lysates were collected at 24 h post infection. The expression of K3 orthologs was determined by immunoblotting with the indicated antibodies. (B). Plaque formation of VC-R4 expressing K3 orthologs. A variety of primate-derived cells and RK13+E3L+K3L were infected with 50 pfu/mL of the indicated viruses, and plaque formation was visualized with crystal violet at 48 h after infection.
Viruses 16 01095 g005
Viruses 16 01095 g006
Figure 6. Virus titer of chimeric vaccinia viruses in primate-derived cell lines. A variety of primate-derived cells and RK13+E3L+K3L were infected with indicated viruses at a MOI of 0.01. Cell lysates were harvested at 30 h post infection and virus titers were determined in RK13+E3L+K3L cells. Error bars represent the standard deviations from three independent infections. Fold changes in virus titer between VC-R4+TPV 012 and VC-R4+YMTV 012 are shown. Data were analyzed using one-way ANOVA followed by Tukey’s multiple comparisons test. p-values are marked as follows: ns = p > 0.05, **** p < 0.0001.
Figure 6. Virus titer of chimeric vaccinia viruses in primate-derived cell lines. A variety of primate-derived cells and RK13+E3L+K3L were infected with indicated viruses at a MOI of 0.01. Cell lysates were harvested at 30 h post infection and virus titers were determined in RK13+E3L+K3L cells. Error bars represent the standard deviations from three independent infections. Fold changes in virus titer between VC-R4+TPV 012 and VC-R4+YMTV 012 are shown. Data were analyzed using one-way ANOVA followed by Tukey’s multiple comparisons test. p-values are marked as follows: ns = p > 0.05, **** p < 0.0001.
Viruses 16 01095 g006
Table 1. Accession numbers of genes used for the phylogenetic analysis.
Table 1. Accession numbers of genes used for the phylogenetic analysis.
Common NameScientific NameGenBank Accession
HumanHomo sapiensNM_002759
GorillaGorilla gorillaEU733258.1
BonoboPan paniscusEU733255.1
ChimpanzeePan troglodytesEU733256.1
OrangutanPongo pygmaeus pgmaeusEU733259.1
White-cheeked gibbonNomascus leucogenysEU733257.1
White-bearded gibbonHylobates albibarbisEU733270.1
SiamangSymphalangus syndactylusEU733271.1
Colobus monkeyColobus guerezaEU733267.1
François’ leaf monkeyTrachypithecus francoisiEU733268.1 *
BaboonPapio anubisXM_009184002.4 *
Sooty mangabeyCercocebus torquatus atysEU733262.1
Rhesus macaque Macaca mulattaEU733261.1
Talapoin monkeyMiopithecus talapoin talapoinEU733269.1
African green monkeyClorocebus aethiopsEU733254.1
Patas monkeyErythrocebus patasEU733260.1
Black-handed spider monkeyAteles geoffroyiEU733263.1
Woolly monkeyLagothrix lagotrichaEU733266.1
TamarinSaguinus labiatusEU733264.1
Dusky titiCallicebus molochEU733265.1
Owl monkeyAotus trivirgatusFJ374685.1
RabbitOryctolagus cuniculusNM_001082213.1
* Nucleotide differences between the PKR sequences we used and the GenBank sequences: François leaf PKR and EU733268.1: all synonymous: T208 > C (L70L), C861 > T (I287I), G1326 > A (G442G), T1446C (H482H). Baboon PKR and XM_009184002.4: one synonymous A465G (Q155Q) and four non-synonymous G1489A (E497K), A1498C (K500Q), G1511A (G504T), and A1558G (K520E).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Megawati, D.; Stroup, J.N.; Park, C.; Clarkson, T.; Tazi, L.; Brennan, G.; Rothenburg, S. Tanapox Virus and Yaba Monkey Tumor Virus K3 Orthologs Inhibit Primate Protein Kinase R in a Species-Specific Fashion. Viruses 2024, 16, 1095. https://doi.org/10.3390/v16071095

AMA Style

Megawati D, Stroup JN, Park C, Clarkson T, Tazi L, Brennan G, Rothenburg S. Tanapox Virus and Yaba Monkey Tumor Virus K3 Orthologs Inhibit Primate Protein Kinase R in a Species-Specific Fashion. Viruses. 2024; 16(7):1095. https://doi.org/10.3390/v16071095

Chicago/Turabian Style

Megawati, Dewi, Jeannine N. Stroup, Chorong Park, Taylor Clarkson, Loubna Tazi, Greg Brennan, and Stefan Rothenburg. 2024. "Tanapox Virus and Yaba Monkey Tumor Virus K3 Orthologs Inhibit Primate Protein Kinase R in a Species-Specific Fashion" Viruses 16, no. 7: 1095. https://doi.org/10.3390/v16071095

APA Style

Megawati, D., Stroup, J. N., Park, C., Clarkson, T., Tazi, L., Brennan, G., & Rothenburg, S. (2024). Tanapox Virus and Yaba Monkey Tumor Virus K3 Orthologs Inhibit Primate Protein Kinase R in a Species-Specific Fashion. Viruses, 16(7), 1095. https://doi.org/10.3390/v16071095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop