key: cord-311204-fc12f845
authors: Zhou, Ling; Tang, Qinghai; Shi, Lijun; Kong, Miaomiao; Liang, Lin; Mao, Qianqian; Bu, Bin; Yao, Lunguang; Zhao, Kai; Cui, Shangjin; Leal, Élcio
title: Full-length genomic characterization and molecular evolution of canine parvovirus in China
date: 2016-04-02
journal: Virus Genes
DOI: 10.1007/s11262-016-1309-y
sha: 
doc_id: 311204
cord_uid: fc12f845

Canine parvovirus type 2 (CPV-2) can cause acute haemorrhagic enteritis in dogs and myocarditis in puppies. This disease has become one of the most serious infectious diseases of dogs. During 2014 in China, there were many cases of acute infectious diarrhoea in dogs. Some faecal samples were negative for the CPV-2 antigen based on a colloidal gold test strip but were positive based on PCR, and a viral strain was isolated from one such sample. The cytopathic effect on susceptible cells and the results of the immunoperoxidase monolayer assay, PCR, and sequencing indicated that the pathogen was CPV-2. The strain was named CPV-NY-14, and the full-length genome was sequenced and analysed. A maximum likelihood tree was constructed using the full-length genome and all available CPV-2 genomes. New strains have replaced the original strain in Taiwan and Italy, although the CPV-2a strain is still predominant there. However, CPV-2a still causes many cases of acute infectious diarrhoea in dogs in China.

Canine parvovirus type 2 (CPV-2) belongs to the genus Protoparvovirus and the Parvoviridae family and was first observed by electron microscopy in 1977 [1] . CPV-2 has three antigenic variants: types 2a, 2b, and 2c. Two antigenic variants, CPV-2a and CPV-2b, are now distributed worldwide [2] . A third CPV-2 variant, which was initially believed to be a Glu-426 mutant and subsequently renamed CPV-2c, was detected in Italy in 2000 [3] and is now circulating there together with types 2a and 2b [4] [5] [6] . The new type 2c has also been reported in Vietnam by Nakamura et al. [7] , who developed monoclonal antibodies specific for that type [7] . Type 2c has also been reported in the United States and South America [8] .

The genome of CPV-2 is a single-stranded DNA molecule that contains two open reading frames (ORF1 and ORF2). ORF1 is located in the 5 0 region of the genome and encodes non-structural protein 1 (NS1) and 2 (NS2). ORF2 is located in the right part of the genome and encodes the structural proteins VP1 and VP2. The middle region of the genome also contains a 500 bp sequence [9, 10] . The non-structural proteins NS1 and NS2 are associated with viral replication. NS1 functions as a nickase and helicase and can form a covalent Edited by William Dundon.

Ling Zhou, Qinghai Tang, Lijun Shi, and Miaomiao Kong have contributed equally to this work. bond at the 5 0 end of DNA. Structural proteins VP1 and VP2 compose the nucleocapsid of CPV-2 [11] .

In 2014 in China, many acute cases of infectious diarrhoea were reported in dogs. To identify the pathogen, we collected and examined faecal samples from symptomatic dogs. We report here that the disease is associated with a strain of parvovirus, and phylogenetic characterization of the causal isolate showed that isolates from China have long branch lengths; this may indicate an extensive process of accumulation of mutations and substitutions in Chinese lineages. This study provides the basis for further exploration of the CPV-2 variation, the selection of vaccines, and the effective prevention and control of CPV-2 infections.

Faeces were collected from dogs exhibiting acute haemorrhagic enteritis. The faecal samples (10 g) were immersed in 10 mL of sterile phosphate-buffered saline (PBS) with antibiotics, and the suspension was centrifuged at 10,0009g for 15 min. The supernatant was filtered and frozen at -80°C.

To determine whether CPV-2 was present in the samples, the filtered supernatants were added to the colloidal gold test strips using the BioNote Rapid Test Kit (Korea, Bio-Note). The samples were mixed with buffer and then tested for the antigen at room temperature according to the manufacturer's instructions. The negative sample was prepared for PCR as follows.

