Genomic analysis of canine pneumoviruses and canine respiratory coronavirus from New Zealand

Abstract

ABSTRACT Aims To isolate canine respiratory coronavirus (CRCoV) and canine pneumovirus (CnPnV) in cell culture and to compare partial genomic sequences of CRCoV and CnPnV from New Zealand with those from other countries. Methods Oropharyngeal swab samples from dogs affected by canine infectious respiratory disease syndrome that were positive for CnPnV (n = 15) or CRCoV (n = 1) by virus-specific reverse transcriptase quantitative PCR (RT-qPCR) in a previous study comprised the starting material. Virus isolation was performed in HRT-18 cells for CRCoV and RAW 264.7 and Vero cells for CnPnV. The entire sequence of CnPnV G protein (1,266 nucleotides) and most (8,063/9,707 nucleotides) of the 3’ region of CRCoV that codes for 10 structural and accessory proteins were amplified and sequenced. The sequences were analysed and compared with other sequences available in GenBank using standard molecular tools including phylogenetic analysis. Results Virus isolation was unsuccessful for both CRCoV and CnPnV. Pneumovirus G protein was amplified from 3/15 (20%) samples that were positive for CnPnV RNA by RT-qPCR. Two of these (NZ-048 and NZ-049) were 100% identical to each other, and 90.9% identical to the third one (NZ-007). Based on phylogenetic analysis of the G protein gene, CnPnV NZ-048 and NZ-049 clustered with sequences from the USA, Thailand and Italy in group A, and CnPnV NZ-007 clustered with sequences from the USA in group B. The characteristics of the predicted genes (length, position) and their putative protein products (size, predicted structure, presence of N- and O-glycosylation sites) of the New Zealand CRCoV sequence were consistent with those reported previously, except for the region located between open reading frame (ORF)3 (coding for S protein) and ORF6 (coding for E protein). The New Zealand virus was predicted to encode 5.9 kDa, 27 kDa and 12.7 kDa proteins, which differed from the putative coding capacity of this region reported for CRCoV from other countries. Conclusions This report represents the first characterisation of partial genomic sequences of CRCoV and CnPnV from New Zealand. Our results suggest that the population of CnPnV circulating in New Zealand is not homogeneous, and that the viruses from two clades described overseas are also present here. Limited conclusions can be made based on only one CRCoV sequence, but the putative differences in the coding capacity of New Zealand CRCoV support the previously reported variability of this region. The reasons for such variability and its biological implications need to be further elucidated.