Amino Acid Substitutions in NS5 Contribute Differentially to Tembusu Virus Attenuation in Ducklings and Cell Cultures

Tembusu virus (TMUV), a highly infectious pathogenic flavivirus, causes severe egg-drop and encephalitis in domestic waterfowl, while the determinants responsible for viral pathogenicity are largely unknown. In our previous studies, virulent strain JXSP2-4 had been completely attenuated by successive passages in BHK-21 cells and the avirulent strain was designated as JXSP-310. Based on the backbone of JXSP2-4, a series of chimeric viruses were generated according to the amino acid substitutions in NS5 and their infectivities were also analyzed in cell cultures and ducklings. The results showed that the viral titers of RNA-dependent RNA polymerase (RdRp) domain-swapped cheimeric mutant (JXSP-310RdRp) in cells and ducklings were both markedly decreased compared with JXSP2-4, indicating that mutations in the RdRp domain affected viral replication. There are R543K and V711A two amino acid substitutions in the RdRp domain. Further site-directed mutagenesis showed that single-point R543K mutant (JXSP-R543K) exhibited similar replication efficacy compared with JXSP2-4 in cells, but the viral loads in JXSP-R543K-infected ducklings were significantly lower than that of JXSP2-4 and higher than JXSP-310RdRp. Surprisingly, the single-point V711A mutation we introduced rapidly reverted. In addition, qRT-PCR and Western blot confirmed that the mutations in the RdRp domain significantly affected the replication of the virus. Taken together, these results show that R543K substitution in the RdRp domain impairs the in vivo growth of TMUV, but sustaining its attenuated infectivity requires the concurrent presence of the V711A mutation.


