key: cord-0719333-17bjnhrp
authors: Xiong, Junyao; Cui, Xingyang; Zhao, Kuan; Wang, Qian; Huang, Xinyi; Li, Dongyan; Yu, Fang; Yang, Yongbo; Liu, Di; Tian, Zhijun; Cai, Xuehui; An, Tongqing
title: A Novel Motif in the 3′-UTR of PRRSV-2 Is Critical for Viral Multiplication and Contributes to Enhanced Replication Ability of Highly Pathogenic or L1 PRRSV
date: 2022-01-18
journal: Viruses
DOI: 10.3390/v14020166
sha: 2a87b929161baed4795f0248cec536c59e9dd75a
doc_id: 719333
cord_uid: 17bjnhrp

Highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) with enhanced replication capability emerged in China and has become dominant epidemic strain since 2006. Up to now, the replication-regulated genes of PRRSV have not been fully clarified. Here, by swapping the genes or elements between HP-PRRSV and classical PRRSV based on infectious clones, NSP1, NSP2, NSP7, NSP9 and 3′-UTR are found to contribute to the high replication efficiency of HP-PRRSV. Further study revealed that mutations at positions 117th or 119th in the 3′-UTR are significantly related to replication efficiency, and the nucleotide at position 120th is critical for viral rescue. The motif composed by 117–120th nucleotides was quite conservative within each lineage of PRRSV; mutations in the motif of HP-PRRSV and currently epidemic lineage 1 (L1) PRRSV showed higher synthesis ability of viral negative genomic RNA, suggesting that those mutations were beneficial for viral replication. RNA structure analysis revealed that this motif maybe involved into a pseudoknot in the 3′-UTR. The results discovered a novel motif, 117–120th nucleotide in the 3′-UTR, that is critical for replication of PRRSV-2, and mutations in the motif contribute to the enhanced replicative ability of HP-PRRSV or L1 PRRSV. Our findings will help to understand the molecular basis of PRRSV replication and find the potential factors resulting in an epidemic strain of PRRSV.

Porcine reproductive and respiratory syndrome (PRRS) is one of the most important infectious diseases in the pig industry, characterized by severe reproductive failure in sows and respiratory diseases in all-age pigs [1] . It was first reported in the United States in 1987 [2] , and it is currently found in all major pig-raising countries [3] [4] [5] [6] . The etiological agent, PRRS virus (PRRSV), is an enveloped, single-stranded and positive-sense RNA virus. The PRRSV genome encodes at least 11 open reading frames (ORFs), flanked by 5 -UTR (untranslated region, UTR) and 3 -UTR [7] [8] [9] . ORF1a and ORF1b encode polyproteins of PP1a and PP1ab, in which at least 16 nonstructural proteins (NSPs) are produced after hydrolysis and cleavage [10] [11] [12] [13] [14] [15] [16] [17] . Meanwhile, ORF2a, ORF2b, ORF3, ORF4, ORF5, ORF5a, ORF6 and ORF7 encode the structural proteins of GP2a, E, GP3, GP4, GP5, ORF5a, M and N, respectively [9, 18] .

PRRSV is a member of the Arteriviridae family in the Nidovirales order (including Coronaviridae family, Arteriviridae family and Roniviridae family). Viruses in these families have

