key: cord-0914300-vtj0ylxa
authors: Zhu, Jie; Qi, Mingpu; Jiang, Chuanwen; Peng, Yongchong; Peng, Qingjie; Chen, Yingyu; Hu, Changmin; Chen, Jianguo; Chen, Xi; Chen, Huanchun; Guo, Aizhen
title: Prevalence of bovine astroviruses and their genotypes in sampled Chinese calves with and without diarrhoea
date: 2021-08-23
journal: J Gen Virol
DOI: 10.1099/jgv.0.001640
sha: 45346c7d44446ad4a9d253aafae119621c51452b
doc_id: 914300
cord_uid: vtj0ylxa

Bovine astrovirus (BoAstV) belongs to genus Mamastravirus (MAstV). It can be detected in the faeces of both diarrhoeal and healthy calves. However, its prevalence, genetic diversity, and association with cattle diarrhoea are poorly understood. In this study, faecal samples of 87 diarrhoeal and 77 asymptomatic calves from 20 farms in 12 provinces were collected, and BoAstV was detected with reverse transcription-polymerase chain reaction (RT-PCR). The overall prevalence rate of this virus in diarrhoeal and asymptomatic calves was 55.17 % (95 % CI: 44.13, 65.85 %) and 36.36 % (95 % CI: 25.70, 48.12 %), respectively, indicating a correlation between BoAstV infection and calf diarrhoea (OR=2.15, P=0.024). BoAstV existed mainly in the form of co-infection (85.53 %) with one to five of nine viruses, and there was a strong positive correlation between BoAstV co-infection and calf diarrhoea (OR=2.83, P=0.004). Binary logistic regression analysis confirmed this correlation between BoAstV co-infection and calf diarrhoea (OR=2.41, P=0.038). The co-infection of BoAstV and bovine rotavirus (BRV) with or without other viruses accounted for 70.77 % of all the co-infection cases. The diarrhoea risk for the calves co-infected with BoAstV and BRV was 8.14-fold higher than that for the calves co-infected with BoAstV and other viruses (OR=8.14, P=0.001). Further, the co-infection of BoAstV/BRV/bovine kobuvirus (BKoV) might increase the risk of calf diarrhoea by 14.82-fold, compared with that of BoAstV and other viruses (OR=14.82, P <0.001). Then, nearly complete genomic sequences of nine BoAstV strains were assembled by using next-generation sequencing (NGS) method. Sequence alignment against known astrovirus (AstV) strains at the levels of both amino acids and nucleotides showed a high genetic diversity. Four genotypes were identified, including two known genotypes MAstV-28 (n=3) and MAstV-33 (n=2) and two novel genotypes designated tentatively as MAstV-34 (n=1) and MAstV-35 (n=3). In addition, seven out of nine BoAstV strains showed possible inter-genotype recombination and cross-species recombination. Therefore, our results increase the knowledge about the prevalence and the genetic evolution of BoAstV and provide evidence for the association between BoAstV infection and calf diarrhoea.

