Analysis of the N Protein Sequence Variability in 13 Isolated PRRSV Strains from China ()
1. Introduction
Caused by the porcine reproductive and respiratory syndrome virus (PRRSV), porcine reproductive and respiratory syndrome (PRRS) is one of the major infectious swine diseases. This disease was reported to have initially emerged in the United States [1]; thereafter, it appeared in many countries, such as Canada, France, German, Britain, the Netherlands, Belgium, Spain, Japan, the Philippines, and Australia [2]. In 1996, Guo Baoqing et al. initially isolated the porcine reproductive and respiratory syndrome virus (PRRSV) from an aborted fetus and hence confirmed the prevalence of porcine reproductive and respiratory syndrome (PRRS) in China [3]. In 2006, an epidemic characterized by high fever, rubefaction and dyspnea emerged in swine farms in most areas of China. Later, this epidemic was confirmed to have been caused by a PRRSV variant strain and was denominated as highly pathogenic porcine reproductive and respiratory syndrome (HP-PRRS) because of its high morbidity [4]. The outbreak of HP-PRRS caused mass mortality of the swine and huge economic losses for swine farms. Considering that the economic loss caused by HP-PRRS is far greater than that of type A infectious disease, PRRS was classified as type B infectious disease by the Office International Des Epizooties (OIE) in 1992. In addition, in 2008, in the revised “first, second or third category of animal disease species list” in China, HP-PRRS was listed in the first category of animal diseases. Although HP-PRRS was effectively controlled with the successful development and application of inactivated and attenuated vaccines, due to an excessive fear of HP-PRRS, vaccines, especially attenuated vaccines, are abused. As Yang Hanchun, a professor of the China Agricultural University, says, “it is quite hard to find PRRS-negative swine farms in China.” This study shows that since the outbreak of HP-PRRS in 2006, domestic PRRSV isolates belong to highly pathogenic porcine reproductive and respiratory syndrome virus branches with a very high homology [5]. Whether these isolates are from the vaccine virus is unknown, but it can be speculated that the field of PRRSV may be facing a new round of selection pressure, so maintaining epidemiological tracking of the field PRRSV-isolated strains is quite significant for the prevention and control of the disease.
PRRSV belongs to the arterial inflammation virus division and is one of the Arteritis viruses. It is spherical or oval under the electron microscope with a particle diameter of 50 - 65 nm and has a cystic membrane structure with a superficial surface. PRRSV is a positive-sense, single-stranded RNA virus with no segment and cystic membrane. With a length of 15 kb, the PRRSV genome contains 10 open reading frames (ORF), ORF1A, ORF1B, ORF2, ORF2A, ORF3-ORF5, ORF5a, ORF6 and ORF7. Each reading code frame of the viral genome partially overlaps with the adjacent reading code frame [6]. The nucleocapsid protein is the only protein that makes up all the capsid of the virus. Encoded by ORF7, the PRRSV N protein is approximately 372 bp in length and approximately 15 kD in width, and without N-glycosylation, it is a nonglycosylated protein. The genetic evolution tree was established on the basis of the whole gene, and the genetic evolution tree based on the ORF7 gene is the most closely inherited; thus, ORF7 is considered to be a marker gene capable of genetic expression of PRRSV information [7]. Rich in antigenic epitopes and with the greatest expression in all PRRSV proteins, the N protein can induce the expression of the body IL-10, and thereby, it can enhance the immune evasion ability of the virus [8]. Compared with the classical PRRSV strains, HP-PRRSV causes more CD4 + CD25 + FOXP3 + T cell-specific response, which is closely related to mutations in the 46th amino acid of the HP-PRRSV N protein [9]. Based on the significance of the N protein in the epidemiological investigation of PRRSV and the immunological significance of antigenic epitopes in resisting the defense from the body, this study analyzed the sequence of the N protein and its variability in 13 strains of PRRSV isolated from fields in Shandong and Liaoning provinces in China. The results provide the basis for the epidemiological analysis of PRRSV and establish the basis for the functional study of N protein antigen epitopes, which will be beneficial for the prevention and control of PRRSV.
2. Materials and Methods
2.1. Isolation and Identification of Virus
All 13 strains of the virus were inoculated in Marc-145 cells in this laboratory (our experiment was approved by the Animal Ethics Committee of Zunyi Medical University); this method of isolates was described in the paper which was published in virus genes by our research group [10]. By Reverse Transcription-Polymerase Chain Reaction (RT-PCR) identification, the 13 strains were identified as PRRSV. Details of the strains are shown in Table 1.
2.2. Primer
Referring to a series of primers published by this project group, the synthesis of one pair of primers amplified the ORF7 gene. The primers used in this study are shown in Table 2.
2.3. Target Gene cDNA Acquisition
Viral RNA was extracted by the acid phenol-thiocyanate-chloroform method, and the target gene was amplified by RT-PCR. Then, the acquired target gene was connected with the PMD18-T carrier and transformed into JM109 sensory E. coli cells. The bacterium liquid identified by Polymerase Chain Reaction (PCR) was sent to the Dalian Bao Biological Company for sequencing.
Table 2. The sequence of primers in the research [10].
2.4. Sequence Analysis
With the help of DNASTAR analysis software, the ORF7 gene and its deduced amino acid sequence from each strain in this study were compared by the Jotun Hein method. The foreign reference strains in this study include the classical American strain VR-2332, the European classical strain Lelystad, the vaccine strain MLV, and the domestic reference strains are CH-1A, CH-1r, R98, JXA1 and 11 strains of PRRSV isolated by this laboratory in 2008.
3. Experimental Results
3.1. PCR Amplification of Target Gene cDNA
As shown in Figure 1, we amplified the ORF7 gene from 13 strains of PRRSV by RT-PCR, and all of which had a 660 bp or so band, which corresponds to the expected size.
3.2. Nucleotide Sequence Analysis of the ORF7 Gene
Sequencing results (all nucleotide sequences and amino acid sequences were submitted to GenBank (numbers: mg922935 - mg922947). To eliminate redundant sequences, the results of Table 3 were obtained by using DNASTAR software to compare the whole gene sequences of 13 strains of PRRSV isolates with
Figure 1. Results of amplification of the ORF7 gene in 13 isolates of PRRSV strains with RT-PCR. M: Marker DL2000; 1 - 13: Results of the amplification of the ORF7 gene from 13 isolates.
21 strains of American PRRSV N protein genes selected in GenBank. The sequencing results showed that the ORF7 gene from 13 strains of PRRSV did not have a base deletion or insertion in comparison with the traditional strains, and the gene mutations were in the partial sites. The homology of the ORF7 genes with the highly pathogenic PRRSV strains isolated from the domestic was 88.7% - 100% with the JXA1 strain, 87.1% - 91.7% with the FJFS strain, 89.3% - 97.1% with the NADC30 and JL58 strains, 89.2% - 96.5% with the classical strains CH-1a and CH-1R strains isolated from China, 89.5% - 94.4% with the vaccine MLV strain, 89.5% - 94.4% with the American standard VR-2332 strain, 89.2% - 94.4% with the classical R98 strain, 55.6% - 58.3% with the European classical Lelystad strain, and 88.2% - 100% with the 11 previously isolated PRRSV strains. The internal homology of 13 isolated strains isolated was 88.4% - 100%. DNASTAR software was used to draw the genetic evolution tree (Figure 2) of the ORF7 gene from the isolated strains. The figure shows that the strains of PRRSV isolated from this experiment belong to the American strain (red mark), and the isolates are divided into two branches. Among them, sdcl03, sdlq08, and sdsg10 and a pedigree represented by the NADC30 strain make up a branch, and the remaining 10 strains belong to a pedigree branch represented by the JXA1 strain.
3.3. Analysis of Amino Acid Sequences Derived from the ORF7 Gene
After the sequence of nucleic acid sequencing was removed, the DNA sequencing results were translated by Primer 5.0 software (amino acid sequence was submitted to GenBank as numbers mg922935 - mg922947). The results showed that the ORF7 gene of 13 isolated strains of PRRSV contains a complete open reading frame, and the coding area consisted of 372 bp and encoded 123 amino acids. The DNASTAR software was used to compare the homology of the 13
Figure 2. Phylogenetic tree of the ORF7 gene from 13 isolates of PRRS.
isolated strains of PRRSV with the amino acid sequence of domestic and foreign strains (see Table 4). The homology of the ORF7 gene with the highly pathogenic PRRSV strains isolated from the domestic JXA1 strain was 88.7% - 100%, and 87.1% - 91.7% with the FJFS strain, 89.3% - 97.1% with the NADC30 and JL58 strains, 89.2% - 96.5% with the classical CH-1a and CH-1R strains isolated from China, 89.5% - 94.4% with the MLV vaccine strain, 89.5% - 94.4% with the American standard VR-2332 strain, 89.2% - 94.4% with the classical R98 strain, 55.6% - 58.3% with the European classical Lelystad strain, 88.2% - 100% with the 11 previously isolated PRRSV strains. The internal homology of 13 isolated viral strains isolated was 88.4% - 100%. Through DNASTAR software, the 13 isolated strains of PRRSV were compared with the amino acid sequences of the domestic and foreign strains (as shown in Table 4); as a whole, there were no deletions or insertions of amino acids in the isolates, and there was a point mutation in some sites. Compared with the classical domestic and foreign strains of CH-1a, CH-1R, VR-2332, MLV and R98 domestic and foreign (Figure 3), the 15th amino acid of the isolates (10/13) was mutated from D (aspartic acid) to N (asparagine), the 46th amino acid of the isolates (11/13) was mutated from K (lysine) to R (arginine), the 91st amino acid of the isolates (10/13) was mutated from T (threonine)
Figure 3. Comparison of the variability of the N protein from 13 PRRSV strains with domestic and foreign PRRSV strains.
to A (alanine), the 109th amino acid of the isolates (10/13) was mutated from H (histidine) to Q (glutamine), and the 117th amino acid of the isolates (12/13) was mutated from V (valine) to A (alanine).
4. Discussion
4.1. The Significance of Homology and Epidemiology Analyses in Disease Prevention and Control
According to the classification system established by M. shi et al., the epidemic strains of the American PRRSV, which was spread worldwide were divided into 9 lineages and 37 sublineages [11]. From the genetic evolution tree (Figure 2), 10 strains of the 13 strains isolated from this study belong to lineage 8 (sublineage 8.7), represented by CH-1A and JXA1 strains. In addition, 11 strains of the virus we isolated in 2008 all belong to sublineage 8.7. In the past decade, sublineage 8.7 was the dominant pedigree in the epidemic field PRRSV strains in China, and this is consistent with the results observed by Zhang Chanping et al. [12]. Three strains of the virus isolated in this experiment belong to lineage one, which was represented by JL580 and NADC30 strains. As the PRRSV pedigree in China appeared and became a trend after 2013, lineage one occupies a certain proportion of the field. The results show that the prevalent strains of PRRSV in the field are still HP-PRRSV, and the field isolates are facing tremendous selective pressure. Whether lineage one represented by the NADC30 strain will present the trend of rising and become the dominant strains in the future is unknown. This virus needs to be illustrated by more epidemiological studies and data. Studies have confirmed that PRRSV can be transmitted over long distances through commercial pigs carrying a virus and with infectious semen [13]. In this study, SDCL05 and LNZD03 strains are from Shandong Province and Liaoning Province, respectively. The homology between these two strains is as high as 100%. Although there is nearly a 1000 kilometer-distance between Liaoning Shandong provinces, viruses isolated from these two places are highly homologous. After further analysis of the background of the isolates, it was found that the semen of the pigs in these two places originated from the same place. This suggests that we should cautiously screen pathogens such as PRRSV while doing breed introducing. Sdsg09 and sdsg10 are two strains from the same city of Shandong Province, but their homology is only 89%. This shows that different epidemic strains exist in the same region, so different vaccines should be selected during the immunization process.
4.2. The Significance of the N Protein Gene and Its Encoded Amino Acid Sequence Variation in the Pathogenesis of Disease
There is a point mutation in the nucleic acid sequence of the N protein gene. Comparing the results of homology of the derived amino acid sequences with the homology of nucleic acid sequences, the variation coefficient of the former is much lower than that of the latter. Therefore, nucleic acid mutations at individual sites did not cause changes in the corresponding amino acid sequences. The results showed that there are three highly mutated regions (3 - 16, 37 - 52, 118 - 123) and three highly conserved regions (16 - 27, 51 - 69, 81 - 90, 102 - 106) in N protein amino acid sequences [14]. By analyzing the variability of domestic and foreign classical PRRSV strains (Figure 3), this study found that the high-frequency mutation sites of 13 strains of PRRSV isolates were located in at the 15th amino acid (10/13), 46th amino acid (11/13), 91st amino acid (10/13), 109th amino acid (10/13), and 117th amino acid (10/13). It is interesting to find that, in addition to the 46th amino acid in the highly mutated region, the remaining four high-frequency mutations sites were not in the highly mutated region. For reference, the N protein of the American PRRS strains contains 5 important epitopes, which are located in the regions of 30 - 52, 37 - 52, 69 - 123 [15]. The high frequency mutations of the 4 sites in this study were all located in the antigenic epitopes, but the effects on the immunogenicity of the PRRSV need further experimental verification. This study confirmed that the PRRSV N protein also contained two typical nuclear location signals (NLS); the single type NLS is located in the N protein at amino acids 10 - 13, and the double type is located in the N protein at amino acids 41 - 47 [16]. The study revealed that the NLS-mediated protein nucleation was an important part of the information exchange of nuclear materials in intracellular signaling networks, and it plays a key role in cell proliferation, differentiation, regulation and transformation [17]. Most viral proteins entering the nucleus require nuclear transporter receptor protein importin α or importin β to identify and bind to the NLS sequence carried by the target protein. The pathogenicity, virulence and replication ability of the virus are significantly affected by the mutation in the NLS structure and sequence carried by virus that enters into the nucleoprotein. For instance, Venezuelan equine encephalitis virus capsid protein helps to infect itself by inhibiting the synthesis of beta-interferon in host cells after it enters into the nucleus protein; therefore, when the nuclear protein NLS sequence is mutated, the virulence of the virus weakens significantly [18]. The dengue virus nonstructural protein 5 (NS5) plays a role in resisting the antiviral effects of host cells, including inhibition of the synthesis of target cells IL-8, and when the NS5 protein NLS sequence is mutated, virus replication capability decreases by nearly 1000 times [19]. The duration of virusemia in infected animals becomes shorter and the ability of neutralizing antibody increases markedly [20]. It was found in a previous study that, compared with the classical PRRSV strains, the N protein dual-type NLS structure of HP-PRRSV, which emerged in China in 2006, had a consistent mutation in the 46th amino acid (K46 was replaced by R46) [10]. This study showed that among the 13 strains of PRRSV field isolates, there were 11 strains whose 46th amino acid mutations in their N protein are arginine. In that previous study, our subjects were strains isolated from 2007-2009, and this research isolated PRRSV strains from 2015-2016. Combining the PRRSV strains published in 2010-2014 by GenBank, we speculate that the mutation at the 46th amino acid (K46→R46) in the N protein dual-type NLS structure has been stably inherited since the outbreak of HP-PRRSV. Studies have been conducted on the mutation of this locus from the perspective of immunology, and the mutation of this locus can cause the host cell to produce more CD4 + CD25 + FOXP3 + T-cell specific responses during the process of virus infection, thus benefiting the immune evasion of the virus [9]. However, there is no research conducted from the perspective of the nucleoprotein (the “duty” of NLS structure) to explain the response to this mutation at this locus.
Author’s Contributions
YX and BH carried out most of the experiments and wrote the manuscript. MY, BW, YZ participated in the planning of the project. JL and YL conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.
Acknowledgements
Financial support for this study was obtained from the National Natural Science Foundation of China (Grant number 31860692), Scientific and Technological Project of Guizhou Province (Grant number Qian Ke He Ji Chu [2016]1166).
NOTES
*These authors contributed equally to this work and should be considered as first authors.
#Corresponding authors.