The RNA secondary structure was predicted by the mFold website server for nucleic acid folding and hybridization prediction (http://www.unafold.org/, accessed on 11 November 2020) under default folding parameters (37 • C, 1 M NaCl, no divalent ions and no limit on the distance between paired bases) [37] . The predicted secondary structure was modified by the RNAviz 2.0 software (http://rnaviz.sourceforge.net/, accessed on 11 November 2020) [38] . According to previous studies, the predicted pseudoknot in the 3 -UTR was depicted. [22, 23, [39] [40] [41] [42] [43] .

Based on the infectious clones of HuN4 or CH-1a, a series of chimeric plasmids of pHC 1 , pCH 1 , pHC 2 , pCH 2 , pHC 3 , pCH 3 , pHC 4 , pCH 4 , pHC 5 , pCH 5 , pHC 7 , pCH 7 , pHC 9 , pCH 9 , pHC 10 , pCH 10 , pHC 11 , pCH 11 , pHC 12 , pCH 12 , pHC GP2 , pCH GP2 , pHC GP3 , pCH GP3 , pHC GP4 , pCH GP4 , pHC GP5 , pCH GP5 , pHC M , pCH M , pHC 3 -UTR and pCH 3 -UTR were constructed, which swapped the corresponding regions of NSP1, NSP2, NSP3, NSP4, NSP5, NSP7, NSP9, NSP10, NSP11, NSP12, GP2, GP3, GP4, GP5, M and 3 -UTR, respectively. Moreover, in order to find the key nucleotides of 3 -UTR, a series of mutant plasmids of pHuN4-3 -UTR-C9T, pHuN4-3 -UTR-17+G, pHuN4-3 -UTR-A117G, pHuN4-3 -UTR-A119G, pHuN4-3 -UTR-G120A, pHuN4-3 -UTR-117+120, pHuN4-3 -UTR-119+120 and pHuN4-3 -UTR-117+119+120 was constructed by site-directed mutagenesis. In brief, all mutant plasmids were introduced into HuN4 infectious clone by site-directed PCR-based mutagenesis. Firstly, PCR products were amplified by using different mutant primers with the template of HuN4. Secondly, the varied PCR products of the same mutation were fused by fusion PCR. Then, HuN4 plasmid and mutant DNA fragments were digested with the same restriction enzyme of Mlu I and Pac I. Finally, the HuN4 vector and the mutant DNA fragments were ligated by DNA ligase, and the desired mutant plasmids were constructed successfully. All plasmids were subject to nucleotide sequencing. All primers were listed in Table 1 . Table 1 . Primers used in this study.

For viral rescue, Marc-145 cells with 80% confluence seeded in 6-well plates were transfected with the chimeric, mutant or backbone plasmids using X-tremeGENE HP DNA transfection reagent (Roche, Switzerland) according to the manufacturer's instruction. The rescued chimeric viruses were individually designated HC 1 , CH 1 , HC 2 , CH 2 , HC 3 , CH 3 , HC 4 , CH 4 , HC 5 , CH 5 , HC 7 , CH 7 , HC 9 , CH 9 , HC 10 , CH 10 , HC 11 , CH 11 , HC 12 , CH 12 , HC GP2 , , HC GP3 , CH GP3 , HC GP4 , CH GP4 , HC GP5 , CH GP5 , HC M , CH M , HC 3 -UTR and CH 3 -UTR ,  respectively. The rescued 3 -UTR mutant viruses were individually designated HuN4-3 -UTR-C9T, HuN4-3 -UTR-17+G, HuN4-3 -UTR-A117G and HuN4-3 -UTR-A119G , respectively. The supernatants of the transfected cells were harvested at 96 h post-transfection, designated as passage 1 (P1) of the rescued viruses. P1 was continuously passaged to the third generation (P3), followed by viral titration on Marc-145 cells, and the viruses were stored at −80 • C for later use. All viruses (P3) were subject to nucleotide sequencing.

For viral titration, confluent Marc-145 cells seeded in 96-well plates were incubated with 10-fold serially diluted viral suspensions, and the viral titers were measured according to the method of Reed and Muench [44] . Briefly, the confluent Marc-145 cells seeded in 96well plates were incubated with 10-fold serially diluted viral suspensions. After absorption for 1 h at 37 • C, the supernatants were removed, and the fresh DMEM with 2% FBS was added. Following this, the plates were incubated for an additional six days, and viral titers were determined by the presence of visible cytopathic effect (CPE). Finally, the viral titers were calculated according to the method of Reed and Muench.

In order to investigate the replication difference of rescued mutant viruses and their parental viruses (rescued HuN4 or CH-1a), confluent Marc-145 cells or PAMs were incubated with the mutant and parental viruses at a multiplicity of infection (MOI) of 0.01. After inoculation for 1 h at 37 • C, the supernatants were discarded, and cells were washed with PBS three times. Following this, fresh DMEM with 2% FBS was added. At last, cell supernatants were harvested at 0, 24, 48, 72 and 96 h post-infection (hpi). The viral titers of different time points were measured by a microtitration assay using Marc-145 cells in 96-well plates and calculated as 50% tissue culture infective doses (TCID 50 ) per milliliter according to the method of Reed and Muench. Each time point was independently repeated three times.

In order to further characterize chimeric or mutant viruses, Marc-145 cells were infected by those viruses and immobilized at 96 hpi; then, they were confirmed with IFA using a monoclonal antibody directed against the PRRSV N protein as described previously [45] .

Parental and mutant viruses were infected with Marc-145 cells in 6-well plates. After seven days, the viral plaques were observed after gentian violet staining. In brief, Marc-145 cells were infected with viruses at a dose of MOI = 0.01. After incubating for 1 h, the supernatants were discarded and washed three times with phosphate-buffered saline (PBS). Then, the pre-prepared DMEM medium containing 2% FBS and 1% low-melting agarose was added to the cells. Subsequently, plates were incubated at 37 • C for 5 to 7 days. Eventually, plaque morphology was observed after paraformaldehyde fixation and crystal violet staining. Ten homologous plaques were selected randomly (avoid selecting large plaques formed by fusion of plaques), for which their diameters were measured by a ruler, and finally analyzed by Graphpad Prism 6 (Graphpad, San Diego, CA, USA).

3 -UTR has several mutant patterns in different PRRSV lineages or branches. In order to detect the transcription activity of 3 -UTR in different mutant patterns, a PRRSV minigenome system harboring luciferase reporter gene was employed according to a previous study [46] . The PRRSV mini-genome system was a kind gift from Dr. Yan-Dong Tang at the Harbin Veterinary Research Institute. Briefly, 3 -UTR was site-directed mutated according to different lineages and replaced the original 3 -UTR in the PRRSV mini-genome system. Followed by DNA sequencing, the mini-genome system plasmid (0.6 µg) and

Renilla luciferase plasmid (0.02 µg) (Promega, Madison, WI, USA) were transfected into Marc-145 cells in 96-well plates. Renilla luciferase activity was used as an internal control for the normalization of luciferase values obtained from cells transfected with the firefly luciferase mini-genome system. At 24 h post-transfection, each well with Marc-145 cells was infected with PRRSV HuN4 at a dose of MOI = 0.1. Twenty-four hours later, the cells were lysed and assayed for firefly luciferase activity using a dual-luciferase reporter assay system (Promega, USA) according to the manufacturer's instruction.

All experiments were performed at least three independent replicates. All data were analyzed by Graphpad Prism 6. The measured values were expressed as the mean with standard deviation (SD). If no virus was detected by plaque assay, the number representing the LOD (limit of detection) was used. Differences were analyzed for statistical significance using two-tailed unpaired t test for two groups or multiple comparison one-way variance (ANOVA) for more than two groups. Differences were considered statistically significant at a value of p < 0.05.

In order to explore the key replication-related viral genes or elements, a series of chimeric viruses by swapping the corresponding genes or elements of HP-PRRSV HuN4 and classical PRRSV CH-1a was constructed based on full-length infectious clones ( Figure 1 ). The replication efficiency of chimeric viruses was compared with their parental viruses in porcine alveolar macrophages (PAMs). The results showed that HC 1 , HC 2 , HC 7 , HC 9 and HC 3 -UTR had significantly lower viral titers than parental HuN4 ( Figure 2A ). Conversely, CH 7 and CH 9 had markedly higher viral titers than the parental CH-1a ( Figure 2B ). The replication efficiency of other chimeric viruses was not different from their parental viruses ( Figure S1 ). and CH-1a and HC9 and CH9, respectively. (φ, p < 0.05; φφ, p < 0.01; φφφ, p < 0.001). Delta (δ) indicates a significant difference between HuN4 and CH-1a and HC3′-UTR and CH3′-UTR, respectively (δ, p < 0.05; δδδ, p < 0.001).

The alignment of 3′-UTR between HuN4 and CH-1a showed that there was a 1-nt deletion at position 17th, which was deleted in all HP-PRRSV strains. In addition, there were four mutations (at positions 9, 117, 119, and 120, respectively) in the 3′-UTR of the HP-PRRSV ( Figure 3A ). In order to investigate the effect of deletion or mutations in the 3′-UTR on viral replication, eight mutant plasmids of pHuN4-3′-UTR-17+G, pHuN4-3′-UTR-C9T, pHuN4-3′-UTR-A117G, pHuN4-3′-UTR-A119G, pHuN4-3′-UTR-G120A, pHuN4-3′-UTR-A117G-G120A, pHuN4-3′-UTR-A119G-G120A and pHuN4-3′-UTR-A117G-A119G-G120A were constructed based on HuN4 infectious clones by replacing the corresponding nucleotide of CH-1a. Only four of them, HuN4-3′-UTR-C9T, HuN4-3′- The data were presented as the mean standard deviation of three independent experiments. LOD: limit of detection. Asterisk (*) indicates a significant difference between HuN4 and CH-1a and HC 1 and CH 1 , respectively. (***, p < 0.001). Pound (#) indicates a significant difference between HuN4 and CH-1a and HC 2 and CH 2 , respectively. (###, p < 0.001). And (&) indicates a significant difference between HuN4 and CH-1a and HC 7 and CH 7 , respectively. (&, p < 0.05; &&&, p < 0.001). Phi (ϕ) indicates a significant difference between HuN4 and CH-1a and HC 9 and CH 9 , respectively. (ϕ, p < 0.05; ϕϕ, p < 0.01; ϕϕϕ, p < 0.001). Delta (δ) indicates a significant difference between HuN4 and CH-1a and HC 3 -UTR and CH 3 -UTR , respectively (δ, p < 0.05; δδδ, p < 0.001).

The alignment of 3 -UTR between HuN4 and CH-1a showed that there was a 1-nt deletion at position 17th, which was deleted in all HP-PRRSV strains. In addition, there were four mutations (at positions 9, 117, 119, and 120, respectively) in the 3 -UTR of the HP-PRRSV ( Figure 3A ). In order to investigate the effect of deletion or mutations in the 3 -UTR on viral replication, eight mutant plasmids of pHuN4-3 -UTR-17+G, pHuN4-3 -UTR-C9T, pHuN4-3 -UTR-A117G, pHuN4-3 -UTR-A119G, pHuN4-3 -UTR-G120A, pHuN4-3 -UTR-A117G-G120A, pHuN4-3 -UTR-A119G-G120A and pHuN4-3 -UTR-A117G-A119G-G120A were constructed based on HuN4 infectious clones by replacing the corresponding nucleotide of CH-1a. Only four of them, HuN4-3 -UTR-C9T, HuN4-3 -UTR-17+G, HuN4-3 -UTR-A117G and HuN4-3 -UTR-A119G, were rescued successfully ( Figure S2 ). It may be surprising that a single-point-mutation of 120th or multiple-point-mutation of 120th with other nucleotides were not rescued, indicating the 120th in the 3 -UTR is critical for rescuing HP-PRRSV.

3 -UTR secondary structures of HuN4-3 -UTR-C9T, HuN4-3 -UTR-17+G, HuN4-3 -UTR-A117G and HuN4-3 -UTR-A119G were similar to that of parental HuN4, suggesting that the deletion or mutations at positions 9, 17, 117 and 119 in the 3 -UTR did not change the overall conformation of the 3 -UTR secondary structure ( Figure 3B ). However, the 3 -UTR secondary structures of HuN4-3 -UTR-G120A, HuN4-3 -UTR-A117G-G120A, HuN4-3 -UTR-A119G-G120A and HuN4-3 -UTR-A117G-A119G-G120A showed obvious differences with that of parental HuN4 ( Figure 3C ), suggesting the 120th plays a critical role in maintaining the steady conformation of 3 -UTR.

After replacement with the nucleotides of CH-1a strain, the replicative level of four mutant viruses was determined and compared with the HuN4 strain in Marc-145 cells. The growth kinetics showed that HuN4-3 -UTR-17+G or HuN4-3 -UTR-C9T displayed a slight decrease ( Figure 4A ). As for HuN4-3 -UTR-A117G and HuN4-3 -UTR-A119G, their titers were remarkably lower than HuN4 at 48, 72 and 96 hpi (p < 0.001) ( Figure 4B ). UTR-A117G and HuN4-3′-UTR-A119G were similar to that of parental HuN4, suggesting that the deletion or mutations at positions 9, 17, 117 and 119 in the 3′-UTR did not change the overall conformation of the 3′-UTR secondary structure ( Figure 3B ). However, the 3′-UTR secondary structures of HuN4-3′-UTR-G120A, HuN4-3′-UTR-A117G-G120A, HuN4-3′-UTR-A119G-G120A and HuN4-3′-UTR-A117G-A119G-G120A showed obvious differences with that of parental HuN4 (Figure 3C ), suggesting the 120th plays a critical role in maintaining the steady conformation of 3′-UTR. 

After replacement with the nucleotides of CH-1a strain, the replicative level of four mutant viruses was determined and compared with the HuN4 strain in Marc-145 cells. The growth kinetics showed that HuN4-3′-UTR-17+G or HuN4-3′-UTR-C9T displayed a slight decrease ( Figure 4A ). As for HuN4-3′-UTR-A117G and HuN4-3′-UTR-A119G, their titers were remarkably lower than HuN4 at 48, 72 and 96 hpi (p < 0.001) ( Figure 4B ). 

The average diameters of plaques caused by HuN4-3 -UTR-C9T, HuN4-3 -UTR-17+G, HuN4-3 -UTR-A117G, HuN4-3 -UTR-A119G and HuN4 were 0.305 cm, 0.285 cm, 0.140 cm, 0.075 cm and 0.320 cm, respectively ( Figure 4C ). The plaques of HuN4-3 -UTR-C9T and HuN4-3 -UTR-17+G were slightly smaller than that of parental HuN4 ( Figure 4D ). In contrast, the average diameters of plaques caused by HuN4-3 -UTR-A117G were about 33% of the parental HuN4, and HuN4-3 -UTR-A119G was approximately 25% of the parental HuN4 ( Figure 4D ), suggesting that both mutations at position 117th or 119th closely related to the replication level of HP-PRRSV HuN4 strain.

Phylogenetic analysis of all the available full-length PRRSV-2 (n = 765) strains from 1992 to 2019 in GenBank showed that those PRRSV-2 were classified into eight lineages ( Figure 5A ). The results also showed that L8 and L1 were the two most abundant PRRSV-2 strains, both in the world and in China. In the 765 PRRSV-2 strains, the number of strains was as follows: L8 (n = 411, 53.73%), L1 (n = 243, 31.76%, 208 out of them were detected in 2013-2019), L5 (n = 79, 10.33%), L3 (n = 21, 2.75%), L9 (n = 6, 0.78%), L7 (n = 3, 0.39%), L4 (n = 1, 0.13%) and L6 (n = 1, 0.13%). It also showed that the full-length PRRSV sequences in China contained four different lineages: L8 (n = 379, 74.31%), L1 (n = 90, 17.65%), L3 (n = 21, 4.12%) and L5 (n = 20, 3.92%) ( Table S1 ). The alignment of 3 -UTR showed that the 117-120th motif was quite conservative within each lineage ( Figure S3 ). The motif in all L9 PRRSVs was AAAG; L8 was AAAG (in HP-PRRSVs) or GAGA (in classical PRRSVs); L7 was AAGA; L6 was AAAA; L5 was GAGA; L4 was AGGA; and L3 was AAGA or AAAA ( Table 2) . No L2 full-length sequence has been made available in GenBank yet. In L1 PRRSV, a newly emerging strain that has recently been spreading rapidly, the motif was complex and diverse, such as AACA (in NADC30-like viruses), AGAA (in NADC31-like viruses) and AATA (in MN184A-or NADC34-like viruses) ( Figure 5B ). 

The average diameters of plaques caused by HuN4-3′-UTR-C9T, HuN4-3′-UTR-17+G, HuN4-3′-UTR-A117G, HuN4-3′-UTR-A119G and HuN4 were 0.305 cm, 0.285 cm, 0.140 cm, 0.075 cm and 0.320 cm, respectively ( Figure 4C ). The plaques of HuN4-3′-UTR-C9T and HuN4-3′-UTR-17+G were slightly smaller than that of parental HuN4 ( Figure  4D ). In contrast, the average diameters of plaques caused by HuN4-3′-UTR-A117G were about 33% of the parental HuN4, and HuN4-3′-UTR-A119G was approximately 25% of the parental HuN4 ( Figure 4D ), suggesting that both mutations at position 117th or 119th closely related to the replication level of HP-PRRSV HuN4 strain.

Phylogenetic analysis of all the available full-length PRRSV-2 (n = 765) strains from 1992 to 2019 in GenBank showed that those PRRSV-2 were classified into eight lineages ( Figure 5A ). The results also showed that L8 and L1 were the two most abundant PRRSV-2 strains, both in the world and in China. In the 765 PRRSV-2 strains, the number of strains was as follows: L8 (n = 411, 53.73%), L1 (n = 243, 31.76%, 208 out of them were detected in 2013-2019), L5 (n = 79, 10.33%), L3 (n = 21, 2.75%), L9 (n = 6, 0.78%), L7 (n = 3, 0.39%), L4 (n = 1, 0.13%) and L6 (n = 1, 0.13%). It also showed that the full-length PRRSV sequences in China contained four different lineages: L8 (n = 379, 74.31%), L1 (n = 90, 17.65%), L3 (n = 21, 4.12%) and L5 (n = 20, 3.92%) ( Table S1 ). The alignment of 3′-UTR showed that the 117-120th motif was quite conservative within each lineage ( Figure S3 ). The motif in all L9 PRRSVs was AAAG; L8 was AAAG (in HP-PRRSVs) or GAGA (in classical PRRSVs); L7 was AAGA; L6 was AAAA; L5 was GAGA; L4 was AGGA; and L3 was AAGA or AAAA ( Table 2) . No L2 full-length sequence has been made available in GenBank yet. In L1 PRRSV, a newly emerging strain that has recently been spreading rapidly, the motif was complex and diverse, such as AACA (in NADC30-like viruses), AGAA (in NADC31-like viruses) and AATA (in MN184A-or NADC34-like viruses) ( Figure 5B ). CH-1a) . 

In order to investigate whether the mutations in the motif of L1 were related to viral replication, a PRRSV mini-genome system harboring the 3 -UTR and double luciferase report system was used to evaluate the synthesis level of viral negative genomic RNA, which is essential for genome replication and synthesis. A series of site-directed mutations was performed according to the nucleotides in the motif of different PRRSV lineages ( Figure 6A ). Firefly luciferase activity was normalized with respect to a co-transfected plasmid encoding Renilla luciferase. The luciferase reporter assay showed that the 3 -UTR of HP-PRRSV HuN4 was 1.2 times than that of classical PRRSV CH-1a, with no significant difference ( Figure 6B ). All different 3 -UTR of L1 PRRSV exhibited higher replication efficiency than that of classical PRRSV (p < 0.05) ( Figure 6B ), suggesting the RNA synthesis level of L1 PRRSV is higher than that of classical PRRSV. was performed according to the nucleotides in the motif of different PRRSV lineages (Figure 6A ). Firefly luciferase activity was normalized with respect to a co-transfected plasmid encoding Renilla luciferase. The luciferase reporter assay showed that the 3′-UTR of HP-PRRSV HuN4 was 1.2 times than that of classical PRRSV CH-1a, with no significant difference ( Figure 6B ). All different 3′-UTR of L1 PRRSV exhibited higher replication efficiency than that of classical PRRSV (p < 0.05) ( Figure 6B ), suggesting the RNA synthesis level of L1 PRRSV is higher than that of classical PRRSV. 

3′-UTR secondary structures of representative PRRSVs were predicted in the mFold software ( Figure 7A) . Meanwhile, the 117-120th motif was located into a pseudoknot (Figure 7B ) in the 3′-UTR, according to previous studies [23, 43] . The rescued A117G or A119G mutant viruses showed obviously lower replication ability than that of HuN4, while 

3 -UTR secondary structures of representative PRRSVs were predicted in the mFold software ( Figure 7A) . Meanwhile, the 117-120th motif was located into a pseudoknot ( Figure 7B ) in the 3 -UTR, according to previous studies [23, 43] . The rescued A117G or A119G mutant viruses showed obviously lower replication ability than that of HuN4, while G120A mutant viruses were not successfully rescued, indicating that the 117-120th motif in the 3 -UTR is critical for viral replication or rescue. The mutations at positions 117 and 119 did not change the overall conformation of 3 -UTR secondary structure, but the mutation at position 120 changed the overall secondary structure of 3 -UTR ( Figure 3B,C) . Nucleotide mutations probably changed the minimum free energy of 3 -UTR secondary structure, which in turn results in alterations in the structure or stability of the pseudoknot. 

PRRS has caused severe economic losses to the pig industry since it was reported in 1987. The first PRRSV-2 strain, VR-2332 (belonging to L5), was isolated in the US in 1992. In China, the first PRRSV-2 strain, CH-1a (classical PRRSV, belonging to L8), was isolated in 1996. With the gradual evolution of PRRSV, HP-PPRSV (belonging to L8) emerged and 

PRRS has caused severe economic losses to the pig industry since it was reported in 1987. The first PRRSV-2 strain, VR-2332 (belonging to L5), was isolated in the US in 1992. In China, the first PRRSV-2 strain, CH-1a (classical PRRSV, belonging to L8), was isolated in 1996. With the gradual evolution of PRRSV, HP-PPRSV (belonging to L8) emerged and spread rapidly, and they quickly became the dominant epidemic strains [5] . HP-PRRSV was well documented, and it has higher replicative ability than classical PRRSV [32, 33] , which helps in understanding why HP-PRRSV spread rapidly. However, the key replicationregulated genes or elements still have not been elucidated.

Our study showed that NSP1, NSP2, NSP7, NSP9 and 3 -UTR are closely related to the replication efficiency of HP-PRRSV. Moreover, previous studies demonstrate that NSP1 has an effect on the synthesis of TNF-α and IFN-β [47, 48] ; NSP2, NSP3 and NSP5 interact with each other and participate in RTC formation and viral replication [19, 49] ; NSP7 interacts with NSP9, and NSP9 contributes to the replication ability and pathogenicity of PRRSV [33, 50] . Meanwhile, the NSP2, NSP3, NSP5, NSP9 and NSP10 were reported to assemble in viral replication and transcription complex (RTC); then, RTC interacts with 3 -UTR, which initiates the transcription of sgmRNA and the replication of viral genome [19, 20] . However, as a critical element of RTC, the function and key regions or sites of 3 -UTR in viral replication are still unclear. Therefore, a follow-up study on 3 -UTR was performed. By swapping the corresponding genes between CH-1a and HuN4, we demonstrated that the 3 -UTR contributes to the enhanced replicative ability of HP-PRRSV. Due to the low replication ability of CH-1a, the chimeric or mutant viruses based on the backbone of CH-1a were not replicated well in Marc-145 or PAM cells. Thus, the site-direct mutations of 3 -UTR were based on the HuN4 strain. Growth kinetics and plaque morphology assays show that the nucleotides at positions 117th and 119th in the 3 -UTR significantly affect viral replication or plaque sizes, suggesting that two nucleotides are critical sites for viral replication. Interestingly, the 120th nucleotide in the 3 -UTR is important for viral rescue, which may be because the 120th nucleotide and those nearby are involved in the formation of a senior RNA structure. Therefore, the 120th and its nearby nucleotides (117th to 119th) consist of a novel motif, 117-120th, which is a key motif related to PRRSV replication.

Alignment of all available entire sequences of PRRSV-2 in 1992-2019 revealed that the 117-120th nucleotides in the 3 -UTR were conservative in individual lineages. For each lineage, there were more or less mutations in the motif when compared with the prototype PRRSV-2, VR-2332 strain. The mutations contribute to the higher replication efficiency of HP-PRRSV HuN4 compared with the CH-1a strain. A highly efficient replication in host cells is one of the important features of dominant epidemic strain. There are many factors that result in HP-PRRSV becoming the dominant epidemic strain, such as increased viral tropism, high level of viral replication, moderate pathogenicity and so on [33, 50, 51] .

In China, L1 PRRSV emerged in 2013, then it spread quickly in main pig-raising provinces since 2015, and L1 PRRSVs mainly include NADC30-like PRRSV and NADC34like PRRSV [6, 31] . At present, HP-PRRSV and L1 PRRSVs are the dominant epidemic strains, and the proportion of L1 PRRSVs is higher than that of HP-PRRSV in some provinces [5, 6, 30] . There is a trend that the L1 PRRSVs surpass HP-PRRSV to become the most dominant epidemic strain in the next few years. In the present study, we found that mutations in the motif of 3 -UTR enhanced the replication of HP-PRRSV or L1 PRRSVs. At present, it is suitable for the most prevalent PRRSV lineages (L1/L3/L4/L5/L6/L7/L8/L9, 96.20%), except the FJ1402-like strain with deletion at positions 118-120th (L1, 3.80%). FJ1402-like strains with deletion at positions 118-120th in 3 -UTR have increased in recent years, and the detailed relationship between deletion and viral replication remains to be further explored.

Pseudoknot in the 3 -UTR of PRRSV can serve as pathogen-associated molecular patterns (PAMPs) and bind to the host cell pattern recognition receptors (PRRs) TLR3 and RIG-I to activate innate immune signaling and produced IFNs [43] . A previous study revealed a hairpin structure in the ORF7 together with the 3 -UTR form a kissing structure, which is essential for PRRSV replication [42] . In Mengovirus, three stem-loop (SL) are formed in its 3 -UTR, in which SL1 is a non-essential region for viral replication, and the absence of SL2 reduces the viral replication capability, and SL3 is necessary for viral replication [52] . In Dengue virus, SLs (A2 and A3) in the D2 region of 3 -UTR are essential for viral replication [53] . Consequently, we speculate that the 117-120th motif plays an important role in PRRSV replication. A pseudoknot is formed in the 3 -UTR according to the previous RNA structure analysis method [23] , and the 117-120th motif was located in the pseudoknot.

In summary, our study has revealed that (i) NSP1, NSP2, NSP7, NSP9 and 3 -UTR are closely related to replication efficiency of PRRSV-2; (ii) the nucleotides at positions 117th and 119th in the 3 -UTR are critical for viral replication; (iii) the 117-120th is an important motif in the 3 -UTR of PRRSV-2, and mutations in the motif contribute to the higher replication level of HP-PRRSV or L1 PRRSVs. Our findings discovered the key motif or nucleotide involved in viral replication and would help to understand the mechanism of viral replication. The motif is potentially related to the formation of dominant epidemic of HP-PRRSV strains or L1 PRRSVs, which will require further investigation.

Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/v14020166/s1, Figure S1 : The comparison of replication efficiency between parental and chimeric viruses in PAMs; Figure S2 : Identification of mutant viruses; Figure S3 : The 3 -UTR alignment of 765 PRRSV-2 strains. Table S1 

Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States

Epidemiology of porcine reproductive and respiratory syndrome (PRRS): An overview

Mystery swine disease in the Netherlands: The isolation of Lelystad virus

Lelystad virus, the cause of porcine epidemic abortion and respiratory syndrome: A review of mystery swine disease research at Lelystad

Genotypic and geographical distribution of porcine reproductive and respiratory syndrome viruses in mainland China in 1996-2016

Phylogenetics, genomic recombination, and NSP2 polymorphic patterns of porcine reproductive and respiratory syndrome virus in China and the United States in 2014-2018

Molecular characterization of porcine reproductive and respiratory syndrome virus, a member of the arterivirus group

A 10-kDa structural protein of porcine reproductive and respiratory syndrome virus encoded by ORF2b

Novel structural protein in porcine reproductive and respiratory syndrome virus encoded by an alternative ORF5 present in all arteriviruses

Non-structural proteins 2 and 3 interact to modify host cell membranes during the formation of the arterivirus replication complex

Heterogeneity in Nsp2 of European-like porcine reproductive and respiratory syndrome viruses isolated in the United States

Processing and evolution of the N-terminal region of the arterivirus replicase ORF1a protein: Identification of two papainlike cysteine proteases

Processing of the equine arteritis virus replicase ORF1b protein: Identification of cleavage products containing the putative viral polymerase and helicase domains

Efficient −2 frameshifting by mammalian ribosomes to synthesize an additional arterivirus protein

Proteolytic processing of the replicase ORF1a protein of equine arteritis virus

The Nsp12-coding region of type 2 PRRSV is required for viral subgenomic mRNA synthesis

Alternative proteolytic processing of the arterivirus replicase ORF1a polyprotein: Evidence that NSP2 acts as a cofactor for the NSP4 serine protease

Full-length sequence of a Canadian porcine reproductive and respiratory syndrome virus (PRRSV) isolate

The network of interactions among porcine reproductive and respiratory syndrome virus non-structural proteins

The PRRSV replicase: Exploring the multifunctionality of an intriguing set of nonstructural proteins

Inhibition of porcine reproductive and respiratory syndrome virus replication with exosome-transferred artificial microRNA targeting the 3 untranslated region

Characterization of the RNA components of a putative molecular switch in the 3 untranslated region of the murine coronavirus genome

An RNA pseudoknot in the 3 end of the arterivirus genome has a critical role in regulating viral RNA synthesis

The S gene is necessary but not sufficient for the virulence of porcine epidemic diarrhea virus novel variant strain BJ2011C

RNA recombination at Chikungunya virus 3 UTR as an evolutionary mechanism that provides adaptability

Porcine reproductive and respiratory syndrome virus comparison: Divergent evolution on two continents

Porcine reproductive and respiratory syndrome virus was isolated from abortive fetus of suspected PRRS

Origin of highly pathogenic porcine reproductive and respiratory syndrome virus

Highly pathogenic porcine reproductive and respiratory syndrome virus

Pathogenesis and control of the Chinese highly pathogenic porcine reproductive and respiratory syndrome virus

Importation and recombination are responsible for the latest emergence of highly pathogenic porcine reproductive and respiratory syndrome virus in China

Quantitative estimation of the replication kinetics of genotype 2 PRRSV strains with different levels of virulence in vitro

Nsp9 and Nsp10 contribute to the fatal virulence of highly pathogenic porcine reproductive and respiratory syndrome virus emerging in China

MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets

Phylogenybased evolutionary, demographical, and geographical dissection of North American type 2 porcine reproductive and respiratory syndrome viruses

Comparative genomic analysis of five pairs of virulent parental/attenuated vaccine strains of PRRSV

Mfold web server for nucleic acid folding and hybridization prediction

RnaViz 2: An improved representation of RNA secondary structure

A conformational switch at the 3 end of a plant virus RNA regulates viral replication

A phylogenetically conserved hairpin-type 3 untranslated region pseudoknot functions in coronavirus RNA replication

Characterization of an essential RNA secondary structure in the 3 untranslated region of the murine coronavirus genome

Kissing interaction between 3 noncoding and coding sequences is essential for porcine arterivirus RNA replication

Identification of the RNA pseudoknot within the 3 end of the porcine reproductive and respiratory syndrome virus genome as a pathogen-associated molecular pattern to activate antiviral signaling via RIG-I and toll-like receptor 3

A simple method of estimating fifty percent endpoints

Conserved nucleotides in the terminus of the 3 UTR region are important for the replication and infectivity of porcine reproductive and respiratory syndrome virus

Open reading frames 1a and 1b of the porcine reproductive and respiratory syndrome virus (PRRSV) collaboratively initiate viral minus-strand RNA synthesis

Identification of two autocleavage products of nonstructural protein 1 (nsp1) in porcine reproductive and respiratory syndrome virus infected cells: Nsp1 function as interferon antagonist

Porcine reproductive and respiratory syndrome virus non-structural protein 1 suppresses tumor necrosis factor-α promoter activation by inhibiting NF-kB and Sp1

Structural analysis of porcine reproductive and respiratory syndrome virus non-structural protein 7α (NSP7α) and identification of its interaction with NSP9

Two residues in NSP9 contribute to the enhanced replication and pathogenicity of highly pathogenic porcine reproductive and respiratory syndrome virus

Adaptions of field PRRSVs in Marc-145 cells were determined by variations in the minor envelope proteins GP2a-GP3

Phenotypic characterization of three phylogenetically conserved stem-loop motifs in the Mengovirus 3 untranslated region

Long-range RNA-RNA interactions circularize the dengue virus genome

Institutional Review Board Statement: Not applicable.

Data Availability Statement: All available data are presented in the article.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.