Astroviruses (AstV) belong to the Astroviridae family, including two genera Mamastravirus (MAstV) and Avastrovirus (AAstV). MAstV can cause diarrhoea and neurological symptoms in mammals, while AAstV can lead to hepatitis, nephritis, and diarrhoea in birds [1] . Among MAstV, human AstV was first discovered to be OPEN ACCESS an etiological agent causing infantile diarrhoea in 1975, while bovine astrovirus (BoAstV) was first reported in 1978 [2] . In 1984, two American bovine strains (US1 and US2) antigenically related to a UK strain (L107) were demonstrated to cause infection and cytopathology of M cells in the dome epithelium [3] . The virulence of BoAstV related to calf diarrhoea was clinically supported in Hong Kong, China, in 2011 [4] , the mainland of China in 2013 [5] , and Brazil in 2015 [6] . However, some other evidence from calf experiments [2] and natural infection [7] has indicated that BoAstV is avirulent. Up to now, many studies have revealed that BoAstV may cause or exacerbate calf diarrhoea by co-infection with other enteric viruses such as bovine rotavirus (BRV), bovine torovirus (BToV) [3] , bovine viral diarrhoea virus (BVDV), bovine coronavirus (BCoV), and bovine kobuvirus (BKoV) [8] . However, the prevalence of BoAstV in diarrhoeal cases varied in different reports, such as 87.5 % in China [5] , 85.7 % in Japan [9] , 66.7 % in South Korea [8] , 64.1 % in Brazil in 2015 [6] , and 32.0 % in Egypt [10] . Since 2010, astroviruses have been expanded to at least 25 animal species including various terrestrial domestic animals, wild animals, and aquatic animals. This list includes cattle, pigs, rabbits, dogs, cats, turkeys, chicken, ducks, deer, minks, bats, sea lions, marine fish, etc. [1, [11] [12] [13] .

At the molecular level, AstVs are a group of small nonenveloped RNA viruses. Their genomes start from a 5′ untranslated region (5′UTR), followed by three open reading frames (ORF1a, ORF1b, and ORF2), a 3′ untranslated region (3′UTR), and a poly (A) tail [11] . According to the latest AstV classification criteria proposed by the International Committee on the Taxonomy of Viruses (ICTV) in 2012, the two genera MAstV and AAstV are further classified into different genotypes based on the mean amino acid genetic distance (p-dist) between 0.338 and 0.783 in ORF2 [14] . Specifically, MAstV includes 19 established genotypes (MAstV-1 to 19), 14 proposed genotypes (MAstV-20 to 33), and other undefined genotypes [15] .

Although BoAstV is one of the earliest discovered AstV, its genomes remain uncharacterized until four bovine strains from adult cattle without diarrhoea in Hong Kong, China, in 2011 [4] and one strain from yak in the Qinghai-Tibetan Plateau, China, in 2013 [16] were sequenced. So far, a total of 45 genome sequences from seven countries have been publicized in GenBank. Six genotypes (MAstV-13, MAstV-24, MAstV-28, MAstV-29, MAstV-30, and MAstV-33) have been identified from BoAstV isolates mostly based on reverse transcription-polymerase chain reaction (RT-PCR) and sequencing of the ORF2 [17] .

However, few studies have addressed the prevalence and genetic diversity of BoAstV and its association with calf diarrhoea, which might be due to the difficulties in isolating this virus from cell and tissue culture and high proportion of co-infection in both diarrhoeal and asymptomatic animals. Fortunately, the current next-generation sequencing (NGS) method can effectively replace the conventional cell culture to identify BoAstV and other co-infected viruses. Considering this, this study is aimed to determine the correlation between BoAstV and calf diarrhoea by investigating the prevalence and genomic characteristics of BoAstVs in Chinese calves with and without diarrhoea based on recent samples collected from a wide range of regions using NGS method. As a result, a strong correlation was observed between BoAstV infection and calf diarrhoea, and two novel genotypes were identified and designated as MAstV-34 and MAstV-35, in addition to the known MAstV-28 and MAstV-33. Our findings would be of significance for developing control measures against BoAstV related to calf diarrhoea.

The 12 provinces were randomly selected from 32 provincial regions in the mainland of China in 2019. Three intensive farms with a scale of over 100 cattle were randomly chosen from each province. First, the occurrence of calf diarrhoea was inquired by phone. In cases of calf diarrhoea occurrence, the farms were visited with the consent of the owners, and faecal samples were randomly collected from both the diarrhoeal and asymptomatic calves based on the principle of case-control study. A total of 164 faecal samples were collected from calves (under 6 weeks old) from 20 farms located in 12 provinces, including 87 from diarrhoeal calves and 77 from asymptomatic calves. The information about the samples was shown in Table 1 . All faecal samples were freshly collected from diarrhoeal calves and their asymptomatic controls in the same farms and stored at about 4 °C, and then they were immediately transported to the laboratory where they were stored at −80 °C until use.

One hundred and sixty-four faecal samples from diarrhoea and asymptomatic calves were tested for BoAstV through RT-PCR by previously reported method [7] . BoAstV-positive samples were divided into five regions based on the geographic locations of China (Table 1) . Each region contained two sample pools: the asymptomatic faecal sample pool and the diarrhoea faecal sample pool. Each pool contained ten samples with two exceptions (five samples in one pool and 15 samples in the other pool). One gram of faeces was taken from each sample and completely mixed by vortexing. Then, 0.5 g of each sample pool was taken, mixed with 5 ml of sterile phosphate buffer saline (PBS), and fully vortexed. Subsequently, the resultant samples were subjected to two freeze-thaw cycles at −80 °C/25 °C to release virus particles and centrifuged at 8000 g at 4 °C for 10 min. The supernatant was collected and filtered through 0.22 µm column filters to remove bacteria and other contaminants. According to the manufacturer's instructions, viral RNA was extracted with a QIAamp Viral RNA/DNA Mini Kit (QIAGEN, Germany), and the rRNA was removed with a QIAseq Fast Select-rRNA HMR Kit (QIAGEN, Germany). Afterwards, the cDNA library was prepared using a QIAseq FX DNA Library Kit (Qiagen, Germany) and commercially sequenced using a MiSeq bench-top sequencer (Illumina) with 151 paired-end reads in Wuhan We Find Biology Co., Ltd (Wuhan, China). The generated raw reads were subjected to quality control by trimming the reads with Phred quality score <10 and by filtering the reads containing adapters. The remaining raw reads were assembled with IDBA [18] and aligned against the Nr (NCBI non-redundant protein sequences) database using an E-value ≤0.001. The virome output was visualized and analysed with MEGAN 6 [19] .

To validate the genome sequences of nine BoAstV strains obtained from NGS, we performed the extra RT-PCR to amplify the whole genomes containing 5-9 fragments and sequencing the RT-PCR products. Briefly, 5-9 specific primer pairs were designed with the overlapping of 100-400 bp between adjacent fragments based on genome sequences of each virus obtained from NGS (Table S1 , available in the online version of this article). Nine calf faecal samples were preliminarily processed, as described above, and then viral RNA was extracted with a QIAamp Viral RNA/DNA Mini Kit (QIAGEN, Germany). Reverse transcription (RT) was performed using the reagents of the PrimeScript II Reverse Transcriptase (Takara, Otsu, Shiga, Japan), following the manufacturer's protocol. Subsequently, RT-PCR was performed to amplify the genome of each BoAstV strain with specific primers (Table S1 ). The RT-PCR mixture contained 25 µl of 2×PrimeSTAR Max Premix (Takara, Otsu, Shiga, Japan), 2 µl of cDNA, 19 µl of ddH 2 O, and 2 µl each of the primers (10 mM). The RT-PCR products were directly sequenced in Tsingke Biotechnology Co., Ltd. (Beijing, China). The fragments with weak bands were cloned into the pMD19-T vector and sequenced in Tsingke Biotech. The nearly complete genome sequences of nine BoAstV strains were obtained by assembling the fragment sequences acquired above. Finally, the sequences obtained by NGS were aligned against those obtained by Sanger sequencing to determine the complete genome sequences of nine BoAstV strains. These nine BoAstV genome sequences were submitted to GenBank with accession numbers of MW373712-MW373720. [20] . The tertiary structure of the RdRp was predicted using I-TASSER (https:// swissmodel. expasy. org/ interactive) based on previously published sequences. The tertiary structure obtained from I-TASSER was visualized with the PyMOL (http://www. pymol. org/).

Multiple sequences were aligned using the Clustal W method with mega 7 [21] . Neighbour-joining analysis of the nucleotide and amino acid sequences in different genomic regions of AstVs strains was performed. The unrooted phylogenetic trees were constructed by mega 7 with bootstrap values calculated for 1000 replicates [22] . The identities between nine BoAstV strains and other AstV strains at levels of both nucleotide and amino acid sequences were examined using the Clustal W algorithm in DNAstar MegAlign software. The information on reference strains, including the strains, species, GenBank accession numbers, isolation time and locations, was shown in Table S2 . Mean amino acid genetic distances of ORF2 among known MAstV 1-33, novel MAstV-34, and novel MAstV-35 were calculated using mega 7 with the pairwise deletion option and 1000 bootstrap replicates. Information on the representative strains of each genotype was displayed in Table S3 .

Similarity plots were drawn using the SimPlot software package (version 3.5.1) [23] . Bootscan analysis was performed using the neighbour-joining tree model by Kimura two-parameter distance algorithm with a window of 200 bp and a step of 20 bp. To determine possible recombination events, the recombination signals in the whole genome of the BoAstV isolates were measured by the Recombinant Detection Program (RDP4, v4.46) with seven methods (RDP, GENECONV, MaxChi, Bootscan, Chimaera, SiScan, and 3Seq) [24] . Briefly, the complete genome sequences of all BoAstV strains (including the isolates in the current study and representative strains from cattle and other species previously published) were inputted to RDP4 to search for recombination signals. The genomes with significant recombination event (P <0.05) derived from at least six methods were presented.

Based on NGS data, the sequences of other viral fragments from vertebrate viruses, plant viruses, and insect viruses were predicted. The possible viruses associated with calf diarrhoea was evaluated by reviewing previous reports [25] [26] [27] [28] , and RT-PCR was used to verify the existence of nine potential enteric viral pathogens including BVDV [29] , BRV [30] , BCoV [31] , bovine enterovirus (BEV) [32] , BToV [33] , BKoV [34] , bovine norovirus (BNoV) [35] , bovine nebovirus (BNeV) [36] , and mammalian orthoreovirus (MRV) [37] according to the methods described in the above-mentioned literatures. 

A chi-square test was used to assess the correlation between BoAstV infection and calf diarrhoea. The P <0.05 were considered statistically significant. Odds ratio (OR) and their 95 % confidence interval (CI) were also used to assess the degree of correlation between BoAstV and calf diarrhoea. Using SPSS software, the binary logistic regression method was used to analyse the effects of BoAstV co-infection with other nine viral pathogens (BVDV, BRV, BCoV, BToV, BNoV, BNeV, BKoV, BEV, and MRV) on calf diarrhoea.

The faecal samples from 20 farms in 12 provinces were examined with RT-PCR. Table 2 ).

