key: cord-0795087-yofqa4zj authors: Hu, Wen; Liu, Qiang; Chen, Qinxi; Ji, Jun title: Molecular characterization of Cachavirus firstly detected in dogs in China date: 2020-09-03 journal: Infect Genet Evol DOI: 10.1016/j.meegid.2020.104529 sha: ec7e999beedcdc713da3e116d868cd4d4d4030a0 doc_id: 795087 cord_uid: yofqa4zj Canine Cachavirus was novel parvovirus species has been firstly identified in dogs in USA and was classified within the proposed Chaphamaparvovirus genus. To investigate Cachavirus infection in dogs in China, 408 rectal swabs from healthy and diarrheic dogs obtained during 2018–2019 were screened. The rate of Cachavirus positivity was 0% and 1.55% in healthy or diarrheic dogs, respectively. However, statistical analysis suggested no association between the presence of the virus and clinical signs (p > 0.05). Nucleotide identity was 98.2%–98.9% for NS1 and 98.6%–99.1% for VP1, and amino acid identity was 97.9%–98.7% for NS1 and 98.8%–99.6% for VP1 between the five Chinese strains and Cachavirus-1A and Cachavirus-1B detected in the United States. Phylogenetic analysis also indicated that these Cachavirus strains are genetically related to Cachavirus-1A and Cachavirus-1B. This study confirms the presence of Cachavirus in pet dogs in China and provides novel findings on its molecular characteristics. Viruses belonging to the Parvoviridae family include 4-6-kb-long linear, single-stranded DNA genomes (Pénzes et al., 2019) . Recently, Parvoviridae family has been re-organized by International Committee on Taxonomy of Viruses (ICTV). The novel subfamily Hamaparvovirinae comprises the current genera Hepandensoviruses, Penstyldensovirus, and Brevidensovirus, together with the yet unclassified chapparvoviruses (ChPV) (Cotmore et al., 2019) . Phylogeny evidence suggested the ChPV detected in many kinds of vertebrates, harbored close relation to members of other genera of subfamily Hamaparvovirinae, which evidently derive from ancient members of this lineage, have been identified in several arthropod genomes. ( Pénzes et al., 2019; Souza et al., 2017) . These reports suggested that ChPV lineage infects both vertebrates and invertebrates, and transmission of ChPVs between distantly related host species may have occurred in the past. Recently, cases of ChPV infection in various organisms have been reported. In addition to rat parvovirus 2 (RPV2), ChPV include fruit bat (Eidolon helvum) parvovirus 1 (Yinda et al., 2018) , E. helvum parvovirus 2 (Baker et al., 2013) , common vampire bat (Desmodus rotundus) parvovirus (Souza et al., 2017) , murine kidney parvovirus (Roediger et al., 2018) , murine ChPV (Williams et al., 2018) , simian parvo-like virus 3 (Kapusinszky et al., 2017) , turkey parvovirus 1 (Reuter et al., 2014) , porcine parvovirus 7 (Palinski et al., 2016) , Tasmanian devil-associated ChPV strains (Chong et al., 2019) , red-crowned crane-associated parvovirus (Wang et al., 2019) , chicken ChPV 1 and 2 (Lima et al., 2019) , and peafowl parvovirus (PePV) 1 and PePV2 (Liu et al., 2020) . J o u r n a l P r e -p r o o f Journal Pre-proof Meanwhile, Cachavirus, which is very similar to ChPV, was first detected in dog feces in the United States in 2019 (Fahsbender et al., 2019) , raising concerns regarding the emergence and coinfection of Cachavirus in dogs. The Cachavirus was shown to harbor two major open reading frames (ORFs) encoding NS1 and VP1. Therefore, Cachavirus may show an extensive geographic distribution and may infect diverse hosts (Souza et al., 2017) . To obtain detailed understanding of the prevalence and evolution of Cachavirus in dogs, we performed a large-scale investigation using rectal swabs from healthy and diarrheic dogs. Rectal swabs from 85 healthy and 323 diarrheic dogs were collected from veterinary hospitals in Henan, Hubei, and Zhejiang Provinces, China, during 2018-2019. Viral DNA and RNA were extracted using the EasyPure Viral DNA/RNA Kit (TransGen Biotechnology, Inc., Beijing, China) in accordance with the manufacturer's instructions. The extracted DNA and RNA sample were stored at −80°C until use. The extracted DNA and RNA samples were tested for Cachavirus (Fahsbender et al., 2019) , canine coronavirus (CCV) (Nguyen et al., 2017) , canine distemper virus (CDV) (Wang et al., 2017) , and canine parvovirus (CPV) (Han et al., 2015) using RT-PCR as previously described. Table S1 . Gene amplification was performed via PCR using a 20-μL reaction mixture containing DNA template (>100 ng/μL), 6-pmol upstream/downstream primers, PrimeSTAR HS DNA polymerase, and supporting reaction buffer (TaKaRa Biotechnology Co., Ltd., Dalian, China). The following cycling conditions were used: initial denaturation at 95°C for 3 min; 30 cycles of denaturation at 95°C for 30 s, annealing at 55°C for 30 s, and extension at 72°C for 1 min; and final extension at 72°C for 10 min. The obtained J o u r n a l P r e -p r o o f Journal Pre-proof amplicons were cloned into pMD18-T Easy Vector (TaKaRa Biotechnology Co., Ltd.) for subsequent sequencing (Syn-Biotechnology, Suzhou, China). PCR and partial-genome sequencing were performed at least three times. The ClustalW method was employed using MEGA 7.0 to analyze the identity between the obtained and reference strains. Moreover, phylogenetic trees were constructed based on NS1 and VP1 sequences to analyze the relationships between the tested Chinese strains and reference strains. Specific information of each viral strain is presented in Figure 1A and 1B. Specifically, the maximum likelihood method was used. MEGA 7.0 was used with the settings of pairwise deletion option and 1000 bootstrap replicates. Illustrator for Biological Sequences was used to obtain a visual overview of NS1 and VP1 from the five Chinese strains and Cachavirus-1A. Based on the sequence alignment results, mutated amino acid sequences were marked on SWISS-MODEL (https://swissmodel.expasy.org/interactive), and post-modeling Pdb files were constructed using PyMOL for collation and preservation. Chi square or fisher's exact test was used to compare the frequency of Cachavirus between healthy and diarrheic dogs. Statistical analyses were performed using GraphPad Prism 8.0 (San Diego, CA, USA). A p < 0.05 was considered statistically significant. J o u r n a l P r e -p r o o f Journal Pre-proof 3. Results Upon viral screening, five Cachavirus-positive samples were collected from diarrheic dogs in Henan and Hubei Provinces and no infected samples were identified from healthy dogs. Moreover, CDV and CPV coinfection with the five Chinese strains was identified ( Table 1 ). The positivity rate in healthy or diarrheic dogs was 0% and 1.55%, respectively. Statistical analysis showed no association between the presence of the virus and clinical signs (p > 0.05). In the phylogenetic tree based on NS1 and VP1 sequences (Figure 1 ), the reference strains and the five Chinese strains were generally divided into two major branches. Intuitively, the five Chinese strains were more closely related to Cachavirus-1A and Cachavirus-1B. Cachavirus-1A was used as a reference to display the putative structure of proteins in the five Chinese strains (Figure 2 ). The genomic structure of the five strains was consistent with that of Cachavirus-1A, with two major ORFs encoding Table 2 . Due to these mutations, predict changes based on tertiary structure modeling with NS1 structure in Cachavirus-1A as a reference were shown in Figure 3 . Viral enteritis is a major threat to pet dogs that can lead to life-threating illness. Canine parvovirus (CPV-2), canine distemper virus (CDV), canine coronavirus (CCoV) and canine adenovirus (CAV) are the main pathogens that cause canine viral enteritis (Cardillo et al., 2020; Deng et al., 2018) . In this study, Cachavirus, belonging to (Nuñez et al., 2020) . However, report about pathogenicity of the ChPVs in vertebrate hosts is still very limited. However, a ChPV called mouse kidney parvovirus (MKPV) circulates among laboratory mice populations has been demonstrated to cause a kidney disease known as inclusion body nephropathy (Roediger et al., 2018) . Meanwhile, murine ChPV, has been detected at a very high prevalence in murine liver tissue, suggesting it is a gastrointestinal agent (Williams et al., 2018) . Herein, whether ChPV infection leads to gastroenteritis in dogs needs warrants further investigation. Primers targeting the whole genome have been designed based on the genome sequences of Cachavirus-1A and Cachavirus-1B; however, these primers failed to amplify the 5′ untranslated region (1-493); this may be because of the special DNA structure or the Chinese strains may markedly differs from that of strains reported from the United States. Therefore, future research should explore the mechanism of virus propagation using second-generation sequencing to obtain accurate whole genome sequences. Unfortunately, Cachavirus replication was not detected in CRFK and MDCK cell line. The reason for the non-cultivatable nature of the viruses remains unclear. According to identity and phylogenetic analyses based on NS1 and VP1 of the five Chinese strains, all strains are closely related to the two strains detected in the United States. Besides these two Cachavirus strains, the detected Chinese strains shared the highest NS1 and VP1 identity with bat parvovirus from Cameroon and simian parvo-like virus 3 from diarrheic rhesus macaques in the United States (Kapusinszky et J o u r n a l P r e -p r o o f Yinda et al., 2018) . Frequent recombination between mouse parvoviruses (MPVs), Minute viruses of mice (MVMs) and hamster parvovirus (HaPV) has been predicted (Shackelton et al., 2007) . Phylogenetic analysis also revealed the closer relationship between the five Chinese Cachavirus strains with bat parvovirus from Cameroon than with other ChPVs; however, no recombination events occurred based on recombination analysis (data not shown). Hence, the role of bats and primates in the transmission and evolution of ChPVs remains unknown. The five Chinese strains harbored increased number of mutations compared with the two strains from the United States. Previous report have highlighted the importance of obtaining the NS1 coding sequence in molecular epidemiology investigations of CPV-2 (Cotmore et al., 2014; Mira et al., 2019) . NS1 of Cachavirus-1A, Cachavirus-CN20181128, and Cachavirus-CN20190714 shows the SGGY252-255CLTF mutation, whereas NS1 of the three Chinese strains showed the GGY253-255VTF mutation. However, these continuous mutations did not lead to remarkable changes in protein structure. In contrast, mutation at the 442 site led to significant structural changes in Cachavirus-CN20190917 and Cachavirus-CN20191013 (Glu442Lys). Meanwhile, Cachavirus-CN20181128, Cachavirus-CN20190714, and Cachavirus-CN20191013 showed different structural changes of the 406-416 loci without mutations, whether these changes are caused by the neighboring effect of adjacent mutations should be studied in future. Furthermore, these mutation sites may show different origins due to regional variations. To sum up, this is the first report about the presence of Cachavirus in dogs in China. Phylogenetic analysis of NS1 and VP1 of the five Chinese strains revealed a close relationship with the two Cachavirus strains from the United States. These findings provide a reference for studying the transmission and evolution of ChPV. Nonetheless, large-scale studies are warranted to confirm the pathogenicity of this virus. Sampling was not harmful to the dogs. The research protocol was approved by the Animal Welfare and Ethics Committee of Nanyang Normal University (approval number: 18032, year: 2018). An ethical statement was not needed sample collection. J o u r n a l P r e -p r o o f Metagenomic study of the viruses of African straw-coloured fruit bats: detection of a chiropteran poxvirus and isolation of a novel adenovirus Lifestyle as Risk Factor for Infectious Causes of Death in Young Dogs: A Retrospective Study in Southern Italy Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics ICTV Virus Taxonomy Profile: Parvoviridae The family Parvoviridae A multiplex PCR method for the simultaneous detection of three viruses associated with canine viral enteric infections Chapparvovirus DNA Found in 4% of Dogs with Diarrhea Full-length genomic characterizations of two canine parvoviruses prevalent in Northwest China Case-Control Comparison of Enteric Viromes in Captive Rhesus Macaques with Acute or Idiopathic Chronic Diarrhea The intestinal virome of malabsorption syndrome-affected and unaffected broilers through shotgun metagenomics Genomic and transcriptional analyses of novel parvoviruses identified from dead peafowl Molecular Characterization and Evolutionary Analyses of Carnivore Protoparvovirus 1 NS1 Gene Molecular Characterization and Pathogenicity of Chicken Parvovirus (ChPV) in Specific Pathogen Free Chicks Infected Experimentally Discovery of a novel Parvovirinae virus, porcine parvovirus 7, by metagenomic sequencing of porcine rectal swabs An Ancient Lineage of Highly Divergent Parvoviruses Infects both Vertebrate and Invertebrate Hosts Novel circular single-stranded DNA virus from turkey faeces An Atypical Parvovirus Drives Chronic Tubulointerstitial Nephropathy and Kidney Fibrosis Comparative analysis reveals frequent recombination in the parvoviruses Chapparvoviruses occur in at least three vertebrate classes and have a broad biogeographic distribution Characterization of canine coronavirus spread among domestic dogs in Vietnam Rapid and sensitive detection of canine distemper virus by real-time reverse transcription recombinase polymerase amplification The fecal virome of red-crowned cranes Viral Diversity of House Mice Cameroonian fruit bats harbor divergent viruses, including rotavirus H, bastroviruses, and picobirnaviruses using an alternative genetic code Wen Hu Data curation: Jun Ji Formal analysis Qinxi Chen Project administration Jun Ji Visualization Qiang Liu; Qinxi Chen Roles/Writing -original draft Additional Supporting Information may be found online in the supporting information tab for this article.