key: cord-0790391-3p0ozn6k
authors: Ji, Jun; Li, Wanyu; Hu, Wen; Xu, Xin; Kan, Yunchao; Yao, Lunguang; Bi, Yingzuo; Xie, Qingmei
title: Novel Genotype Definition and the First Epidemiological Investigation of Canine Adenovirus Type 2 in Dogs in Central China
date: 2020-08-19
journal: Front Vet Sci
DOI: 10.3389/fvets.2020.00534
sha: c22a49e1d059d78328442b45d8a8d342e6859cb8
doc_id: 790391
cord_uid: 3p0ozn6k

Infections caused by canine adenovirus (CAdV) type 1 have been reported worldwide in the past two decades. However, only few studies have specifically reported the prevalence of CAdV type 2 (CAdV-2). The present study investigated the persistent circulation of CAdV-2 in dogs with diarrhea in the Henan, Hubei, and Jiangsu provinces in central China from 2017 to 2019. We conducted polymerase chain reaction for detecting CAdV-2 and other related pathogens in 224 rectal swabs of pet dogs and the co-infection of canine diseases was also analyzed. In addition, the structural protein genes—Fiber, Hexon, and Penton—of the isolated CAdV-2 strains were sequenced and analyzed. The similarity between Hexon and Penton among the 19 strains was 97.4%, as revealed by sequence alignment. Multiple sequence alignment results showed that the Fiber gene sequences of these CAdV-2 strains shared 97.4–99.8% nucleotide and 94.1–99.3% amino acid identity with reference sequences and shared only 79.0–80.5% nucleotide and 77.3–80.5% amino acid identity with the vaccine strain CLL, indicating that Fiber harbored most of the variant sites. Furthermore, pairwise sequence comparisons of Hexon of CH-JS-1901 and CH-HN-1801 with that of India2006 revealed a novel genotype. Furthermore, protein model prediction showed that the amino acid mutation of fiber protein in 19 strains was located in the head region, that may cause structural changes on the surface of the fiber protein. These findings are of significance for monitoring the epidemiology of CAdV-2 infection and developing a novel vaccine which contribute to understanding genetic evolution of CAdV-2 in China.

Canine adenovirus (CAdV), a member of the Mastadenovirus genus of the Adenoviridae family, is one of the pathogens causing severe viral diseases in dogs worldwide (1) . There are three major capsid proteins encoded by CAdV-Fiber, Hexon, and Penton. Adenoviral Fiber protein plays a crucial role in infection by attaching the virus to a specific cell surface receptor, whereas the Penton base is responsible for virus internalization (2) . The capsid is icosahedral: and the faces are composed of 240 Hexons, each comprising three identical [Hexon] proteins, and the trimer constitutes the outer sheath of the capsid. Because it contains a specific antigenic determinant, it is also the main protein targeted by the host neutralizing antibody (3) .

To date, epidemic reports of CAdVs have been mainly conducted in Italy, France, Britain, Turkey, and North America (4) (5) (6) (7) , where there are wild animals such as wolves, foxes, skunks, wild dogs, and brown bears (8) (9) (10) (11) (12) (13) . On contact with the carrier, the virus is transmitted via air, body fluids, and excreta (14) .

CAdVs are classified into two types: CAdV type 1 (CAdV-1), which causes canine infectious hepatitis and fox encephalitis, and CAdV type 2 (CAdV-2), which causes infectious laryngotracheitis and enteritis (15) . CAdV-2 differs from CAdV-1 in virulence, antigen structure, infection ability, and erythrocyte aggregation range. Dogs infected with CAdV-2 show persistent high fever, cough, and other clinical manifestations, such as mucinous rhinorrhea, tonsillitis, laryngotracheitis, pneumonia, vomiting, and diarrhea (16) . CAdV-2 is one of the main pathogens that cause canine viral enteritis. Recently, more studies have reported diarrhea in dogs, and the number of cases of CAdV-2 infection has gradually increased. However, only few studies have focused on the association between CAdV-2 and diarrhea in dogs in clinical diagnosis (17) (18) (19) (20) .