In addition, the overall prevalence rate of BoAstV varied greatly with the individual provinces (each including 1-3 farms), ranging from 14.29 % to 75.00 % in all the sampled calves; while ranging from 25.00-100.00 % in the diarrhoeal calves, and from 0.00-75.00 % in asymptomatic calves (Table S4) .

Nearly complete genome sequences of nine BoAstV strains (including five strains from diarrhoeal calves and four strains from asymptomatic calves) from NGS data were assembled, and these nine BoAstV strains were distributed in seven provinces in China (Fig. 1) . The genome length of these nine strains ranged from 6048 to 6316 nucleotides (nt) excluding the poly(A) tail, and the guanine cytosine (GC) contents ranged from 48-53 % ( Table 4 ). The typical structures of genomes of all the nine BoAstV strains contained one 5′-UTR, three coding region ORFs (ORF1a, ORF 1b, and ORF2), and one 3′-UTR. Like other known astroviruses, these nine BoAstV strains possessed the overlapped region between ORF 1a and ORF 1b with a highly conserved frame-shift ribosomal heptameric sequence (5′-AAAAAAC-3′). In addition, strains BoAstV Henan-1 and Henan-2 exhibited two more amino acids in the highly conserved ORF1b region than seven other BoAstV strains (Table 4) .

