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

Surveillance of Human Rotaviruses in Wuhan, China (2019-2022): Whole-Genome Analysis of Emerging DS-1-like G8P[8] Rotavirus

Group A rotaviruses (RVAs) are major etiologic agents of gastroenteritis in infants and young children worldwide. To study the prevalence and genetic characteristics of RVAs, a hospital-based surveillance study was conducted in Wuhan, China from June 2019 through May 2022. The detection rates of RVAs were 19.40% (142/732) and 3.51% (8/228) in children and adults, respectively. G9P[8] was the predominant genotype, followed by G8P[8] and G3P[8]. G8P[8] emerged and was dominant in the 2021-2022 epidemic season. The genome constellation of six G8P[8] strains was assigned to G8-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2. Phylogenetic analysis revealed that the VP7, VP4, VP2, VP3, NSP1, NSP2, NSP3, and NSP5 genes of these G8P[8] strains clustered closely with those of the G8P[8] strains in Asia and were distant from those of the P[8] and G2P[4] strains simultaneously detected in Wuhan. In contrast, the VP1, VP6, and NSP4 genes were closely related to the typical G2P[4] rotavirus, including those of G2P[4] strains simultaneously detected in Wuhan. The detection rate of RVAs decreased in the COVID-19 pandemic era. It was deduced that the G8P[8] rotaviruses that emerged in China may be reassortants, carrying the VP6, VP1, and NSP4 genes derived from the G2P[4] rotavirus in the backbone of the neighboring DS-1-like G8P[8] strains represented by CAU17L-103.

Isolation and Characteristics of a Novel Aichivirus D from Yak

Aichivirus D (AiV-D) is a newly emerging Kobuvirus detected in bovine and sheep, and information is limited regarding its biological significance and prevalence. This study aimed to explore both the prevalence and characteristics of AiV-D in yaks. From May to August 2021, 117 fecal samples were collected from yaks with diarrhea in three provinces of China's Qinghai-Tibet Plateau, 15 of which were selected and pooled for metagenomic analysis. A high abundance of AiV-D sequences was obtained. Of the 117 diarrhea samples, 29 (24.8%) tested AiV-D-positive, including 33.3% (14/42) from Sichuan, 21.1% (8/38) from Qinghai, and 18.9% (7/37) from Tibet, respectively, suggesting a wide geographical distribution of the AiV-D in yaks in the Qinghai-Tibet Plateau. Furthermore, three AiV-D strains were successfully isolated using Vero cells. Significantly, the AiV-D strain could cause diarrhea, intestinal bleeding, and inflammation in yak calves via oral inoculation. The virus was distributed in the ileum, jejunum, duodenum, colon, cecum, and rectum. Based on phylogenetic analysis of the genome and capsid protein P1 (VP0, VP3, and VP1 genes), the yak AiV-D strains likely represent a novel genotype of AiV-D. On the whole, this study identified a novel genotype of AiV-D from yaks, which was successfully isolated, and confirmed that this virus is a diarrhea pathogen in yaks and has a wide geographical distribution in the Qinghai-Tibet Plateau. Our results expand the host range of AiV-D and the pathogen spectrum of yaks and have significant implications for diagnosing and controlling diarrhea in yaks. IMPORTANCE In this study, we identified and successfully isolated a novel genotype of AiV-D from yaks. Animal infection confirmed that this virus can cause diarrhea, intestinal bleeding, and inflammation in yak calves via oral inoculation. The virus was distributed in the ileum, jejunum, cecum, duodenum, colon, and rectum. All of these results have significant implications for diagnosing and controlling diarrhea in yaks. These novel AiV-D strains have a wide geographical distribution in yaks from the Qinghai-Tibet Plateau in China. In addition to expanding the host range of AiV-D and the pathogen spectrum of yaks, these findings can increase knowledge of the prevalence and diversity of AiV-D.

Genetic Diversity of Bovine Group A Rotavirus Strains Circulating in Korean Calves during 2014 and 2018