In 1962, Dithfield isolated the A26 strain (Toronto A26/61) of adenoviruses that caused respiratory tract lesions (laryngotracheitis) but not hepatitis (21) . Subsequently, CAdV-2 has become popular worldwide. In 1985, Hamelin confirmed that CAdV-2 can cause not only upper respiratory tract diseases but can also proliferate in the digestive tract epithelium and lead to diarrhea in dogs, which may also be associated with the variation of CAdV-2 (22) . In China, CAdV-2 infection was previously reported only in Yunnan, Jilin, and other regions (23) . At present, most vaccines used to prevent CAdV infection in China are imported, and the epidemic strains may not be completely homologous with the vaccine strains isolated in foreign countries. Moreover, only few studies have reported the prevalence of adenovirus in most parts of China. Therefore, it is essential to investigate the prevalence, variation, and regularity of CAdV-2 and determine the possible association between CAdV-2 and structural changes in major capsid proteins.

Here we studied the prevalence of CAdV-2 in some provinces in central China, amplified and sequenced the structural protein genes of CAdV-2 strains, and compared the obtained sequences with reference sequences of CAdV-2 strains reported in the literature worldwide. We believe that our study will contribute reliable data for research on developing a canine adenoviral vaccine, which is beneficial for the targeted prevention and treatment of CAdV infection.

Rectal swabs from 224 dogs with evidence of diarrhea and 164 health dogs were collected from December 2017 to March 2019 in 

Supernatants from the diluted CAdV-2-positive samples were filtered through a 0.22-µm membrane and inoculated onto a monolayer of Madin-Darby canine kidney cells at 37 • C in 5% CO 2 , and this process repeated until cytopathic effects were apparent (27) . To detect the presence of other viruses mentioned above, total DNA/RNA from cell culture isolates was extracted using the Viral DNA/RNA kit (TaKaRa, Beijing, China).

To amplify the sequences of Fiber, Hexon, and Penton, three pairs of oligonucleotide primers were used ( Table 1) . Amplifications were performed in 50 µL of reaction mixture comprising 2.5 mM dNTPs (4 µL), 10 µM of each primer (2 µL), 500 ng DNA template, and PrimerSTAR HotStart DNA polymerase (TaKaRa Biotechnology Co., Ltd., Dalian, China). Sequence amplification was performed under the following conditions: initial denaturation at 94 • C for 3 min, followed by 30 cycles of 94 • C for 30 s, 56 • C for 1 min, 72 • C for 1-3 min for different length of gene targets, and a final extension at 72 • C for 10 min in an automated thermal cycler. All the purified PCR products were cloned into a pMD18-T easy vector (TaKaRa Biotechnology Co., Ltd.) for Sanger sequencing (Syn-Biotechnology, Suzhou, China).