To verify the genome sequences obtained by NGS method from the sample pools, we performed the extra RT-PCR to amplify the whole genomes containing 5-9 fragments (Fig. S1 ) and sequencing the RT-PCR products of each fragment. We further assembled the fragments of each strain to obtain a nearly full-length genome and found that the re-assembled genome sequences were the same with those generated by NGS (Fig. S2) , indicating the reliability of NGS results of pooled samples.

ORFs' amino acid and nucleotide sequences and genomewide sequences of nine BoAstV strains were compared with those of 15 BoAstV reference strains retrieved from GenBank (Table 5) . Among the three ORFs, the ORF1b exhibited a higher sequence identity than the other two, indicating that ORF1b was conserved in astroviruses. The genetic distances in the amino acids of ORF2 between nine strains and 43 reference MAstV strains ranged from 0.023 to 0.751 (Table 5) .

Likewise, the nine strains also showed significant difference in the sequence identity at the levels of both amino acids and nucleotides, and the genetic distance in the amino acids of ORF2 among these nine strains ranged from 0.049 to 0.557 (Table 6) . Fig. 2 . Neighbour-joining phylogenetic tree of AstVs based on ORF 2 amino acid sequences. Unrooted trees were generated using the neighbour-joining method with 1000 bootstrap replicates and sequence alignments were performed by using ClustalW in mega 7.0 software. The strains isolated in this study and the proposed two novel genotypes, MAstV-34 and MAstV-35 were marked in red. Strains were coloured by hosts and indicated by the host and strain names. Information on the reference strains was shown in Table S2 .

