Ning Li
Jing Liu1
Lei Xu
Kangkang Guo- 1College of Veterinary Medicine, Northwest A&F University, Yangling, China
- 2College of Life Sciences, Northwest A&F University, Yangling, China
As the major pathogen for porcine circovirus-associated disease (PCVAD), porcine circovirus type 2 (PCV2) is no longer treated as an emerging virus anymore. The wide distribution of PCV2 infection in China causes huge economic losses in the swine industry. Currently, it is generally believed that PCV2 has eight genotypes (PCV2a to PCV2h), with PCV2a, PCV2b, and PCV2d being widely distributed. To comprehensively explore the genetic diversity and prevalence of PCV2 in China, PCV-2 sequences submitted from China in the GenBank database were retrieved. With a total of 714 PCV2 strains were retrieved, we found that early-submitted PCV2 sequences were mainly collected from coastal provinces in the southeast part of China, which may indicate PCV2 was initially circulating in those regions. From 2002 to 2008, PCV2b was the dominant prevalent genotype in those retrieved sequences. From 2009, PCV2d became the dominant genotype in those sequences, dropping a hint that a potential shift of PCV2b to PCV2d might occur in 2009, which is similar to the patterns at the global level. In addition to the PCV2a, PCV2b, and PCV2d genotypes, novel strains were also characterized. We further revealed that the amino acid sequences consistency of PCV2a Cap is higher than those in other genotypes. Together, this study provided clues for the possible prevalent genotypes and dynamics of genetic diversity in China from 2000 to 2019.
Introduction
As the smallest non-enveloped virus, porcine circoviruses (PCVs) have a circular, single-stranded DNA genome. They are members of the genus Circovirus in the Circoviridae family (1). First Porcine circovirus-associated disease (PCVAD) was reported in 1998 and the pathogen was characterized as porcine circovirus type 2 (PCV2) (2–4). In addition to PCV2, other species of PCVs including porcine circovirus type 1 (PCV1) and porcine circovirus type 3 (PCV3) were also identified (1, 5). Recently, a novel circovirus PCV4 was first reported in Hunan province, China (6) and later in other provinces of China (7, 8).
Currently, eight genotypes (PCV2a-PCV2h) have been identified (9). A nucleotide diversity cut-off of 3.5% was used to define different PCV2 genotypes. PCV2a and PCV2b were identified in 2008 (10) and retrospectively study characterized PCV2c from Denmark (11). Recent studies indicated this genotype also circulating in feral pigs from Brazil and in China (12, 13). Subsequently, a fourth genotype named PCV2d has emerged (14) and this genotype is prevalent worldwide. Increasingly new genotypes have been identified like PCV2e was identified from the USA and Mexico (15, 16). Furthermore, PCV2f, which has a novel viral sequence clustering that does not match the existing genotypes was identified in China (17). The distribution of PCV2g was reported in Austrian pigs in a sample from 2009 and also in the India pig population (18, 19). PCV2h was detected in China and Korean (20, 21).
It is commonly believed the length of the PCV2 genome is about 1,767 nucleotides (nt), with six open reading frames (ORFs) have been characterized (22). The ORF1 and ORF2 encode replicase (Rep) and capsid (Cap) proteins are major ORFs (23). ORF2 encoded Cap protein is the major structural protein that responsible for viral entry and triggering neutralization antibodies production. Compared with other ORFs, the ORF2 gene can be used as a phylogenetic and epidemiological marker due to its high genetic variations. For instance, compared with PCV2a-2d, the ORF2 sequence of PCV2e has 12 or 15 additional nts (CTT/TCTTATATGTAA) (24).
In China, the first identification of PCV2 in pig farms was reported in 2000 and increasing prevalence was reported after that (14, 25–29). Although remarkable knowledge of the epidemiology of PCV2 infection in China has been obtained (20, 30), there are few studies focusing on the genetic diversity and prevalence of PCV2 from the emerging of this important virus. In this study, we collected and analyzed the sequences from GenBank that were submitted from China between 2001 and 2019. Based on those sequences, we suggested that PCV2 genotypes in China are undergoing constant changes and the predominant strains were shifted to PCV2d in recent years.
Materials and Methods
Sequence Collection
Seven hundred and fourteen complete PCV2 genomes, collected between 2001 and 2019 from different provinces in China, were retrieved from GenBank at the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/). ORF2 of these sequences was analyzed by ORFfinder in NCBI.
GenBank Accession Numbers
The GenBank accession numbers of 714 PCV2 genomes retrieved from NCBI are listed in Supplementary Table 1. The reference sequences of PCV2a (Accession: KY810321.1 and KY940532.1), PCV2b (Accession: KX831482.1 and MF150189.1), PCV2d (Accession: MF169693.1 and MF169747.1) and genotypes were obtained.
Sequence Alignment and Phylogenetic Analysis
All sequences were aligned by ClustalW using MEGAX. The maximum likelihood (ML) fits of 24 different nucleotide substitution models were used to find the most suitable parameters of the phylogenetic tree. Models with the lowest BIC scores (Bayesian Information Criterion) are considered to describe the substitution pattern the best. Phylogenetic trees were constructed using ORF2 sequences of each complete sequence, by using the maximum-likelihood method in MEGAX. The alignment sequence was used to construct the phylogenetic tree with MEGAX, and the parameters were as follows: the statistical method is ML, the test of phylogenetic is the bootstrap method, the number of bootstrap replications is 1,000, the gaps or missing data treatment is partial deletion and the ML heuristic method is Nearest-Neighbor-Interchange (NNI).
Amino Acid Sequence Alignment of Cap Protein
15 sequences of PCV2a, PCV2b and PCV2d were selected randomly for each genotype and two PCV2e sequences found in 2017 were also included and amino acid sequences of Cap protein were generated by using Lasergene software (DNASTAR, Madison, WI, USA). Next, the amino acid sequence difference between those sequences was compared Based on the Molecular Evolutionary Genetics Analysis X (MEGA X).
Possible Spatial Structure Analysis
To analyze the possible spatial structure analysis of Cap protein of different genotypes, the consensus amino acid sequences of different genotypes were obtained by Megalign in Lasergene software (DNASTAR, Madison, WI, USA). The protein sequences of different genotypes to DNAMAN 8.0 (Lynnon Biosoft, Foster City, CA, USA) and then the possibility of different secondary structures at different locations were analyzed.
Results
The Prevalence of PCV2 in China During the Period 2000-2019
To explore the transmission of PCV2 in China, PCV2 sequences deposited in the GenBank database were retrieved. A total of 714 PCV2 strains were collected from the GenBank that were submitted from China between 2001 and 2019 (Supplementary Table 1). In 2001, four sequences of PCV2 strains were collected from the GenBank. However, the location that collected these strains was not mentioned. In 2002, sequences for PCV2 strain collected from four coastal provinces including Guangdong, Shandong, Tianjin and Zhejiang provinces were submitted to GenBank. In 2003, PCV2 infection was mentioned in Hubei, Shanghai, Fujian, Guangxi, Hainan and Hunan provinces, which are coastal provinces or next to the provinces that already have PCV2 outbreaks. Next year, PCV2 infection spread to Jiangxi, Beijing and Henan provinces. After that, PCV2 gene sequences in GenBank suggested PCV2 infection might spread in Jiangsu, Hebei provinces (2005), Hong Kong (2006) and the northeast part of China, Heilongjiang (2007) (Figure 1). Thenceforth, an enormous amount of PCV2 sequences were submitted from China (Supplementary Table 1). Based on this chronological information, we speculated that PCV2 may initially circulate in coastal provinces in the southeast part of China and then spread to other provinces.
Figure 1. The detection PCV2 in China. The year on the map indicate the time when PCV2 sequences were submitted for the first time in these provinces.
Phylogenetic Analysis of PCV2 Strains
To more comprehensively understand the extent of genetic diversity of PCV2 in China, 714 PCV2 strains sequences in GenBank were analyzed by phylogenetic analysis. Phylogenetic construction was conducted using ORF2 sequences. As shown in Figures 2A,B, from 2001 to 2006, the PCV2a was the most prevalent genotype while PCV2d may become the second prevalent genotype from 2007 to 2009 in those retrieved sequences. Of note, after 2009, PCV2d became the dominant strain while PCV2a is the minority in those sequences (Figures 2C,D). From 2014 to 2016, two strains that were recombined by PCV2b and PCV2d were identified (Figure 2E) and from 2017 to 2019, novel PCV2e was identified in those sequences (Figure 2F).
Figure 2. The phylogenetic ML tree of ORF2 of PCV2 strains. The phylogenetic ML tree of ORF2 from PCV2 strains deposited in GenBank from China between 2001 and 2019. (A) ML tree of ORF2 from PCV2 strains from 2001 to 2006. (B) ML tree of ORF2 from PCV2 strains from 2007 to 2009. (C) ML tree of ORF2 from PCV2 strains from 2010 to 2011. (D) ML tree of ORF2 from PCV2 strains from 2012 to 2013. (E) ML tree of ORF2 from PCV2 strains from 2014 to 2016. (F) ML tree of ORF2 from PCV2 strains from 2017 to 2019. Different colors indicate different genotypes. Red, PCV2a; Yellow, PCV2b; Blue, PCV2d; Purple, PCV2e; Green, PCV2b+2d. The tree was constructed by the maximum likelihood method with bootstrap values calculated for 1,000 replicates.
The Prevalence of PCV2 Genotypes in China
To further explore the prevalence of PCV2 in China, the genotype of 714 PCV2 strains collected that from GenBank were analyzed. Between 2001 and 2019, the main prevalent PCV2 genotypes in those sequences were PCV2a, PCV2b and PCV2d (Figure 3 and Supplementary Table 2). By analyzing the complete genome sequences, we showed that all sequences summited in 2001 were belong to genotype PCV2a. However, this did not suggest no PCV2b and PCV2d was circulating at the time since only four sequences were collected. From 2002 to 2008, PCV2b was the dominant prevalent genotype in those sequences, while the prevalence of PCV2a was lower than PCV2d. During this time, the proportion of PCV2d witnessed a fluctuation but its overall tendency was increased. From 2009, PCV2d become the predominant genotype in those sequences, while the proportion of the PCV2b in all sequences was substantially decreased. Of note, PCV2a was still the lowest prevalent strain in those sequences during this period. The analysis of those sequences may suggest PCV2d is still the main epidemic strain currently and the prevalence of PCV2a, PCV2b, and PCV2d are rather stable from 2009. Specifically, the proportion of PCV2a in all retrieved sequences at that year was around 10%, while the proportion of PCV2b was around 20-30% and PCV2d was <70%.
Figure 3. The proportion of different PCV2 genotypes from 2001 to 2019. The proportion of different PCV2 genotypes in each year was calculated. Red, PCV2a; Blue, PCV2b; Yellow, PCV2d; Black, Others.
Although these three genotypes (PCV2a, PCV2b, and PCV2d) were the most widely detected strains in those retrieved sequences from the report of PCV2 infection, novel strains were also characterized. In 2014 and 2017, some sequences that did not belong to those major genotypes were also identified. The new genotype identified in 2014 was the recombinant of PCV2b and PCV2d, and in 2017, the emerged genotype was PCV2e (Figure 3 and Supplementary Table 2), which was first reported in the United States in 2015 with a 702 bp ORF2 (15, 31).
Amino Acid Sequence Differences of PCV2 Cap Protein
To further characterize differences of Cap protein among PCV2a, PCV2b, PCV2d, and PCV2e, amino acid sequence alignment was performed. We showed that the length of Cap protein varies in different genotypes (Figure 4). Most of the PCV2a and PCV2b strains have a Cap protein with 233 amino acids, while the majority of PCV2d strains have a Cap protein with 234 amino acids. Of note, many PCV2e strains have a larger Cap protein with 238 amino acids. Amino acid sequences similarity of Cap within specific genotype are high. Particularly, amino acid sequences consistency of PCV2a Cap is greater than that in other genotypes. However, significant differences in Cap protein amino acid sequence were observed between PCV2e and other genotypes.
Figure 4. Amino acid sequence differences of PCV2 Cap protein. 15 PCV2a, 15 PCV2b, 15 PCV2d, and 2 PCV2e were selected respectively, and their amino acid sequences were compared by using the Alignment. The Alignment was performed using the software MegAlign, the major amino acid mutations are displayed in a box. The overlapping antigen sites (47-63, 165-200, and 230-233) of Cap protein among different genotypes are indicated by red dashed boxes.
To explore whether those essential residues are changed in different PCV2 genotypes, we compared the amino acid residues of PCV2a, PCV2b, and PCV2d. As shown in Table 1, amino acid sequences of three spatially overlapping antigen sites 47-63, 165-200, and 230-233 of Cap protein are largely different among different genotypes. Finally, analysis of the possible spatial structures of different genotypes revealed PCV2a and PCV2b showed high similarity while significantly different were observed PCV2d and PCV2e (Figure 5).
Table 1. Amino acid sequences difference among PCV2a, PCV2b, and PCV2d.
Figure 5. Predicated spatial structure of PCV2 ORF2 of different genotypes. The possible spatial structure of PCV2 ORF2 of different genotypes as helices, strands and coils.
Discussion
Porcine circovirus 2 (PCV2) is no longer treated as an emerging virus anymore, due to its widespread worldwide. The first report in 1998 (2–4), tremendous attention has been drawn by this virus, since its infection causes huge economic losses in the swine industry, including China. The first report of the postweaning multisystemic wasting syndrome (PMWS) related to PCV2 infection in China was published in 2000 (25). Then, PCV2 was detected in many tissue samples of the diseased pig (32). In this study, we retrieved PCV2 sequences from GenBank that were submitted from China between 2001 and 2019. Atm the early period, PCV2 sequences were mainly collected from Guangdong, Shandong, Tianjin, and Zhejiang provinces, all of which are coastal provinces in China. The submission of PCV2 sequences from inner provinces in China happened at later times, which may suggest the PCV2 was spread from the coastal provinces, perhaps following major swine commercial routes.
By analyzing those retrieved sequences, we noticed that between 2002 and 2008, the dominant prevalent genotype is PCV2b, which is consistent with the previous report (29). Afterward, although the PCV2 vaccines were widely used, the prevalent genotype may shift from PCV2b to PCV2d. Such genotype shift seems occurred globally and the worldwide use of PCV-2 vaccines may be the driven force (33, 34). Until now, this genotype is widely circulated and has become the predominant genotype worldwide (35–39). In our previous study, we showed that 38 PCV2 positive samples were detected from 300 blood samples collected from the slaughterhouse in Shaanxi province from 2018 to 2019. Of note, genotyping of those PCV2-positive samples showed that 11 samples were PCV2b positive out of 20 sequenced samples, while nine samples were PCV2d positive (40). However, this may only partially reflect the prevalence of PCV2 in Shaanxi due to a small number of samples being obtained. Recently, a study collected clinical samples from Henan province of China and the whole genome sequences of 20 PCV2 strains were analyzed (41). The results showed that among those 20 strains, 30% (6 strains) belonged to PCV2a, 30% (6 strains) belonged to PCV2b and 40% (8 strains) belonged to PCV2d. In another study, they complete genome sequences of 34 PCV2 strains that were collected from 55 different pig farms in central China from 2018 to 2020 (42). Consistently, among those 34 strains, 10 strains belonged to PCV2b and 16 strains belonged to PCV2d and only 8 strains belonged to PCV2a. Taken together, those observations implied that current dominant prevalent genotypes in China are PCV2b and PCV2d. Furthermore, genotype PCV2e was detected in 2017 in the Fujian province in China (31). Although the prevalence of this genotype was rare in our analysis (Figure 2) and was not detected in 300 blood samples collected in Shaanxi province, attention should be paid to this genotype. Importantly, the genotypes shift we showed in this study may not precisely reflect the reality that happened in the field and further clinical samples should be sequenced to determine and predict the circulation of different PCV2 genotypes.
As the principal strategies to control PCV2 infection, vaccination is broadly employed in the field. A commercial vaccine against PCV2a was first developed in the USA with efficient effects in alleviating the viral loads and clinical signs (43, 44). Despite different PCV2 genotypes has been circulating, most commercial vaccines are still generated based on PCV2a (44, 45). In particular, although one study has been reported that the PCV2a-based vaccine is capable to reduce PCV2d infection (46), PCV2d genotype was emerging reported worldwide (47–51). Compared with PCV2a and 2b, PCV2d are more virulent with severe clinical signs, viremia and lesions (52, 53). In line with these studies, we also revealed a potential shift of PCV2b to PCV2d from 2009 (Figure 3), which is similar to the patterns at the global level. Therefore, efforts should be taken to develop new commercially PCV2 vaccines based on the PCV2d genotype, but not PCV2a, since this genotype is probably no longer the predominant circulating genotype in China. The current PCV2a based vaccines should be evaluated for their ability to prevent the infection or transmission of other genotypes.
It has been reported that position 1-41 of Cap protein, a special arginine-rich N-terminus, the nuclear localization signal (NLS) that responds for its nucleus accumulation (54, 55). In particular, positions 12-18 and 34-41 of the amino acid residues are important for the nuclear localization of Cap protein (54). Amino acids 143-145 (N143YS) is the N-linked glycosylation site PCV2 Cap protein (56). The mutation of those sites leads to a higher level of immune responses and improved immunogenicity of Cap-based vaccination and stronger T lymphocyte proliferative activity and cytokine production (56), indicting those conserved sites might responsible for the failure of PCV2 vaccination. In addition, amino acid residues 47-63, 165-200, and 230-233 are characterized as essential for epitope recognition of PCV2 (57). Intriguingly, positions 169, 185, 187, 190, 191, and 232 in PCV2a were different from PCV2b and 2d and such substitutions may relate to the increased circulating of PCV2b and PCV2d in the PCV2a-vaccinated fields. Particularly, it has been proved that the Alanine (A) in Cap protein at position 59 of PCV2a genotype is an essential conformational neutralizing antibody epitope (58). Notably, this site was changed to R/K in PCV2b and K in PCV2d (Table 1), which might lead to a reduced recognition of this epitope in these two emerging genotypes. Consequently, animals that were vaccinated by PCV2a-based vaccine may have decreased antibody neutralizing ability against PCV2b infection.
PCV2a and PCV2b strains encode a Cap protein with 233 amino acids and PCV2d strains have a Cap protein with 233 or 234 amino acids. Of note, many PCV2e strains have extra12 nucleotides at the carboxyl terminus of ORF, which encodes Cap protein with 238 amino acids. A possible explanation of the extra length of the PCV2e capsid is that it is the prototype and a 12 nt deletion generated the progenitor of PCV2a to PCV2d (16). Importantly, such extra carboxyl terminus seems not to affect the predicted protein structure of PCV2e capsid protein (16). The different length between the PCV2a/PCV2b Cap protein and some of the PCV2d Cap protein is due to one extra Lysine (K) at the carboxyl terminus (Figure 4).
Cap protein is the main determinant of the PCV2 antigenicity and immunogenicity and the mutation of Cap protein may facilitate PCV2 to evade the host immune response. A single amino acid mutation in Cap protein led to a strain that can not be neutralized by currently available antibodies (58). The PCV2 nucleotide substitution was estimated to be 1.2 × 10−3 sites per year, which is extremely high for a single-stranded DNA virus (59). With such high evolution rates, prevalent genotypes are constantly varying and difficulties for vaccine development were substantially increased. In this study, we found that amino acid residues 47-63, 165-200, and 230-233, which are essential for epitope recognition of PCV2, are largely different among different genotypes (57). Further investigations regarding the antigens and protein structure variations of different PCV2 genotypes could pave the road for developing the new vaccine and deepen the understanding of PCV2 pathogenesis. In addition, co-infections of two different PCV2 genotypes in the field were also reported (60). In China, co-infection of different PCV2 genotypes (61) and even co-infection of PCV2 and PCV3 were also reported (62). Therefore, a better understanding of the molecular epidemiology of PCV2 is vital for controlling further infections.
In conclusion, by analyzing PCV2 strains collected from GenBank between 2001 and 2019 in China, we may speculate that PCV was initially circulating in coastal regions and spread to other provinces soon afterward. Furthermore, genotyping analysis of these sequences provides clues that PCV2b might be the dominant prevalent genotype from 2002 to 2008 and a potential shift to PCV2d might occur in 2009. Finally, we found amino acid sequences consistency of PCV2a Cap is greater than those in other genotypes. These results highlighted prevalent genotypes and dynamics of genetic diversity in China from 2000 to 2019, which is essential for the understanding of molecular epidemiology of PCV2-related diseases and the development of new PCV2 vaccines.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
KG had the idea for the article. NL and LX analyzed the data. JL, JQ, and FH helped to analyze the data. KG and LX wrote the draft of the manuscript and critically revised the work. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the National Natural Science Foundation of China (Project No. 31672580).
Author Disclaimer
Frontiers Media SA remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Conflict of Interest
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.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.788172/full#supplementary-material
References
1. Rosario K, Breitbart M, Harrach B, Segales J, Delwart E, Biagini P, et al. Revisiting the taxonomy of the family Circoviridae: establishment of the genus Cyclovirus and removal of the genus Gyrovirus. Arch Virol. (2017) 162:1447–63. doi: 10.1007/s00705-017-3247-y
2. Allan GM, Mcneilly F, Kennedy S, Daft B, Clarke EG, Ellis JA, et al. Isolation of porcine circovirus-like viruses from pigs with a wasting disease in the USA and Europe. J Vet Diagn Invest. (1998) 10:3–10. doi: 10.1177/104063879801000102
3. Ellis J, Hassard L, Clark E, Harding J, Allan G, Willson P, et al. Isolation of circovirus from lesions of pigs with postweaning multisystemic wasting syndrome. Can Vet J. (1998) 39:44–51.
4. Meehan BM, Mcneilly F, Todd D, Kennedy S, Jewhurst VA, Ellis JA, et al. Characterization of novel circovirus DNAs associated with wasting syndromes in pigs. J Gen Virol. (1998) 79:2171–9. doi: 10.1099/0022-1317-79-9-2171
5. Palinski R, Pineyro P, Shang P, Yuan F, Guo R, Fang Y, et al. A novel porcine circovirus distantly related to known circoviruses is associated with porcine dermatitis and nephropathy syndrome and reproductive failure. J Virol. (2017) 91:e01879. doi: 10.1128/JVI.01879-16
6. Zhang HH, Hu WQ, Li JY, Liu TN, Zhou JY, Opriessnig T, et al. Novel circovirus species identified in farmed pigs designated as Porcine circovirus 4, Hunan province, China. Transbound Emerg Dis. (2020) 67:1057–61. doi: 10.1111/tbed.13446
7. Chen N, Xiao Y, Li X, Li S, Xie N, Yan X, et al. Development and application of a quadruplex real-time PCR assay for differential detection of porcine circoviruses (PCV1 to PCV4) in Jiangsu province of China from 2016 to 2020. Transbound Emerg Dis. (2021) 68:1615–24. doi: 10.1111/tbed.13833
8. Sun W, Du Q, Han Z, Bi J, Lan T, Wang W, et al. Detection and genetic characterization of porcine circovirus 4 (PCV4) in Guangxi, China. Gene. (2021) 773:145384. doi: 10.1016/j.gene.2020.145384
9. Franzo G, Segales J. Porcine circovirus 2 (PCV-2) genotype update and proposal of a new genotyping methodology. PLoS ONE. (2018) 13:e0208585. doi: 10.1371/journal.pone.0208585
10. Segales J, Olvera A, Grau-Roma L, Charreyre C, Nauwynck H, Larsen L, et al. PCV-2 genotype definition and nomenclature. Vet Rec. (2008) 162:867–8. doi: 10.1136/vr.162.26.867
11. Dupont K, Nielsen EO, Baekbo P, Larsen LE. Genomic analysis of PCV2 isolates from Danish archives and a current PMWS case-control study supports a shift in genotypes with time. Vet Microbiol. (2008) 128:56–64. doi: 10.1016/j.vetmic.2007.09.016
12. Franzo G, Cortey M, De Castro AM, Piovezan U, Szabo MP, Drigo M, et al. Genetic characterisation of Porcine circovirus type 2 (PCV2) strains from feral pigs in the Brazilian Pantanal: an opportunity to reconstruct the history of PCV2 evolution. Vet Microbiol. (2015) 178:158–62. doi: 10.1016/j.vetmic.2015.05.003
13. Liu X, Wang FX, Zhu HW, Sun N, Wu H. Phylogenetic analysis of porcine circovirus type 2 (PCV2) isolates from China with high homology to PCV2c. Arch Virol. (2016) 161:1591–9. doi: 10.1007/s00705-016-2823-x
14. Guo LJ, Lu YH, Wei YW, Huang LP, Liu CM. Porcine circovirus type 2 (PCV2): genetic variation and newly emerging genotypes in China. Virol J. (2010) 7:273. doi: 10.1186/1743-422X-7-273
15. Harmon KM, Gauger PC, Zhang J, Pineyro PE, Dunn DD, Chriswell AJ. Whole-genome sequences of novel porcine circovirus type 2 viruses detected in swine from Mexico and the United States. Genome Announc. (2015) 3:e01315–15. doi: 10.1128/genomeA.01315-15
16. Davies B, Wang X, Dvorak CM, Marthaler D, Murtaugh MP. Diagnostic phylogenetics reveals a new Porcine circovirus 2 cluster. Virus Res. (2016) 217:32–7. doi: 10.1016/j.virusres.2016.02.010
17. Bao F, Mi S, Luo Q, Guo H, Tu C, Zhu G, et al. Retrospective study of porcine circovirus type 2 infection reveals a novel genotype PCV2f. Transbound Emerg Dis. (2018) 65:432–40. doi: 10.1111/tbed.12721
18. Weissenbacher-Lang C, Kristen T, Mendel V, Brunthaler R, Schwarz L, Weissenbock H. Porcine circovirus type 2 (PCV2) genotyping in Austrian pigs in the years 2002 to 2017. BMC Vet Res. (2020) 16:198. doi: 10.1186/s12917-020-02413-4
19. Rajkhowa TK, Lalnunthanga P, Rao PL, Subbiah M, Lalrohlua B. Emergence of porcine circovirus 2g (PCV2g) and evidence for recombination between genotypes 2g, 2b and 2d among field isolates from non-vaccinated pigs in Mizoram, India. Infect Genet Evol. (2021) 90:104775. doi: 10.1016/j.meegid.2021.104775
20. Lv Q, Wang T, Deng J, Chen Y, Yan Q, Wang D, et al. Genomic analysis of porcine circovirus type 2 from southern China. Vet Med Sci. (2020) 6:875–89. doi: 10.1002/vms3.288
21. Song S, Park GN, Choe S, Cha RM, Kim SY, Hyun BH, et al. Genetic diversity of porcine circovirus isolated from Korean Wild Boars. Pathogens. (2020) 9:457. doi: 10.3390/pathogens9060457
22. Li D, Wang J, Xu S, Cai S, Ao C, Fang L, et al. Identification and functional analysis of the novel ORF6 protein of porcine circovirus type 2 in vitro. Vet Res Commun. (2018) 42:1–10. doi: 10.1007/s11259-017-9702-0
23. Zhai SL, Zhou X, Zhang H, Hause BM, Lin T, Liu R, et al. Comparative epidemiology of porcine circovirus type 3 in pigs with different clinical presentations. Virol J. (2017) 14:222. doi: 10.1186/s12985-017-0892-4
24. Park KH, Oh T, Cho H, Yang S, Chae C. The first isolation of porcine circovirus type 2e from a Korean pig. Arch Virol. (2020) 165:2927–30. doi: 10.1007/s00705-020-04827-9
25. Lang H, Zhang G, Wu F, Zhang C. Detection of serum antibody against postweaning multisystemic wasting syndrome in pigs. Chin J Vet Sci Technol. (2000) 30:3–5. doi: 10.16656/j.issn.1673-4696.2000.03.001
26. Larochelle R, Magar R, D'allaire S. Genetic characterization and phylogenetic analysis of porcine circovirus type 2 (PCV2) strains from cases presenting various clinical conditions. Virus Res. (2002) 90:101–12. doi: 10.1016/S0168-1702(02)00141-7
27. Wang F, Guo X, Ge X, Wang Z, Chen Y, Cha Z, et al. Genetic variation analysis of Chinese strains of porcine circovirus type 2. Virus Res. (2009) 145:151–6. doi: 10.1016/j.virusres.2009.05.015
28. Yang S, Yin S, Shang Y, Liu B, Yuan L, Zafar Khan MU, et al. Phylogenetic and genetic variation analyses of porcine circovirus type 2 isolated from China. Transbound Emerg Dis. (2018) 65:e383–92. doi: 10.1111/tbed.12768
29. Wang H, Gu J, Xing G, Qiu X, An S, Wang Y, et al. Genetic diversity of porcine circovirus type 2 in China between 1999-2017. Transbound Emerg Dis. (2019) 66:599–605. doi: 10.1111/tbed.13040
30. Lv N, Zhu L, Li W, Li Z, Qian Q, Zhang T, et al. Molecular epidemiology and genetic variation analyses of porcine circovirus type 2 isolated from Yunnan Province in China from 2016-2019. BMC Vet Res. (2020) 16:96. doi: 10.1186/s12917-020-02304-8
31. Liu J, Wei C, Dai A, Lin Z, Fan K, Fan J, et al. Detection of PCV2e strains in Southeast China. PeerJ. (2018) 6:e4476. doi: 10.7717/peerj.4476
32. Wen L, Guo X, Yang H. Genotyping of porcine circovirus type 2 from a variety of clinical conditions in China. Vet Microbiol. (2005) 110:141–6. doi: 10.1016/j.vetmic.2005.07.003
33. Xiao CT, Halbur PG, Opriessnig T. Global molecular genetic analysis of porcine circovirus type 2 (PCV2) sequences confirms the presence of four main PCV2 genotypes and reveals a rapid increase of PCV2d. J Gen Virol. (2015) 96:1830–41. doi: 10.1099/vir.0.000100
34. Franzo G, Cortey M, Segales J, Hughes J, Drigo M. Phylodynamic analysis of porcine circovirus type 2 reveals global waves of emerging genotypes and the circulation of recombinant forms. Mol Phylogenet Evol. (2016) 100:269–80. doi: 10.1016/j.ympev.2016.04.028
35. Kwon T, Lee DU, Yoo SJ, Je SH, Shin JY, Lyoo YS. Genotypic diversity of porcine circovirus type 2 (PCV2) and genotype shift to PCV2d in Korean pig population. Virus Res. (2017) 228:24–9. doi: 10.1016/j.virusres.2016.11.015
36. Dei Giudici S, Lo Presti A, Bonelli P, Angioi PP, Sanna G, Zinellu S, et al. Phylogenetic analysis of porcine circovirus type 2 in Sardinia, Italy, shows genotype 2d circulation among domestic pigs and wild boars. Infect Genet Evol. (2019) 71:189–96. doi: 10.1016/j.meegid.2019.03.013
37. Hou Z, Wang H, Feng Y, Song M, Li Q, Li J. Genetic variation and phylogenetic analysis of Porcine circovirus type 2 in China from 2016 to 2018. Acta Virol. (2019) 63:459–68. doi: 10.4149/av_2019_413
38. Wang Y, Noll L, Lu N, Porter E, Stoy C, Zheng W, et al. Genetic diversity and prevalence of porcine circovirus type 3 (PCV3) and type 2 (PCV2) in the Midwest of the USA during 2016-2018. Transbound Emerg Dis. (2020) 67:1284–94. doi: 10.1111/tbed.13467
39. Raev S, Yuzhakov A, Alekseev K, Kostina L, Gulyukin M, Stepanova T, et al. Transmission of porcine circovirus genotype type 2 (PCV2) in Russia and genotype association (PCV2d) with porcine dermatitis and nephropathy syndrome (PDNS). In: IOP Conference Series: Earth and Environmental Science. Krasnoyarsk: IOP Publishing. (2019).
40. Luo B, Xu H, Luo L, Fan Y, Yang Z, Wang Y, et al. Detection and genetic variation of PCV2 in slaughtered pigs in Shaanxi Province. Anim Husband Vet Med.(2020) 53:100–5.
41. Xu T, Hou CY, Zhang YH, Li HX, Chen XM, Pan JJ, et al. Simultaneous detection and genetic characterization of porcine circovirus 2 and 4 in Henan province of China. Gene. (2022) 808:145991. doi: 10.1016/j.gene.2021.145991
42. Xu T, Zhang YH, Tian RB, Hou CY, Li XS, Zheng LL, et al. Prevalence and genetic analysis of porcine circovirus type 2 (PCV2) and type 3 (PCV3) between 2018 and 2020 in central China. Infect Genet Evol. (2021) 94:105016. doi: 10.1016/j.meegid.2021.105016
43. Dvorak CM, Yang Y, Haley C, Sharma N, Murtaugh MP. National reduction in porcine circovirus type 2 prevalence following introduction of vaccination. Vet Microbiol. (2016) 189:86–90. doi: 10.1016/j.vetmic.2016.05.002
44. Afghah Z, Webb B, Meng XJ, Ramamoorthy S. Ten years of PCV2 vaccines and vaccination: Is eradication a possibility? Vet Microbiol. (2017) 206:21–8. doi: 10.1016/j.vetmic.2016.10.002
45. Segales J. Best practice and future challenges for vaccination against porcine circovirus type 2. Expert Rev Vaccines. (2015) 14:473–87. doi: 10.1586/14760584.2015.983084
46. Opriessnig T, Xiao CT, Halbur PG, Gerber PF, Matzinger SR, Meng XJ. A commercial porcine circovirus (PCV) type 2a-based vaccine reduces PCV2d viremia and shedding and prevents PCV2d transmission to naive pigs under experimental conditions. Vaccine. (2017) 35:248–54. doi: 10.1016/j.vaccine.2016.11.085
47. Xiao CT, Halbur PG, Opriessnig T. Complete genome sequence of a novel porcine circovirus type 2b variant present in cases of vaccine failures in the United States. J Virol. (2012) 86:12469. doi: 10.1128/JVI.02345-12
48. Salgado RL, Vidigal PM, De Souza LF, Onofre TS, Gonzaga NF, Eller MR, et al. Identification of an emergent porcine circovirus-2 in vaccinated pigs from a Brazilian farm during a postweaning multisystemic wasting syndrome outbreak. Genome Announc. (2014) 2:e00163–14. doi: 10.1128/genomeA.00163-14
49. Seo HW, Park C, Kang I, Choi K, Jeong J, Park SJ, et al. Genetic and antigenic characterization of a newly emerging porcine circovirus type 2b mutant first isolated in cases of vaccine failure in Korea. Arch Virol. (2014) 159:3107–11. doi: 10.1007/s00705-014-2164-6
50. Eddicks M, Fux R, Szikora F, Eddicks L, Majzoub-Altweck M, Hermanns W, et al. Detection of a new cluster of porcine circovirus type 2b strains in domestic pigs in Germany. Vet Microbiol. (2015) 176:337–43. doi: 10.1016/j.vetmic.2015.01.013
51. Xiao C-T, Harmon KM, Halbur PG, Opriessnig T. PCV2d-2 is the predominant type of PCV2 DNA in pig samples collected in the US during 2014–2016. Vet Microbiol. (2016) 197:72–7. doi: 10.1016/j.vetmic.2016.11.009
52. Guo L, Fu Y, Wang Y, Lu Y, Wei Y, Tang Q, et al. A porcine circovirus type 2 (PCV2) mutant with 234 amino acids in capsid protein showed more virulence in vivo, compared with classical PCV2a/b strain. PLoS ONE. (2012) 7:e41463. doi: 10.1371/journal.pone.0041463
53. Opriessnig T, Xiao CT, Gerber PF, Halbur PG, Matzinger SR, Meng XJ. Mutant USA strain of porcine circovirus type 2 (mPCV2) exhibits similar virulence to the classical PCV2a and PCV2b strains in caesarean-derived, colostrum-deprived pigs. J Gen Virol. (2014) 95:2495–503. doi: 10.1099/vir.0.066423-0
54. Liu Q, Tikoo SK, Babiuk LA. Nuclear localization of the ORF2 protein encoded by porcine circovirus type 2. Virology. (2001) 285:91–9. doi: 10.1006/viro.2001.0922
55. Hou Q, Hou S, Chen Q, Jia H, Xin T, Jiang Y, et al. Nuclear localization signal regulates porcine circovirus type 2 capsid protein nuclear export through phosphorylation. Virus Res. (2018) 246:12–22. doi: 10.1016/j.virusres.2017.12.012
56. Gu J, Cao R, Zhang Y, Lian X, Ishag H, Chen P. Deletion of the single putative N-glycosylation site of the porcine circovirus type 2 Cap protein enhances specific immune responses by DNA immunisation in mice. Vet J. (2012) 192:385–9. doi: 10.1016/j.tvjl.2011.08.005
57. Lekcharoensuk P, Morozov I, Paul PS, Thangthumniyom N, Wajjawalku W, Meng XJ. Epitope mapping of the major capsid protein of type 2 porcine circovirus (PCV2) by using chimeric PCV1 and PCV2. J Virol. (2004) 78:8135–45. doi: 10.1128/JVI.78.15.8135-8145.2004
58. Huang LP, Lu YH, Wei YW, Guo LJ, Liu CM. Identification of one critical amino acid that determines a conformational neutralizing epitope in the capsid protein of porcine circovirus type 2. BMC Microbiol. (2011) 11:188. doi: 10.1186/1471-2180-11-188
59. Firth C, Charleston MA, Duffy S, Shapiro B, Holmes EC. Insights into the evolutionary history of an emerging livestock pathogen: porcine circovirus 2. J Virol. (2009) 83:12813–21. doi: 10.1128/JVI.01719-09
60. Hesse R, Kerrigan M, Rowland RR. Evidence for recombination between PCV2a and PCV2b in the field. Virus Res. (2008) 132:201–7. doi: 10.1016/j.virusres.2007.10.013
61. Zhai SL, Chen SN, Wei ZZ, Zhang JW, Huang L, Lin T, et al. Co-existence of multiple strains of porcine circovirus type 2 in the same pig from China. Virol J. (2011) 8:517. doi: 10.1186/1743-422X-8-517
Keywords: PCV2, nucleotide sequence, genetic analysis, cap protein, vaccine
Citation: Li N, Liu J, Qi J, Hao F, Xu L and Guo K (2021) Genetic Diversity and Prevalence of Porcine Circovirus Type 2 in China During 2000-2019. Front. Vet. Sci. 8:788172. doi: 10.3389/fvets.2021.788172
Received: 01 October 2021; Accepted: 16 November 2021;
Published: 15 December 2021.
Edited by:
Vassilis Papatsiros, University of Thessaly, GreeceReviewed by:
Stefan Vilcek, University of Veterinary Medicine and Pharmacy in Košice, SlovakiaFernando Bauermann, Oklahoma State University, United States
Copyright © 2021 Li, Liu, Qi, Hao, Xu and Guo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kangkang Guo, guokk2007@nwsuaf.edu.cn; Lei Xu, xulei@nwafu.edu.cn