- Research article
- Open Access
Genetic characterization and mutation analysis of Qihe547 Aujeszky’s disease virus in China
- Cun Liu†1,
- Yanhan Liu†2,
- Ye Tian2,
- Xuehua Wei1,
- Yue Zhang1 and
- Fulin Tian1Email authorView ORCID ID profile
© The Author(s). 2018
Received: 3 January 2017
Accepted: 11 May 2018
Published: 6 July 2018
Aujeszky’s disease virus (ADV) can cause neurologic disease in young pigs, respiratory disease in older pigs and abortion or birth of mummified fetuses or stillborn neonates. The re-emergence of Aujeszky’s disease (AD) in pig farms vaccinated with live vaccine (Bartha-K61) caused substantial economic losses to Chinese pig industry since late 2011. A field ADV, named Qihe547, was isolated from pigs that exhibited suspected AD clinical symptoms. To better understand the genetic characteristics and mutations of Qihe547 ADV, the whole genome was sequenced and analyzed.
The genomic length of Qihe547 ADV was 143,404 bp, with 73.59% G + C contents. Phylogenetic analysis based on the whole genome of ADV strains revealed that Chinese ADV strains were located to one group with three subgroups. Qihe547 ADV was closely related to these novel ADV strains isolated in China since 2012. Qihe547 presented numerous hypervariable regions compared with oversea ADV strains. In 34 genes of Qihe547 ADV, amino acid (AA) insertion or deletion were observed. In addition, numerous AA mutations were found in the main protective antigen genes (gB, gC and gD genes). The differences of potential antigenic peptides in the main protective antigens between Qihe547 ADV and ADV Bartha were discovered in the dominant antigenic regions of gB (AA59-AA126, AA507-AA734),the extracellular region of gC and gD.
High diversity was observed between Qihe547 and foreign ADV isolates. The AA variations and the differences of potential antigenic peptides in the important functional regions of the main protective antigen (gB, gC and gD) of ADV Qihe547 may contribute to immune evasion of the virus and may be partial reason that the virus escapes from the vaccination of Bartha-K61 vaccine. In a word, the effect of the variations obviously requires further research.
Aujeszky’s disease (AD), a devastating swine disease, also known as Pseudorabies (PR), is an economically important viral disease of swine worldwide [1–3]. AD is characterized by high morbidity in piglets, severe respiratory disorders in older pigs and sow’s reproductive failure [2, 4]. AD occurs all over the world and caused great economic losses to the pig industry in many countries .
The causative agent responsible for AD is Aujeszky’s disease virus (ADV) which is a member of the genus Varicellovirus of the subfamily Alphaherpesvirinae within the family Herpesviridae . The genome of ADV is double-stranded DNA, approximately 150 kb, with G + C contents more than 70%. Swine are the natural host for ADV and the reservoir of the virus in nature , but a wide range of mammals can be infected with ADV .
Vaccination is widely used to control the AD. Attenuated live vaccine has been widely used for many years in many countries to diminish the economic losses caused by AD . Due to large-scale compulsory vaccination with gE-deleted vaccines and the DIVA (differentiating infected from vaccinated animals) strategy, AD has been eradicated from domesticated pigs in many countries [3, 10]. In the 1970s, China imported the Bartha-K61 vaccine from Hungary . Then, it was widely used in swine farmers, from the 1990s until late 2011, AD well controlled. Since late 2011, however, AD broke out in China and has caused huge economic losses [11, 12]. It is great necessary to understand the genetical mutation of re-emergenced ADV strain.
In this study, we genetically characterized Qihe547 ADV strains from swine herd vaccinated with live vaccine by the Sanger dideoxy-chain termination method and phylogenetic analysis to get a better insight into the genetic relatedness between strains in both populations and to compare Qihe547 strains with strains circulating in neighboring countries. To understand the mutations of Qihe547 ADV and the causes of immune failure, comparative genomic analyses and phylogenetic analysis were performed.
Ethics and consents
Sampling was performed under the permission from the owner of the swine herd in Qihe, Shandong province, China and the ethical approval was granted by the Administration Bureau of Animal Husbandry and Veterinary Medicine of Shandong. Animal manipulation procedures adhered to the guidelines under the Animal Ethics Committee of the Administration Bureau of Animal Husbandry and Veterinary Medicine of Shandong.
In 2014, outbreak of AD occurred in a pig farm in Qihe, Shangdong province, China. In the farm, the pregnant sows showed high rates of abortion and some fattening pigs showed severe respiratory symptoms along with high fever and then died. The tissue samples, including brain, lung, kidney, liver, and lymph node, were collected from the clinical dead pigs for diagnostic tests and the corpses were harmlessly treated after anatomical examination. The clinically dead pigs were from the swine herd in Qihe area, Shandong province, China. The tissue samples were homogenated. The homologized samples were centrifuged at 8000 g for 5 min. Then, the supernatant was used to extracted DNA and virus isolation.
Virus isolation and identification
Nucleic acid was extracted from the supernatant using AxyPrep Body Fluid viral DNA/RNA Mini-prep kit according to the manufacturer’s protocol (Corning Co., China). The nucleic acid was detected for Aujeszky’s disease virus (ADV), porcine reproductive and respiratory syndrome virus (PRRSV), porcine circovirus type 2 (PCV2), swine fever virus (SFV) by the real-time PCR. The supernatant, ADV positive, was filtrated through 0.22 μm filter (Millipore, Milford, MA) and inoculated to monolayer culture of Monolayer culture of porcine kidney 15 (PK-15) cells cultured in dulbecco’s modified eagle mediun (DMEM) (Hyclone) supplemented with 8% fetal bovine serum (Shanghai ExCell Biology Co.). PK-15 cells were maintained at 37 °C in an atmosphere of 5% CO2. The cytopathic effect (CPE) was examined daily.
Then, the virus was identified by indirect immuofluorescence assay (IFA). Briefly, inoculated cells were fixed with cold absolute ethanol for 1 h, blocked with 1% bovine serum albumin (Sigma, St. Louis, MO) for 30 min, and incubated with 200 times diluted Rabbit anti-ADV antibody (Abcam Ab3534) for 1 h, followed by a 2000 dilution of goat anti-rabbit IgG-FITC-conjugated secondary antibody (Abcam Ab6717) for 1 h. Cell staining was examined under a fluorescence microscope. Finally, the ADV was plaque-purified three times and propagated on PK-15 cells. Virus titer was determined by measuring the 50% tissue culture infective dose (TCID50).
PCR and sequencing
Primers for the amplification of Qihe547 ADV whole genome were designed based on ADV Becker (JF797219). The plaque-purified virus was used to extract DNA with the AxyPrep Body Fluid viral DNA/RNA Mini-prep kit according to the manufacturer’s protocol (Corning Co., China). PCR was performed with the Takara LA Taq with GC buffer kit (Takara, Dalian), following the manufacturers’ instructions. Cycling conditions for PCR was as follows: denaturation at 95 °C for 10 min, followed by 38 cycles of 95 °C for 1 min, gradiently annealing temperature 50 °C~ 70 °C for 50 s, and 72 °C for 1 min, and final extension of 72 °C for 10 min. The products were cloned into pMD18-T vector (Takara, Dalian). The recombinant plasmid, containing the target DNA fragments, was sent to BioSune Biotech Co. Ltd. (Jinan, China) for DNA sequencing by the Sanger dideoxy-chain termination method after restriction endonuclease digestion with EcoRI and Hind III.
Sequence assembling and alignments
Sequences of the target DNA fragments were assembled using SeqMan (DNASTAR, Madison, WI). Gene’s annotation of Qihe547 ADV was based on ADV Becker and ADV JS-2012(KP257591). The alignment based on full-length genome sequences between Qihe547 ADV and reference ADV strains from GenBank was performed using the mVISTA genomics analysis tool with global LAGAN alignment .
Each gene’s protein-coding sequence was filtered out from the whole genome sequence of Qihe547 ADV and foreign ADV strains (ADV32751/Italy2014, Kolchis, Bartha, Becker, and Kaplan). To study the deduced amino acid (AA) sequence mutations in each genes of Qihe547 ADV relative to foreign ADV isolates (ADV32751/Italy2014, Kolchis, Bartha, Becker, Kaplan), AA sequence alignments were performed. Alignments of the deduced amino acid sequence of each gene’s protein-coding sequence were performed with MegAlign (DNASTAR, Madison, WI). The potential antigenic peptides in glycoprotein (g) B, gC, gD of Qihe547 ADV and Bartha were predicted by Antigenic Peptide Prediction online tool (http://imed.med.ucm.es/Tools/antigenic.html).
The phylogenetic analysis was performed by the neighbor-joining method. The phylogenetic tree based on full-length ADV genome sequences and gB, gC and gD genes sequences were constructed by MEGA 6.0 software with 1000 bootstrap replication.
Virus isolation and identification
Genomic characterization of ADV Qihe547
The genome of Qihe547 ADV was sequenced and approximately 143 kb long. Its G + C contents was 73.59%. Its genome organization was similar to other ADV strains submitted to GenBank. The genomic annotation was performed according to the sequences information of ADV Becker and JS-2012. Genome sequence of Qihe547 ADV submitted to Genbank, the accession No. is KU056477. Interestingly, compared to other ADV isolates from GenBank, there are special polynucleotide sequences (5’ CCCAGCTCTCCCCCGAGGGCCCAGCTCTCCCCC 3′) at the genome 5′- terminal. Upstream gB gene, the unique repeats region, possible repeating unit (5’ GGGGAGAGGGGAGACGAGA3’) different from ADV Becker (5’GGGACGGAGGGGAGA3’), exists a gap. Multi-genome alignments based on ADV isolates indicated that numerous hypervariable regions were found between Qihe547 ADV and oversea ADV isolates (ADV32751/Italy2014, Batha, Becker), especially in intergenic sequences (Fig. 2a). However, between Qihe547 and other China ADV isolates (JS–2012, TJ, Ea, SC), the genetic diversity showed little differences.
Multiple sequences comparison
Amino acid Indels in Qihe547 ADV genes, compared to other ADV strains
428-432(+QAHSQ), 596-597(−AL), 622(+V), 663-665(+SSS)
71-75(−(E/R)(E/G)ER(M/T)), 128(−E), 224(+L)
502-503(-PP), 568-570(+EEE), 588-592(+GNAAD), 618-621(-GSFR)
345(+G), 1123(+G), 1140(+G)
121(+V), 134(+K), 142-144(+DED)
168(+A), 189-190(-GP), 200-203(+RGES), 214-218(+GSGAK), 648-649(+PG)
249-257(+GAPAVAAVG), 562-563(−GA), 2256-2269(+APPAEAAPAAKPAP), 2339-2347(+RAPPPPQPQ), 2368-2392(+PPPQQQQQQQQQRQEPPAATASKKA), 2449-2456(+TQPAAPAE), 2476-2478(+TAK), 2736(+S), 2911(+V), 2973-2974(+PA), 3040-3045(+RAEPAR), 3056(+E), 3110-3111(+PF), 3126-3128(+SPL), 3142(+E), 3211(+P)
The differences of the potential antigenic peptide in gB, gC and gD between Qihe547 and Bartha ADV. In the N terminus of ADV gB, there are antigenically important regions, such as (AA59-AA126, AA214-AA279, AA507-AA734) . The extracellular region (AA23-AA453) of ADV gC plays an important role in ADV adsorption, neutralizing antibody generation and virus neutralization. gD also is the main target of ADV antibody. Mutation occurred in these regions may result antigen drift. Numerous AA mutations and potential antigenic peptides differences were found in gB, gC and gD of Qihe547, this may lead to immune failure
AD caused substantial economic losses to the pig industry in many countries, especially in regions with dense pig populations . Although it has been eradicated from domesticated pigs in many countries, AD remains one of the most important diseases of swine. Mounting evidences showed that ADV infections are more widespread in wild animals and it is a threat to domestic pigs [2, 14–16]. In China, tremendous progress has been made in controlling and eliminating ADV in herds, including large-scale vaccination with live vaccine from the 1990s until late 2011. Since late 2011, the re-emergence of AD in pig farms has attracted the attention of researchers. In this study, Qihe547 strain, a field Aujeszky’s disease virus, was isolated and identified from herd vaccinated with live vaccine (Bartha-K61). For a better understanding of Qihe547 ADV genome characteristic, the full-length genome was sequenced.
Similar genome organization, with 73.59% G+C contents, were observed. The differences between Qihe547 and foreign ADV isolates distributed whole-genome-widely. Some interesting regions were found in Qihe547 ADV genome sequence, such as the unique upstream repeat region of gB gene which the possible repeat unit was different with Becker strain. Further study was needed to research the function of the repeat region for the virus. According to the phylogenetic analysis, Qihe547 ADV within subgroup 1a was closely related to these novel ADV strains isolated since 2012. Compared with foreign ADV strains (ADV32751/Italy2014, Kolchis, Bartha, Becker, Kaplan), unique AA indels were observed in 34 genes of Qihe547 ADV, mainly in proteins associated with viral egress, viral replication and genetic regulation. The highest divergence (more than 20%) was found in US1 protein acting as regulatory protein between Qihe547 ADV and Bartha. UL49.5 protein (gN) associated with immune evasion and UL36 and UL 51 proteins associated with viral egress also showed higher variation rates. Previous studies also draw similar conclusions . Varying degrees of variation were also found in other proteins.
Mutations in antigen contribute to immune failure. Glycoprotein gB, gC and gD, the main protective antigen of ADV, could induce protective immune responses. The variability of gB, gC and gD was beneficial to ADV evading the host immune defense mechanism . The differences displayed on glycoproteins of gB, gC and gD between ADV vaccine strain and prevalent ADV strains might contribute to immune failure. Vaccination or challenge experiments in sheep and pigs indicated that Bartha-K61 vaccine could not provide complete protection against the current prevalent ADV in China [19, 20]. In this study, numerous AA mutations were found in gB, gC and gD of Qihe547 ADV, especially in the dominant antigenic regions of gB (AA59-AA126 and AA507-AA734) .The heparan binding domains of gC (AA44-AA290) were also the antibody-binding domains . As a result of AA variations, the potential antigenic peptides in gB, gC and gD were inconsistent between Qihe547 ADV and ADV Bartha. The differences of potential antigenic peptides in Qihe547 gB were mainly found in the dominant antigenic region (AA59-AA126, AA507-AA734). The extracellular region (AA23-AA453) in gC of ADV plays an important role in ADV adsorption, neutralizing antibody generation and virus neutralization. Many AA variations and the differences of potential antigenic peptides were observed in extracellular region in gC of ADV Qihe547. The differences of potential antigenic peptide were discovered in gD of ADV Qihe547. These antigenic peptide differences were found in the dominant antigenic region (AA59-AA126, AA507-AA734) of gB, the extracellular region (AA23-AA453) of gC, the dominant antigenic region of gD. The main protective antigen of Qihe547 ADV may be closely related with antigenicity of the virus, probably resulting neutralizing antibody differences between field strains and Bartha. And this may be the partial reason that the virus escapes from the vaccination of Bartha-K61 vaccine.
In this study, high diversity was observed between Qihe547 and other ADV isolates (ADV32751/Italy2014, Kolchis, Bartha, Becker, and Kaplan). AA mutations were found in nearly each gene and unique AA indels were observed in 34 genes of ADV Qihe547. The AA variations and the differences of potential antigenic peptides in the important functional regions of the main protective antigen (gB, gC and gD) of ADV Qihe547 may contribute to immune evasion of the virus. The effect of the variations obviously requires further research.
This study was funded by Agricultural Technology Innovation Projects of Livestock and Poultry Breeding.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
CL carried out the studies and drafted the manuscript. YL helped to draft the manuscript and modified the manuscript. FT designed the study. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ao JQ, Wang JW, Chen XH, Wang XZ, Long QX. Expression of pseudorabies virus gE epitopes in Pichia pastoris and its utilization in an indirect PRV gE-ELISA. J Virol Methods. 2003;114:145–50.View ArticlePubMedGoogle Scholar
- Cramer SD, Campbell GA, Njaa BL, Morgan SE, Smith SK, McLin WR, Brodersen BW, Wise AG, Scherba G, Langohr IM, et al. Pseudorabies virus infection in Oklahoma hunting dogs. J Vet Diagn Investig. 2011;23:915–23.View ArticleGoogle Scholar
- Pomeranz LE, Reynolds AE, Hengartner CJ. Molecular biology of pseudorabies virus: impact on neurovirology and veterinary medicine. Microbiol Mol Biol Rev. 2005;69:462–500.View ArticlePubMedPubMed CentralGoogle Scholar
- Nauwynck HJ. Functional aspects of Aujeszky's disease (pseudorabies) viral proteins with relation to invasion, virulence and immunogenicity. Vet Microbiol. 1997;55:3–11.View ArticlePubMedGoogle Scholar
- Klupp BG, Hengartner CJ, Mettenleiter TC, Enquist LW. Complete, annotated sequence of the pseudorabies virus genome. J Virol. 2003;78:424–40.View ArticleGoogle Scholar
- Crandell RA. Pseudorabies (Aujeszky's disease). Vet Clin North Am Large Anim Pract. 1982;4:221–9.Google Scholar
- Marcaccini A, Lopez Pena M, Quiroga MI, Bermudez R, Nieto JM, Aleman N. Pseudorabies virus infection in mink: a host-specific pathogenesis. Vet Immunol Immunopathol. 2008;124:264–73.View ArticlePubMedGoogle Scholar
- Müller TF, Teuffert J, Zellmer R, Conraths FJ. Experimental infection of European wild boars and domestic pigs with pseudorabies viruses with differing virulence. Am J Vet Res. 2005;62:252–8.View ArticleGoogle Scholar
- Ficinska J, Bienkowska-Szewczyk K, Jacobs L, Plucienniczak G, Plucienniczak A, Szewczyk B. Characterization of changes in the short unique segment of pseudorabies virus BUK-TK900 (Suivac a) vaccine strain. Arch Virol. 2003;148:1593–612.View ArticlePubMedGoogle Scholar
- Muller T, Hahn EC, Tottewitz F, Kramer M, Klupp BG, Mettenleiter TC, Freuling C. Pseudorabies virus in wild swine: a global perspective. Arch Virol. 2011;156(10):1691–705.View ArticlePubMedGoogle Scholar
- An TQ, Peng JM, Tian ZJ, Zhao HY, Li N, Liu YM, Chen JZ, Leng CL, Sun Y, Chang D, et al. Pseudorabies virus variant in Bartha-K61-vaccinated pigs, China, 2012. Emerg Infect Dis. 2013;19(11):1749–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu R, Bai C, Sun J, Chang S, Zhang X. Emergence of virulent pseudorabies virus infection in northern China. J Vet Sci. 2013;14(3):363.View ArticlePubMedPubMed CentralGoogle Scholar
- Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, Sidow A, Batzoglou S. LAGAN and multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003;13(4):721–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Verin R, Varuzza P, Mazzei M, Poli A. Serologic, molecular, and pathologic survey of pseudorabies virus infection in hunted wild boars (Sus scrofa) in Italy. J wildlife Dis. 2014;50(3):559–65.View ArticleGoogle Scholar
- Verpoest S, Cay AB, Bertrand O, Saulmont M, De Regge N. Isolation and characterization of pseudorabies virus from a wolf (Canis lupus) from Belgium. Eur J Wildlife Res. 2013;60(1):149–53.View ArticleGoogle Scholar
- Hahn EC, Fadl-Alla B, Lichtensteiger CA. Variation of Aujeszky's disease viruses in wild swine in USA. Vet Microbiol. 2010;143(1):45–51.View ArticlePubMedGoogle Scholar
- Yu T, Chen F, Ku X, Fan J, Zhu Y, Ma H, Li S, Wu B, He Q. Growth characteristics and complete genomic sequence analysis of a novel pseudorabies virus in China. Virus Genes. 2016;52(4):474–83.View ArticlePubMedGoogle Scholar
- Takada A, Kida H. Induction of protective antibody responses against pseudorabies virus by intranasal vaccination with glycoprotein B in mice. Arch Virol. 1995;140(9):1629–35.View ArticlePubMedGoogle Scholar
- Wang CH, Yuan J, Qin HY, Luo Y, Cong X, Li Y, Chen J, Li S, Sun Y, Qiu HJ. A novel gE-deleted pseudorabies virus (ADV) provides rapid and complete protection from lethal challenge with the ADV variant emerging in Bartha-K61-vaccinated swine population in China. Vaccine. 2014;32(27):3379–85.View ArticlePubMedGoogle Scholar
- Luo Y, Li N, Cong X, Wang CH, Du M, Li L, Zhao B, Yuan J, Liu DD, Li S, et al. Pathogenicity and genomic characterization of a pseudorabies virus variant isolated from Bartha-K61-vaccinated swine population in China. Vet Microbiol. 2014;174(1-2):107–15.View ArticlePubMedGoogle Scholar
- Zaripov MM, Morenkov OS, Fodor N, Braun A, Schmatchenko VV, Fodor I. Distribution of B-cell epitopes on the pseudorabies virus glycoprotein B. J Gen Virol. 1999;80(3):537–41.View ArticlePubMedGoogle Scholar
- Ober BT, Teufel B, Wiesmüller KH, Jung G, Pfaff E, Saalmüller A, Rziha HJ. The porcine humoral immune response against pseudorabies virus specifically targets attachment sites on glycoprotein gC. J Virol. 2000;74(4):1752–60.View ArticlePubMedPubMed CentralGoogle Scholar