To reveal the evolutionary relationship between nine BoAstV strains and 96 other representative AstVs strains retrieved from the GenBank, a phylogenetic tree was constructed based on the nucleotide sequences of complete genomes, and amino acid sequences of RdRp, ORF1a, and ORF2. The results indicated that the nine BoAstV strains fell into four different clusters within the Mamastrovirus genus with ORF2 amino acid genetic distance (ranging from 0.338 to 0.783) as the criterion to define different genotypes (Fig. 2) [14] . Strains BoAstV Henan-2, HLJ-1, and Hebei-1 were clustered with Japanese and Italian BoAstV strains, and this cluster was reported as MAstV-28 [1] . Strains BoAstV HLJ-2 and Hunan-1 were clustered with Chinese and Japanese BoAstV strains, and this cluster was previously defined as MAstV-33 [1] . One BoAstV Henan-1 strain and two Japanese BoAstV strains exhibited the closest genetic distance, falling into one cluster, but the genetic distances between these three strains and 33 other known genotype strains ranged from 0.399 to 0.762 (Table 7) , thus we identified this cluster as a tentative novel genotype and referred to it as MastV-34. In addition, our data showed that three other strains (BoAstV JL-1, Hubei-1, and SD-1) were clustered together with two BoAstV strains isolated from Uruguay and eight BoAstV strains isolated from Japan. The genetic distance between these strains in this cluster and 33 other known strains ranged from 0.399 to 0.716 (Table 7) . Therefore, we identified this cluster as another tentative novel genotype and named it MAstV-35. The ORF2 amino acid genetic distances within these two novel genotypes MAstV-34 and MAstV-35 were 0.128-0.227 and 0.001-0.219, respectively. Furthermore, the genetic distances between the members of MAstV-34 and those of MAstV-35 ranged from 0.404 to 0.507 (Table 7) , indicating that our classification of two tentative novel genotypes was highly compliant with ICTV criteria for the determination of novel genotypes. Overall, the nine BoAstV strains in this study were clustered into five groups based on the complete genome nucleotide sequences (Fig. 3a) , ORF1a amino acid sequences (Fig. 3b) , and RdRp amino acid sequences (Fig. 3c ).

The complete genomes of the nine BoAstV strains were aligned against other published AstVs genomes by the ClustalW program in mega 7 for similarity plotting (Fig. 4) and recombination analyses (Table 8 ). All the nine BoAstV strains were used to identify potential genetic recombination sites with representative strains of other species retrieved from GenBank (until June 2020) as the reference sequences (Table S2 ). All the other seven BoAstV strains except BoAstV JL-1 and Henan-2 strains were identified to have obvious recombination events by at least six methods in RDP 4 program.

The similarity plot and RDP 4 program analysis indicated that both BoAstV Hubei-1 (MAstV-35) and SD-1 (MAstV-35) were identified as recombinants, and that all their parental strains belonged to MAstV-35 (Fig. 4a, b, Table 8 ), which was a common recombination within the same genotype to produce progeny viruses. In contrast, although both BoAstV Hebei-1 and HLJ-1 strains were identified as MAstV-28 and detected to have significant recombination events, their parental strains belonged to MAstV-33 (Fig. 4c , d, Table 8 ), suggesting that AstVs parental strains with the same genotype produced the progeny viruses with different genotypes through recombination. In addition, a recombination event was identified in BoAstV Henan-1 strain (MAstV-34), and its parental strains were identified as MAstV-34 and MAstV-28, respectively (Fig. 4e, Table 8 ). BoAstV Henan-1 strain recombination event supported inter-genotype recombination of AstV. To the best of our knowledge, this is the first report on inter-genotype recombination of BoAstV. In addition, BoAstV Hunan-1 and HLJ-2 strains were identified as MAstV-33, and they were observed to have significant recombination events (Fig. 4f , g, Table 8 ). The parental strains of BoAstV Hunan-1 belonged to BoAstV (MAstV-28) and roe deer AstV (MAstV- 33) , and that of BoAstV HLJ-2 belonged to roe deer AstV (unclassified genotype) and Sichuan takin AstV (MAstV-33). The recombination of BoAstV Hunan-1 and HLJ-2 strains illustrated that recombination events in AstV can occur across species and genotypes. The phylogenetic analyses confirmed the recombination events of BoAstV and their parental strains (Fig. S3 ).

We further investigated the differences in VPg and RdRp protein sequences (Fig. 5) . Although the VPg sequences at the end of ORF 1a were conserved ( 689 TEEEY 693 ) in MAstV, a few mutations were still observed. Compared with the other seven viruses, BoAstV Henan-1 and Henan-2 had unique mutation sites in VPg, which was similar to mink and dromedary AstV strains. Additionally, the strains were obtained from both diarrhoeal and asymptomatic calves. No correlation was observed between strain mutations and calf health conditions (diarrhoeal vs asymptomatic) (Fig. 5a) . The RdRp of BoAstV strains exhibited three types of tertiary structures with two similar to each other and one completely different across the nine strains. No correlation was observed between BoAstV strains' RdRp tertiary structures and calf health conditions (diarrhoeal vs asymptomatic) (Fig. 5b) . In other words, the nine investigated BoAstV strains exhibited genetic diversity in the amino acid sequences of VPg and protein structures of RdRp. *The reference strain information for each genotype was shown in Table S3 .

