PIV3 is an enveloped, single-stranded negative sense RNA virus within the respirovirus genus of the Paramyxoviridae family, which causes respiratory diseases in many host species (Maidana et al., 2012) . The respirovirus genus includes human parainfluenza virus types 1 and 3 (HPIV1 and HPIV3, respectively), bovine parainfluenza virus type 3 (BPIV3) and Sendai virus (Fauquet et al., 2005) . BPIV3 is an important pathogen associated with bovine respiratory disease complex (BRDC) which is considered the most significant disease associated with feedlot cattle worldwide (Horwood et al., 2008) . Parainfluenza virus type 3 (PIV3) is one of the most important viral respiratory pathogens for humans and for many animals, but goat infection has been rarely reported. Starting in Aug 2013, goats in the Jiangsu and Anhui provinces of eastern China suffered severe respiratory diseases. In order to identify the causative agent, numerous related pathogens were tested with RT-PCR or PCR. A unique PIV3 strain was detected in most of the clinical nasal swabs or serum samples. The virus was isolated on MDBK cells and characterized by RT-PCR, nucleotide sequence analysis and hemagglutination test. The entire M and F gene coding regions, HN, 5 0 -UTR-N and L gene fragments were amplified using pairs of degenerate primers. Nucleotide, amino acid sequence alignments and phylogenetic analyses based on these genes indicated that the goat-derived PIV3 strain was distinct from previously reported BPIV3 genotypes and HPIV3 strains. The novel isolate, named JS2013, might be a potentially new member of the respirovirus genus. Goats were experimentally infected with JS2013 culture. The virus-inoculated goats displayed coughing and nasal discharges that were related to respiratory diseases. Viremia and virus shedding were detected during 4-10 days post-inoculation (dpi). Virus-specific HI antibodies became positive from 14 dpi. This is the first report of the detection of PIV3 from Chinese goat herds and genetic and pathogenetic characterization of the novel goat-derived PIV3.
ß 2014 Published by Elsevier B.V.
BPIV3 infection causes severe bronchopneumonia and secondary bacterial infections in instances of high stress, such as transportation, weather change and feedlot situations (Haanes et al., 1997) . The main clinical signs of BPIV3 infection are coughing, anorexia, pyrexia, nasal and ocular discharges, dyspnea and sometimes diarrhea (Zhu et al., 2011) . Three BPIV3 genotypes, A (BPIV3a), B (BPIV3b) and C (BPIV3c) have been described and BPIV3a and BPIV3c have been found in Chinese cattle (Zhu et al., 2011) . However, reports of PIV3 infection in goat/sheep are limited and the genetic characteristics of the virus have not been determined (Elazhary et al., 1984; Lamontagne et al., 1985; Lyon et al., 1997; Taylor et al., 1975) . Furthermore, nothing is known about the circulation of PIV3 in goats/ sheep in China and the relationship between virus infections and clinical diseases, as well as the relationship of small ruminants' PIV3 with BPIV3 and HPIV3.
In 2013, goats in Jiangsu and Anhui provinces of China suffered severe respiratory diseases. We used RT-PCR or PCR to test for various potential agents from clinical samples and PIV3 viral RNA was detected in most samples. Virus isolation was performed and was further confirmed by RT-PCR, sequencing and HA assay. Genetic and phylogenetic analysis was performed and suggested that the virus might be a new member of PIV3.
From Aug 2013 to Mar 2014, goats of many farms in Jiangsu and Anhui provinces in eastern China suffered severe respiratory diseases, with high morbidity and varied mortality. These farms were fattening farms, which purchased weaning goats from several small backyard farms for fattening until slaughter time. Most of the diseased goats (3-6 months) showed respiratory syndrome (coughing, nasal discharges, dyspnea), anorexia emaciation and some of them suffered diarrhea. The diseased goats usually showed no obvious body temperature changes. Antibiotic treatment had no apparent effect, although a few diseased animals recovered within 10-15 days. Necropsy of the diseased goats showed mainly seroperitoneum, pathological lesions in the lung and pancreas, white punctiform tuberculum in epiploica or mesenterium and enlarged liver ( Fig. 1) . 
Seventy-seven nasal swab/serum samples were collected from eleven farms in Jiangsu and Anhui provinces during Oct 2013 and Mar 2014 (Table 1 ). Samples were stored at À70 8C for RT-PCR and virus isolation. All sampled goats showed clinical signs of pneumonia when samples were collected.
Total RNA was extracted from nasal swabs/serum samples using TRIzol reagent (Invitrogen) according to the manufacturer's instruction. DNA was extracted from nasal swab/serum samples using phenol/chloroform. PIV3-specific RT-PCR amplification was carried out with Easyscript one-step RT-PCR supermix (Transgen, Bio, Inc.) in a 20 ml reaction mixture containing 2Â R-Mix buffer, 20 pM of each primer (MF1/MR1, Table 2), 0.5 ml of E-Mix and 4 ml extracted RNA. The reaction was run in a thermocycler (Mjmini, BIO-RAD) with the following program: 45 8C 40 min for reverse transcription; denaturation at 94 8C for 5 min, 35 cycles composed of denaturation at 94 8C for 30 s, annealing at 50 8C for 30 s and extension at 72 8C for 45 s; and was terminated with a final extension of 10 min at 72 8C. Amplification products were detected by electrophoresis in 1.2% agarose gels.
Other related pathogens, ovine respiratory syncytial virus (ORSV) (Eleraky et al., 2001) , peste des petits ruminants virus (PPRV) (Mao et al., 2010) , pestivirus (Vilcek et al., 1994) , foot and mouth disease virus (FMDV) (Xu et al., 2013) , bovine coronavirus (Amer et al., 2013) and Mycoplasma (Bascunana et al., 1994; McAuliffe et al., 2003) were also tested by previously reported methods.
Positive nasal swab samples were centrifuged at 12,000 rpm for 20 min, the supernatants were filtered through 0.22 mm filter (Millipore) and inoculated onto Madin-Darby bovine kidney (MDBK) cell monolayers cultured in 24-well cell culture plates. Positive serum samples were used for inoculation directly after centrifugation at 12,000 rpm for 20 min. The MDBK cultures were observed for 5-7 days. Cell cultures were harvested and passaged 4 more times. Each viral stock was stored at À70 8C for RT-PCR detection. Virus isolation was confirmed by RT-PCR using primers MF1/MR1 and HNF/HNR, and the isolates were further identified by hemagglutination test and genetic analysis.
HA test was performed with 0.25% guinea pig erythrocyte suspension in 96-well V-bottom plates. Equal volumes (0.05 mL) of serial twofold dilutions of the virus isolates and erythrocytes were mixed, shaken and allowed to sediment for 45 min at 37 8C. The end point was the last well in which the guinea pig erythrocytes did not agglutinate.
To confirm the specificity of JS2013, the HN, 5 0 -UTR-N and L gene fragments and entire M and F gene coding regions were amplified from positive samples using the degenerate primers listed in Table 2 , purified with a purification kit (Axygen Bio, Inc.), cloned into pMD18-T vector (Takara Bio, Inc.), and transformed into Escherichia Table 2 Primers used for PIV3 detection and gene amplification.
Sequences ( coli DH5a competent cells. Three positive clones for each PCR products were sequenced. Each sequence was determined in both directions. The nucleotide sequences and corresponding predicted amino acid sequences were edited by Editseq (DNASTAR Inc., Madison, WI) and subjected to MegAlign (DNASTAR Inc., Madison, WI) and Blast analyses. Multiple sequence alignment was done by using Clustal X 1.83 (Thompson et al., 1997) , together with other reference sequences of previously identified BPIV3 and human parainfluenza virus type 3 (HPIV3). After determining the percentages of sequence identity among different PIV 3 strains, phylogenetic trees were generated with the distance-based neighbor-joining (NJ) method by using MEGA 4.0.2 software (Tamura et al., 2007) . The robustness of the phylogenetic trees was determined by bootstrap resampling analysis carried out on 1000 replicates.
Eight goats were randomly divided into two groups of four each and housed separately. Group CC was intranasally inoculated with JS2013 (3 mL/goat, HA titer of 256). Group NC was inoculated with PBS and used as negative control. After challenge, all animals were monitored for 35 days. Rectal temperatures and clinical signs were observed daily and clinical observations were recorded. The serum samples were collected from all animals at 0, 4, 7, 10, 14, 21, 28 and 35 dpi. Nasal swabs were taken at 0, 4, 7, 10, 14, 21 dpi. RNA was extracted from serum and nasal swab samples at 0, 4, 7, 10, 14, 21 dpi; RT-PCR was performed as mentioned above (Section 2.3) to detect viremia and virus shedding. Hemagglutination inhibition (HI) assay was performed to determine PIV3-specific antibodies in serum samples of 0, 14, 21, 28 and 35 dpi.
For the HI assay, serial two-fold dilutions of sera were allowed to react with 4 HA units of PIV3 for 30 min at 37 8C, followed by addition of guinea pig erythrocytes. After 40 min incubation at 37 8C, HI titers were calculated as the highest serum dilution that inhibited 4HA unit antigen. If HI titer was less than or equal to 3 log 2 , it was considered negative; if HI titer was more than or equal to 4 log 2 , it was considered positive.
The collection of nasal and serum samples and the artificial infection experiment were performed in strict accordance with the guidelines of Jiangsu Province Animal Regulations (Government Decree No. 45).
M gene fragments were detected in 43 of 77 (56%) nasal swab/serum samples by RT-PCR, and all tested farms had positive samples (Table 1) . The amplified products were purified, cloned and sequenced. Blast search revealed that the sequences shared the highest identity of 82% with Japanese BN-CE strain of bovine PIV3. All amplified M gene products from the different farms we tested had 100% identity with each other, indicating they might be the same PIV3 strain.
In addition, ovine respiratory syncytial virus was detected in 5/77 samples (farm A: 2, C: 1 and G: 2). Mycoplasma ovipneumoniae was detected in 7/77 samples (farm D: 3, E: 2 and I: 2). PPRV, bovine coronavirus, pestivirus and FMDV were negative.
After culture and passage in MDBK cells, one virus isolate was obtained from a nasal swab sample (isolation from serum samples was not successful). The isolate uniformly produced cytopathic effect from the 2nd passage. The 329 bp M gene and 674 bp HN gene fragments were amplified from virus cultures by RT-PCR (data not shown). In addition, these two genes of different passages did not change (100% identity). The 3rd, 4th and 5th passages of isolated viral stock showed titers of 32, 128 and 256 by hemagglutination test, respectively. These results suggested successful viral isolation and amplification; the novel virus was named JS2013.
To confirm the specificity of JS2013, the entire M and F gene coding regions, HN, 5 0 -UTR-N and L gene fragments were amplified and sequenced (GenBank accession No. M: KJ850331, HN: KJ767525; N: KJ843145; F: KM279938; L: KJ843144). Blast research and MegAlign analysis indicated that the M, F, HN, N and L sequences shared the highest nucleotide identity of 80.2% with SF strain, 77% with SF strain, 80% with BN-CE strain, 82.5% with Q5592 strain and 83.5% with Q5592 strain, respectively. They possessed lower identity with HPIV3 sequence. In brief, they showed nucleotide identity of 77.5-80. 2%, 72.5-76.9%, 74.1-80%, 78.4-82.5% and 80.7-83 .5% with the selected representative sequences of different BPIV3 genotypes and HPIV3 strains, respectively (Table 3 ). The F gene possessed the lowest identity and L gene was the highest with other compared sequences (Table 3) .
The corresponding amino acid sequences had higher identities than their nucleotide sequences (Table 3) . Similarly, they had lower identity with HPIV3 sequences. Among them, the L protein was the highest (93.9-95.6%), while F had the lowest (75.2-79.1%), indicating that F was more variable amongst different strains.
Phylogenetic reconstructions were compiled using a 1056 nucleotide region of M, a 674 nucleotide region of HN, a 848 nucleotide region of N and a 1623/1626 nucleotide region of F, respectively. Phylogenetic analysis based on M and F gene sequences showed similar grouping results, while HN-and N-based phylogenetic trees had similar structures, in which JS2013 located between BPIV3 and HPIV3 clusters. All of them showed that JS2013 formed a unique branch, which clearly differed from HPIV3 and other reported BPIV3 genotypes (Fig. 2) . The results, together with nucleotide and amino acid alignment results, indicated that JS2013 might represent a novel member of the respirovirus genus in the Paramyxoviridae family.
Goats were intranasally infected with JS2013 to preliminarily assess the pathogenicity of JS2013. After challenge, rectal temperatures of all animals remained normal, with no significant differences observed between the inoculated and non-treated groups (data not shown). NC  D1  -----------2  2  3  3  2  D2  -----------2  2  2  3  2  D3  -----------3  2  3  2  3  D4  -----------2  3  2  3 As shown in Table 4 , three out of the four virusinoculated goats (D5, D7 and D8) displayed coughing and mild to severe nasal discharges during 4-10 dpi. Viremia was detected in all goats of group CC during 4-10 dpi. Virus shedding was confirmed by RT-PCR from nasal swabs and M gene sequencing indicated that it was same as the originally infected virus. In addition, virus specific HI antibodies were positive from 14 dpi; 3/4 goats of CC group showed !2 fold change of HI titers between 0 and 14 dpi, and all goats showed !2 fold HI titers increase in 28 dpi.
Respiratory diseases have been causing major economic loss to feedlot goat farms. In feedlot cattle, BRDC is considered the most significant disease, and BPIV3 plays an important role in BRDC development. Viral infection plus stress conditions and other secondary or co-infections result in severe clinical observations (Haanes et al., 1997; Maidana et al., 2012) . However, there are few reports about related goat viral pathogens. Most previous reports regarding PIV3 in goat/sheep were published in the 1970s-1980s and only serological tests were done; thus, related virus genome characteristics and virus isolations had not been identified (Afshar and Terlecki, 1979; Elazhary et al., 1984; Kennedy-Stoskopf et al., 1983; Lamontagne et al., 1985; Lyon et al., 1997; Taylor et al., 1975) . In the present cases, the goats were imported from different places and transport stress might trigger disease development. The diseased goats were treated with different antibiotics and did not show significant changes in clinical signs. The diseased and dead animals showed similar pathological characteristics. In the primary laboratory examination, several potential pathogens were tested. Most of the samples (56%) showed PIV3 positive results, and sequencing results indicated that viruses from different farms shared 100% identity, indicating the contribution of this novel PIV3 to this clinical disease. ORSV and M. ovipneumoniae were also detected from some samples, indicating different co-infection situations in different farms and this might result in the difference of recovery rate between farms; however, ORSV and M. ovipneumoniae were not the major causes of the goat infections. PIV3 was successfully isolated from nasal swabs and was intranasally infected into goats to assess its pathogenicity. From 4 to 10 dpi, viremia and virus shedding could be detected in all virus-inoculated goats; and 3/4 goats displayed coughing, with mild to moderate nasal discharges. Virus-specific HI antibodies were positive from 14 dpi (Table 4 ). These results further confirmed the goat-derived PIV3 strain had pathogenicity in goats. Of course, other potential unknown pathogens might also be present, requiring further studies.
Three BPIV3 genotypes have been described. BPIV3a was relatively older and widespread. BPIV3b and BPIV3c were found in Australia (Horwood et al., 2008) and Chinese cattle (Zhu et al., 2011) , respectively. Blast searches of the 329 bp M region of JS2013 showed the highest identity of only 82% with Japanese BN-CE strain, indicating that JS2013 was distinct from the previously reported BPIV3 strains. Phylogenetic analysis based on this region gave the same result (data not shown). It was hypothesized that JS2013 might represent a potentially new group of PIV3. To confirm this hypothesis, several pairs of degenerate primers were designed to amplify the entire M and F coding region, partial N, HN and L gene fragments. The sequences were used for genetic and phylogenetic analysis. Similarly, the highest nucleotide identities of M, F, N, HN and L genes with reported strains in Genbank were only 76.9-83.5% (Table 3 ). The new PIV3 strain shared slightly higher identity with BPIV3 than HPIV3. As shown in Table 3 , the F and HN genes and proteins possessed relatively lower identities while the L gene and protein had the highest identity with other sequences, suggesting that F and HN were highly varied among different strains. Considering that HN protein was an immunodominant protein of PIV3, HN-based serological methods and vaccine should be developed in the future. Zhu et al. had proposed the use of PIV3 M and HN genes as appropriate targets for genotyping (Zhu et al., 2011) . F gene was also used for other paramyxovirus genotyping. Therefore, M, F, HN and N genes were used for phylogenetic analysis in this study and they showed similar grouping results (Fig. 2) : JS2013 was clearly differed from HPIV3 and other reported BPIV3 genotypes. The results presented here strongly indicate that JS2013 might represent a novel member of PIV3 in the respirovirus genus and the name caprine parainfluenza virus type 3 is proposed (Fig. 2) . Full-length genome sequencing is currently underway for detailed understanding the genomic characteristics of this new virus.
To our knowledge, this is the first report about the detection and isolation of PIV3 in Chinese goat flocks with respiratory diseases. A broader epidemiological investigation of PIV3 in goats should be carried out and serological methods should be established for detecting specific antibodies against PIV3 in goat as well as sheep herds. The pathogenicity of this new isolate needs further detailed studies and related preventative techniques should be developed. This would be a good beginning for detailed research on respiratory diseases of goat/sheep in China.