The PCR assay was performed in a 20-lL reaction volume that included 2.4 lL of template, 1.8 lL each of VP2-F and VP2-R (Table 1) , 2.0 lL of dNTPs, 1.0 lL of KOD FX Neo polymerase (1 U/lL) (TOYOBO Biotechnology Company, Shanghai, China), 5.0 lL of 2 9 PCR buffer for KOD FX Neo (TOYOBO Biotechnology Company, Shanghai, China), and sufficient ddH 2 O to increase the volume to 20 lL. Amplification was carried out in a pre-heated thermocycler (Applied Biosystems 2720 Thermal Cycler) as follows: one cycle at 94°C for 5 min; followed by 35 cycles at 94°C for 45 s, 50°C for 60 s, and 72°C for 60 s; and a final extension at 72°C for 10 min. Amplicons were detected by electrophoresing 10-lL aliquots in 2 % agarose gels in 1 9 TAE [40 mM Tris-acetate (pH 8.0), 1 mM EDTA].

Identification of CPV-2 with cell culture and IPMA Feline kidney cell line F81 was obtained from the American Type Culture Collection, USA and was used to isolate viruses from clinical samples and to observe cytopathic effects associated with viral replication. The treated samples were inoculated at a concentration of 1 mL per 25 cm 2 flask in a cell monolayer. After adsorbing for 1 h at 37°C, the inoculum was removed, DMEM with 2 % foetal bovine serum was added, and the cells were again incubated at 37°C. Cell cultures were observed daily for 4-5 days to monitor the appearance of cytopathic effects (CPE). Subsequently, an IPMA was conducted, and the other flask was frozen at -20°C and submitted to further passages following the same procedure after freezing and thawing for three times, until the eventual appearance of CPE [3, 12] .

Cloning the full-length genomic sequence of CPV-2 DNA and RNA were extracted from the homogenized samples (faeces from infected dogs) with the TIANamp Virus DNA/RNA Kit (Beijing TIANGEN Biotech Company, Beijing, China) according to the manufacturer's protocol. DNA and RNA were analysed by PCR or RT-PCR, and the presence of other potential pathogens, such as canine distemper virus (CDV), canine adenovirus (CAV), canine coronavirus (CCV), and canine rotavirus (CRV), was investigated. Table 1 summarizes the sequences of primers used to amplify the virus, and the PCR was conducted as previously described [13] [14] [15] . These primers were the same set used in Calderon et al. [13] , Chaturvedi et al. [14] , and Pratelli et al. [15] DNA fragments were then cloned into the pMD18-T sample vector (TaKaRa Biotechnology Co., Ltd. Japan). The same set of PCR primers was used to sequence the full-length genome of the isolate NY-14. Sequencing reactions were performed by a commercial corporation (Huada, Beijing).