Diarrhoea in animals and humans has been reported to be associated with intestinal AstV infection [11, 39] , but whether BoAstV is the causative agent of bovine diarrhoea is controversial [2, 40] . Compared with previous studies of BoAstV, this research focused on the prevalence rate of faecal BoAstV from diarrhoeal and asymptomatic calves with the samples derived from 20 farms in 12 provinces of China based on a case-control study. We reported that BoAstV infection including single infection and co-infection was positively correlated with calf diarrhoea, which was in agreement with previous investigations [6, 16] . Although BoAstV co-infection was previously reported in different countries, such as Korea [8] , China [5, 16] , and Italy [41] , this study first quantified the correlation between different BoAstV co-infection patterns and calf diarrhoea. Our results revealed that the BoAstV/ BRV co-infection represented the most frequent co-infection pattern and increased the diarrhoea risk by 8.14-fold; the BoAstV/BRV/BKoV co-infection further increased the diarrhoea risk by 14.82-fold. In addition, BRV has been known as one of the common and important causative agent of calf diarrhoea [27, 38] , and its co-infection with BoAstV or even Fig. 3 . BoAstV strains were separated into five groups based on phylogenetic analysis. Neighbour-joining phylogenetic trees of AstV complete genome nucleotide (a), ORF 1a amino acid (b) and RdRp amino acid sequences (c) were shown respectively. Unrooted trees were generated using the neighbour-joining method with 1000 bootstrap replicates and sequence alignments were performed by using ClustalW in mega 7.0 software. The strains in this study were marked red. BoAstVs were coloured by groups and indicated by the hosts and strain names. Information on the reference strains was shown in Table S2 .

BoAstV/BKoV greatly increased the risk of calf diarrhoea.

These findings indicate that the combined vaccine against BoAstV/BRV/BKoV infection might be a potential control measure over calf viral diarrhoea in China.

In addition, the prevalence rate of BoAstV varied greatly with individual provinces. The main factors affecting the prevalence rate in individual provinces lay in farm management measures after the birth of calves such as the proper pen temperature, dry bedding materials, sufficient timely colostrum, and biosafety levels determining the source of the viruses. These farm management measures differed greatly between individual provinces. In addition, provincial difference in prevalence rate might be due to the different sample sizes from one to three farms in each province (Table 1) . *Nucleotide localization of the genomic segments implicated in recombination. In some instances, the region is shown in place of a nucleotide position. The 99 % confidence interval (CI) is

shown as the nucleotide segment from X to Y. †The seven tests for recombination were implemented in the RDP4 program. ‡Asterisks indicate the number of methods used to determine statistical evidence of a recombination event: *, six methods; **, seven methods. Genomes in which statistical evidence of a recombination event was obtained with fewer than six methods were not indicated.

Currently, BoAstVs are composed of five genotypes, namely, MAstV-13 [42] [43] [44] , MAstV 28-30, and MAstV-33 [4, 5, 7, 9, 41, 45, 46] . In this study, two known genotypes (MAstV-28 and MAstV-33) were determined by directly sequencing the genomes of the faecal samples and two novel genotypes (MAstV-34 and MAstV-35) were identified in terms of the criterion of ORF2 amino acid genetic distance [14] . Consistently, multiple genotypes of AstVs have been identified in other hosts, such as humans, bats, and pigs [1, 47] . However, to date, most BoAstVs have not been assigned to known genotypes (MAstV 1-33) [1, 47] . Therefore, more novel genotypes of BoAstVs might appear in the future. To facilitate the accurate classification of BoAstVs, we attempted to propose a method for classifying BoAstVs in this study based on the phylogenetic analyses of RdRp, ORF1a, and complete genome sequences. This method could cluster most current BoAstV strains (96.08 %, 49/51) into five independent branches, except two strains BoAstV-VC34/338 (accession number: MK987099) and BoAstV/JPN/Hokkaido12-25/2009 (accession number: LC047793). These two strains might be derived from cross-species transmission [9] which was similar to human AstV [48] and porcine AstV [49, 50] . By using our proposed method, seven AstV strains from roe deer, sheep, ovibos moschatus, and Sichuan takin were included and clustered into BoAstV group 1 and group 4. The above-mentioned ruminants were similar to cattle in species or breeding mode, indicating that these strains were closely related to BoAstVs. Since the Astroviridae family is an emerging family with significant genetic diversity, long evolution history, and a wide range of host species, the nomenclature and taxonomy of this virus family should be more fully discussed, agreed upon, and updated [47] . Abundant BoAstV gene sequences have been published in GenBank, and most sequences are relatively conservative amplicons with the length of 300-400 bp.