The purpose of this study was to investigate annual changes in BoRVA strains by examining the VP4 and VP7 genes of rotaviruses in Korean calves. Between 2014 and 2018, 35 out of 138 samples of calf diarrhea feces collected nationwide were positive for BoRVA. Further genetic characterization of the VP7 and VP4 genes of 35 BoRVA isolates identified three different G-genotypes (G6, G8, and G10) and two different P genotypes (P[5] and P[11]). The G6 genotype was most common (94.3%) in BoRVA-positive calves, followed by the P[5] genotype (82.9%). Four genotypes comprised combinations of VP4 and VP7: 80% were G6P[5], 14.2% were G6P[11], 2.9% were G8P[5], and 2.9% were G10P[11]. Susceptibility to infection was highest in calves aged < 10 days (35%) and lowest in calves aged 30−50 days (15.4%). The data presented herein suggest that the G6P[5] genotype is the main causative agent of diarrhea in Korean calves. In addition, it is predicted that G6P[5] will continue to act as a major cause of diarrhea in Korean calves.

Rapid detection of bovine rotavirus a by isothermal reverse transcription recombinase polymerase amplification assays

Background

Bovine rotavirus A (BRVA) is considered to be the most common pathogen of severe diarrhea in cattle worldwide, which could lead to the death of newborn calves and cause the significant economic losses to the cattle industry. As a novel isothermal nucleic acid amplification technique, recombinase polymerase amplification (RPA) has been applied widely for the rapid detection of different important pathogens in human and animals.

Results

An RT-RPA assay based on the real time fluorescence monitoring (real-time RT-RPA) and an RT-RPA assay combined with a lateral flow strip (LFS RT-RPA) were successfully developed by targeting the VP6 gene of BRVA. The RT-RPA assays allowed the exponential amplification of the target fragment in 20 min. After incubation of the LFS RT-RPA on a metal bath at 40 °C, the results were displayed on the lateral flow strip within 5 min, while real-time RT-RPA allowed the real-time observation of the results in Genie III at 42 °C. Both of the two assays showed high specificity for BRVA without any cross-reaction with the other tested pathogens causing diarrhea in cattle. With the standard RNA of BRVA serving as a template, the limit of detection for real-time RT-RPA and LFS RT-RPA were 1.4 × 10² copies per reaction and 1.4 × 10¹ copies per reaction, respectively. In the 134 fecal samples collected from cattle with diarrhea, the BRVA positive rate were 45.52% (61/134) and 46.27% (62/134) in real-time RT-RPA and LFS RT-RPA, respectively. Compared to a previously published real-time PCR, the real-time RT-RPA and LFS RT-RPA showed a diagnostic specificity of 100%, diagnostic sensitivity of 98.39% and 100%, and a kappa coefficient of 0.985 and 1.0, respectively.

Conclusions

In this study, BRVA was successfully detected in cattle fecal samples by the developed real-time RT-RPA and LFS RT-RPA assays. The developed RT-RPA assays had great potential for the rapid detection of BRVA in under-equipped diagnostic laboratory and the point-of-need diagnosis at quarantine stations and farms, which is of great importance to control BRVA-associated diarrhea in cattle herds.

Prevalence of bovine rotavirus among Bovidae in China during 1984-2021: A systematic review and meta-analysis

Bovine rotavirus (BRV) is a potential zoonotic intestinal pathogen that brings a serious threat to calf health, and has resulted in huge economic losses to China's breeding industry. Here, a systematic review and meta-analysis was conducted to estimate the prevalence of BRV among Bovidae from 1984 to 2021 in China. A total of 64 publications on BRV investigation in China were screened from the databases Chinese National Knowledge Infrastructure (CNKI), Wan Fang Database, Technology Periodical Database (VIP), PubMed, and ScienceDirect. The random-effect model was used to calculate the pooled prevalence of BRV, and the analyzed data were derived from 25 provinces in China. The estimated pooled prevalence of BRV in China was 35.7% (8176/17,292). In addition, the prevalence of BRV in Southwestern China (77.1%; 2924/3600) was significantly higher than that in other regions of China. Regarding geographic and climatic factors, the prevalence of BRV in the subgroup of latitude 30-35° (76.8%; 3303/4659) was significantly higher than that in the subgroup of latitude less than 30° (37.0%; 485/1275) or more than 35° (32.6%; 1703/5722), while the prevalence of BRV in the subgroup of longitude 100-105° (75.4%; 2513/3849) was significantly higher than that in the subgroup of longitude less than 100° (32.6%; 619/2255) or more than 105° (48.9%; 2359/5552). Rainfall was positively correlated with the prevalence of BRV, whereas temperature was negatively correlated with the positive rate of BRV (P < 0.05). Our data showed that the prevalence of BRV was strongly correlated with geographical and climatic conditions. Thus, we recommend that the corresponding prevention and control programs should be formulated according to different geographical conditions. The strengthening of BRV surveillance in areas with high altitude, low temperature, and heavy rainfall may contribute to the decrease of the incidence of BRV infection among Bovidae herds in China.

