{'fraction_non_alphanumeric': 0.03526046171246092, 'fraction_numerical': 0.0445412813518728, 'mean_word_length': 4.739404352806415, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Nonstructural protein 9 residues 586 and 592 are critical sites in determining the replication efficiency and fatal virulence of the Chinese highly pathogenic porcine reproductive and respiratory syndrome virus

The highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) has caused huge economic losses to the swine industry in China. Understanding the molecular basis in relation to the virulence of HP-PRRSV is essential for effectively controlling clinical infection and disease. In the current study, we constructed and rescued a serial of mutant viruses in nsp9 and nsp10 based on the differential amino acid sites between HP-PRRSV JXwn06 and LP-PRRSV HB-1/3.9. The replication efficiency in pulmonary alveolar macrophages (PAMs) and the pathogenicity of the mutant viruses for piglets were analyzed. Our results showed that the mutation of Thr to Ala in 586 and Ser to Thr in 592 of nsp9 decreased the replication efficiency of HP-PRRSV in PAMs, and could attenuate its virulence for piglets, suggesting that the residues 586 and 592 of nsp9 are critical sites natively in determining the fatal virulence of the Chinese HP-PRRSV for piglets.

Tremendous molecular epidemiological data have shown that either type 1 or type 2 PRRSV exhibits broad genetic variation and diversified strains (Shi et al., 2010). The persistent evolution nature of PRRSV has led to the emergence of novel and variant strains with higher pathogenicity or virulence (Han et al., 2006;Karniychuk et al., 2010;Tian et al., 2007;Wang et al., 2015;Zhou et al., 2015). The Chinese highly pathogenic PRRSV (HP-PRRSV) with a unique hallmark of 30-amino acid deletion in the viral nsp2 emerged and prevailed in 2006 (Tian et al., 2007), causing colossal economic losses to the swine production (Zhou and Yang, 2010). A number of studies have provided essential evidence for understanding the pathogenesis of the Chinese HP-PRRSV (Han et al., 2017). HP-PRRSV infection is shown to display expanded tissue tropism and increased virus loads in a variety of organs (Hu et al., 2013;Li et al., 2012), and cause more severe lung injuries and histopathological changes (Han et al., 2014;Hu et al., 2013), and thymus atrophy, depletion of thymocytes, and apoptotic cell death of thymic epithelial cells (Guo et al., 2013;He et al., 2012;Li et al., 2014b). Importantly, much attention has been paid to the molecular basis in relation to its virulence. Our previous study indicated that the 30-amino acid deletion in its nsp2 is not related to the fatal virulence for piglets (Zhou et al., 2009). Subsequent works clearly revealed that the nsp9 and nsp10 together are closely related to the replication efficiency both in vitro and in vivo, and contribute to the fatal virulence of the Chinese HP-PRRSV (Li et al., 2014c). However, it remains unclear whether the Chinese HP-PRRSV shares the amino acid sites in nsp9 or/and nsp10 determining its fatal virulence.
In arterivirus, the nsp9 possesses RNA-dependent RNA polymerase (RdRp) activity, and the nsp10 contains metal binding region with RNA helicase activity, which are considered to be assembled into the replication and transcription complex (RTC), and to play a crucial role in virus replication (Beerens et al., 2007;van Dinten et al., 2000). In terms of the structural model of coronavirus nsp12 (Gorbalenya et al., 1989;Xu et al., 2003), the nsp9 of PRRSV is predicated to consist of at least two domains, an N-terminal domain with RdRp-associated nucleotidyltransferase (NiRAN) activity and a canonical RdRp occupying its Cterminal domain (Lehmann et al., 2015a). A conserved SDD motif and Asp residue within the RdRp domain of nsp9 are recognized to be critical for RNA polymerase activity and RNA synthesis (Lehmann et al., 2016;Snijder et al., 1990;Subissi et al., 2014). The nsp10 of arterivirus belongs to superfamily 1 helicase (HEL1), and its helicase activity domain has similar structure to a cellular helicase Upf1 that participates in nonsense-mediated mRNA decay pathway, implying that the nsp10 might also play a role in post-transcriptional quality control of viral RNA (Deng et al., 2014;Lehmann et al., 2015b). In the present study, we constructed a serial of mutant viruses with the mutated amino acid residues in the nsp9 and nsp10 of PRRSV by using reverse genetic approach, based on the differential amino acid sites in nsp9 and nsp10 between HP-PRRSV and low pathogenic PRRSV (LP-PRRSV), and analyzed the replication efficiency and pathogenicity of mutant viruses for piglets, in an attempt to identify the critical amino acids in relation to the fatal virulence in nsp9 or/and nsp10 of HP-PRRSV.