The contig assemblies of Fiber, Hexon, and Penton of CAdV-2 were arranged using the DNASTAR 7.0 software (DNASTAR, Inc., USA) and were compared with those of reference CAdV strains obtained from GenBank. Information of the reference strains is presented in Table 2 . The amino acid mutation points of three gene segments were displayed using the MegAlign software. A divergence analysis of the Fiber protein of CAdV-2 was performed with WebLogo (http://weblogo.threeplusone. com/) for sequence logo generator. A phylogenetic tree of the CAdVs identified in this study and other reference canine adenovirus strains was constructed based on nucleotide alignments using the maximum likelihood (ML) method and the Molecular Evolutionary Genetics Analysis software (MEGA X). A bootstrap analysis was performed using 1,000 replications, and the phylogenetic tree was visualized using Tree-view (http://www.evolgenius.info/evolview/).

A divergence analysis of the Fiber protein of CAdV-2 was performed with WebLogo (http://weblogo.threeplusone.com/). For the protein encoded by Fiber that is partly associated with the receptor binding of virus-infected cells, it is necessary to analyze the mutation sites in the model. The Fiber protein sequence of the CAdV-2 reference strain CC0710 in China was downloaded from GenBank, and the model was predicted using the Fiber amino acid sequence as a template. Amino acid site mutation was identified in 19 isolates, and the tertiary structure model of Fiber was established using the online software SWISS-MODEL (https://swissmodel.expasy.org/interactive). After prediction, the protein analysis software PyMOL was used for image processing and data analysis of the prediction model.

The association among clinical signs and presence of the virus in the samples of healthy dogs and dogs with evidence of diarrhea were evaluated using the chi-squared test. Logistic regression was used to identify possible bivariate associations between the presence of CAdV-2 DNA and the presence of other pathogens in the samples. Data were processed using SPSS 19.0 software for Windows (SPSS Inc.; Chicago, IL, USA). A p < 0.05 was considered statistically significant. Figure 1 . Furthermore, none of the pathogen was detected in the healthy dogs.

A statistically significant difference (p < 0.05) could be detected compared the dogs with evidence of diarrhea with health dogs. Logistic regression analysis was performed to evaluate possible bivariate association between the presence of co-infection of CadV-2 and other pathogens with the diarrhea signs, and showed no association between the presence of the CadV-2 and clinical signs (p > 0.05).

The Penton, Hexon, and Fiber sequences of the CAdV-2 strains identified in this study are deposited in GenBank, and the accession numbers of each gene fragments are shown in Table 3 . Sequence analysis based on the CAdV-2 strains and other reference CAdV strains from GenBank demonstrated that the Penton of the 19 CAdV-2 strains identified herein was 1,434 nucleotides (bp), in length, whereas Fiber was 1,629 bp and Hexon is 2,718 bp in length.

Sequence comparisons revealed that the Penton of the 19 CAdV-2 strains identified in this study showed 99.5-99.9% nucleotide identity and 99.5-99.9% amino acid identity. Amino acid mutations were the same as those of amino acid sites such as Val18Ala, Ser20Pro, Leu22Pro, Asn38His, and Ile57Thr. Mutations at Gly32Glu were found in two CAdV-2 strains identified here (CH-HB-1801 and CH-JS-1701). At the 471 locus, the reference strain A26/61 was found to be different from all the isolated strains, which had Val at the 471 locus; however, the A26/61 strain had Ala.

Hexon showed 98. Compared with that observed for the vaccine strain CLL, the nucleotide and amino acid identity of 19 isolated strains with Fiber was only 79.0-80.5% and 77.3-80.5%, respectively. Among them, the CH-JS-1801 strain shared a high degree of amino acid identity with the Chinese reference strain CAdV-2 CC0710QB. Conversely, alignment of the amino acid sequence of Fiber with the sequences of the reference strains indicated that the variation of the 19 isolated strains is complex. A total of 20 mutations occurred at the CH-HN-1901 amino acid sites, among which 6 were the same as those at the CH-HN-1701 amino acid sites; however, only 12 amino acid sites showed mutations in the CH-HN-1701 strain. The mutation sites of only five tested strains were the same as those of the reference strain CC0710 isolated in Shanghai, China, indicating a large difference in mutation sites. Except for the fact that the amino acids of CH-HN-1903 and CH-HN-1904 are basically the same, all the additional strains have 1-to 20-point mutations. Furthermore, both CH-HN-1701 and CH-HN-1801 had a mutation that occurred at site 25, where the mutation rate for Thr-25 mutation to Ala-25 was 10.5% (two strains in the Henan Province), and it was the same as that observed for CC0710. The main amino acid mutations of fiber proteins were displayed in Figure 2 (Amino acid mutation sites of each strain were displayed detailedly in Table S1 ).

Phylogenetic tree and cluster information for the Penton, Hexon, and Fiber nucleotide sequences are shown in Figure 3 .

According to the phylogenetic analysis of Penton, 19 CAdV-2 mainly belonged to the type 2 clade and were further divided into three subsets. Among them, the strains aligned with the Canadian A26/61 reference strain belonged to a small branch mainly from the Henan and Jiangsu provinces. However, the CH-JS-1701 and CH-HB-1801 strains were located on a small branch alone, suggesting that these two strains have a similar origin. There was no significant regional difference among the strains.

As shown in Figure 3B , the phylogenetic tree constructed according to the Hexon sequence showed that the strain was divided into three branches, that 19 isolated strains belonged to CAdV-2, and that the other reference CAdV-1 strains were located on an independent branch. Among them, 18 strains were closely related to the reference strains isolated from Jilin Province, whereas the CH-HN-1801 strain was similar to the strains isolated from Canada and India. Hexon Following this definition, the CAdV-2 strains were subdivided into CAdV-2.1 and CAdV-2.2 genotypes. Overall, the two strains CH-JS-1901, and CH-HN-1801 were distant from the other strains in the phylogenetic tree. Considering the geographical distribution, the strains isolated from the Henan and Hebei provinces were clustered on the same branch, whereas those isolated from Jiangsu Province were scattered among the branches. The phylogenetic tree constructed according to Fiber sequences showed that the CH-HN-1801 strain was clustered with the SH01 and A26-61 strains on a small branch, and the remaining strains were located in another small branch. Isolated strains of Henan Province CH-HN-1802, CH-HN-1902, and Jiangsu Province isolate strains in the evolutionary tree distribution distance is relatively close, CH-JS-1701, and Henan Province isolate strains clustered on another small branch of the evolutionary tree. However, most of the strains were close to the reference strain CC0710 and far from the SH01 and A26-61 strains. The specific distributions are shown in Figure 3C .

CAdV-2 binds to cell surface receptors through the Fiber protein.

As shown in Figure 4 , most of the amino acid mutation sites of the 19 isolated strains are concentrated on the knob; the shaft, which comprises proline and hydrophobic amino acid residues, was not conservative, which is associated with the antigen specificity of distinct adenoviruses serotypes. Based on the composition of the fibrillar structural functional regions of human adenoviruses, we found that the structural functional regions of CAdV-2 are composed of amino acids 1-42, 43-365, and 366-542, which are located in the tail, shaft, and head regions, respectively (28) (29) (30) . The mutant amino acids in the 19 CAdV-2 strains were located at sites 18, 25, 26, 30, 36, 37, 40 , and 41, which were located in the tail of the fibrillar process; sites 44, 45, 50, etc., which were located on the basal axis of the fibrillar process, also known as the shaft; and sites 370, 374, 393, etc., which were located in the ganglion region of Fiber, also known as the knob (Figure 4A ). There are four differences between Figures 4B,C at sites 37, 77, 408 , 520, 521, 523-532 of Fiber. 

The present study reports the findings of investigations on the presence of CAdV-2 in 224 dogs from the Henan, Hubei, and Jiangsu provinces, China, from 2017 to 2019. All the dogs presented with diarrheal symptoms, with no obvious respiratory problems, and the prevalence of CAdV-2 infection (8.5%, 19/224 dogs) was verified; CAdV-2 infection was often concomitant with CPV-2 and other pathogenic infections. Even though the CAdV-2 was only detected in dogs with diarrheal symptoms, logistic regression analysis showed no association between the presence of the CAdV-2 and clinical signs. Previous studies have reported the infection of CAdV-2 in dogs with diarrhea, and the number of cases of CAdV-2 infection has gradually increase (17, 19) . We speculated the poor association between the presence of the CAdV-2 and diarrhea may derive from small quantity of cases singly positive for CAdV-2. Therefore, whether variant CAdV-2 caused canine enteritis required larger investigation studies and animal regression with CAdV-2 challenge experiments.

Penton, Hexon, and Fiber are the three main structural proteins of CAdV-2. Penton serves as a pedestal, which helps the virus adhere to the cell surface (31, 32) , and its sequence is relatively conservative compared with the sequences of reference strains. In the present study, we found that only few amino acid sites on Penton differed from each other, suggesting an important role of Penton in the stability of virus survival function.

Previous study have suggested that the specific epitopes present on the Hexon protein are the target of neutralizing antibodies and the most sensitive part of adenoviruses to immune selection pressure and that its specificity is likely determined by the amino acid sequence of the primary structure of the Hexon protein (33) . A comparison between the amino acid sequences of Hexon showed that the continuous site mutations of CH-HB-1701 and the discontinuous site mutations of CH-HN-1901 were first identified in the CAdV-2 isolates, which were different compared with the CAdV-1 reference strains isolated in Europe (8, 9, 11) . Whether these continuous site mutations cause substantial changes in antigenic epitopes needs to be further analyzed. The phylogenetic tree based on Hexon gene sequences indicated new genotypes for CAdV-2.1 and CAdV-2.2. Even though the 19 CAdV-2 strains could also been sub-classed as new genotypes based on phylogenetic trees about Penton and Fiber genes, phylogenetic analysis of Hexon gene was choosed to be the genotyping criteria for Hexon was the main protein targeted by the host neutralizing antibody. The Penton and Fiber genes for the proposed Clade 2.2 strains do not also separate into a distinct sub-clade might cause by the independent evolution of each gene. The relationship between CAdV-2 genotyping and evolution of Penton and Fiber genes needs deeply discussion in future studies.

Analysis of the Fiber gene sequences revealed the same pattern of amino acid substitutions observed in three strains isolated in this study (CH-HN-1701, CH-HN-1801, and CH-HN-1802) and two previously reported strains isolated from bats (Mm32 and TJM) (33, 34) . In particular, these amino acid substitutions occurred at positions 51 and 210, which are different from those observed in other proteins. The properties of 50 amino acids sites in the fiber section of CAdV-2 were altered. Seven mutant amino acids (Pro27Arg, Pro40Arg, Pro50His, Pro66Arg, Pro113His, Pro171Arg, Pro278His) changed their properties from uncharged to acidic; two mutant amino acids (Arg155Cys and Asp290Gly) changed from acidic to polar; one mutant amino acid (Asp519Ala) changed from acidic to medium; and three mutant amino acids (Glu147Lys, Glu198Lys, and Asp531His) changed from basic to acidic, and two mutant amino acids (Lys505Glu and Lys529Glu) changed from acidic to basic. The length of the basal axis of adenoviruses has been shown to play a very important role in viral infection and affinity (35) . The conserved amino acid region of the fibrillar process is the tail, which is the subgroup specific antigen of adenoviruses and is connected to the five adjacent Penton base proteins (28) . In the present study, a total of 49 sporadic mutant amino acids were found to be located on the basal axis. The amino acid properties of 22 sites did not change and only two amino acid sites showed acid to alkaline changes, which had a limited effect on the variation of the properties of the isolated strains.

The RDP 4.36 software was used to predict the possibility of recombination events in three major capsid proteins of the CAdV-2 but no recombination event was predicted (data was not shown). Analyzing the mutations in the CAdV-2 strain is helpful for screening the vaccine strain used in the present study. At present, most vaccines used for preventing CAdV infection in China are replicas of foreign vaccines (36) (37) (38) . However, over time, some mutations occur in CAdV-2 and the amino acid sites of these mutations are the same as those of the mutations in the CAdV-1 strain. For instance, a mutation from Asn-71 to Arg-71 occurred in the CH-HN-1901 strain and a mutation from Pro-279 to Thr-279 occurred in the CH-HN-1902 strain. The effect of these site mutations on the genetic evolution of the strain needs to be verified in further experiments.

In conclusion, this study describes the prevalence and coinfection of CAdV-2 in some regions of China, revealed that CAdV-2 may be one of the main pathogens causing canine diarrhea. Novel genotypes of CAdV-2 were analyzed. Among 19 CAdV-2 strains, there were much more mutation sites in Fiber than in Penton and Hexon. Thus, new vaccines should be screened to effectively control the dominant strains and expand the scope of immunization.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material. 

JJ and XX designed the study and wrote the manuscript. WL and WH performed the sampling and data collection. WL performed the clinical investigations, necropsies, and molecular tests. YK, YB, QX, and LY performed the molecular genetic studies and helped to draft the manuscript. All authors have read and approved the final manuscript. 

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Ji, Li, Hu, Xu, Kan, Yao, Bi and Xie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Genetic content and evolution of adenoviruses

Integrins alpha v beta 3 and alpha v beta 5 promote adenovirus internalization but not virus attachment

Exploration of new sites in adenovirus hexon for foreign peptides insertion

A retrospective investigation of canine adenovirus (CAV) infection in adult dogs in Turkey

Infectious canine hepatitis: an "old" disease reemerging in Italy

Infectious canine hepatitis in a brown bear (Ursus arctos horribilis) from Alaska, USA

Investigation of the presence of canine adenovirus (CAdV) in owned dogs in Northern Italy

Molecular epidemiology of canine adenovirus type 1 and type 2 in freeranging red foxes (Vulpes vulpes) in Italy

An update on canine adenovirus type 2 and its vectors

Sequential circulation of canine adenoviruses 1 and 2 in captive wild carnivores, France

Characterization of a novel adenovirus isolated from a skunk

Serological and molecular epidemiology of canine adenovirus type 1 in red foxes (Vulpes vulpes) in the United Kingdom

Longitudinal study of viruses associated with canine infectious respiratory disease

Isolation and identification of canine adenovirus type-2 from the upper respiratory tract of a dog

Canine adenovirus type 2 infection in four puppies with neurological signs

Identification of enteric viruses circulating in a dog population with low vaccine coverage

Canine morbillivirus (canine distemper virus) with concomitant canine adenovirus, canine parvovirus-2, and Neospora caninum in puppies: a retrospective immunohistochemical study

Characterisation of the canine faecal virome in healthy dogs and dogs with acute diarrhoea using shotgun metagenomics

Severe enteric disease in an animal shelter associated with dual infections by canine adenovirus type 1 and canine coronavirus

Association of canine Adenovirus (Toronto A 26/61) with an outbreak of Laryngotracheitis

Genotypic characterization of type-2 variants of canine adenovirus

Canine circovirus 1 (CaCV-1) and canine parvovirus 2 (CPV-2): recurrent dual infections in a papillon breeding colony

Detection and differentiation of CAV-1 and CAV-2 by polymerase chain reaction

A fast and simple one-step duplex PCR assay for canine distemper virus (CDV) and canine coronavirus (CCoV) detection

Epidemiological investigation reveals genetic diversity and high co-infection rate of canine bocavirus strains circulating in Heilongjiang province, Northeast China

Bioprocess development for canine adenovirus type 2 vectors

Sequence analysis of the fiber genomic region of a porcine adenovirus predicts a novel fiber protein

Threedimensional structure of canine adenovirus serotype 2 capsid

Structural and mutational analysis of human Ad37 and canine adenovirus 2 fiber heads in complex with the D1 domain of coxsackie and adenovirus receptor

Sequence analysis of the canine adenovirus 2 fiber-encoding gene

Nucleotide sequence analysis of the leader segments in a cloned copy of adenovirus 2 fiber mRNA

Comparative sequence analysis of the hexon gene in the entire spectrum of human adenovirus serotypes: phylogenetic, taxonomic, and clinical implications

Genome analysis of bat adenovirus 2: indications of interspecies transmission

Phylogeny and primary structure analysis of fiber shafts of all human adenovirus types for rational design of adenoviral gene-therapy vectors

Evaluation of the efficacy and duration of immunity of a canine

Special thanks to members of the Agriculture and Engineering College, Nanyang Normal University, for supporting our work.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets. 2020.00534/full#supplementary-material