Introduction
Tembusu virus (TMUV) was first isolated from Culex tritaeniorhynchus mosquitoes in Malaysia in 1955 [1]. TMUV, a member of the Ntaya group in the Flaviviridae family, has a 10,990-nucleotide-long single-stranded, positive-sense RNA genome. The large, single precursor polyprotein encoded by the viral genome is processed into three structural proteins: the capsid (C), the premembrane/membrane protein (prM), the envelope (E) and seven nonstructural proteins: NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5 [1][2][3]. TMUV was associated with the explosive epizootics that occurred in commercial duck flocks in China in 2010 and in Southeast Asian countries later on [4][5][6]. The virus is transmitted among susceptible duck flocks by multiple routes (e.g., mosquito bites, infectious aerosol inhalation and eating contaminated feed) and results in severe egg drop and fatal encephalitis [7][8][9]. Because of the widespread duck farms and diversified duck housing systems in China and Southeastern Asia, vaccination is the preferred prophylactic countermeasure against highly contagious TMUV disease. Live-attenuated TMUV vaccines are expected to have several desirable attributes (e.g., stimulation of durable systemic immunity, production BHK-21 cells with Lipofectamine 3000 Transfection Reagent (Invitrogen, Carlsbad, C USA). The supernatants from the RNA-transfected cells were harvested at 72 h post-tra fection when typical cytopathic effects were observed. The supernatants were then u to infect fresh BHK-21 cells to amplify the rescued viruses. The sequences of the rescu viruses were confirmed by full-genome sequencing. Viral titers were determined standard plaque assays in BHK-21 cells. The amplified viral stocks were aliquoted a stored at −80 °C.   The culture supernatants were collected every 12 h and viral titers were quantified using plaque assay. All experiments were performed in triplicate and the mean ± the SD values are shown.

Immunofluorescence Assays
BHK-21 cells grown in 12-well plates were infected with the rescued viruses or were mock infected. At 48 h post-infection, the cells were fixed in cold (−20 • C) acetone and methanol (1:1) for 20 min at room temperature. The fixed cells were washed with PBS three times and then incubated with anti-E protein mAb (1:1000 dilution in PBS) at 37 • C for 1 h. After three washes with PBS, the cells were incubated with DyLight488-conjugated goat anti-mouse IgG (Canlifesci Inc., Co., Ltd., Beijing, China) (1:1000 dilution in PBS) at 37 • C for another 1 h. Finally, the plates were washed three times and examined by fluorescence microscopy.

Virus Replication Kinetics
To investigate replication capacity of the chimeric viruses, mutant viruses and their parental strain, BHK-21 cells and DEF cells were grown in 6-well plates to 80% confluence and then infected with each virus at an multiplicity of infection (MOI) of 0.1. After inoculation at 37 • C for 1 h, the supernatants were discarded and the cells were washed three times with PBS. Fresh DMEM (2 mL) and 2% FBS were added. The cell culture supernatants were collected at various time points post-infection. Progeny virus yields were quantitated by plaque assay in BHK-21 cells.

Plaque Assays
Serial 10-fold dilutions of each viral sample were prepared in DMEM (plus 1% FBS) and inoculated onto BHK-21 monolayers in 6-well plates. The plates were incubated at 37 • C for 1 h. After inoculum removal, the cells were washed twice with PBS and overlaid with DMEM containing 1% (wt/vol) low-melting point agarose (Amresco, Houston, TX, USA) and 2% FBS. After further incubation at 37 • C for 72 h, the cells were stained with 0.02% neutral red to visualize the plaques.

Duck Infection Experiments
One-day-old Pekin ducklings, purchased from Beijing Golden Star Duck Co., Ltd., (Beijing, China) were used to examine the virulence of the chimeric and mutant viruses. Before challenge, TMUV negativity was confirmed in the ducklings using a blocking enzyme-linked immunosorbent assay [20]. For the challenge groups, 7-day-old healthy Pekin ducklings were infected subcutaneously with 10 5 plaque forming units (PFUs) of the parental virus, the chimeric viruses or the mutant viruses. The control group was mock infected with 0.5 mL DMEM. Viral loads were detected in the tissues of five euthanized ducks from each group at various time points. The blood and tissue samples from the heart, liver, spleen, lungs, kidneys, thymus and bursa of Fabricius were collected for viral titer determination. Tissue samples were homogenized and diluted in DMEM to a final concentration of 20% (wt/vol). After three cycles of freeze-thawing, the tissue suspensions were centrifuged (5000× g at 4 • C for 20 min). Viral supernatant titers were determined using plaque assay in BHK-21 cells.

RNA Extraction and qRT-PCR Analysis
RNAs in the supernatants from the infected DEF cells and tissue samples were extracted using a SV Total RNA Isolation System (Promega, Fitchburg, WI, USA). cDNAs were synthesized using a reverse transcription system (Promega, Fitchburg, WI, USA) with Viruses 2021, 13, 921 5 of 14 random primers. Real-time PCR was conducted in a Bio-Rad CFX Connect real-time system using SYBR green SuperReal PreMix Plus (Tiangen, Beijing, China) with primers (forward, 5 -TCATTGATAGAATTTGAGGAG-3 ; reverse, 5 -TTCCAATTTGCTTCCAGAGTA-3 ) targeting the TMUV E region. A plasmid fragment containing the E region's genome sequence was used to prepare a standard curve for quantifying viral genome copy numbers. Relative qRT-PCR was performed to detect the transcription level of IFN-α, IFN-β, IFN-γ, IL-6 and IL-1β and the primers used in this study were listed in Table 2.

Western Blot Analysis
DEF cells in 6-well plates, which were separately infected with the parental virus and mutant viruses, were collected at various time points post-infection. Cells were lysed with Cell Lysis Buffer for Western and IP (IP: immunoprecipitation; Beyotime, Shanghai, China). Protein concentrations were determined by the BCA method (Solarbio, Beijing, China). Equal protein amounts were separated by 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis and then transferred from the gel to polyvinylidene difluoride membranes by wet transfer. The membranes were blocked with 5% skimmed milk in PBST (PBS/Tween) at 37 • C for 2 h. The membranes were incubated at 37 • C for 1 h with anti-E protein mAb or anti-β-actin mAb (1:5000 dilution in PBST). After three washes with PBST, the membranes were incubated with secondary antibodies for 1 h at 37 • C. Finally, the protein bands were detected using the BeyoECL Plus Kit after three washes (Beyotime, Shanghai, China).

Statistical Analysis
All data were processed with GraphPad Prism 8 (GraphPad Software Inc., La Jolla, CA, USA). The student's t-test or two-way multiple ANOVA comparisons test was used to analyze differences between the values from the two groups. A value of p < 0.05 was considered statistically significant.

Replacing the NS5 RdRp Domain Decreases the Growth of the Virus in Cell Cultures
The data from our previous study suggested that amino acid substitutions in NS5 contributed to virulence attenuation of TMUV JXSP 2-4 [13]. To identify the sub-regions within TMUV NS5 responsible for virulence attenuation, we first constructed three chimeric mutants by individually exchanging the MTase, Linker and RdRp domains in attenuated JXSP-310 with the corresponding sites in the genome of the JXSP 2-4 -derived infectious clone by reverse genetic manipulation. The overall cloning strategy is shown in Figure 1A. The recombinant viruses were detected by immunofluorescence staining with a monoclonal antibody against the TMUV E protein ( Figure 1B Linker and RdRp domains, respectively. We noted that the JXSP-310 RdRp chimeric virus produced larger plaque sizes in BHK-21 cell monolayers than its parental JXSP 2-4 virus (p < 0.001) ( Figure 1C), suggesting that the fitness of JXSP-310 RdRp might have been altered. We then compared the in vitro replication kinetics of the chimeric viruses with that of parental JXSP 2-4 using multi-step growth curves. BHK-21 and DEF cells were infected with either chimeric virus or the parental virus at an MOI of 0.1. The JXSP-310 RdRp chimeric virus grew significantly slower than the parental virus JXSP 2-4 in both cell types with 10to 25-fold lower infectious virus titers and its peak titer lagged by approximately 12 h in BHK-21 cells ( Figure 1D). In contrast, JXSP-310 MTase and JXSP-310 Linker exhibited growth kinetics comparable to those of the parental virus. These results suggested that mutations in the NS5 RdRp domain decreased the replication ability of TMUV in cultured cells. However, the increased plaque size of JXSP-310 RdRp was not correlated positively with virus multiplication in cell cultures.

The In Vivo Infectivity of Chimeric JXSP-310 RdRp Is Attenuated
To examine the effects of mutations in the aforementioned domains on viral virulence, we initially compared the in vivo infectivity of the chimeric viruses by measuring viral loads in the tissues of 7-day-old ducklings after subcutaneous infection with 10 5 PFU of virus per duckling (n = 5 per group). Plaque assay quantification showed that the infectious virus loads of JXSP-310 RdRp in representative tissues from the inoculated ducklings were 30to 100-fold lower than those of the parental virus on day 2 post-infection (pi) ( Figure 2). No significant difference was observed between the JXSP-310 MTase -and JXSP-310 Linker -infected ducklings. To further monitor the in vivo growth kinetics of the JXSP-310 RdRp chimeric virus, a group of 7-day-old ducklings (n = 20 per group) were subcutaneously infected with 10 5 PFU of JXSP-310 RdRp and another group of 20 with the JXSP 2-4 parental virus. Ducklings were sacrificed on days 1, 2, 4 and 6 pi and the viral loads in representative tissues were measured by plaque assays. The infectivity titers in the tested tissues of the JXSP-310 RdRp -infected ducklings were all significantly lower than those of JXSP 2-4 ( Figure 3). On day 1 pi the viremia and viral titers in the spleen peaked and these values in JXSP-310 RdRp -infected ducklings were 10-fold and 150-fold lower, respectively, than in the JXSP 2-4 -infected ducklings. On day 2 pi, the viral titers in the other organs peaked and JXSP-310 RdRp was significantly attenuated (p < 0.01). By day 4 pi, the viral titers in the heart tissues from the JXSP-310 RdRp -infected ducklings were 200-fold lower than in the JXSP 2-4 -infected ducklings (p < 0.001). These results indicate that replacing the RdRp domain in TMUV attenuates its in vivo infectivity.

The NS5 Val711 to Ala Mutation Reverts Rapidly when the R543K Substitution Is Absent
Because the chimeric JXSP-310 RdRp virus consistently exhibited attenuated phenotypes, we determined whether specific amino acid substitutions were associated with these phenotypic alterations. A single arginine to lysine (R543K) or valine to alanine (V711A) in NS5 was introduced into the JXSP 2-4 genome by site-directed mutagenesis ( Figure 4A). Sequence analysis indicated that the R543K mutation site appeared to be genetically stable after 10 serial passages in BHK-21 cells ( Figure 4B). When the multi-step growth curves were compared with those from the parental virus, the single-point R543K mutant we constructed exhibited growth kinetics comparable to that of the parental virus in BHK-21 or DEF cells ( Figure 4D). This indicated that the R543K substitution alone did not influence viral growth in vitro. However, the Val to Ala mutation (GTA to GCA) reverted to Val (GTA) via a single nucleotide substitution in the V711A mutant after one passage in BHK-21 cells ( Figure 4C). Additional attempts at virus rescue resulted in the same reversion.      mutant we constructed exhibited growth kinetics comparable to that of the parent in BHK-21 or DEF cells ( Figure 4D). This indicated that the R543K substitution al not influence viral growth in vitro. However, the Val to Ala mutation (GTA to G verted to Val (GTA) via a single nucleotide substitution in the V711A mutant af passage in BHK-21 cells ( Figure 4C). Additional attempts at virus rescue resulted same reversion.

The NS5 R543K Mutant Is Less Attenuated In Vivo when the V711A Mutation Is A
To investigate whether the single R543K change in NS5 would affect the abilit virus to replicate in vivo, individual groups of 7-day-old ducklings (n = 15) were s neously infected with JXSP2-4, JXSP-310 RdRp , or JXSP-R543K and the infectivity titer

The NS5 R543K Mutant Is Less Attenuated In Vivo when the V711A Mutation Is Absent
To investigate whether the single R543K change in NS5 would affect the ability of the virus to replicate in vivo, individual groups of 7-day-old ducklings (n = 15) were subcutaneously infected with JXSP 2-4 , JXSP-310 RdRp , or JXSP-R543K and the infectivity titers in the representative tissues were quantified by plaque assays. Ducklings infected with JXSP-310 RdRp had significantly lower infectious virus titers in their sera and in various organs than those infected with parental JXSP 2-4 ( Figure 5), which were consistent with the abovementioned results. By comparison, the single-point R543K mutant appeared to have a less attenuated replication capacity when compared with the parental virus and JXSP-310 RdRp . However, the peak viral titers in the serum, spleen, heart and kidneys were significantly lower, although the differences in the lung and thymus were not statistically significant when compared with the parental virus. This suggests that the R543K and V711A mutations in combination confer better in vivo infectivity attenuation.
To compare the innate immune responses to these viruses, we examined the expression of interferons and early pro-inflammatory cytokines in the spleen by qRT-PCR. As shown in Figure 6, the innate immune responses in ducklings were strongly elicited by infection with TMUV JXSP2-4 and the mutant viruses. The increase of transcript abundance in the spleen was positively correlated with the replication levels of these viruses. The magnitude of spleen-tissue interferon (IFN)-α, IFN-β, IFN-γ, interleukin (IL)-1β and IL-6 expression induced by infection with the R543K mutant was significantly lower than those of JXSP2-4, but higher than those of the JXSP-310 RdRp -infected ducklings. Overall, the early induction of interferons and pro-inflammatory cytokines reflected the infectivity profiles of the viruses. To compare the innate immune responses to these viruses, we examined the expression of interferons and early pro-inflammatory cytokines in the spleen by qRT-PCR. As shown in Figure 6, the innate immune responses in ducklings were strongly elicited by infection with TMUV JXSP 2-4 and the mutant viruses. The increase of transcript abundance in the spleen was positively correlated with the replication levels of these viruses. The magnitude of spleen-tissue interferon (IFN)-α, IFN-β, IFN-γ, interleukin (IL)-1β and IL-6 expression induced by infection with the R543K mutant was significantly lower than those of JXSP 2-4 , but higher than those of the JXSP-310 RdRp -infected ducklings. Overall, the early induction of interferons and pro-inflammatory cytokines reflected the infectivity profiles of the viruses.

Coupled Mutations of NS5 R543K and V711A Decrease Viral RNA Synthesis
To further analyze the effect of RdRp mutations on virus replication, DEF cells were separately infected with either of the mutated viruses or the parental virus at an MOI of 0.1 and RNA copy numbers were measured by qRT-PCR. The JXSP-310 RdRp RNA copy numbers were 5-to 15-fold lower than those of the parental virus ( Figure 7A), whereas the RNA copy numbers for chimeric JXSP-310 MTase and JXSP-310 Linker were similar to those of the parental virus-infected DEF cells, ranging from equal to 2-fold differences. Consistent with the in vitro infectivity observed herein and in common with other chimeric viruses, no significant RNA replication deficiency was observed in the single-point R543K mutant. The RNA copy numbers in the virus-infected DEFs correlated with the infectious virus titers detected by the plaque assays. In agreement with the RNA replication results, the western blot showed that E protein expression level of the JXSP-310 MTase , JXSP-310 Linker and JXSP-R543K was similar to that of JXSP [2][3][4] in DEFs at 24 and 36 hpi. However, the E protein production in JXSP-310 RdRp -infected cells was 2.7-fold and 2.3-fold lower, respectively, than that of JXSP 2-4 ( Figure 7B,C).

Coupled Mutations of NS5 R543K and V711A Decrease Viral RNA Sy
To further analyze the effect of RdRp mutations on virus replica separately infected with either of the mutated viruses or the parenta 0.1 and RNA copy numbers were measured by qRT-PCR. The JXSP numbers were 5-to 15-fold lower than those of the parental virus ( the RNA copy numbers for chimeric JXSP-310 MTase and JXSP-310 Linker w of the parental virus-infected DEF cells, ranging from equal to 2-fo sistent with the in vitro infectivity observed herein and in common viruses, no significant RNA replication deficiency was observed in the mutant. The RNA copy numbers in the virus-infected DEFs correlate virus titers detected by the plaque assays. In agreement with the RNA the western blot showed that E protein expression level of the JXSP-3 and JXSP-R543K was similar to that of JXSP2-4 in DEFs at 24 and 36 protein production in JXSP-310 RdRp -infected cells was 2.7-fold and 2.3 tively, than that of JXSP2-4 ( Figure 7B,C).   IFN-α, IFN-β, IFN-γ, IL-1β and IL-6 in the spleens were analyzed by qRT-PCR and normalized to GADPH. Asterisks indicate significant differences between JXSP2 and JXSP-R543K, (*, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001). Number signs indicate significant differences between JXSP-310 RdRp and JXSP-R543K (##, p < 0.01; ####, p < 0.0001). The data were tested for statistical significance by two-way multiple ANOVA comparisons.

Coupled Mutations of NS5 R543K and V711A Decrease Viral RNA Synthesis
To further analyze the effect of RdRp mutations on virus replication, DEF cells wer separately infected with either of the mutated viruses or the parental virus at an MOI o 0.1 and RNA copy numbers were measured by qRT-PCR. The JXSP-310 RdRp RNA cop numbers were 5-to 15-fold lower than those of the parental virus ( Figure 7A), wherea the RNA copy numbers for chimeric JXSP-310 MTase and JXSP-310 Linker were similar to thos of the parental virus-infected DEF cells, ranging from equal to 2-fold differences. Con sistent with the in vitro infectivity observed herein and in common with other chimer viruses, no significant RNA replication deficiency was observed in the single-point R543 mutant. The RNA copy numbers in the virus-infected DEFs correlated with the infectiou virus titers detected by the plaque assays. In agreement with the RNA replication result the western blot showed that E protein expression level of the JXSP-310 MTase , JXSP-310 Link and JXSP-R543K was similar to that of JXSP2-4 in DEFs at 24 and 36 hpi. However, the protein production in JXSP-310 RdRp -infected cells was 2.7-fold and 2.3-fold lower, respe tively, than that of JXSP2-4 ( Figure 7B,C).

Discussion
Our previous studies showed that 23 amino acid substitutions in multiple genes cocontributed to the full attenuation of this virus [13]. Here, the effects of mutations in NS5 protein on viral attenuation were analyzed. We demonstrated, for the first time, that a single R543K substitution in the RdRp domain of NS5 impaired the infectivity of TMUV in vivo and that the virus was more strongly attenuated when R543K substitution was coupled with the V711A mutation.
NS5, the largest (905 amino acids) and most conserved flavivirus protein, contains three independent domains: N-terminal methyltransferase (MTase) domain, which is involved in RNA capping and prevents host recognition as well as promotes translation [21][22][23]; C-terminal RdRp domain, which is responsible for de novo RNA synthesis [24,25]; and the short Linker in between, which regulates the cross talk between RdRp and MTase dmoains [26]. The critical roles played by NS5 in viral replication and interferon antagonism make it a prime target for antiviral therapy and vaccine development [27,28]. To map the specific amino acid substitutions in the NS5 that contributed to duck TMUV attenuation, we generated a series of chimeric viruses by replacing the individual domains of the wild-type JXSP 2-4 with those of the attenuated JXSP-310. The in vitro cell cultures and in vivo ducklings infection results showed that replacing JXSP 2-4 with either the mutated MTase or Linker domain did not significantly alter the chimeric viruses' phenotypes, whereas the R543K and V711A mutations in the RdRp domain resulted in significant lower infectivity titers when compared with its parental virus. These results demonstrated that mutations in MTase and Linker had no influence on the activity of MTase and the cross-talk between RdRp and MTase. Thus, the adaptive mutations in the RdRp domain selected from laboratory cell culture passages changed the fitness of TMUV.
Generally, the plaque size is determined by virus replication rates and cell-to-cell spread. Previous study showed that PDK53 strain of DENV replicated rapidly, but giving rise to small plaques due to its inability to evade antiviral responses that constrained its spread in primate or rodent cells [29]. Figure 6 showed that JXSP 2-4 -infected ducklings produced higher inflammatory cytokines compared with JXSP-310 RdRp ; thus, the enhanced immune responses might account for the smaller plaque size.
Flavivirus RdRp, whose overall shape resembles a right hand with fingers, palm and thumb subdomains, contains seven conserved motifs (A to G) [30,31]. Studies showed that amino acid substitutions in the RdRp domain affected genome replication and disease development [16,32]. In foot-and-mouth disease virus, replacing the conserved residues in the motif A of RdRp resulted in virulence attenuation in mice by altering the polymerase's fidelity [33]. Predictive structural models of this protein suggest that residue Arg543 is highly conserved within the motif A (YADDTAGWDTR/KIT). The motif A forms the catalytic site in the RdRp with the motif C in the structural palm domain. The single-point R543K mutation in the RdRp described herein did not affect TMUV replication in vitro, as determined by plaque assays and RNA copy quantification. However, significantly lower viremia and viral titers of the mutant were observed in most of the tested parenchymal organs, suggesting that factors other than the physiological condition of the ducklings restricted virus proliferation. One possible explanation is that replacing arginine with lysine partially retains the properties of arginine through conservation of structure; however, the smaller lysine residue may affect interactions with other viral components and possibly host-factor proteins, too. Interchanging lysine and arginine at residue 185 of the PA polymerase of avian H5N1 virus TY165 significantly altered its virulence in mice [34]. Homology modelling analysis predicts that R543 is located in the TMUV RdRp module's α helix, which forms three hydrogen bonds with residues Trp540, Ser690 and Lys691. Contrastingly, the R543K substitution is predicted to break the hydrogen bond between Arg543 and Ser690, which is the first residue of motif D (Figure 8). Interestingly, the NS5 R543K mutation was also observed in a recent laboratory attenuated TMUV strain [35]. The structure of NS5 RdRp of TMUV shows the fingers, palm and thumb domains colored in blue, grey and orange, respectively. The TMUV NS5 protein model was created online using Swiss-Model with reference to the crystal structure of the JEV NS5 protein (PDB identifier [ID] 4k6m). (B) R543 is located in motif A and possibly forms three solid hydrogen bonds with W540, S690 and K691 spatially. V711 is close to motif E and forms one solid hydrogen bond with W708. (C) An amino acid mutation at site 543 changes arginine to lysine, which breaks the hydrogen bond between sites 543 and 690. A valine to alanine change at site 711 slightly changes the distance between motif A and residue W708. This figure was generated in PyMol.
The second mutated residue Val711 is just before motif E in the thumb subdomain of RdRp. However, this mutation rapidly reverted when the rescued single-point mutant was cultured in BHK-21 cells, implying that the V711A mutation was deleterious to viral proliferation and fitness in cell cultures in the absence of the R543K substitution. We noted that a reverse mutation of alanine to valine in the corresponding site of the duck TMUV Du/CH/LSD/110128 strain was reported to occur during serial passaging in chicken embryos [10]; this Val711 to Ala mutation occurred after 30 passages and the reversion was identified in the 50th-passage virus. This also suggested that there was a strong preference for alanine at this position. This amino acid substitution, however, was stable in the R543K/V711A double mutant. In our previous study, the R543K mutation in JXSP2-4 was identified after the 150th passage during BHK-21 cell passaging. The mutation occurred far earlier than the V711A change, which was detected from the 280th passage [13]. The stability of the R543K/V711A double mutant suggests that the R543K mutation somewhat counterbalances the negative effect. It is not clear whether this residue change affects the conformation and/or connection with the palm subdomain because the correct conformation of motif E is important for de novo initiation and elongation efficiency [36]. Although further biochemical studies are needed, we observed that the double mutations caused significantly lower RNA replication and infectivity reduction in infected DEFs. Moreover, the more dramatic attenuation of the R543K/V711A double mutant suggests that the two mutations act in concert to significantly alter the biological properties of TMUV, possibly through the additive effects of these mutations on the functionality of the viral proteins and/or genes. Interestingly, the amiloride-induced coxsackie virus B3 variant with an Ala to Val change at position 372 in the RdRp (3Dpol-A372V), which corresponds to residue 711 in TMUV RdRp, was found to increase enzyme fidelity [37]. In the study by Harrison et al., the A372V mutants of coxsackie B3 virus that were screened by amiloride derivative treatment contained the aspartic acid to glycine substitution in residue 48 of the 2A protein (2A-D48G) and double mutations were found to play a synergistic effect in replication, resulting in higher viral titers in the presence of antiviral compounds [38]. RdRp of TMUV shows the fingers, palm and thumb domains colored in blue, grey and orange, respectively. The TMUV NS5 protein model was created online using Swiss-Model with reference to the crystal structure of the JEV NS5 protein (PDB identifier [ID] 4k6m). (B) R543 is located in motif A and possibly forms three solid hydrogen bonds with W540, S690 and K691 spatially. V711 is close to motif E and forms one solid hydrogen bond with W708. (C) An amino acid mutation at site 543 changes arginine to lysine, which breaks the hydrogen bond between sites 543 and 690. A valine to alanine change at site 711 slightly changes the distance between motif A and residue W708. This figure was generated in PyMol.
The second mutated residue Val711 is just before motif E in the thumb subdomain of RdRp. However, this mutation rapidly reverted when the rescued single-point mutant was cultured in BHK-21 cells, implying that the V711A mutation was deleterious to viral proliferation and fitness in cell cultures in the absence of the R543K substitution. We noted that a reverse mutation of alanine to valine in the corresponding site of the duck TMUV Du/CH/LSD/110128 strain was reported to occur during serial passaging in chicken embryos [10]; this Val711 to Ala mutation occurred after 30 passages and the reversion was identified in the 50th-passage virus. This also suggested that there was a strong preference for alanine at this position. This amino acid substitution, however, was stable in the R543K/V711A double mutant. In our previous study, the R543K mutation in JXSP 2-4 was identified after the 150th passage during BHK-21 cell passaging. The mutation occurred far earlier than the V711A change, which was detected from the 280th passage [13]. The stability of the R543K/V711A double mutant suggests that the R543K mutation somewhat counterbalances the negative effect. It is not clear whether this residue change affects the conformation and/or connection with the palm subdomain because the correct conformation of motif E is important for de novo initiation and elongation efficiency [36]. Although further biochemical studies are needed, we observed that the double mutations caused significantly lower RNA replication and infectivity reduction in infected DEFs. Moreover, the more dramatic attenuation of the R543K/V711A double mutant suggests that the two mutations act in concert to significantly alter the biological properties of TMUV, possibly through the additive effects of these mutations on the functionality of the viral proteins and/or genes. Interestingly, the amiloride-induced coxsackie virus B3 variant with an Ala to Val change at position 372 in the RdRp (3Dpol-A372V), which corresponds to residue 711 in TMUV RdRp, was found to increase enzyme fidelity [37]. In the study by Harrison et al., the A372V mutants of coxsackie B3 virus that were screened by amiloride derivative treatment contained the aspartic acid to glycine substitution in residue 48 of the 2A protein (2A-D48G) and double mutations were found to play a synergistic effect in replication, resulting in higher viral titers in the presence of antiviral compounds [38].
In combination with the phenotypic differences observed between the single-point R543K mutant and the fact that the R543K/V711A double mutant consistently displayed reduced viral replication in vitro and in vivo, it appears that both mutations significantly