Only small fragment sequencing is not conducive to AstVs classification. Therefore, whole genome sequencing is highly suggested when possible. The subsequent comprehensive analysis of multiple regional sequences will be beneficial to classification, like the classification of the rotavirus and influenza virus [51, 52] .

Previous studies have shown that recombination events play an important role in the evolution of AstV [53, 54] . In this study, to identify significant recombination signals in the AstV genomes, the recombination events were identified by at least six methods in RDP 4 program coupled with a similarity plotting. In this study, BoAstV Hunan-1 and HLJ-2 strains were both observed to have very obvious recombination signals with their parental strains identified as BoAstV and roe deer AstV (CcAstV), Sichuan takin AstV and CcAstV, respectively. Our results were in line with the previous report on the cross-species recombination between human and California sea lion AstV strains [55] , between porcine AstV and human AstV strains [56] . Previously, it was speculated that there was a cross-species recombination event between BoAstV and CcAstV [4] . Our data confirmed that the recombination of BoAstVs occurred between cattle and deer. Recombination has been reported to require co-infection with different viruses in the same host cells [57] . It was worth noting that the roe deer AstV (CcAstV) strain (AstV SLO/D12-14) from Slovenia [58] and the Sichuan takin AstV strain from China [59] were the major and minor parents of the recombinant BoAstV HLJ-2, respectively. Since Sichuan takin and roe deer belonged to completely different species, and they were geographically separately, recombination of these two viruses hinted the occurrence of cross-species crossgeographical regions transmission event, which could be explained as follows. AstVs, as common water pollutants, might spread through the faecal-oral route [60] , and water used for breeding cattle could be contaminated by infected wild animals.

