key: cord-0005409-8qja4j9h authors: Li, Weike; Li, Tiansong; Liu, Yuxiu; Gao, Yuwei; Yang, Songtao; Feng, Na; Sun, Heting; Wang, Shengle; Wang, Lei; Bu, Zhigao; Xia, Xianzhu title: Genetic characterization of an isolate of canine distemper virus from a Tibetan Mastiff in China date: 2014-04-02 journal: Virus Genes DOI: 10.1007/s11262-014-1062-z sha: 8b83b7a705bdd5012b017cfe1c070da3e130e872 doc_id: 5409 cord_uid: 8qja4j9h Canine distemper (CD) is a highly contagious, often fatal, multisystemic, and incurable disease in dogs and other carnivores, which is caused by canine distemper virus (CDV). Although vaccines have been used as the principal means of controlling the disease, CD has been reported in vaccinated animals. The hemoagglutinin (H) protein is one of the most important antigens for inducing protective immunity against CD, and antigenic variation of recent CDV strains may explain vaccination failure. In this study, a new CDV isolate (TM-CC) was obtained from a Tibetan Mastiff that died of distemper, and its genome was characterized. Phylogenetic analysis of the H gene revealed that the CDV-TM-CC strain is unique among 20 other CDV strains and can be classified into the Asia-1 group with the Chinese strains, Hebei and HLJ1-06, and the Japanese strain, CYN07-hV. The H gene of CDV-TM-CC shows low identity (90.4 % nt and 88.9 % aa) with the H gene of the classical Onderstepoort vaccine strain, which may explain the inability of the Tibetan Mastiff to mount a protective immune response. We also performed a comprehensive phylogenetic analysis of the N, P, and F protein sequences, as well as potential N-glycosylation sites and cysteine residues. This analysis shows that an N-glycosylation site at aa 108-110 within the F protein of CDV-TM-CC is specific for the wild-type strains (5804P, A75/17, and 164071) and the Asia-1 group strains, and may be another important factor for the poor immune response. These results provide important information for the design of CD vaccines in the China region and elsewhere. Analysis of CDV strains from various animal samples has demonstrated an important relationship with the H gene/ glycoprotein, which has changed by genetic/antigenic drift. As the key protein for CDV, H is used for attachment to cell receptors as the first step of infection and mediates adequate host immune response [9] . The H protein is considered to have the highest antigenic variation and can reflect genetic changes in comparative studies of CDV strains [10] [11] [12] [13] . This variation may affect neutralization-related sites with disruption of important epitopes. Analysis of CDV strains from different animal species and geographical settings has revealed that the geographic pattern is an important factor in the genetic/antigenic drift affecting the H gene/glycoprotein of CDV [14] [15] [16] [17] [18] [19] [20] . Therefore, the H gene may be used for identification and phylogenetic classification of CDV strains, which have been identified into seven major genetic lineages, namely America-1 and -2, Asia-1 and -2, Arcticlike, Europe, and wild-life [21, 22] , as well as an indication of the antigenic response of the virus. Three other proteins, the nucleocapsid (N) protein, the phosphoprotein (P) protein, and the fusion (F) protein, also have important roles for CDV and could provide additional sources of antigenic variability among strains. The N protein has immunosuppressive properties and is the major component of the CDV virion. The N-terminal domain of the N protein is generally well conserved, while the C-terminal end is poorly conserved and is considered hypervariable. The C-terminal tail of the N protein also contains the majority of its phosphorylation sites and antigenic sites [23, 24] . During active infection, antibodies made against the N protein in the host are predominant and account for most of the complement-fixing antibody [25, 26] . The P protein is relatively well conserved and plays a vital role in transcription and replication [27] . This protein is an essential component of the viral RNA phosphoprotein complex (vRNAP) [28] and also function as a chaperone for the N protein. The F protein is a type I integral membrane protein that mediates viral penetration by fusion between the virion envelope and the host cell plasma membrane at neutral pH. It is synthesized as an inactive precursor, F0, and must be proteolytically cleaved to produce the functionally active fusion protein, which consists of disulfide-linked F1 and F2 polypeptides [29] . Like the H protein, the F protein has high antigenic variation. In this study, the wild-type CDV-TM-CC strain was isolated from the spleen of a 1-year-old Tibetan Mastiff that developed clinical signs of CD after having received all standard vaccines. To determine whether this occurrence may be explained by variations in specific nucleotide or amino acid residues of the CDV circulating in China, we sought to genetically characterize the CDV-TM-CC strain. VerodogSLAM cells constitutively expressing the CDV receptor dog signaling lymphocyte activation molecule (SLAM) were cultured in Dulbecco's modified Eagle medium (DMEM; Gibco) supplemented with 10 % heatinactivated fetal bovine serum (FBS) with an additional 8 lg of G418 per ml. The wild-type CDV-TM-CC strain was originally isolated from spleen homogenate (10 % w/v suspension) from a Tibetan Mastiff that succumbed to naturally infection. Virus was propagated in VerodogSLAM cells and stored at -80°C. Total RNA was prepared from VerodogSLAM cells infected with CDV-TM-CC according to the manufacturer's instructions (Total RNA Kit I, OMEGA). The reverse transcription reactions were performed using M-MLV Reverse Transcriptase (Invitrogen) with oligo d(T) and random primers. According to the complete consensus genomic sequence of CDV (GenBank), two sets of primers were designed to amplify the entire genome (Oligo6.0 design software), as shown in Table 1 . Sequences were assembled and compared using DNA sequence analysis software (DNAStar), and the complete consensus genomic sequence was determined. PCR amplification was carried out using Phusion High-Fidelity DNA Polymerase (New England BioLabs). Clones (amplicons emcompassing the full-length CDV-TM-CC genome) were obtained A GenBank number is provided for each of the strains of CDV that were compared with CDV-TM-CC in this study. The geographical location of strain isolation and the species/organ of isolation are also indicated, as well as the clade into which the strains are categorized Virus Genes (2014) 49: 45-57 47 from thirty RT-PCR reactions using CDV-specific oligonucleotides. To genetically characterize the CDV-TM-CC strain, the deduced amino acid sequence was compared to F and H gene fragments of the variant field isolates shown in Table 2 . A phylogenetic tree was constructed based on the deduced amino acid sequences in supplementary Table 1 using MEGA 5.0, and multiple sequence alignment was carried out using ClustalW. Statistical significance of the phylogeny was estimated by bootstrap analysis over a 1,000 pseudoreplicate data set. The wild-type CDV-TM-CC strain was isolated from the spleen of a 1-year-old Tibetan Mastiff in Jilin province that had succumbed to CD after having received all standard vaccines (6 weeks first immunization, 8 weeks second immunization, 10 weeks third immunization with Distemper, adenovirus type 2, parvovirus, parainfluenza quadruple vaccine; Canine coronavirus disease killed virus vaccine portion, USA). The virus was propagated in VerodogSLAM cells and the virulence of the strain was confirmed (data not shown). To identify sequence features that may explain the failure of the vaccine strain to protect the dog against CD, we sequenced the entire genome, using two sets of overlapping primers (Table 1) . Within the CDV genome, the H gene is a major causative disease determinant and also has one of the highest rates of mutation. Consequently, the phylogenetic relationship of CDV strains is often based on the deduced amino acid sequence of the H protein. The H gene of the CVT-TM-CC strain has 1,824 nucleotides and the inferred protein sequence has 607 amino acids, similar to the other CDV strains. Amino acid analysis of the H protein from CDV-TM-CC and 20 other CDV strains in GenBank (Table 2) identified seven clades of CDV strains (America-1, America-2, Asia-1, Asia-2, Europe, Arctic-like, and Europe wildlife). CDV-TM-CC was classified into the Asia-1 group with the strains CYN07-hV (Japan), Hebei (China), and HLJ-06 (China) (Fig. 1 (Fig. 2a, b) . Glycosylation is an important factor in determining the antigenicity of many proteins [30] . Prediction of the glycosylation sites of the H gene (http://www.cbs.dtu.dk/ser vices/NetNGlyc/, NetNGlyc 1.0 Servera) identified a total of eight potential glycosylation sites at positions 19 Notably, the 309-311 N-glycosylation site is specific for virulent strains [14, 18] with the exception of A75/17. The 584 N-glycosylation site is specific for the Asian-1 strains, suggesting that it was acquired later [18, 20] . CDV-TM-CC has both of these predicted glycosylation sites, which could explain its virulence properties. Phylogenetic analyses of the amino acid sequence of the N and P proteins To determine whether the conservation of CDV-TM-CC also extends to other proteins within the virus, we assessed the similarity of the N and P proteins. Consistent with the results for the H protein, the homology of the deduced CDV-TM-CC amino acid sequence of the N protein to the Asia-1 strains (CYN07-hV, HLJ1-06, and Hebei) was high with 98.7-98.9 % identity, as shown in Fig. 4 . The N protein sequence of CDV-TM-CC also showed 98.1 % identity with the Asia-2 group (strains M25CR, 007Lm, 011C, 50Con, and 55L), and 97.5 % identify with the Onderstepoort strain. Moreover, CDV-TM-CC had high similarity (98.5, 98.7, and 97.9 % identity) with wild-type strains 164071, A75/17, and 5804P. The lowest homology of the CDV-TM-CC N protein sequence (96.6-96.8 % aa identity) was found with Arctic-like strains CDV3, Shuskiy, and Phoca-Caspian-2007. This relatively high similarity between the N protein of CDV-TM-CC and other CDV strains is consistent with the generally high conservation among N proteins. The phylogenetic relationship of CDV-TM-CC based on the deduced amino acid sequence of the P protein was also analyzed (Fig. 5 ). Similar to the results for the H protein, CDV-TM-CC classified into the Asia-1 group, but was in a separate branch from the classical Onderstepoort vaccine strains. These results verify the classification of CDV-TM-CC as an Asia-1 group virus. The signal peptide is a short amino acid sequence at the N-terminus of the majority of newly synthesized proteins that are destined towards the secretory pathway and is a highly divergent region [31] . Analysis of the 1-135 aa signal peptide region of the F protein of CDV-TM-CC demonstrated the same set of amino acid variations in comparison with the Onderstepoort strain as for the other Asia-1 strains (CYN07-hV, HLJ1-06, and Hebei): 8 S/ 8 K, 11 T/ 11 P, 19 (Fig. 6) . Among the Asia-2 strains (M25CR, 007Lm, 011C, 50Con, and 55L), variations in comparison with the Onderstepoort strain were found in 30 T/ 30 S, 53 S/ 53 A, 55 R/ 55 W, 59 S/ 59 Y, 62 N/ 62 K, 99 R/ 99 K, 110 I/ 110 V, and 111 N/ 111 K. Additionally, both the Asia-1 and Asia-2 strains had clade-specific amino acid variation in 21 P/ 21 Q. Moreover, the CDV-TM-CC strain had characteristic additional variations in 107 P/ 107 Y and 116 C/ 116 Y. Therefore, the signal peptide region of CDV-TM-CC has both Asia group-specific and individual variations. Among the CDV strains, amino acid variation was also found in 208 K/ 208 N in the F2 region (aa 136-224) for the Asia-1 group. Generally, there was high conservation within the hydrophobic fusion peptide (FP) domain at the N-terminus of the membrane anchored F1 subunit, with the exception of 233 A/ 233 V in the 98-2654 and 98-2646 strains (Fig. 6) . Amino acid variations between the Asia-1 and Asia-2 groups were also found in a region between the helical bundles (HB) and heptad repeats B (HRB) at 394 V/ 394 S, 429 R/ 429 K, and 466 L/ 466 I; within the trans-membrane (TM) domain at 627 C/ 627 Y, 634 Q/ 634 R, and 637 H/ 637 F; and within the cytoplasmic tail (CT) domain, at 656 R/ 656 K. Among the Asia group strains, the HRA (aa 250-307) and HB (aa 328-374) domains were highly conserved, with the exception of a 280 Q/ 280 A variation in the HRA domain. Likewise, the amino acids were highly conserved in the HRB (aa 557-601) domain in all CDV strains except for Hebei ( 583 D/ 583 N) and 5804P ( 587 V/ 587 I). Common amino acid changes in other regions of CDV strains in comparison to the Onderstepoort strain were found at 317 K/ 317 R and 556 S/ 556 G. The potential N-glycosylation sites (N-X-S/T) of the F protein were highly conserved at 141 NLS, 173 NVS, 179 NCT, and 517 NQS in the F1 region among all CDV strains as reported previously [32] [33] [34] (Fig. 6) . Moreover, the Asia-1 group (strains CYN07-hV, HLJ1-06 Hebei, and CDV-TM-CC) had specific potential N-glycosylation sites at 62 NRT and 108 NAT in the signal peptide region, with the exception of the CDV-TM-CC strain, which had the sequence 62 NKT. Five of these six potential glycosylation sites of the CDV-TM-CC strain were at the same positions within the known virulent CDV strains (A75/17, 5804P and 164071) at aa 108-110, 141-143, 173-175, 179-181, and 517-519, whereas 62 NKT was unique for CDV-TM-CC, and 62 NRT and 38 NIT were unique for 5804P. Cysteine is an a-amino acid that plays an important role in intramolecular disulfide bond formation and the steric structure of proteins. In the F protein of CDV-TM-CC, a total of 16 cysteine residues were detected. Among them, 14 residues (aa 123, 132, 180, 307, 446, 455, 470, 478, 502 , 507, 509, 531, 628, and 629) were located at identical positions in all CDV strains; however, several amino acid(s) were characteristic to individual strain(s), such as 67 R/ 67 C in the America group (strains Onderstepoort, 98-2654, 98-2646, Snyder Hill, CDV3, Shuskiy and Phoca-Caspian-2007) and 116 Y/ 116 C in CDV-TM-CC. The presence of amino acid variations, as well as specific N-glycosylation sites and cysteine residues within the F protein, could affect the immune response to CDV-TM-CC. Improved vaccination has reduced the frequency and magnitude of CD [35] . Distemper vaccination failures are uncommon, but outbreaks of CD continue to occur among vaccinated individuals and populations [4, 5, 36, 37] . The most common factor in CD occurrence is a lack of The H protein, a major structural protein of CDV, mediates host selection and pathogenicity, and the rate of genetic variation for its gene is greater than for other genes. With geographically distinct lineages, many studies have demonstrated that phylogenetic analysis can be carried out in accordance with the deduced amino acid sequences of the H protein [14, 18, 21, 38] . In this study, phylogenetic analysis based on the H protein identified seven clades of CDV strains (America-1, America-2, Asia-1, Asia-2, Europe, Arctic-like, and Europe wild-life), and CDV-TM-CC was classified into the Asia-1 group, with the highest identity to the Chinese strains, HLJ1-6 and Hebei, and the Japanese strain, CYN07-hV. Potential N-glycosylation sites may differ for the H protein of the wild-type and vaccine strains of CDV. Usually, only 4-7 potential sites are found within vaccine strains (such as Onderstepoort), in comparison with 8-9 sites in wild-type CDV strains (for example, 5804P). In particular, the 309-311 N-glycosylation site, which is specific for the wild-type strain [14, 18] , is suggestive of the pathogenicity of the CDV-TM-CC strain. Furthermore, the 584-586 N-glycosylation site has been acquired in the Asian-1 strains [18, 20] . Further study may determine whether these differences in glycosylation The N protein is a highly conserved immunogenic protein that can elicit cellular and humoral immunity [39] . Based on sequence differences between the gene of the wild strains and vaccine strain, the N protein may affect the seroprotection rate of the host and lead to immune failure. Like the H protein, the N protein of CDV-TM-CC showed the highest homology with the Asia-1 group. High homology was also observed with the Asia-2 group (strains M25CR, 007Lm, 011C, 50Con, and 55L) and wild-type strains (164071, A75/17, and 5804P). Moreover, the lowest homology was found between CDV-TM-CC and the Onderstepoort strain. Variation in the immunodominant epitope of the virus may change the structure, and therefore, we can speculate that the T cell-mediated immune response may be altered by variations in this protein. The P gene is extremely well conserved and, therefore, is particularly important in the phylogenetic classification. Based on the phylogenetic relationship of the deduced amino acid sequence of the P protein, CDV-TM-CC was also classified into the Asia-1 group. These results highlight the importance of considering the geographical setting to control the occurrence of the disease in a more efficient manner. The F protein is a surface glycoprotein that mediates viral entry into the host cell by fusion of the virion envelope and the host cell plasma membrane at a neutral pH. Within the F protein, the signal peptide region (aa 1-135) has the lowest amino acid homology, especially at positions 13-37 and 72-112. However, our analysis shows that the signal peptide region is relatively well conserved among the Asia-1 group, except for specific individual amino acids, indicating that the signal peptide of the F protein is geographically distinct. In addition, three amino acids specific to the CDV-TM-CC strain ( 62 K, 107 Y, and 116 Y) are located in the signal peptide region. The previous study reported that the amino acids 208 K and 216 L are specific for the CDV vaccine strains; however, we also found 208 K in the wild-type strains in the America group (A75/17, 164071, and 5804P) and Asia-2 group (011C, M25CR, 55L,50Con, and 007Lm). The F protein of the CDV-TM-CC strain has six potential glycosylation sites. Among them, differences were found to reside mainly in the signal peptide region, but no clear rule could obviously explain the differences in the wide-type and vaccine strains or the geographical variation, including the occurrence of a strain-specific site (62-64 NKT) for CDV-TM-CC. Four additional potential glycosylation sites were recognized at positions 141-143, 173-175, 179-181 in the F2 region and 517-519 in the F1 region, as reported previously [32] [33] [34] . The 108-110 N-glycosylation site is specific for the wildtype strains (5804P, A75/17, and 164071) and the Asia-1 group (Hebei, HLJ1-06, and CYN07-hV), and may be another important factor in vaccination failure. The fusion peptide (FP) domain also was found to be highly conserved among all CDV strains, except for 233 A/ 233 V in 98-2654 and 98-2646. In short, the genetic/antigenic drift observed in the currently circulating CDV strains should be considered as a possible factor leading to the resurgence of CD cases. Analysis of CDV strains detected globally and from a variety of host species will provide a more in-depth understanding of the global ecology of CDV and will provide the basis for the improvement of current CDV vaccines. The wild-type CDV-TM-CC strain, originally isolated from spleen homogenate from a fully vaccinated Tibetan Mastiff in China, was classified into the Asia-1 group cluster of CDV strains based on the sequence of its H protein and verified by the sequence of its P protein. Variations in specific amino acid residues, N-glycosylation sites, and cysteine residues throughout the CDV-TM-CC genome may explain the failure of the dog to mount vaccine-mediated protection against CD. These results provide the foundations for the global improvement in current CDV vaccines. Virus Infections of Carnivores Role of Glycosylation of Notch in Development Acknowledgments This work was supported by Ecology of Zoonoses and Research of Infection and Immunity mechanisms (2012CB722501).