Effects of Short-Chain Fatty Acid Modulation on Potentially Diarrhea-Causing Pathogens in Yaks Through Metagenomic Sequencing

Short-chain fatty acids (SCFA) are principal nutrient substrates of intestinal epithelial cells that regulate the epithelial barrier in yaks. Until now, metagenomics sequencing has not been reported in diarrheal yaks. Scarce information is available regarding the levels of fecal SCFA and diarrhea in yaks. So, our study aims to identify the potential pathogens that cause the emerging diarrhea and explore the potential relationship of short-chain fatty acids in this issue. We estimated diarrhea rate in yaks after collecting an equal number of fecal samples from affected animals. Metagenomics sequencing and quantitative analysis of SCFA were performed, which revealed 15%–25% and 5%–10% prevalence of diarrhea in yak’s calves and adults, respectively. Violin box plot also showed a higher degree of dispersion in gene abundance distribution of diarrheal yaks, as compared to normal yaks. We found 366,163 significant differential abundance genes in diarrheal yaks, with 141,305 upregulated and 224,858 downregulated genes compared with normal yaks via DESeq analysis. Metagenomics binning analysis indicated the higher significance of bin 33 (Bacteroidales) (p < 0.05) in diarrheal animals, while bin 10 (p < 0.0001), bin 30 (Clostridiales) (p < 0.05), bin 51 (Lactobacillales) (p < 0.05), bin 8 (Lachnospiraceae) (p < 0.05), and bin 47 (Bacteria) (p < 0.05) were significantly higher in normal yaks. At different levels, a significant difference in phylum (n = 4), class (n = 8), oder (n = 8), family (n = 16), genus (n = 17), and species (n = 30) was noticed, respectively. Compared with healthy yaks, acetic acid (p < 0.01), propionic acid (p < 0.01), butyric acid (p < 0.01), isobutyric acid (p < 0.01), isovaleric acid (p < 0.05), and caproic acid (p < 0.01) were all observed significantly at a lower rate in diarrheal yaks. In conclusion, besides the increased Staphylococcus aureus, Babesia ovata, Anaplasma phagocytophilum, Bacteroides fluxus, viruses, Klebsiella pneumonia, and inflammation-related bacteria, the decrease of SCFA caused by the imbalance of intestinal microbiota was potentially observed in diarrheal yaks.

First Detection and Genomic Characterization of Bovine Norovirus from Yak

Yak are a unique free-grazing bovine species in high-altitude areas. The objective of this study was to investigate the presence and molecular characteristics of BNoV in yak. A total of 205 diarrheal samples of yak (aged ≤ 3 months) were collected from 10 farms in Sichuan Province, China, from May 2018 to October 2020, and four samples were detected as BNoV-positive with RT-PCR. Moreover, a nearly full-length genome of SMU-YAK-J1 containing three complete ORFs was successfully sequenced. Sequence analysis with only nine genome sequences of the GIII genogroup showed that SMU-YAK-J1 was most closely related with GIII.P2 GIII.4, sharing 90.9% gnomic nucleotide identity, but only shared 71.6–85.9% with other genotypes, which confirmed that SMU-YAK-J1 belongs to genotype GIII.P2 GIII.4. However, compared with the sole genome of GIII.4 in GenBank, the BNoV in this study also exhibited many unique amino acid changes among all the three ORFs, which may represent the unique genetic evolution of BNoV in yak. This study first determined the presence of BNoV in yak, contributing to a better understanding of the prevalence and genetic evolution of BNoV.

Isolation and Characterization of Bovine RVA from Northeast China, 2017-2020