Sequences were aligned using Clustal X software version 2.1 [16] . After the first alignment, the sequences were manually edited to maintain the reading frames using the SE-Al programme version 2.0 (http://evolve.zoo.ox.ac.uk/ software/). The following CPV-2 references were used for the phylogenetic analyses: Maximum likelihood trees and bootstrap values were obtained using PhyML software [17] . The KHY model and gamma distribution (C) was selected according to the likelihood ratio test (LRT) using jModelTest software [18, 19] . The resulting trees were visualized and edited using FigTree (http://tree.bio.ed.ac.uk/software/figtree/).

To determine the extent of recombination in the CPV-2 sequences, RDP v.4 software [20] was used for analyses, which utilizes a collection of methods. An excellent and detailed explanation of each method implemented in the RDP programme can be found in the user's manual (http:// darwin.uvigo.es/rdp/rdp.html).

The ill dogs showed acute infectious diarrhoea, and the faecal samples were bloody stools (Fig. 1a) . The colloidal gold test strip showed positive (Fig. 1b1) and negative ( Fig. 1b2) results, while the PCR results were both positive (Fig. 1c) .

One viral strain was isolated from the faeces of dogs that tested negative in the rapid test but positive using PCR. CPE was observed in F81 cells after three passages. CPV-2-infected cell cultures showed vast regions of cell detachment (Fig. 2) . The IPMA was conducted with positive sera (Fig. 3) .

PCR amplification showed a fragment of 1755 bp corresponding to the VP2 gene of CPV-2. PCR amplification detected the presence of CPV-2 but was negative for other pathogens in the samples analysed. The sequencing results showed that the VP2 gene of CPV-NY-14 is 1755 bp long with no insertions or deletions in the coding region. The results were once again positive for CPV-2. Three experimental animals were inoculated with NY-14 and exhibited clinical symptoms at 14 days post inoculation. CPV-2 was re-isolated from the faeces of the experimental dogs. The whole genome was sequenced, and the CPV-NY-14 complete genome structure was analysed and is shown in Fig. 4 .

A maximum likelihood tree was constructed using the VP2 sequence and full-length genomes of CPV-NY-14 and all available CPV-2 genomes. Figure 4 shows that all the CPV-2 sequences from China clustered in a group, and the CPV-NY-14 isolate is also located in the Chinese clade. Interestingly, located at the base of the Chinese clade is one isolate from Russia (JN033694) sampled in 1993. All sequences used to construct the above-mentioned tree were absent for recombination.

To further examine CPV-2 isolates, we used a larger dataset composed only of VP2 protein sequences. The maximum likelihood tree again showed a monophyletic cluster containing all CPV-2 isolates from China. CPV-NY-14 belongs to the CPV-2a type, and it was inferred that CPV-2a is still circulating in China and is also the main agent of acute haemorrhagic enteritis in dogs.

In addition, both trees showed that isolates from China have long branch lengths; this may indicate an extensive process of accumulation of mutations and substitutions in Chinese lineages. Since it emerged in 1978, canine parvovirus has spread in the domestic and wild canine population, where it is continuously evolving into new adaptive viral variants. The variability and the intrinsically high mutation rate of the CPV-2 genome allowed the diversity of CPV-2 to rapidly increase as the virus spread through canine populations [21] . Because some viral variants replicate more successfully than others, the virus population changes over time [22] . CPV-2a and CPV-2b are the best examples of variants with a fast spread and replacement capacity since the emergence of CPV-2 [9, 23] . During 2014, there were many cases of acute infectious diarrhoea in dogs. Some samples tested negative for the CPV-2 antigen by a colloidal gold test strip, while they were positive for CPV-2 by PCR. One viral strain was isolated from the faeces of dogs that tested negative in the colloidal gold test strip but positive with PCR. CPV-NY-14, which was described in this study, was found in dogs that exhibited acute haemorrhagic enteritis. Based on the PCR results, the only virus consistently detected in the samples was CPV-2. Some samples tested with the serology method (colloidal gold test strip) failed to detect the virus, while PCR was able to detect CPV-2, which is a valid and important objective. This may be due to differences in the test sensitivity (PCR is more sensitive than the serological test) or due to a change in the viral genome that leads to a major change in the type of the antigen (viral protein) and the failure of the antibodies to recognize the viral antigen in the serology test. This question needs to be addressed by future experiments. The isolated virus was then cultivated in F81 cells, and we next inoculated adult dogs with the isolated CPV-2. These animals presented clinical conditions characteristic of CPV-2 infection. PCR was used to demonstrate that CPV-2 was the agent causing the disease in these experimental animals.

CPV-2 is pandemic, and the frequencies of the different antigenic types of CPV-2 vary in different countries. In Uruguay, for example, CPV-2c is the major epidemic strain [24] . In the USA and southern Africa, CPV-2b is the main viral type that causes most outbreaks of CPV-2 infection [23, 25] . In the UK, both CPV-2a and CPV-2b are present, and Germany and Spain have similar frequencies of isolation [26, 27] . In India, the major epidemic strains are CPV-2a and CPV-2c. The new strains have also replaced the original one in Taiwan and Italy, although the CPV-2a strain is still predominant there [28, 29] . However, CPV-2 still causes many cases of acute infectious diarrhoea in dogs in Japan [30] .

In this study, molecular phylogenetic analysis of CPV-NY-14 and other CPV-2 isolates in GenBank revealed that CPV-NY-14 is closely related to isolates s5, SC02-11, LZ1, LZ2, and nj01-06, especially isolate LZ2. CPV-2a still causes many cases of acute infectious diarrhoea in dogs in China. The sequence analysis showed that there is little variation, and this is important in choosing a vaccine to prevent this disease. Given the results obtained in the current study, the continuous surveillance of CPV-2 in China is imperative for determining whether CPV-2a will colonize and spread into new territories.

Diarrhea in puppies: parvovirus-like particles demonstrated in their feces

Clustal W and Clustal X version 2.0

Proc. Natl. Acad. Sci. USA

Acknowledgments This work was partly supported by the Agricultural Science and Technology Innovation Programme of China (ASTIP-IAS15).

Conflict of interest The authors declare that they have no competing interests.