Abstract
Group A porcine rotavirus (RVA) is a serious threat to the breeding industry worldwide, which was associated with severe diarrhea in piglets. However, the prevalence and molecular characterizations of RVA circulating in farms of East China remains largely unknown. Five hundred and ninety-four samples were collected from 35 farms in East China from September 2017 to December 2019. The results showed that 16.8% was positive for RVA of all samples. Among different types of samples, the highest positive rate of RVA was intestinal samples (19.5%), and among pigs at different growth stages, the highest detection rate of RVA in piglets was 18.5%. Furthermore, the VP7 and VP4 genes of nine positive samples were sequenced for alignment and phylogenetic analysis. Phylogenetic analysis revealed that the nine isolates belong to four kinds of genotype combinations correspondingly: G9P[7](5/9), G5P[13](2/9), G9P[13](1/9), and G5P[7](1/9).The data suggested that multiple genotypes combinations of RVA were circulating in pigs in East China. Thus, it's necessary to continuously survey the prevalence of RVA in pigs, aiding the rational application of vaccines or other measures for the prevention and control of RVA spread.
1. Introduction
Rotaviruses (RVs) are one of the main pathogens implicated in acute diarrhea in children and young animals, including nursing and weaned piglets (–). RVs mainly infect and disrupt small intestinal mature enterocytes and enteroendocrine cells leading to acute gastroenteritis, which are transmitted via the fecal-oral route (). Porcine rotavirus (PoRV) was first isolated from infected pigs in 1976 (). Subsequently, many studies documented the widespread prevalence of PoRV infections all over the world (). PoRV had resulted in enormous economic loss in the worldwide swine industry since it was discovered (, ).
RVs are non-enveloped, segmented dsRNA viruses, which belong to the member of the family Reoviridae. Its genome was consisted of 11 segments encoding 6 structural proteins (VP1-VP4, VP6, and VP7) and 5 non-structural proteins (NSP1-NSP5/6). RVs are classified into 10 groups (A–J) based on the antigenic characteristics of VP6 gene (). Besides, the G and P dual typing system was established, based on the outer capsid proteins VP7 and VP4 (). To date, group A rotavirus (RVA), RVB, RVC, RVE, and RVH have been described in pigs (, ). However, among these groups of RVs, RVA represents the main cause of acute diarrhea in piglets because of its high prevalence and pathogenicity (). Previous study showed that the prevalence rates of RVA in pigs vary from 3.3% to 67.3% (). Furthermore, 12 G genotypes (G1-G6, G8-G12, and G26) and 18 P genotypes (P[1]-P[8], P[11], P[13], P[19], P[23], P[25], P[26], P[27], P[32], P[34], and P[49]) of RVA have been detected in pigs (–).
Until now, the surveillance of pigs RVA remains rare in China. Previous studies showed that the prevalence rate of porcine RVA was 28.76% in Shandong province and the dominant genotypes were G3, G5 and G9 (). However, the prevalence genotypes and molecular characterization of RVA circulating in pigs in East China remains unknown. In this study, we investigated the prevalence and genetic characterization of porcine RVA in pig herds from East China (Jiangsu, Anhui, Shanghai, Zhejiang, and Shandong) from September 2017 to December 2019.
2. Materials and methods
2.1. Sample collection
From September 2017 to December 2019, a total of 594 samples (fecal swabs and small intestine tissue samples) were collected from 35 different pig farms in five provinces/cities (Jiangsu, Anhui, Shanghai, Zhejiang, and Shandong) in East China (details are summarized in Supplementary Table S1). Samples were collected from both healthy and sick pigs in different growing stages (within 30 days of piglets, 35–50 days of nursery pigs, 60–120 days of fattening pigs). There were no repeated samples collected from those farms. All collected samples were diluted in sterilized phosphate buffered saline (PBS), and centrifuged at 8,000 × g for 10 min. The supernatants were collected for further analysis.
2.2. RNA extraction and RT-PCR
Total RNA was extracted from pretreated samples using FastPure Cell/Tissue Total RNA Isolation Kit according to manufacturer's instructions (Vazyme, China). One-step RT-PCR was performed by using HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, China) following the manufacturer's protocol. The partial VP6 gene and full-length VP7 and VP4 genes were amplified by using 2× Taq Master Mix (Vazyme, China). The primers sequences were presented in Supplementary Table S2. Briefly, a total of 12.5 μl 2× Taq Master Mix were mixed with the specific primer pair for the individual genes, 5 μl cDNA of each sample used as template, and nuclease-free water to 25 μl. The reaction conditions were as follows: initial denaturation at 95°C for 5 min, followed by 35 cycles of denaturation at 95°C for 30 s, annealing for 30 s, and extension at 72°C for 3 min, followed by extension fully at 72°C for 10 min. All RT-PCR products were analyzed on 1% agarose gels containing 0.5 μg/ml ethidium bromide and visualized under UV transilluminator.
2.3. VP7 and VP4 gene sequencing and phylogenetic analysis
Each PCR product was purified with the E.Z.N.A.® Gel Extraction Kit according to manufacturer's instruction (Omega, China). The purified PCR products of VP7 and VP4 gene were cloned to pMD19-T vector (Takara, China), transformed to DH5α competent cells, and then cultured overnight at 37°C on Ampicillin agar plates. Recombinant DNA clones identified as positive for the respective gene by using PCR and then were sequenced by Sangon Biotech (Shanghai, China). The sequences were determined using Sanger sequencing and deposited in GenBank with accession numbers OP454313–OP454330, respectively (Supplementary Table S3). Sequences data were aligned with reference strains from GenBank using Clustal W and were analyzed using the software MEGA 11 (version 11.0.13) for phylogenetic analysis (the reference strains information was summarized in Supplementary Tables S4, S5).
3. Results
3.1. The prevalence of RVA in East China
A total of 594 samples were collected for detection and 100 out of 594 (16.83%) were positive for RVA by established RT-PCR. Furthermore, it was demonstrated that the positive rates of 2017, 2018, and 2019 were 11.8% (15/127), 20.8% (79/380) and 6.9% (6/87), respectively (Figure 1A). In regard to regional analysis, it was observed that the highest RVA prevalence was Anhui province (82.4%, 42/51), followed by Shanghai (50.0%, 1/2), Jiangsu (12.6%, 31/247), Shandong (11.0%, 25/228), and Zhejiang (1.5%, 1/66), respectively (Figure 1B). Besides, we also analyzed the positive rates from different kinds of samples. As shown in Figure 1C, the detection rates of intestinal samples and anal swab samples were 19.5% (68/349) and 12.7% (31/245), respectively. Also, the detection rates of RVA in relation to the age of the infected pigs were analyzed. It was demonstrated that the highest proportion was detected among piglets (18.5%, 70/379), followed by fattening pigs (16.2%, 21/130), nursery pigs (10.0%, 1/10), and sows (9.3%, 7/75), respectively (Figure 1D).
Figure 1
3.2. Phylogenetic analysis of the VP7 gene
To further investigate the genotypes of RVA strains circulating in East China, nine positive samples from different provinces or cities were chosen to analyze VP7 and VP4 genes by using RT-PCR. The VP7 genes of the nine different positive samples were sequenced for genotyping. Subsequently, a phylogenetic tree was established based on the VP7 gene sequences of the nine positive samples and publicly available reference strains from GenBank. Phylogenetic analysis revealed that the nine isolates from different areas and years belong to two distinct genetic groups (Figure 2). Six of nine rotavirus isolates clustered in G9 group with nucleotide p-distances ranging from 0.012 to 0.020 to the reference strain, while the remaining three isolates belonged to G5 group with nucleotide p-distances ranging from 0.127 to 0.147 to the reference strain. The VP7 sequences of CHN-SH1701, CHN-AH1702, CHN-SD1703, CHN-JS1801, CHN-ZJ1802, and CHN-JS1902 strains were most similar to NJ2012 strain (98.0%−99.9%). CHN-JS1901, CHN-AH1903, and CHN-ZJ1904 shared high homology with LNCY1 strain of 85.3%−87.3%.
Figure 2
3.3. Phylogenetic analysis of the VP4 gene
At the same time, the genotypes of VP4 gene of the nine isolates were also analyzed based on sequencing and phylogenetic analysis. As shown in Figure 3, the phylogenetic tree indicated that the nine isolates belonged to two distinct genetic groups. Six of the nine rotavirus isolates were classified into P7 group with nucleotide p-distances ranging from 0.007 to 0.010 to the reference strain, and the rest three isolates belonged to P13 group with nucleotide p-distances ranging from 0.064 to 0.073 to the reference strain. CHN-SH1701, CHN-AH1702, CHN-SD1703, CHN-JS1801, CHN-ZJ1802, and CHN-ZJ1904 strains showed higher similarities (98.9%−99.3%) with NJ2012 strain. CHN-JS1901, CHN-JS1902, and CHN-AH1903 showed more closed relationship with ET8B strain (92.5%−93.5%). In summary, it was concluded that G9P[7] (five of nine) was the outstanding dominant genotype combination among the sequenced positive samples, followed by G5P[13] (two of nine), G9P[13] (one of nine) and G5P[7] (one of nine) (Table 1).
Figure 3
Table 1
| Year | G and P genotype | Number | Area |
|---|---|---|---|
| 2017 | G9P[7] | 1 | Shanghai |
| 1 | Anhui | ||
| 1 | Shandong | ||
| 2018 | G9P[7] | 1 | Jiangsu |
| 1 | Zhejiang | ||
| 2019 | G5P[13] | 1 | Jiangsu |
| G5P[13] | 1 | Anhui | |
| G9P[13] | 1 | Jiangsu | |
| G5P[7] | 1 | Zhejiang |
Summary of genotype combinations of the sequenced samples.
4. Discussion
Porcine diarrhea caused by viral infections continuously damages herd enteric health and leads to enormous economic loss in the global swine industry (). However, as one of the most important pathogens causing swine acute gastroenteritis, it is the economic impact of RVA infection that cannot be ignored (, , ). Considering the impact of RVA in the swine industry, it is essential to investigate the prevalence and genetic characterization of porcine RVA in pig herds.
During 2017–2019, we investigated the presence of RVA from 35 farms collected from five provinces/cities in East China, and then carried out the phylogenetic analysis of some isolates from different areas. Our results showed that the detection rate of RVA from these farms was 16.83%. The positive rate of RVA in this study was consistent with the results presented in other studies performed in other areas in China during the same period, which suggested that PoRV seemed to be one of the most important pathogens causing diarrhea in pigs in China (). At the same time, we also analyzed the detection rate in different years. It was noted that there was a lower detection rate (6.9%) in 2019 compared with other years. This was probably caused by the fewer samples collected from pig farms since African swine fever virus (ASFV) circulation in China. Besides, the prevalence rates of RVA in pigs varied greatly from 1.5% to 82.4% among different areas. It cannot be ignored that the impact of RVA in swine industry in Anhui province, for which the positive rate was as high as 82.4%. The samples sources were included intestinal samples and anal swab samples, and the detection rates of intestinal samples were higher than swab samples. Furthermore, RVA infections were more common in piglets and fattening pigs compared with nursery pigs and sows. These observations were also consistent with previous published study ().
Although previous studies had investigated the prevalence rates of PoRV in several areas in China, very limited genotypes data of dominant PoRV isolates were available until now (, , ). The G/P genotype combinations of RVA identified mainly included G9P[7], G3P[13], G5P[13], G9P[13], G9P[23], G1P[7], and G9P[6] based on several isolate strains of publications during past decades in China (, –). Considering the rare data of genotypes, the VP7 gene and VP4 gene from the nine RVA strains were sequenced to demonstrate the genotypic classifications in this study. The G-types and P-types combinations in this study were also common in many areas worldwide (). G9 and G5 types were dominant genotypes of the sequenced samples in our study, which was consistent with the previous study in Shandong province (). Moreover, many public reports showed that positive rates of G9 and G5 types of RVA were highest in many areas of China (data not shown), so we speculated that G9 and G5 types may spread in many regions in China. However, it should be noted that there are some limitations in this study. Only nine positive samples were chosen to sequence and analyze genotypes, which could not represent the prevalent characterizations of RVA infections in the entire regions. The size of collected samples in some regions was too small to reveal the real prevalence of RVA year by year for some reasons (samples collection or samples transportations). In spite of the limitations, it was still revealed that multiple genotype combinations of RVA circulating in East China in this study. The epidemiological analysis data of RVA were several years ago in this study, which were little value for prevention the emergence of new strains. Therefore, it is necessary to maintain surveillance on the prevalent characterizations of RVA in different regions, which facilitated vaccines and other optimal strategies rational application in swine industry to prevent RVA infections.
Statements
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: GenBank, accession numbers OP454313-OP454330.
Ethics statement
The animal study was reviewed and approved by Jiangsu Academy of Agricultural Sciences Experimental Animal Ethics Committee.
Author contributions
Conceptualization: BL, XZha, and XC. Methodology: XZhu, XC, and JZ. Software, data curation, writing—original draft preparation, and visualization: XC and RT. Validation: XC, KL, and LG. Formal analysis: XC, XZhu, and RT. Investigation: XC, XZhu, LG, and KL. Resources: BL, XC, KL, and XZhu. Writing—review and editing: BL, XC, RT, and SY. Supervision: BL, XZha, XC, and RT. Project administration and funding acquisition: BL. All authors have read and agreed to the published version of the manuscript.
Funding
This work was financially supported by National Key Research and Development Program (2022YFD1800601), National Natural Science Foundation of China (31872481 and 32272996), Jiangsu Province Natural Sciences Foundation (BK20210158), Jiangsu Agricultural Science and Technology Innovation Fund [CX(22)3028], Innovation Foundation of Jiangsu Academy of Agricultural Sciences [ZX(21)1217], the Special Project of Northern Jiangsu (SZ-LYG202109), the JBGS Project of Seed Industry Revitalization in Jiangsu Province [JBGS(2021)024], and Jiangsu Funding Program for Excellent Postdoctoral Talent (2022ZB766).
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1138419/full#supplementary-material
References
1.
TateJEBurtonAHBoschi-PintoCParasharUDWorld Health Organization-Coordinated Global Rotavirus Surveillance. Global, regional, and national estimates of rotavirus mortality in children < 5 years of age, 2000–2013. Clin Infect Dis. (2016) 62 Suppl 2:S96–S105. 10.1093/cid/civ1013
2.
PatelMMPitzerVEAlonsoWJVeraDLopmanBTateJet al. Global seasonality of rotavirus disease. Pediatr Infect Dis J. (2013) 32:e134–47. 10.1097/INF.0b013e31827d3b68
3.
VlasovaANAmimoJOSaifLJ. Porcine rotaviruses: epidemiology, immune responses and control strategies. Viruses. (2017) 9:48. 10.3390/v9030048
4.
EngevikMABanksLDEngevikKAChang-GrahamALPerryJLHutchinsonDSet al. Rotavirus infection induces glycan availability to promote ileum-specific changes in the microbiome aiding rotavirus virulence. Gut Microbes. (2020) 11:1324–47. 10.1080/19490976.2020.1754714
5.
WoodeGNBridgerJHallGAJonesJMJacksonG. The isolation of reovirus-like agents (rota-viruses) from acute gastroenteritis of piglets. J Med Microbiol. (1976) 9:203–9. 10.1099/00222615-9-2-203
6.
PappHLaszloBJakabFGaneshBDe GraziaSMatthijnssensJet al. Review of group A rotavirus strains reported in swine and cattle. Vet Microbiol. (2013) 165:190–9. 10.1016/j.vetmic.2013.03.020
7.
MatthijnssensJOttoPHCiarletMDesselbergerUVan RanstMJohneR. Vp6-sequence-based cutoff values as a criterion for rotavirus species demarcation. Arch Virol. (2012) 157:1177–82. 10.1007/s00705-012-1273-3
8.
MatthijnssensJCiarletMRahmanMAttouiHBanyaiKEstesMKet al. Recommendations for the classification of group a rotaviruses using all 11 genomic RNA segments. Arch Virol. (2008) 153:1621–9. 10.1007/s00705-008-0155-1
9.
ChaseyDBridgerJCMcCraeMA. A new type of atypical rotavirus in pigs. Arch Virol. (1986) 89:235–43. 10.1007/BF01309892
10.
WakudaMIdeTSasakiJKomotoSIshiiJSanekataTet al. Porcine rotavirus closely related to novel group of human rotaviruses. Emerg Infect Dis. (2011) 17:1491–3. 10.3201/eid1708.101466
11.
MarthalerDHomwongNRossowKCulhaneMGoyalSCollinsJet al. Rapid detection and high occurrence of porcine rotavirus A, B, and C by Rt-QPCR in diagnostic samples. J Virol Methods. (2014) 209:30–4. 10.1016/j.jviromet.2014.08.018
12.
SarkarSEsonaMDGautamRCastroCJNgTFFHaqueWet al. Outbreak of diarrhoea in piglets caused by novel rotavirus genotype G4p[49] in North-Western District of Bangladesh, February 2014. Transbound Emerg Dis. (2020) 67:442–9. 10.1111/tbed.13343
13.
AmimoJOVlasovaANSaifLJ. Detection and genetic diversity of porcine group a rotaviruses in historic (2004) and recent (2011 and 2012) swine fecal samples in Ohio: predominance of the G9p[13] genotype in nursing piglets. J Clin Microbiol. (2013) 51:1142–51. 10.1128/JCM.03193-12
14.
PhamHACarrique-MasJJNguyenVCNgoTHNguyetLADoTDet al. The prevalence and genetic diversity of group a rotaviruses on pig farms in the Mekong Delta Region of Vietnam. Vet Microbiol. (2014) 170:258–65. 10.1016/j.vetmic.2014.02.030
15.
CollinsPJMartellaVSleatorRDFanningSO'SheaH. Detection and characterisation of group a rotavirus in asymptomatic piglets in Southern Ireland. Arch Virol. (2010) 155:1247–59. 10.1007/s00705-010-0713-1
16.
OkitsuSKhamrinPThongprachumAKongkaewAManeekarnNMizuguchiMet al. Whole-genomic analysis of G3p[23], G9p[23] and G3p[13] rotavirus strains isolated from piglets with diarrhea in Thailand, 2006-2008. Infect Genet Evol. (2013) 18:74–86. 10.1016/j.meegid.2013.05.005
17.
OkitsuSKhamrinPThongprachumAManeekarnNMizuguchiMUshijimaH. Predominance of porcine P[23] genotype rotaviruses in piglets with diarrhea in Northern Thailand. J Clin Microbiol. (2011) 49:442–5. 10.1128/JCM.02263-10
18.
MiyazakiAKugaKSuzukiTKohmotoMKatsudaKTsunemitsuH. Genetic diversity of group A rotaviruses associated with repeated outbreaks of diarrhea in a farrow-to-finish farm: identification of a porcine rotavirus strain bearing a novel Vp7 Genotype, G26. Vet Res. (2011) 42:112. 10.1186/1297-9716-42-112
19.
CollinsPJMartellaVBuonavogliaCO'SheaH. Identification of a G2-like porcine rotavirus bearing a novel Vp4 type, P[32]. Vet Res. (2010) 41:73. 10.1051/vetres/2010045
20.
MoutelikovaRDufkovaLKamlerJDrimajJPlhalRProdelalovaJ. Epidemiological survey of enteric viruses in wild boars in the Czech Republic: first evidence of close relationship between wild boar and human rotavirus a strains. Vet Microbiol. (2016) 193:28–35. 10.1016/j.vetmic.2016.08.003
21.
XueRTianYZhangYZhangMLiZChenSet al. Diversity of group a rotavirus of porcine rotavirus in Shandong Province China. Acta Virol. (2018) 62:229–34. 10.4149/av_2018_216
22.
SaifLJ. Comparative pathogenesis of enteric viral infections of swine. Adv Exp Med Biol. (1999) 473:47–59. 10.1007/978-1-4615-4143-1_4
23.
SvensmarkBNielsenKDalsgaardKWillebergP. Epidemiological studies of piglet diarrhoea in intensively managed danish sow herds. Iii Rotavirus Infection Acta Vet Scand. (1989) 30:63–70. 10.1186/BF03548069
24.
WenFYangJLiAGongZYangLChengQet al. Genetic characterization and phylogenetic analysis of porcine epidemic diarrhea virus in Guangdong, China, between 2018 and 2019. PLoS ONE. (2021) 16:e0253622. 10.1371/journal.pone.0253622
25.
MiyazakiAKugaKSuzukiTKohmotoMKatsudaKTsunemitsuH. Annual changes in predominant genotypes of rotavirus a detected in the feces of pigs in various developmental stages raised on a conventional farm. Vet Microbiol. (2013) 163:162–6. 10.1016/j.vetmic.2012.11.044
26.
ZhangFLuoSGuJLiZLiKYuanWet al. Prevalence and phylogenetic analysis of porcine diarrhea associated viruses in Southern China from 2012 to 2018. BMC Vet Res. (2019) 15:470. 10.1186/s12917-019-2212-2
27.
DingGFuYLiBChenJWangJYinBet al. Development of a multiplex RT-PCR for the detection of major diarrhoeal viruses in pig herds in China. Transbound Emerg Dis. (2020) 67:678–85. 10.1111/tbed.13385
28.
ShiHChenJLiHSunDWangCFengL. Molecular characterization of a rare G9p[23] porcine rotavirus isolate from China. Arch Virol. (2012) 157:1897–903. 10.1007/s00705-012-1363-2
29.
YiJLiuC. Detection of a Porcine Rotavirus Strain with Vp4, Vp7 and Nsp4 Genes of Different Animal Origins. Arch Virol. (2011) 156:2045–52. 10.1007/s00705-011-1064-2
30.
YanNYueHWangYZhangBTangC. Genomic analysis reveals G3p[13] porcine rotavirus a interspecific transmission to human from pigs in a swine farm with diarrhoea outbreak. J Gen Virol. (2021) 102(2). 10.1099/jgv.0.001532
31.
WangZLvCXuXLiXYaoYGaoXet al. The dynamics of a Chinese Porcine G9p[23] rotavirus production in Ma-104 cells and intestines of 3-day-old piglets. J Vet Med Sci. (2018) 80:790–7. 10.1292/jvms.17-0657
32.
ZhangHZhangZWangYWangXXiaMWuH. Isolation, molecular characterization and evaluation of the pathogenicity of a porcine rotavirus isolated from Jiangsu Province, China. Arch Virol. (2015) 160:1333–8. 10.1007/s00705-015-2347-9
33.
ZhouJHuangSFanBNiuBGuoRGuJet al. ITRAQ-based proteome analysis of porcine group a rotavirus-infected porcine IPEC-J2 intestinal epithelial cells. J Proteomics. (2021) 248:104354. 10.1016/j.jprot.2021.104354
Summary
Keywords
porcine rotavirus, prevalence, East China, VP7, VP4
Citation
Tao R, Chang X, Zhou J, Zhu X, Yang S, Li K, Gu L, Zhang X and Li B (2023) Molecular epidemiological investigation of group A porcine rotavirus in East China. Front. Vet. Sci. 10:1138419. doi: 10.3389/fvets.2023.1138419
Received
05 January 2023
Accepted
06 March 2023
Published
21 March 2023
Volume
10 - 2023
Edited by
Shao-Lun Zhai, Guangdong Academy of Agricultural Sciences, China
Reviewed by
Ahsan Naveed, University of Kentucky, United States; Nader Sobhy, Zagazig University, Egypt
Updates
Copyright
© 2023 Tao, Chang, Zhou, Zhu, Yang, Li, Gu, Zhang and Li.
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.
*Correspondence: Bin Li libinana@126.com
†These authors have contributed equally to this work and share first authorship
This article was submitted to Veterinary Infectious Diseases, a section of the journal Frontiers in Veterinary Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.