Group A rotaviruses (RVAs) are major enteric pathogens causing infections in calves. To investigate the epidemiological characteristics and genetic diversity of bovine rotavirus (BRV), 233 fecal samples were collected from calves with diarrhea in northeast China. The samples were analyzed for sequences encoding the inner capsid protein VP6 (subgroup) and the outer capsid proteins VP7 and VP4 (G and P type, respectively) using RT-PCR. Ten of the 233 samples (4.3%) were identified as BRV positive and were used for virus isolation and sequence analysis, revealing that all strains analyzed were of the G6P[1] genotype. The isolates exhibited high VP6 sequence identity to the USA cow RVA NCDV strain (>99% amino acid identity) and were further shown to be closely related to Japanese cow RVA BRV101 and Israelian human RVA G6P[1] strains, with >99% amino acid identity to VP7 and VP4 proteins, respectively. Comparative analyses of genome-predicted amino acid sequences between the isolates and the NCDV strains indicated that the antigenicity and infectivity of the strains isolated had changed. In this study, BRV genotypes and the genetic diversity among vaccinated cattle herds were monitored to provide epidemiological data and references for early diagnosis, allowing for early detection of new, potentially pathogenic RVA strains.

Prevalence of bovine astroviruses and their genotypes in sampled Chinese calves with and without diarrhoea

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.

Detection and molecular characteristics of bovine rotavirus A in dairy calves in China

Background

Bovine group A rotavirus (BoRVA) is a major cause of severe gastroenteritis in newborn dairy calves. Only one study has investigated the G and P genotypes among dairy calves in a few regions of China, which were G6 and P[5]. Therefore, data on the prevalence and molecular characteristics of BoRVA in dairy calves in China remains limited.

Objectives

The purpose of this study was to investigate the prevalence and molecular characteristics of BoRVA in dairy calves in China.

Methods

269 dairy calves diarrheic samples from 23 farms in six provinces in China were collected to detect BoRVA using reverse transcription polymerase chain reaction.

Results

71% of samples were determined to be BoRVA-positive. Two G genotypes (G6, G10) and two P genotypes (P[1], P[5]) were identified, and G6P[1] BoRVA was the predominant strain. Moreover, the VP7 and VP4 gene sequences of these dairy calf BoRVA strains revealed abundant genetic diversity. Interestingly, eight out of 17 complete G6 VP7 sequences were clustered into G6 lineage VI and analysis showed the strains were closely related to Chinese yak BoRVA strains.

Conclusions

The results of this study show that BoRVA circulates widely among dairy calves in China, and the dominant genotype in circulation is G6P[1], first report on molecular characteristics of complete P[5] VP4 genes in chinese dairy calves. These results will help us to further understand the prevalence and genetic evolution of BoRVA among dairy calves in China and, thus, prevent the disease more effectively.

Identification and genomic characterization of bovine parvovirus 1 in yaks

Bovine parvovirus 1 (BPV1) is a causative agent of respiratory, gastrointestinal, and reproductive cattle diseases. We collected 149 yak diarrhea fecal samples from 9 farms in the Qinghai-Tibet Plateau. The samples were screened for BPV1 by PCR, and 2 samples were positive for BPV1. The complete genomes of these BPV1 isolates were sequenced successfully. The sequences of these 2 variants were both 5,515 bp in length and shared 96.5-96.8% identity with 2 previously reported BPV1 genomes (GenBank DQ335247, NC_001540). Twenty-six identical amino acid mutations were found in the 2 yak variants, including 7 amino acid substitutions in receptor-binding regions of the VP2 protein, and 5 amino acid substitutions in the NS1 protein C-terminal region that functions to activate transcription. The new genome sequences contribute to better understanding of the evolution and molecular characteristics of BPV1.

Development and use of a reverse transcription insulated isothermal PCR assay for detection and characterization of bovine torovirus in yaks

Bovine torovirus (BToV) is an important diarrhea-causing pathogen affecting bovines. To facilitate BToV detection, a reverse transcription insulated isothermal PCR (RT-iiPCR) assay was developed that targets the BToV M gene with high specificity and reproducibility. The assay has a limit of detection of 23 copies/μL. Out of 69 diarrheic fecal samples from yaks collected on six farms in Tibet and Sichuan provinces in China, 11.59% (8/69) tested positive for BToV using this assay. The full-length spike (S) and hemagglutinin-esterase (HE) genes of three positive samples were subsequently sequenced. Notably, an identical recombination event was identified in the S1 subunit of the S protein of three isolates. All of the HE genes were found to belong to genotype III and shared the same unique aa variation (P44S) in the esterase domain. This study is the first confirmation of BToV in yaks and the first report of an S gene recombination event in BToV. Our findings will enhance the current understanding of the molecular characteristics and genetic evolution of BToV.