Ethics statements
The animal experiments in this study were approved by the Laboratory Animal Ethical Committee of China Agricultural University with the license number (CAU20160828-2). All animal experiments were performed according to the Chinese Regulations of Laboratory Animals-The Guidelines for the Care of Laboratory Animals (Ministry of Science and Technology of People's Republic of China) and Laboratory Animal-Requirements of Environment and Housing Facilities (National Laboratory Animal Standardization Technical Committee).

Construction of infectious cDNA clones with mutated sites in nsp9 and nsp10 of PRRSV and virus rescue
By aligning the amino acid sequences, totally nine different amino acid residues existed in nsp9, including nsp8 that is considered an Nterminal domain of nsp9 (Lunney et al., 2016), and nsp10 between the HP-PRRSV JXwn06 and LP-PRRSV HB-1/3.9 with an amino acid deletion in nsp9 that is not counted (Zhou et al., 2009). Of them, four amino acids located in the position 427, 586, 592 and 609 of nsp9, and five located in the position 11, 51, 69, 296 and 408 of nsp10 (Table 1). According to the strategies as previously described (Zhou et al., 2009), respective nucleotide mutation targeting the differential residues in the fragment C (containing nsp9-and nsp10-coding region) of infectious cDNA clones was conducted with the unique restriction enzymes NheI and AscI (New England Biolabs, Ipswich, MA). The plasmid Pjet1.2/ blunt (Fisher Scientific) was used to construct the plasmid Pjet1.2-C, which was then inserted the fragment C. Then the target nucleotides were mutated in the plasmid by fast mutagenesis system (TransGen, Beijing, China). The mutant plasmids and the backbones of infectious cDNA clones were digested by the restriction enzymes NheI and AscI, and finally the mutated fragment C was inserted into the backbone to generate the full-length infectious cDNA clone with respective mutated site.
MARC-145 cells were transfected with the mutant full-length cDNA clone by Lipofectamine LTX (Fisher Scientific) and serially passaged for three times in MARC-145 cells. Then the rescued viruses were examined by indirect immunofluorescence assay (IFA) using the PRRSV N-specific monoclonal antibody (mAb) SDOW17 (Rural Technologies, Brookings, SD). The RNAs of third-passage rescued viruses were extracted, and amplified by RT-PCR, and the amplified fragment was sequenced to further confirm the mutated sites.

Animals and animal trials
Healthy, six-week old, landrace piglets were obtained from Beijing Center for SPF Swine Breeding & Management that is known to be free of PRRSV. All piglets were confirmed to be negative for PRRSV infection by commercial IDEXX ELISA kits and RT-PCR detection in our laboratory. The piglets were raised in the animal facilities at China Agricultural University (CAU).

The growth kinetics of rescued viruses in PAMs
Primary PAMs were infected with respective virus at a multiplicity of infection (MOI) of 0.1. The cell cultures were collected at 12 h, 24 h, 36 h, 48 h, 60 h and 72 h post-infection (pi) and the virus titers were determined using a microtitration infectivity assay, as previously described (Li et al., 2014c).

Quantification of PRRSV nsp9 gene by real-time PCR
Primary PAMs were infected with respective virus at a MOI of 1. The PAMs were collected prior to infection and at 12 h pi, and total RNA was extracted from the cells using TRIZOL Reagent according to the manufacturer's instructions (Fisher Scientific). 200 ng of RNAs were used for further RT-PCR. Real-time PCR was carried out using the SYBRGreen (ABI) on an ABI 7500 Fast Real Time PCR system (Life Technologies, Grand Island, NY). The primers used for detecting the copy numbers of PRRSV nsp9 gene were designed as previously described (Spear and Faaberg, 2015). Eight serial dilutions of plasmid with the copy number from 10 1 to 10 8 copies/μl served as template to run reaction on the ABI 7500 quantitative PCR machine to generate a standard curve. Cycling parameters used for all the reactions were as follows: 50°C for 2 min; 95°C for 2 min; and 40 cycles of 95°C for 15 s, 60°C for 60 s. Data collection was performed during the 60°C elongation step. The standard curve was automatically generated and the numbers of nsp9 gene copy were calculated with the 7500 System SDS software (Applied Biosystems, Foster City, CA).

Microscopic pathological changes and immunohistochemistry (IHC) examinations
Microscopic pathological changes and IHC examinations for lung of each piglet were performed as previously described (Halbur et al., , 1996. Briefly, lung tissue samples were taken, and fixed with 4% paraformaldehyde solution at room temperature for 48 h and then processed by conventional histopathological procedures. Two sections for each sample were prepared. One section was stained with hematoxylin and eosin (H&E) for microscopic pathological change examination, and another was stained with the mAb specific for PRRSV N protein (SDOW17) at 1:1000 dilutions for IHC examination of PRRSV antigen. The scores of lung microscopic lesions were blindly evaluated from 0 to 4, which accounted for the distribution and severity of interstitial pneumonia. The IHC scores of PRRSV antigen were conducted through a range score of 0-4 for evaluating the numbers of PRRSVpositive cells.

Viremia analysis
Serum samples of animals were collected at 0, 3, 5, 7, 10, 14 and 21 dpi to analyze the viral loads in the sera of the inoculated animals using a microtitration infectivity assay.

Statistical analysis
The significant differences were analyzed using a one-way or twoway RM ANOVA in the GraphPad Prism (version 5.0) software. A nonparametric test was employed to determine the significant difference in microscopic lung lesion scores and immunohistochemistry scores. Differences were considered statistically significant at a p value of ＜0.05.

Results
3.1. The residues in the position 586 and 592 of nsp9 determined the replication efficiency of HP-PRRSV in PAMs To screen for the potential amino acid sites affecting the replication efficiency of PRRSV, we used the chimeric full-length cDNA clones (pWSK-RvJHn9n10 and pWSK-RvHJn9n10), and constructed a serial of full-length cDNA clones with mutated differential amino acid sites in nsp9 and nsp10 between JXwn06 and HB-1/3.9. Totally, 18 viruses with a single residue mutation were rescued ( Table 2). The replication efficiencies of these mutated viruses in PAMs were examined and compared by drawing growth curves, together with the parental viruses (RvJXwn, RvHB-1/3.9, RvJHn9n10 and RvHJn9n10). As shown in Fig. 1A, compared with the parental virus RvHB-1/3.9, the RvHJn9n10 with swapping nsp9-and nsp10-coding regions of JXwn06 had an increased replication efficiency in PAMs. While the HJn9n10-T586A with the mutation of Thr to Ala in 586 of nsp9 had lower titers than RvHJn9n10, with significant differences at 36 h and 48 h pi (p＜0.01); the HJn9n10-S592T with the mutation of Ser to Thr in 592 of nsp9 also showed the decreased titers, with a similar level to RvHB-1/3.9, and its titers were significantly lower than RvHJn9n10 at 36 h (p＜0.01) and 48 h pi (p＜0.05). The viruses with other mutations in nsp9 and nsp10 showed similar growth kinetics to the RvHJn9n10 (Fig. 1B). On other hand, we measured the growth kinetics of mutant viruses based on RvJHn9n10 with swapping nsp9-and nsp10-coding regions of HB-1/ 3.9. The results showed that RvJHn9n10 had a decreased replication efficiency than RvJXwn, but the viruses with mutated differential amino acid sites did not exhibit an increased growth efficiency than RvJHn9n10, even the viruses with the mutation of Ala to Thr in 586 and Thr to Ser in 592 of nsp9 (JHn9n10-A586T and JHn9n10-T592S) had slightly decreased titers ( Fig. 1C and D). Obviously, these data indicate that it is insufficient to increase the replication efficiency of the chimeric PRRSV RvJHn9n10 by mutating only one amino acid.
To further confirm the influence of the two residues (586 and 592 in nsp9) on PRRSV replication, we used the full-length cDNA clone (pWSK-RvJXwn) of JXwn06 and constructed a serial of full-length cDNA clones with mutated differential amino acid sites in nsp9 and rescued the viruses (JXwn-T586A, JXwn-S592T, JXwn-T586A/S592T, JXwn-A427T/D609G) (Table 2). Meanwhile, considering that the nsp10 of HP-PRRSV also plays a role in its fatal virulence for piglets (Li et al., 2014c), we used the chimeric full-length infectious cDNA clone pWSK-RvJHn10 as backbone and mutated the corresponding amino acid sites, and rescued the viruses (JHn10-T586A, JHn10-S592T, JHn10-T586A/ S592T and JHn10-A427T/D609G) (Table 2). Firstly, we evaluated the replication efficiency of these mutant viruses and their parental viruses in PAMs. As shown in Fig. 2A, the viruses with the mutation of Thr to Ala in 586 or Ser to Thr in 592 or both had lower titers than the parental virus RvJXwn, with significant difference at 24 h pi (p＜0.05), and delayed peak. While the viruses with the mutation of other two residues together showed no change of the replication efficiency. Similar data were obtained in the mutant viruses based on pWSK-RvJHn10. The titers of the viruses with the mutation of Thr to Ala in 586 or Ser to Thr in 592 or both decreased and their titer peaks delayed, while the viruses with the mutation of other two residues together had no significant differences in viral titers, compared with RvJHn10 (Fig. 2B).
Moreover, PAMs were infected with respective virus at a MOI of 1, and the copy numbers of nsp9 gene for the viruses were detected in order to evaluate the RNA synthesis efficiency by real-time PCR. As shown in Fig. 2C, the copy numbers of nsp9 gene of all viruses had a similar level at 0 h pi, while at 12 h pi, the average copy numbers of JXwn-T586A, JXwn-S592T and JXwn-T586A/S592T reached 7.16log10/μl, 7.46log10/μl, and 6.92log10/μl, respectively, which were remarkable lower than that of RvJXwn (7.74log10/μl) (p＜0.01 or 0.001), in addition to the JXwn-A427T/D609G with a similar level to RvJXwn. The RNA synthesis efficiency of the mutant viruses based on pWSK-RnJHn10 was also measured. The mutant viruses (JHn10-T586A, JHn10-S592T and JHn10-T586A/S592T), had the nsp9 copy numbers of 7.31log10/μl, 7.23log10/μl, and 6.95log10/μl, respectively, which were significantly lower than that of RvJHn10 (7.75log10/ul) (p＜ 0.001), while the virus JHn10-A427T/D609G showed no statistically significant differences in nsp9 copy numbers with RvJHn10 (Fig. 2D). Taken together, the above results indicated that the mutations of the  residues 586 and 592 in nsp9 decreased the replication efficiency of HP-PRRSV in PAMs, while the mutations of the residues 427 and 609 in nsp9 did not, suggesting that the residues 586 and 592 in nsp9 contribute to the replication efficiency of HP-PRRSV.

Mutation of Thr to Ala in 586 and Ser to Thr in 592 of nsp9 attenuated the virulence of HP-PRRSV for piglets
The pathogenicity of the mutant viruses based on pWSK-RvJXwn and pWSK-RvJHn10 for piglets was analyzed. The piglets were inoculated with respective virus, and rectal temperatures of the inoculated piglets were daily recorded. JXwn-A427T/D609G-infected piglets displayed a similar temperature reaction to RvJXwn-infected piglets which reached the highest temperature (41.5°C) and maintained above 41°C from 4 to 12 dpi, while the piglets inoculated with the viruses with the mutation of Thr to Ala in 586 and Ser to Thr in 592 of nsp9 or both had a slow rising of body temperature, with the highest temperature of 41.2°C and short lasting of above 41°C (Fig. 3A). JHn10-A427T/D609G-infected piglets presented a similar temperature reaction to RvJHn10-infected group with the highest temperature (41.3°C) and hovered above 41°C from 8 to 14 dpi, while the piglets inoculated with the viruses with the mutation of Thr to Ala in 586 and Ser to Thr in 592 of nsp9 or both had a body temperature of less 41°C (Fig. 3B).
As a whole, the mutation of the residue 586 or 592 and the both in nsp9 of RvJXwn06 or RvJHn10 resulted in the alleviated temperature reaction, increased ADG and decreased mortality of the inoculated pigs, while the mutation of the residue 427 and 609 in nsp9 did not.
The lung microscopic lesions of the inoculated piglets were examined. The lung microscopic lesions of dead piglets post-inoculation and euthanized piglets at the termination of trials were shown in Fig. 4A and B, respectively. All the lungs of dead piglets in the inoculated groups exhibited severe histopathological changes with higher histopathological scores, including large amounts of inflammatory cells and fibrin exudation of the alveolar and bronchioles spaces, necrotic debris and hemorrhage, and the thicken interlobular septal and the destroyed lung structure (Fig. 4A). The euthanized piglets in JXwn-T586A-and JXwn-T586A/S592T-inoculated groups displayed moderate histopathological lesions (Fig. 4B). Compared with RvJHn10-inoculated piglets, the piglets inoculated with the mutant viruses showed the alleviated lung lesions with less inflammatory cells and fibrin exudation, and lower histopathological scores, and the piglets in JHn10-A427T/ D609G-inoculated group exhibited no differences in lung lesions with RvJHn10-inoculated piglets. No microscopic lesions were observed in the lungs of control piglets which survived until the trial termination. These results showed that the mutation of the residue 586 or both the residue 586 and 592 in nsp9 of RvJXwn06 or RvJHn10 alleviated the lung lesions of the survived pigs.
The virus loads in the sera of inoculated piglets were examined using a microtitration assay. The data showed that the virus titers in sera of piglets infected with JXwn-T586A, JXwn-T586A/S592T were lower than those of RvJXwn-infected group with significant differences at 3, 5, 10 dpi (p＜0.01, p＜0.001, p＜0.05), and the virus titers in JXwn-S592T-infected group were also slightly lower than those of RvJXwn-infected group, while JXwn-A427T/D609G-infected group had  L. Xu et al. Virology 517 (2018) 135-147 no statistically significant differences in virus loads with RvJXwn-infected group (Fig. 5A). The virus titers in sera of the piglets infected with JHn10-T586A, JHn10-S592T, and JHn10-T586A/S592T were lower than those of RvJHn10-infected group with significant differences at some time points, and JHn10-T586A/S592T-infected piglets had the lowest virus loads, with significant differences at 3, 5, 7, 10, 14 dpi, while JHn10-A427T/D609G-infected piglets had no differences in virus loads with RvJHn10-infected group (Fig. 5B). The results indicated that the mutation of the residue 586 or 592 and the both in nsp9 of RvJXwn06 or RvJHn10 led to the declined virus load in the sera of the inoculated pigs. Moreover, PRRSV antigen in lungs of the inoculated piglets was examined by immunohistochemistry. The PRRSV-positive signals were distributed in all the lung section in dead piglets, and the scores of PRRSV antigen in lungs of the piglets inoculated with mutant viruses were lower than their respective parental viruses (Fig. 6A). For euthanized piglets, the PRRSV-positive signals in lungs of JHn10-T586A-, JHn10-S592T-and JHn10-T586A/S592T-inoculated piglets were less than those of RvJHn10-inoculated piglets, whereas JHn10-A427T/ D609G-inoculated piglets showed no differences with RvJHn10-inoculated piglets (Fig. 6B). No PRRSV-positive signals were observed in lungs of control piglets. Collectively, the above data revealed that the mutation of Thr to Ala in 586 and Ser to Thr in 592 of HP-PRRSV nsp9 alleviated the temperature reaction and lung lesions of the infected piglets, and declined the replication efficiency of HP-PRRSV in vivo, and delayed the death time of the infected piglets, and resulted in the decreased mortality of piglets, suggesting that the mutation of these two residues in nsp9 of HP-PRRSV attenuates its virulence for piglets.
3.3. Mutation of the residues in 586 and 592 in nsp9 of LP-PRRSV did not enhance its pathogenicity for piglets To address whether the mutation of amino acid sites in 586 and 592 of LP-PRRSV nsp9 enhances its pathogenicity for piglets, we performed the mutation of Ala to Thr in 586 and Thr to Ser in 592 in nsp9 of HB-1/ 3.9 using another chimeric full-length infectious clone pWSK-RvHJn10 and rescued the mutant virus (Table 2). Compared with RvHJn10-infected group, the piglets infected with HJn10-A586T/T592S showed no differences in body temperature reaction, with the highest temperatures of 40.2°C (Fig. 7A), ADG (Fig. 7B), and virus loads in sera (Fig. 7C). Similar to RvHJn10, HJn10-A586T/T592S infection induced mild inflammatory reaction with low scores of microscopic lung lesions (Fig. 8A), and all the piglets survived during the whole period of experiments. Meanwhile the viral titers in sera of HJn10-A586T/T592Sinoculated piglets were lower and only a fewer PRRSV-positive signals distributed in lung (Fig. 8B). These results indicated that the mutation of the residue 586 or 592 and the both in nsp9 of RvHJn10 could not impact its pathogenicity for piglets, suggesting that the mutation of the two residues in nsp9 did not enhance the pathogenicity of LP-PRRSV for piglets.

Discussion
Chinese swine industry has been experiencing the damage of PRRSV for over two decades, and this situation has become even worse since HP-PRRSV emerged in 2006 (Han et al., 2017). Despite numerous efforts on HP-PRRSV pathogenesis and antiviral immunity investigation for several years, to effectively control the clinical diseases caused by this virus remains a great challenge for the Chinese pig production. As a consequence, it is essential to find out the pathogenesis of HP-PRRSV. In our previous work, we utilized a HP-PRRSV strain JXwn06 and a LP-PRRSV strain HB-1/3.9 with genetic similarity to HP-PRRSV and swapped different viral coding regions by infectious cDNA clones, and clearly demonstrated that the nsp9 and nsp10 together were related to viral replication efficiency of HP-PRRSV in vivo and in vitro, contributing its fatal virulence for piglets (Li et al., 2014c). To further dissect the molecular basis of HP-PRRSV virulence, we screened for the critical amino acids in relation to the fatal virulence of HP-PRRSV in nsp9 and nsp10 in the current study.
As well known, many viruses can strengthen their pathogenicity by increasing viral replication efficiency (Brault et al., 2007;Hanley et al., 2002;Watanabe et al., 2009). Similar findings have been documented on PRRSV. In comparison to LP-PRRSV HB-1/3.9, HP-PRRSV JXwn06 exhibits stronger replication efficiency in vivo and in vitro (Zhou et al., 2009). The PRRSV Lena, a strain belonging to the subtype 3 of genotype 1, shares an increased virulence and higher replication efficiency (Karniychuk et al., 2010). Li et al. analyzed the replication efficiency differences of the chimeric viruses (RvJHn9n10 and RvHJn9n10) generated by swapping nsp9-and nsp10-coding region of LP-PRRSV HB-1/ 3.9 and HP-PRRSV JXwn06 and determined that the nsp9-and nsp10coding regions together are closely related to the replication efficiency of the Chinese HP-PRRSV in vitro and in vivo (Li et al., 2014c). Moreover, the changes of replication efficiency were more significant in PAMs than in MARC-145 cells. Thus, we performed the mutation in nsp9 and nsp10 using the chimeric full-length infectious clones of RvJHn9n10 and RvHJn9n10 based on the differential amino acid sites between HP-PRRSV JXwn06 and LP-PRRSV HB-1/3.9. Considering that PAMs are the target cells of PRRSV, we evaluated the growth kinetics of all the rescued mutant viruses in PAMs instead of MARC-145 cells in the L. Xu et al. Virology 517 (2018) 135-147 present study. Our analyses showed that the mutation of Thr to Ala in 586 and Ser to Thr in 592 of RvHJn9n10 nsp9 led to the decreased replication efficiency of the virus, while the mutation of respective amino acid sites of RvJHn9n10 could not result in the increased replication efficiency of the virus. Then we continued to conduct the mutation of these two sites using the infectious cDNA clones of RvJXwn06 or RvJHn10, and confirmed the influences of these two residues mutation for replication efficiency. The data indicated that the respective mutations decreased the replication efficiency and viral RNA synthesis efficiency of the viruses based on RvJXwn and RvJHn10.
These results together suggest that the residues in 586 and 592 of nsp9 are the key sites contributing to the replication efficiency of the Chinese HP-PRRSV, but are insufficient to determine the replication efficiency of PRRSV in general. Animal trials data of serial mutant viruses based on the infectious cDNA clones of RvJXwn and RvJHn10 indicated that the mutation in 586 and 592 of nsp9 could lead to the decreased virulence of the viruses RvJXwn and RvJHn10 for piglets, while the mutation of these two residues in nsp9 based on the chimeric infectious clone of RvHJn10 could not enhance the pathogenicity of the virus RvHJn10 for piglets. In terms of our data, we suggest that the residues in 586 and 592 of nsp9 are critical sites natively conferring the fatal virulence of the Chinese HP-PRRSV for piglets.
Our present study further proved that the replication efficiency of PRRSV is closely related to its pathogenicity or virulence, and the nsp9 can be considered the key replicase directing the replication of PRRSV in host cells and in vivo. To date, the crystal structure of PRRSV nsp9 remains unknown. By using sequence alignments and homology modeling prediction with equine arterivirus (EAV) and coronavirus, the RdRp domain of PRRSV contains at least six canonical conserved motifs (A to F) (Beerens et al., 2007;Poch et al., 1989;te Velthuis, 2014). By amino acid alignment analysis of nsp9 between HP-PRRSV JXwn06 and LP-PRRSV HB-1/3.9, it is found that the differential amino acid in 427 is located between motif F and A, while other three differential amino acids is situated in motif E. Our results revealed that the Thr in 586 and Ser in 592 of HP-PRRSV nsp9 are selected to affect its replication efficiency and pathogenicity, but it is difficult to explain how these amino acids impact the feature of nsp9 due to unknown function of these motifs. What makes us more frustrated is that it is difficult to assess the RdRp activity of nsp9 because of its weak and controversial activity in vitro (Lehmann et al., 2016). Comparative analysis of nsp9 amino acid sequence showed that the amino acid in 586 is a conserved residue in all strains of HP-PRRSV, namely the amino acid in 586 of nsp9 is Thr in HP-PRRSV, whereas Ala in LP-PRRSV. However, the amino acid in 592 is Thr only in HB-1/3.9 and its wild-type HB-1(sh)/2002 (Gao et al., 2004), whereas Ser in other strains. In addition, our animal trails indicated that the amino acid in 586 impacted more obviously the virus pathogenicity than the amino acid in 592. Thus, it is proposed that the amino acid in 586 of nsp9 shares common feature affecting the replication efficiency and pathogenicity of PRRSV, while the effect of the amino acid in 592 on viral replication efficiency and pathogenicity might be strain-specific. Another concern is that, by mutating these two amino acid sites, the virulence of HP-PRRSV JXwn06 could be attenuated, whereas the replication efficiency and pathogenicity of LP-PRRSV HB-1/3.9 could not be enhanced reversely. These findings imply that other unknown factors may participate in this influence for the pathogenicity of PRRSV although the two amino acids in 586 and 592 of nsp9 are recognized as critical sites conferring the fatal virulence of HP-PRRSV in the current study. One example is that HP-PRRSV can be attenuated by de-optimization of codon pair bias in nsp9-coding region without any amino acid changes (Gao et al., 2015).
Many host cellular factors have been proved to interact with nsp9 and regulate the replication or RNA synthesis of PRRSV (Dong et al., 2014;Liu et al., 2016;Zhao et al., 2015). In another aspect, our results also suggested that the virulence of mutant viruses could be more remarkable decreased by the nsp10 of HB-1/3.9, although the replication efficiency changes were not observed. Previous studies have demonstrated that some nonstructural proteins of nidovirus can interact with each other, and indicated that the polymerase and helicase can interact with each other in vitro (von Brunn et al., 2007). We propose that there are two possible mechanisms of affecting virus replication efficiency and pathogenicity by mutating these two amino acids: they might directly affect the normal function of viral RNA amplification, or impact the interaction between nsp9 and host factors or other viral proteins. These are required to be further investigated in the future work.
As a whole, our findings indicate that the amino acids in 586 and 592 of nsp9 contribute to the replication efficiency of the Chinese HP-PRRSV in PAMs, and are critical sites natively in determining its fatal virulence for piglets, providing essential evidence for understanding the molecular basis in relation to the virulence of the Chinese HP-PRRSV for pigs. of Science and Technology (2014CB542700), and the earmarked fund for China Agriculture Research System from the Chinese Ministry of Agriculture (CARS-35).