Our data showed that four out of seven recombination events occurred in the ORF1b-ORF2 junction region which contained conserved sequences and a stable hairpin structure [61, 62] . In this specific junction region, the recombination of AstVs has been reported in several studies [57, [63] [64] [65] [66] . The sequence conservation in this region increased the likelihood of homologous recombination, and the recombination events in the ORF2 region facilitated the swapping of viral capsid protein. Considering the role of capsid protein in defining the cell tropism and host range of the virus, the recombination events in this region might contribute to the separation of virus from its host species [4] . After across-species spreading, the genetic diversity of a virus is greatly increased by adapting to new host species. The genetic diversity of the viruses might further affect their virulence. The findings of virus genetic diversity have great significance for developing a novel vaccine and provide reference for molecular epidemiological studies.

This case-control study demonstrated that there was an association between BoAstV infection and calf diarrhoea, and a very strong positive correlation between co-infection of BoAstV/BRV/BKoV with calf diarrhoea. Further, the genomes of nine BoAstV strains were sequenced and clustered into four genotypes, including two novel genotypes MAstV-34 and MAstV-35 and two known genotypes MAstV-28 and MAstV-33. Finally, our data confirmed recombination of BoAstVs, which further hinted their cross-species transmission. These findings provide significant epidemiological evidence for developing novel measures to prevent and control BoAstV-induced diarrhoea in calves.

Bovine astrovirus B18/ HK

BoAstV/JPN/Kagoshima1-7/2014(MAstV-33)

BoAstV Henna-2(MAstV-28)

The broad host range and genetic diversity of mammalian and avian astroviruses

Isolation of small viruses resembling astroviruses and caliciviruses from acute enteritis of calves

Astrovirus and Breda virus infections of dome cell epithelium of bovine ileum

Rediscovery and genomic characterization of bovine astroviruses

Molecular epidemiology and phylogenetic analysis of diverse bovine astroviruses associated with diarrhea in cattle and water buffalo calves in China

Molecular detection and phylogenetic analysis of bovine astrovirus in Brazil

High prevalence and diversity of bovine astroviruses in the faeces of healthy and diarrhoeic calves in South West Scotland

Phylogenetic analysis of bovine astrovirus in Korean cattle

Full genome analysis of bovine astrovirus from fecal samples of cattle in Japan: identification of possible interspecies transmission of bovine astrovirus

Molecular detection of enteric viruses from diarrheic calves in Egypt

Astrovirus infections in humans and animals -molecular biology, genetic diversity, and interspecies transmissions

Virome composition in marine fish revealed by meta-transcriptomics

Genetic diversity and expanded host range of astroviruses detected in small mammals in Singapore

Ninth Report of the International Committee on Taxonomy of Viruses

The expanding field of mammalian astroviruses: Opportunities and challenges in clinical virology

A novel astrovirus species in the gut of yaks with diarrhoea in the Qinghai-Tibetan Plateau

Bovine astrovirus surveillance in Uruguay reveals high detection rate of a novel Mamastrovirus species

IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth

MEGAN community edition -interactive exploration and analysis of Large-Scale microbiome sequencing data

FoldIndex: a simple tool to predict whether a given protein sequence is intrinsically unfolded

CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

Prospects for inferring very large phylogenies by using the neighbor-joining method

Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination

RDP4: Detection and analysis of recombination patterns in virus genomes

An overview of calf diarrhea -infectious etiology, diagnosis, and intervention

Viral enteritis in calves

Case-control study of microbiological etiology associated with calf diarrhea

Genome sequence of the novel reassortant mammalian orthoreovirus strain MRV00304/13, isolated from a calf with diarrhea from the United States

Prevalence study and genetic typing of bovine viral diarrhea virus (BVDV) in four bovine species in China

High prevalence and genomic characteristics of G6P[1] Bovine Rotavirus A in yak in China

Detection and molecular characteristics of neboviruses in dairy cows in China

Molecular identification and characterization of a new type of bovine enterovirus

Epidemiological analysis of bovine torovirus in Japan

Three clusters of bovine kobuvirus isolated in Korea

Reverse transcription-PCR assays for detection of bovine enteric caliciviruses (BEC) and analysis of the genetic relationships among BEC and human caliciviruses

First detection of Nebovirus and Norovirus from cattle in China

Isolation and identification of bat viruses closely related to human, porcine and mink orthoreoviruses

Pathophysiology of diarrhea in calves

Pathogenesis of astrovirus infection

Serotypes of bovine astrovirus

Identification of astroviruses in bovine and buffalo calves with enteritis

Identification of a second encephalitis-associated astrovirus in cattle

Nonsuppurative (Aseptic) Meningoencephalomyelitis associated with neurovirulent astrovirus infections in humans and animals

The first case of Bovine astrovirus-associated Encephalitis in the Southern hemisphere (Uruguay), uncovers evidence of viral introduction to the Americas from

Fecal shedding of Bovine astrovirus CH13/NeuroS1 in veal calves

Full-genome based molecular characterization of encephalitis-associated bovine astroviruses

Human astroviruses

Epidemiology of classic and novel human Astrovirus: Gastroenteritis and beyond

Multiple novel and prevalent astroviruses in pigs

Whole genome analysis of porcine astroviruses detected in Japanese pigs reveals genetic diversity and possible intra-genotypic recombination

Review of G and P typing results from a global collection of rotavirus strains: implications for vaccine development

Uniformity of rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group (RCWG)

Recombination analysis based on the HAstV-2 and HAstV-4 complete genomes

Next-generation amplicon sequencing identifies genetically diverse human astroviruses, including recombinant strains, in environmental waters

Characterization of phylogenetically diverse astroviruses of marine mammals

Genomic analysis of two ORF2 segments of new porcine astrovirus isolates and their close relationship with human astroviruses

Evidence of a recombinant wild-type human astrovirus strain from a Kenyan child with gastroenteritis

A molecular survey, whole genome sequencing and phylogenetic analysis of astroviruses from roe deer

Metagenomic analysis of Sichuan takin fecal sample viromes reveals novel enterovirus and astrovirus

Astrovirus survival in drinking water

Molecular epidemiology of human astrovirus diarrhea among children from a periurban community of Mexico City

High evolutionary rate of human astrovirus

Genetic heterogeneity and recombination in human type 2 astroviruses

Lineage diversification and recombination in type-4 human astroviruses

Genetic heterogeneity and recombination in type-3 human astroviruses

Molecular characterization of a novel recombinant strain of human astrovirus associated with gastroenteritis in children

The authors declare that there are no conflicts of interest.

The protocol (HZAUMO-2018-027) regarding animal sample collection was approved by the Committee on the Ethics of Animal Experiments at Huazhong Agricultural University and conducted in strict accordance with the Guide for the Care and Use of Laboratory Animals, Hubei Province, China.

Five reasons to publish your next article with a Microbiology Society journal 1 . The Microbiology Society is a not-for-profit organization. 2. We offer fast and rigorous peer review -average time to first decision is 4-6 weeks. 3. Our journals have a global readership with subscriptions held in research institutions around the world. 4. 80% of our authors rate our submission process as 'excellent' or 'very good'. 5. Your article will be published on an interactive journal platform with advanced metrics.Find out more and submit your article at microbiologyresearch.org.