- Department of Veterinary Infectious Diseases and Avian Diseases, Center for Poultry Diseases Control, College of Veterinary Medicine, Jeonbuk National University, Iksan, South Korea
Avian reoviruses (ARVs) cause severe arthritis, tenosynovitis, pericarditis, and depressed growth in chickens, and these conditions have become increasingly frequent in recent years. Studies on the role of wild birds in the epidemiology of ARVs are insufficient. This study provides information about currently circulating ARVs in wild birds by gene detection using diagnostic RT-PCR, virus isolation, and genomic characterization. In this study, we isolated and identified 10 ARV isolates from 7,390 wild birds' fecal samples, including migratory bird species (bean goose, Eurasian teal, Indian spot-billed duck, and mallard duck) from 2015 to 2019 in South Korea. On comparing the amino acid sequences of the σC-encoding gene, most isolates, except A18-13, shared higher sequence similarity with the commercial vaccine isolate S1133 and Chinese isolates. However, the A18-13 isolate is similar to live attenuated vaccine av-S1133 and vaccine break isolates (SD09-1, LN09-1, and GX110116). For the p10- and p17-encoding genes, all isolates have identical fusion associated small transmembrane (FAST) protein and nuclear localization signal (SNL) motif to chicken-origin ARVs. Phylogenetic analysis of the amino acid sequences of the σC-encoding gene revealed that all isolates were belonged to genotypic cluster I. For the p10- and p17-encoding genes, the nucleotide sequences of all isolates indicated close relationship with commercial vaccine isolate S1133 and Chinese isolates. For the σNS-encoding gene, the nucleotide sequences of all isolates indicated close relationship with the Californian chicken-origin isolate K1600657 and belonged to chicken-origin ARV cluster. Our data indicates that wild birds ARVs were derived from the chicken farms. This finding suggests that wild birds serve as natural carriers of such viruses for domestic poultry.
Introduction
Avian reoviruses (ARVs), classified under the Orthoreovirus genus of Reoviridae family, have a non-enveloped icosahedral double capsid with 10 double-stranded RNA (dsRNA) genome segments which are further divided into three size classes based on their electrophoretic mobility: large (L1, L2, L3), medium (M1, M2, M3), and small (S1, S2, S3, S4) (1). ARV genome encodes eight structural (λA, λB, λC, μA, μB, σA, σB, and σC) and four non-structural viral proteins (μNS, σNS, p10, and p17) (2). Among the S-class segments of ARV, the segment S1 contains three open reading frames (ORFs) that are translated into p10, p17, and σC proteins (2). The p10 protein induces cell fusion (3), while the p17 protein induces cellular protein translation shut-off and cell cycle arrest (4, 5). The σC protein encoded by the third ORF of the segment S1 is the cell attachment protein and elicits the production of ARV-specific neutralizing antibodies (6, 7). The σC sequence is used as a genetic marker to characterize and classify ARV isolates into different six or seven genotypic clusters (8, 9). The σNS protein encoded by the segment S4, has been reported for its single-stranded RNA (ssRNA) binding activity, and it is divided into diverse lineages that conserved by host origin of ARVs (10, 11); as well it has been used for diagnostic analysis (10, 12).
ARVs are widespread in nature and affect various commercial and wild avian species. ARVs have been isolated from poultry, such as chickens, ducks, turkeys, ostriches, and wild birds (13–17). Major damages from ARV infections are observed in young chickens, and ARV infections in broiler chickens are associated with up to 10% mortality and up to 40% morbidity. ARV infections in broilers have been associated with various clinical diseases, such as viral arthritis/tenosynovitis, malabsorption syndrome (MAS), runting–stunting syndrome (RSS), respiratory diseases, hepatitis, myocarditis, neurological signs (incoordination, tremors, and twisted necks), and immunosuppression (18, 19). ARV infections are responsible for significant economic losses in the poultry industry (20). There have been ARV infection outbreaks in the USA, Canada, Brazil, Europe, Israel, China, and Korea (8, 21–26).
ARVs have also been identified in wild birds, such as crows, magpies, partridges, black-capped chickadees, brown-eared bulbuls, psittacine bird species, and mallard ducks (27–35). ARV infections in wild birds have been associated with mortality and/or various clinical diseases, such as hepatitis, enteritidis, and neurological signs (29, 31, 33, 34, 36, 37). A previous study on the seroprevalence of ARVs in wild birds reported 25% seroprevalence in bean goose (Anser fabalis) and 34% seroprevalence in white-fronted goose (Anser albifrons frontalis) in Germany in 1998 (38). In addition, antigenic prevalence in wild birds representing 32 species was reportedly ~30% in Poland from 2014 to 2016 (12). ARVs have been isolated from hooded crows in Finland, American crows in the USA, brown-eared bulbuls in Japan, psittacine bird species in Germany, and mallard ducks in China (27, 29, 30, 34, 35). Various hypotheses have been proposed about the role of wild birds in the dissemination and maintenance of ARVs and pathogenic potential of ARVs in wild birds that usually commingle around poultry farms. In previous studies, ARVs isolated from wild birds were somewhat distant from the those from chicken farms (29, 30, 33–35). However, recently research showed that the σA-encoding gene isolated from a healthy ostrich at a domestic farm in Japan showed great similarity to the chicken-origin ARVs, and the σC-encoding gene isolated from magpies was found to be genetically similar to the chicken-origin ARVs (13, 31). This genetic association between isolates from wild birds and poultry can be explained by the potential of wild birds to act as a carrier for the transmission of the ARVs (31). However, studies on the role of wild birds in the epidemiology of ARVs are insufficient.
This study provides information about currently circulating ARVs in wild birds by gene detection using diagnostic RT-PCR, virus isolation, and genomic characterization of S1 and S4 gene. Furthermore, it investigates the possibility of wild birds being potential sources of infection in poultry.
Materials and Methods
Sample
A total of 7,390 fresh fecal samples were obtained from major migratory bird habitats, including near commercial/domestic chicken premises in South Korea by an avian influenza virus national monitoring program (2015–2019). All live migratory birds were trapped and collected feces by the South Korea Animal and Plant Quarantine Agency (QIA) and the Ministry of Agriculture, Food and Rural Affairs (MAFRA) (39). Ten fecal samples obtained from pure or separately standing flocks of the same species were pooled into one test sample. Finally, the 739 resultant test samples were 10-fold diluted in sterile phosphate-buffered saline [PBS, pH7.4; supplemented with 100X antibiotic–antimycotic (Gibco, New York, USA)], thoroughly mixed by vortexing, and centrifugated at 600 x g for 10 min at 4°C. The supernatant was subsequently filtered (0.45-μm pore size; Minisart® NML, Sartorious, Germany). The filtered supernatant was then conserved in aliquots at −70°C for virus isolation and viral RNA extraction.
RNA Extraction and RT-PCR
Viral RNA was extracted from the clarified fecal samples using the MagMAXTM - 96 AI/ND Viral RNA isolation kit (Thermo Fisher Scientific, Vilnius, Lithuania) with KingFisher Duo Prime Purification system (Thermo Fisher Scientific, Waltham, MA) following the manufacturer's protocol. Viral cDNA was generated from RNA samples using GoScriptTM reverse transcriptase (Promega, Madison, WI USA) with random primers (9-mers; TaKaRa Bio. Inc., Otsu, Shiga, Japan). In the reverse transcription (RT) reaction, 8 μL of extracted RNA and 2 μL of dimethyl sulfoxide (DMSO, Tedia, USA) were heated at 100°C for 5 min and then placed in an ice bath for 5 min. Then, the following were added to this reaction mixture: 8 μL of GoScriptTM 5X RT reaction buffer (Promega, Madison, WI USA), 10 μL of 2.0 mM of each dNTP (SolGent, Daejeon, Korea), 4 μL of MgCl2 (Promega, Charbonnie‘re, France), 1 μL of 20 units Recombinant RNasin® Ribonuclease Inhibitor (Promega, Madison, WI, USA), 1 μL of GoScriptTM reverse transcriptase, 1 μL of 50 pmol random primer, and 4 μL of diethylpyrocarbonate-treated water (DEPC water; Biosesang, Seoul, Korea); a final volume of 39 μL was obtained. The RT reaction mixture was incubated in this sequence: 25°C for 5 min, 42°C for 60 min, and 70°C for 15 min to inactivate the enzyme. The cDNA amplification was performed in 50 μL volume, containing 10X Taq buffer (Solgent, Daejeon, Korea), 5U Taq DNA polymerase (Solgent, Daejeon, Korea), 10 pmol of each primer (Table 1). For the detection of ARVs in wild bird feces, partial S1 gene was amplified. Thermal cycling protocols were as follows: Initial denaturation at 94°C for 5 min, 35 cycles (denaturation at 94°C for 1 min, annealing at 50°C for 1 min, extension at 72°C for 1 min) and one final extension at 72°C for 10 min (40). For the sequence analysis, viral RNA was extracted from virus stocks using the Viral Gene-SpinTM Viral DNA/RNA Extraction kit (iNtRON, Daejeon, Korea). Viral cDNA was obtained from RNA samples by the above-mentioned method. For the amplification of the full S1 gene, thermal cycling protocols were as follows: Initial denaturation at 94°C for 5 min, 40 cycles (denaturation at 94°C for 1 min, annealing at 60°C for 1 min, extension at 72°C for 100 s) and one final extension at 72°C for 10 min (2). For the amplification of the full S4 gene, thermal cycling protocols were as follows: Initial denaturation at 94°C for 5 min, 35 cycles (denaturation at 94°C for 30 s, annealing at 58°C for 30 s, extension at 72°C for 90 s) and one final extension at 72°C for 10 min (41).
DNA Barcoding for Species Identification of Fecal Samples
The host species of the PCR-positive fecal samples were identified using the DNA barcoding technique as previously described (42). Briefly, DNA was extracted from fresh fecal samples using an Accuprep® Stool DNA Extraction kit (Bioneer, Daejeon, Korea) according to the manufacturer's protocol. The primers Aves-F (5′-GCATGAGCAGGAATAGTTGG-3′) and Aves-R (5′-AAGATGTAGACTTCTGGGTG-3′) were used to amplify the mitochondrial cytochrome oxidase gene subunit I present in host feces (43). PCR products were sequenced and identified using the information available at the Barcode of Life Data Systems website (Biodiversity Institute of Ontario, University of Guelph, Guelph, Ontario, Canada; http://www.barcodinglife.org/views/login.php) (44).
Virus Isolation
The PCR-positive samples were inoculated into the yolk sac of 6-day-old specific pathogen-free (SPF) embryos (SPAFAS Poultry Company, Jinan, China). Embryonated eggs were candled daily for 5 days. Chorioallantoic membranes (CAMs) of the infected embryos were collected by three cycles of freezing and thawing; then, 1.5 mL PBS was added, followed by centrifugation at 6,000 x g for 10 min. For cell culture passage, chicken embryo liver (CEL) cells were prepared from 14-day-old SPF chicken embryos as per the standard protocol and then dispensed into a 6-well-cell culture plate (SPL life sciences, Pocheon, Korea). The medium for CEL cell culture was Eagle's minimum essential medium supplemented with 8% fetal bovine serum (FBS) and 1% addition of 100X antibiotic–antimycotic. When the cell monolayers were ~80% confluent, the medium was aspirated. Cell monolayers were infected with 0.2 mL of 10-fold diluted CAM fluids from the first passage of chicken embryos and incubated at 37°C for 60 min. Then, a maintenance medium containing 4% FBS was added. The cultures were incubated at 37°C under 5% CO2 and were observed daily under a microscope to check for a cytopathic effect (CPE). Once 70–80% CPE was developed, the cultures were subjected to three cycles of freezing and thawing and then clarified with low centrifugation at 600 x g for 20 min (41).
Cloning and Sequencing
PCR products of full-length S1 and S4 genes were cloned and sequenced. The PCR products of the expected lengths were purified with the QIAquick gel extraction kit (Qiagen, Chatsworth, CA) and then cloned into the pGEM-T easy vector (Promega, Madison, WI, USA) according to the manufacturer's instructions, and nucleotide sequences were determined using an ABI 3730XL DNA Analyzer (Applied Biosystems, Foster City, CA) by SolGent (Daegeon, Korea) (45).
Sequence Analysis and Phylogenetic Analysis
Based on the nucleotide and deduced amino acid sequences of the σC- and σNS-encoding genes, phylogenetic analysis was carried out. Sequence alignments and pairwise sequence comparisons were performed using the GENETYX software (Genetyx Corp., Tokyo, Japan) (46). Phylogenetic analysis was performed using the maximum likelihood method of MEGA-X software package (version 10.0.5) with bootstrap values calculated from 1,000 replicates (47). Recombination detection was accomplished by using Bootscan analysis within the Simplot program version 3.5.1, using the neighbor-joining method, with a Kimura 2-parameter applied and 100 replicates (48).
Results
Virus Isolation and Propagation
A total of 14 PCR-positive samples from 739 samples (1.9%) were identified. Ten ARVs were isolated by using CEL cell from the 14 PCR-positive samples, and all isolates were confirmed by diagnostic RT-PCR as described above in Table 1. Five ARVs were detected from bean goose (Anser fabalis), one from mallard duck (Anas platyrhtchos), one from Eurasian teal (Anas crecca), one from oriental turtle dove (Streptopelia orientalis), and one from Indian spot-billed duck (Anas poecilorhyncha); the source could not be identified for one ARV (Table 2). CPE was observed in CEL cell infected with the 10 isolates. The CPE of isolates manifested as detachment of the monolayer or syncytium formation (Figure 1).
Table 2. ARV obtained in the study, origin of sample, year of isolation, host species, sample type, and GenBank accession number.
Figure 1. Cytopathic effect (CPE) in ARV-infected CEL cells. (A) S1133/Chicken/USA/1971 infected CEL cells forming a syncytial cluster at 1 day post infection. (B) Representative wild bird isolate (A18-13/Wild bird/Korea/2018) infected CEL cells forming a syncytial cluster at 1 day post infection. (C) Non-infected CEL cell control.
Sequence Analysis
The full-length S1 and S4 genes of all ARV isolates were sequenced, and the obtained sequences were submitted to the GenBank database under the accession numbers showed in Table 2. The σC- and σNS-encoding genes of all isolates were analyzed and compared with reference isolates (Supplementary Table 1). The alignment of nucleotide sequences of the σC-encoding gene of 10 isolates showed high similarity with commercial vaccine isolates (av-S1133, S1133, 2408, and 1733) (98.4–99.9%) and Chinese chicken-origin isolates (GuangxiR1, SD09-1, LN09-1, and GX110116) (98.6–99.6%) in nucleotide sequence. Conversely, the 10 isolates showed low similarity with previously identified Korean chicken-origin isolates (SNU0044, SNU0046, K738-14, and iREO0309) (50.2–74.5%) and Tvärminne avian virus (TVAV)-like wild bird isolates (Pycno-1, Corvus corone cornix 2002, and Tvärminne avian virus) (48.5–53.9%) (Supplementary Table 2). An alignment of nucleotide sequences of the σNS-encoding gene of the 10 isolates showed high similarity with the Californian chicken-origin isolate K1600657 (92.7%−92.9%) and low similarity with the TVAV-like wild bird isolates (Pycno-1, SD-12, Tvärminne avian virus, 71-03, Chickadee, and SRK) (58.2%−77.7%) (Supplementary Table 3).
We compared the amino acid sequences of the 10 isolates with those of 15 reference isolates. The results were found to be consistent with the results of pairwise comparisons of the σC-encoding gene (Supplementary Table 4). Compared at six conservative amino acid mutations (24, 71, 106, 113, 134, and 135), most newly isolated viruses, except A18-13 isolate, shared higher sequence similarity with commercial vaccine isolate S1133 and Chinese chicken-origin isolates (GuangxiR1 and GuangxiR2). Whereas 10 isolates shared lower similarity with Korean chicken-origin isolates (SNU0044, SNU0046, iREO0309, and K738-14). Nucleotide changes in the σC ORF unique to A18-13 isolate and the live attenuated vaccine av-S1133 and vaccine break isolates (SD09-1, LN09-1, and GX110116) were identified. A18-13 isolate, live attenuated vaccine, and vaccine break isolates carry four nonsynonymous nucleotide substitutions (A71, G317, T338, and C405), whereas commercial vaccine isolate S1133 have nucleotide substitutions at C71, C317, C338, and A405 in the σC-encoding gene (Supplementary Table 4). These four non-synonymous nucleotide substitutions lead to deduced amino acid substitutions (T24N, T106R, T113I, and V135I) in the σC-encoding gene (Table 3). Nucleotide changes in the σNS ORF unique to all isolates and the Californian chicken-origin isolate K1600657 were identified. All isolates and K1600657 carry 10 non-synonymous nucleotide substitutions (G215, T361, A455, A472, G478, G484, T652, A760, G817, and C1049), whereas commercial vaccine isolates, including live attenuated vaccine have nucleotide substitutions at A215, C361, G455, G472, A478, C484, G652, C760, A817, and T1049 in the σNS-encoding gene (Supplementary Table 4). These 10 non-synonymous nucleotide substitutions lead to deduced amino acid substitutions (D72G, P121S, R152H, A158T, S160D, P162A, A218S, L254T, T273V, and M350T) in the σNS-encoding gene (Table 3). For the p10 protein, all isolates were not observed for deduced amino acid substitutions in putative transmembrane domains. For the p17 protein, all isolates were not observed for deduced amino acid substitutions in nuclear localization signal (NLS) motif (Supplementary Table 5).
Table 3. Deduced amino acid substitutions in the σC- and σNS-encoding genes of ARV isolates from wild bird.
Phylogenetic Analysis
Phylogenetic analysis of the σC- and σNS-encoding genes was demonstrated. For the σC-encoding gene, nucleotide sequences of the 10 isolates showed that they were closely related to each other and clustered with Chinese chicken-origin isolates and commercial vaccine isolates. Conversely, previously identified TVAV-like wild bird isolates (Pycno-1, SD-12, Tvärminne avian virus, and Corvus corone cornix 2002) were distinct from all isolates in this study (Figure 2). To classify the genotypes of our isolates, we constructed phylogenetic trees using deduced amino acid sequences of the σC-encoding gene. This analysis revealed that all isolates were in the same cluster (genotypic cluster I) as the commercial vaccine isolates (Figure 3). For the σNS-encoding gene, nucleotide sequences of the 10 isolates showed that they are grouped together in one branch and clustered with the Californian chicken-origin isolate K1600657. Conversely, Chinese chicken-origin isolates, TVAV-like wild bird isolates (Pycno-1, SD-12, Tvärminne avian virus, 71-03, Chickadee, and SRK) and commercial vaccine isolates were distinct from all isolates in this study (Figure 4). For the p10- and p17-encoding gene, nucleotide sequences of the 10 isolates showed that they are clustered with Chinese chicken-origin isolates and commercial vaccine isolates. Conversely, previously identified TVAV-like wild bird isolates (Pycno-1 and Tvärminne avian virus) were distinct from all isolates in this study (Supplementary Figures 1, 2).
Figure 2. Maximum likelihood phylogenetic tree for the σC-encoding gene based on nucleotide sequences (981 nt). Maximum likelihood phylogenetic analyses were conducted using MEGA-X software with the Kimura 2-parameter model and 1,000 bootstrap replicates. The 43 reference sequences were obtained from GenBank. The black circle (•) indicates our isolates, while the white circle (◦) indicates the vaccine isolates. Additionally, the black diamond (♦) indicates field isolates in Korea, and the black triangle (▴) indicates previously isolated wild bird isolates. Each sequence on the tree is identified by the isolate name, host, country of origin, year of isolation, and GenBank accession number.
Figure 3. Phylogenetic tree of ARV isolates based on the deduced amino acid sequences of the σC-encoding gene. Maximum likelihood phylogenetic analyses were conducted using MEGA-X software with the Jones–Taylor–Thornton (JTT) model and 1,000 bootstrap replicates. The tree shows genetic relationship between the σC protein sequences (326 amino acids) of our 10 isolates and 48 reference isolates that were isolated from around the world. The virus isolates are clustered into six genotypic clusters. The black circle (•) indicates our isolates, and the white circle (◦) indicates the vaccine isolates. Additionally, the black diamond (♦) indicates field isolates in Korea, and the black triangle (▴) indicates previously isolated wild bird isolates. Each sequence on the tree is identified by the isolate name, host, country of origin, year of isolation, and GenBank accession number.
Figure 4. Maximum likelihood phylogenetic tree for the σNS-encoding gene based on nucleotide sequences (1,104 nt). Maximum likelihood phylogenetic analyses were conducted using MEGA-X software with the Kimura 2-parameter model and 1,000 bootstrap replicates. The 52 reference sequences were obtained from GenBank. The black circle (•) indicates our isolates, and the white circle (◦) indicates the vaccine isolates. Additionally, the black diamond (♦) indicates field isolates in Korea, and the black triangle (▴) indicates previously isolated wild bird isolates. Each sequence on the tree is identified by the isolate name, host, country of origin, year of isolation, and GenBank accession number.
In Recombination analysis, representative 5 isolates (A15-113, A15-157, A18-13, A18-19, and A19-106) were selected by five independent recombination patterns in the σC-encoding gene of our 10 isolates. A15-113 has a similar pattern to A15-108, and A15-157 has a similar pattern to A15-19, A15-48, and A18-205. In addition, A18-13 has a similar pattern to A15-71, and A18-19 and A19-106 showed independent patterns. Whereas all isolates have identical recombination patterns in σNS-encoding gene. All reference sequences in Supplementary Table 1 were analyzed. However, for the σC-encoding gene, meaningful recombination with 3 isolates (GX110116, GuangxiR1, and K1600657) were detected. And for the σNS-encoding gene, meaningful recombination with 3 isolates (GX110116, AVS-B, and 01224A-14) were detected. According to a Bootscan analysis, it appears that the σC-encoding gene of A18-13 were recombined with Chinese vaccine break isolates GX110116 and Californian chicken-origin isolate K1600657. However, 4 isolates (A15-113, A15-157, A18-19, and A19-106) did not show clear evidence for recombination with these isolates (GX110116 and K1600657). For the σNS-encoding gene, there is evidence that 5 isolates may have recombined with Chinese vaccine break isolate GX110116 and American chicken-origin isolate AVS-B (Figure 5).
Figure 5. Bootscan analysis of wild bird ARV for detecting recombination was performed using Simplot program version 3.5.1, the neighbor-joining method with a Kimura 2-parameter applied and 100 replicates. (A) σC-encoding gene, (B) σNS-encoding gene.
Discussion
It has been hypothesized that wild birds have a role in the transmission and maintenance of ARVs and pathogenic potential of ARVs in wild birds that usually commingle around poultry farms. However, lack of evidence to support these hypotheses. Molecular surveillance is crucial to strategize control and prevention of endemic diseases. Contrary to the previous study results, recently detected reovirus isolates in ostrich (Stuthio camelus) and free-living magpie were genetically related to chicken-origin ARVs (13, 31). Despite of the need for a reliable and relevant source of genetic information of wild birds, available genetic information is rare. In this study, we conducted surveillance and molecular characterization of ARVs in wild birds, and ARVs were detected in 1.9% (14/739) of the wild bird feces. This result suggests a much lower prevalence of ARVs than reported previously in some studies in Poland (11.5%; 9/78 dead wild birds; intestine) (12) and the USA (25.0%; 4/20 dead American woodcock; intestine) (49). Both studies reported that ARVs in wild bird's intestine had high prevalence than ARVs in wild bird feces. These high prevalence in intestines are probably due to ARVs have enteric tropism (50) and low prevalence in feces may have influenced by several factors, such as environmental factors or feces can act as an anti-viral factor (51). A similar result to our study was reported in Brazil that conducted surveillance of ARVs in fecal samples from broiler chickens (1.9%; 7/378) (52). The differences in prevalence could be attributed to the sample type (feces vs. tissue) or regional differences.
We classified reovirus isolates in this study based on the σC protein, all isolates belonged to genotypic cluster I, which is the predominant cluster in chicken-origin ARVs and includes commercial vaccine isolates (Figure 3) (8, 53). Similar results in wild birds were described in the UK and ARV isolated from free-living magpie (Pica pica) belonged to genotypic cluster I (31). Conversely, in previous studies, wild bird-origin ARVs belonged to the TVAV-like cluster (27–30, 34, 35). Moreover, the phylogenetic analysis of the p10- and p17-encoding genes revealed that all isolates were closely clustered with Chinese chicken-origin isolates. For the σNS-encoding gene, our isolates were more closely clustered with the Californian chicken-origin isolate K1600657 that belonged to chicken-origin ARV cluster (Figure 4). These findings suggest that chicken-origin ARV can infect wild birds and can thus be a potential source of ARV infections that can spillover to chickens and vice versa.
Throughout this study, the 10 ARVs of the σC protein showed high degree of antigenic homogeneity to commercial vaccine isolates (inactivated vaccine) and Chinese chicken-origin isolates. However, these isolates were distinct from Korean chicken-origin isolates. In addition, analysis of the σNS protein showed high similarity to Californian chicken-origin isolate K1600657. Our results indicated that the deduced amino acid substitution patterns of all isolates are identical to those of chicken-origin ARVs in the σC and σNS proteins. Particularly, the A18-13 isolate had identical patterns with avirulent ARV (av-S1133) and virulent ARV (GX110116) in terms of the σC protein, which contributed to ARV virulence (54, 55). For the σNS protein, deduced amino acid substitutions at positions 72, 121, 152, 158, 160, 162, 218, 254, 273, and 350 were observed for all isolates with respect to S1133 (Table 3). These substitutions were also observed for the Californian chicken-origin isolate K1600657 that was associated with severe clinical signs (53). These findings suggested that pathotyping based on the σC protein has limitations. More analyses of polymorphisms and deduced amino acid substitutions of other segments are required for pathotyping of ARVs. For the p10 protein, all isolates had a putative transmembrane domains, which are members of the fusion associated small transmembrane (FAST) protein family that conserved in chicken-origin ARVS (56). For the p17 protein, all isolates had identical nuclear localization signal (NLS) motif with chicken-origin isolates, which is conserved in chicken-origin ARVS (Supplementary Table 5) (57). These finding suggested that 10 isolates have identical FAST protein and NLS to chicken-origin ARVs.
According to recombination analysis, the σC-encoding gene of A18-13 isolates might have recombined with Chinese vaccine break isolate GX110116 at ~300–500 nucleotide sequence positions that include four non-synonymous nucleotide substitutions (G317, T338, A403, and C405) and Californian chicken-origin isolate K1600657 at approximate 700–900 nucleotide sequence positions with no non-synonymous nucleotide substitutions. Whereas other four isolates showed independent patterns in recombination analysis, although these isolates did not show clear evidence of recombination. The σNS-encoding gene of five isolates (A15-113, A15-157, A18-13, A18-19, and A19-106) might have recombined with American chicken-origin isolate AVS-B at ~350–470 nucleotide sequence positions that include three non-synonymous nucleotide substitutions (T361, A455, and A472) and Chinese vaccine break isolate GX110116 at 100–350 and 470–1,000 nucleotide sequence positions that include seven non-synonymous nucleotide substitutions (G215, G478, G484, T652, A760, G817, and C1049). In addition, all isolates showed identical recombination patterns in the σNS-encoding gene (Figure 5 and Supplementary Table 4). These findings suggest that recombination patterns are seems to be associated with the characteristics of the σC and σNS protein. In terms of the σC protein, which is the most variable region in the viral genome (58). For the σNS protein, which is relatively conserved among other genome segments (59). Moreover, our 10 isolates genome reflects genomic reassortment events and intra-segmental recombination between Chinese vaccine break isolate and American chicken-origin isolate in the σNS-encoding gene (Figure 5). These finding indicated that our isolates are intracontinental recombined ARVs and can be potential risk in commercial chicken farms.
Wild bird ARVs have been isolated in sedentary bird species, including magpie, crow, partridge, black-capped chickadee, and psittacine birds (27, 29, 31–33). A previous study reported the isolation of wild bird ARV from a migratory bird species (mallard duck) (35). Eight isolates identified in this study (A15-48, A15-71, A15-108, A15-113, A18-13, A18-19, A18-205, and A19-106) were isolated from migratory bird species (bean goose, Eurasian teal, Indian spot-billed duck, and mallard duck) (Table 2). In the flyways of these birds, South Korea is an important wintering site for wild migratory birds, such as waterfowl which flies across the East Asia–Australian flyway, with Korea and China located along this flyway (60). The A15-157 isolate was from an oriental turtle dove. Notably, the oriental turtle dove is a sedentary bird species commonly found around chicken farms in Korea. These sedentary bird species may play a role in ARV transmission between poultry and wild birds.
Conclusion
We describe genetic characterization of 10 ARVs isolated from wild birds' feces and showed that all isolates were closely related to chicken-origin ARVs. It is possible that wild birds may play a potential role in the epidemiology of chicken-origin ARVs as a carrier. For further study to evaluate risk factors, it is essential to conduct whole-genome sequencing and investigate the potential pathogenicity of ARVs isolated from wild birds for domestic poultry.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
S-YC and MK: contributed to conception and design of experiments. H-KJ, S-WK, Y-RC, KS, BW, J-YP, and J-FZ: contributed to acquisition, analysis, and interpretation of data. S-WK, Y-RC, S-YC, and MK: drafted and/or revised the article. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by Korea Institute of Planning and Evaluation for Technology in Food, Agriculture and Forestry (IPET) through Agriculture, Food and Rural Affairs Convergence Technologies Program for Educating Creative Global Leader (716002-7, 320005-4) funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
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.2022.794934/full#supplementary-material
Supplementary Figure 1. Maximum likely phylogenetic tree for p10-encoding gene based on nucleotide sequences (297 nt). Maximum likelihood phylogenetic analyses were conducted using MEGA X software with the Kimura 2-parameter model and 1,000 bootstrap replicates. The 31 reference sequences were obtained from GenBank. The black circle (•) indicates our isolates while the white circle (◦) indicates the vaccine isolates. Additionally, the black diamond (♦) indicates field isolates in Korea while the black triangle (▴) indicates previously isolated wild bird isolates. Each sequence on the tree is identified by isolated name, host and country of origin, year of isolation, GenBank accession number.
Supplementary Figure 2. Maximum likely phylogenetic tree for p17-encoding gene based on nucleotide sequences (441 nt). Maximum likelihood phylogenetic analyses were conducted using MEGA X software with the Kimura 2-parameter model and 1,000 bootstrap replicates. The 31 reference sequences were obtained from GenBank. The black circle (•) indicates our isolates while the white circle (◦) indicates the vaccine isolates. Additionally, the black diamond (♦) indicates field isolates in Korea while the black triangle (▴) indicates previously isolated wild bird isolates. Each sequence on the tree is identified by isolated name, host and country of origin, year of isolation, GenBank accession number.
Supplementary Table 1. Information of new isolated ARV and reference isolates used in this study.
Supplementary Table 2. Sequence alignment of the σC-encoding gene (981 nt).
Supplementary Table 3. Sequence alignment of the σNS-encoding gene (1,104 nt).
Supplementary Table 4. Single nucleotide polymorphisms in the σC- and σNS-encoding genes of ARV isolates from wild bird.
Supplementary Table 5. Deduced amino acid substitutions in the putative transmembrane domains of the p10 and nuclear localization signal motif in the p17.
References
1. Spandidos Da, Graham A. Physical and chemical characterization of an avian reovirus. J Virol. (1976) 19:968–76. doi: 10.1128/jvi.19.3.968-976.1976
2. Bodelón G, Labrada L, Martinez-Costas J, Benavente J. The avian reovirus genome segment S1 is a functionally tricistronic gene that expresses one structural and two nonstructural proteins in infected cells. Virology. (2001) 290:181–91. doi: 10.1006/viro.2001.1159
3. Liu H-J, Lin P-Y, Wang L-R, Hsu H-Y, Liao M-H, Shih W-L. Activation of small GTPases RhoA and Rac1 is required for avian reovirus p10-induced syncytium formation. Mol Cells. (2008) 26:396–403.
4. Liu H-J, Lin P-Y, Lee J-W, Hsu H-Y, Shih W-L. Retardation of cell growth by avian reovirus p17 through the activation of p53 pathway. Biochem Biophys Res Commun. (2005) 336:709–15. doi: 10.1016/j.bbrc.2005.08.149
5. Chulu JL, Huang WR, Wang L, Shih WL, Liu HJ. Avian reovirus nonstructural protein p17-induced G2/M cell cycle arrest and host cellular protein translation shutoff involve activation of p53-dependent pathways. J Virol. (2010) 84:7683–94. doi: 10.1128/JVI.02604-09
6. Calvo PG, Fox GC, Parrado XLH, Llamas-Saiz AL, Costas C, Martínez-Costas J, et al. Structure of the carboxy-terminal receptor-binding domain of avian reovirus fibre sigmaC. J Mol Biol. (2005) 354:137–49. doi: 10.1016/j.jmb.2005.09.034
7. Wickramasinghe R, Meanger J, Enriquez CE, Wilcox GE. Avian reovirus proteins associated with neutralization of virus infectivity. Virology. (1993) 194:688–96. doi: 10.1006/viro.1993.1309
8. Lu H, Tang Y, Dunn PA, Wallner-Pendleton EA, Lin L, Knoll EA. Isolation and molecular characterization of newly emerging avian reovirus variants and novel strains in Pennsylvania, USA, 2011–2014. Sci Rep. (2015) 5:14727. doi: 10.1038/srep14727
9. De la Torre D, Astolfi-Ferreira CS, Chacón RD, Puga B, Piantino Ferreira A. Emerging new avian reovirus variants from cases of enteric disorders and arthritis/tenosynovitis in Brazilian poultry flocks. Br Poult Sci. (2021) 62:361–72. doi: 10.1080/00071668.2020.1864808
10. Pantin-Jackwood MJ, Day JM, Jackwood MW, Spackman E. Enteric viruses detected by molecular methods in commercial chicken and turkey flocks in the United States between 2005 and 2006. Avian Dis. (2008) 52:235–44. doi: 10.1637/8174-111507-Reg.1
11. Wang D, Shi J, Yuan Y, Zheng L, Zhang D. Complete sequence of a reovirus associated with necrotic focus formation in the liver and spleen of muscovy ducklings. Vet Microbiol. (2013) 166:109–22. doi: 10.1016/j.vetmic.2013.05.022
12. Styś-Fijoł N, Kozdruń W, Czekaj H. Detection of avian reoviruses in wild birds in Poland. J Vet Res. (2017) 61:239. doi: 10.1515/jvetres-2017-0033
13. Sakai K, Ueno Y, Ueda S, Yada K, Fukushi S, Saijo M, et al. Novel reovirus isolation from an Ostrich (Struthio camelus) in Japan. Vet Microbiol. (2009) 134:227–32. doi: 10.1016/j.vetmic.2008.08.022
14. Yan T, Guo L, Jiang X, Wang H, Yao Z, Zhu S, et al. Discovery of a novel recombinant avian orthoreovirus in China. Vet Microbiol. (2021) 260:109094. doi: 10.1016/j.vetmic.2021.109094
15. De Gussem J, Swam H, Lievens K, De Herdt P. Reovirus tenosynovitis in a flock of layer breeders. Avian Pathology. (2010) 39:169–70. doi: 10.1080/03079451003717597
16. Wang H, Gao B, Liu X, Zhang S, Diao Y, Tang Y. Pathogenicity of a variant duck orthoreovirus strain in Cherry Valley Ducklings. Vet Microbiol. (2020) 242:108546. doi: 10.1016/j.vetmic.2019.108546
17. Mor SK, Sharafeldin TA, Porter RE, Ziegler A, Patnayak DP, Goyal SM. Isolation and characterization of a turkey arthritis reovirus. Avian Dis. (2013) 57:97–103. doi: 10.1637/10353-090712-Reg.1
18. Neelima S, Ram G, Kataria J, Goswami T. Avian reovirus induces an inhibitory effect on lymphoproliferation in chickens. Vet Res Commun. (2003) 27:73–85. doi: 10.1023/A:1022014825451
19. Van de Zande S, Kuhn E-M. Central nervous system signs in chickens caused by a new avian reovirus strain: a pathogenesis study. Vet Microbiol. (2007) 120:42–9. doi: 10.1016/j.vetmic.2006.10.024
20. Troxler S, Rigomier P, Bilic I, Liebhart D, Prokofieva I, Robineau B, et al. Identification of a new reovirus causing substantial losses in broiler production in France, despite routine vaccination of breeders. Vet Rec. (2013) 172:556. doi: 10.1136/vr.101262
21. Palomino-Tapia V, Mitevski D, Inglis T, van der Meer F, Abdul-Careem MF. Molecular characterization of emerging avian reovirus variants isolated from viral arthritis cases in Western Canada 2012–2017 based on partial sigma (σ) C gene. Virology. (2018) 522:138–46. doi: 10.1016/j.virol.2018.06.006
22. De Carli S, Wolf JM, Gräf T, Lehmann FK, Fonseca AS, Canal CW, et al. Genotypic characterization and molecular evolution of avian reovirus in poultry flocks from Brazil. Avian Pathology. (2020) 49:611–20. doi: 10.1080/03079457.2020.1804528
23. Kant A, Balk F, Born L, van Roozelaar D, Heijmans J, Gielkens A, et al. Classification of Dutch and German avian reoviruses by sequencing the $\sigma $ C protein. Vet Res. (2003) 34:203–12. doi: 10.1051/vetres:2002067
24. Goldenberg D, Pasmanik-Chor M, Pirak M, Kass N, Lublin A, Yeheskel A, et al. Genetic and antigenic characterization of sigma C protein from avian reovirus. Avian Pathology. (2010) 39:189–99. doi: 10.1080/03079457.2010.480969
25. Teng L, Xie Z, Xie L, Liu J, Pang Y, Deng X, et al. Sequencing and phylogenetic analysis of an avian reovirus genome. Virus Genes. (2014) 48:381–6. doi: 10.1007/s11262-013-1029-5
26. Noh JY, Lee DH, Lim TH, Lee JH, Day JM, Song CS. Isolation and genomic characterization of a novel avian orthoreovirus strain in Korea, 2014. Arch Virol. (2018) 163:1307–16. doi: 10.1007/s00705-017-3667-8
27. de Kloet SR. Sequence analysis of four double-stranded RNA genomic segments reveals an orthoreovirus with a unique genotype infecting psittaciformes. Avian Dis. (2008) 52:480–6. doi: 10.1637/8212-011908-Reg.1
28. Huhtamo E, Uzcátegui NY, Manni T, Munsterhjelm R, Brummer-Korvenkontio M, Vaheri A, et al. Novel orthoreovirus from diseased crow, Finland. Emerg Infect Dis. (2007) 13:1967. doi: 10.3201/eid1312.070394
29. Forzán MJ, Renshaw RW, Bunting EM, Buckles E, Okoniewski J, Hynes K, et al. A novel orthoreovirus associated with epizootic necrotizing enteritis and splenic necrosis in American crows (Corvus brachyrhynchos). J Wildl Dis. (2019) 55:812–22. doi: 10.7589/2019-01-015
30. Dandar E, Huhtamo E, Farkas SL, Oldal M, Jakab F, Vapalahti O, et al. Complete genome analysis identifies Tvärminne avian virus as a candidate new species within the genus Orthoreovirus. J Gen Virol. (2014) 95:898–904. doi: 10.1099/vir.0.060699-0
31. Lawson B, Dastjerdi A, Shah S, Everest D, Núñez A, Pocknell A, et al. Mortality associated with avian reovirus infection in a free-living magpie (Pica pica) in Great Britain. BMC Vet Res. (2015) 11:20. doi: 10.1186/s12917-015-0329-5
32. Kugler R, Dandár E, Fehér E, Jakab F, Mató T, Palya V, et al. Phylogenetic analysis of a novel reassortant orthoreovirus strain detected in partridge (Perdix perdix). Virus Res. (2016) 215:99–103. doi: 10.1016/j.virusres.2015.11.018
33. Mor SK, Armién AG, Reed L, Schott R, Goyal SM. Detection and molecular characterization of a reovirus in black-capped chickadees (Poecile atricapillus) from Minnesota, USA. J Wildl Dis. (2014) 50:928–32. doi: 10.7589/2014-02-045
34. Ogasawara Y, Ueda H, Kikuchi N, Kirisawa R. Isolation and genomic characterization of a novel orthoreovirus from a brown-eared bulbul (Hypsipetes amaurotis) in Japan. J Gen Virol. (2015) 96:1777–86. doi: 10.1099/vir.0.000110
35. Yu K, Li Y, Han H, Song M, Ma X, Liu C, et al. Complete genome sequence of an avian reovirus isolated from wild mallard ducks in China. Genome Announc. (2014) 2:e00813–14. doi: 10.1128/genomeA.00813-14
36. Gough R, Alexander D, Collins M, Lister S, Cox W. Routine virus isolation or detection in the diagnosis of diseases in birds. Avian Pathology. (1988) 17:893–907. doi: 10.1080/03079458808436511
37. Vindevogel H, Meulemans G, Pastoret PP, Schwers A, Calberg-Bacq CM. Reovirus infection in the pigeon. In: Annales de Recherches Vétérinaires. (1982) 13:149–52.
38. Hlinak A, Müller T, Kramer M, Mühle R-U, Liebherr H, Ziedler K. Serological survey of viral pathogens in bean and white-fronted geese from Germany. J Wildl Dis. (1998) 34:479–86. doi: 10.7589/0090-3558-34.3.479
39. Lee Y-N, Cheon S-H, Lee E-K, Heo G-B, Bae Y-C, Joh S-J, et al. Pathogenesis and genetic characteristics of novel reassortant low-pathogenic avian influenza H7 viruses isolated from migratory birds in the Republic of Korea in the winter of 2016–2017. Emerg Microbes Infect. (2018) 7:182. doi: 10.1038/s41426-018-0181-3
40. Caterina KM, Frasca S Jr, Girshick T, Khan MI. Development of a multiplex PCR for detection of avian adenovirus, avian reovirus, infectious bursal disease virus, and chicken anemia virus. Mol Cell Probes. (2004) 18:293–8. doi: 10.1016/j.mcp.2004.04.003
41. Liu Q, Zhang G, Huang Y, Ren G, Chen L, Gao J, et al. Isolation and characterization of a reovirus causing spleen necrosis in Pekin ducklings. Vet Microbiol. (2011) 148:200–6. doi: 10.1016/j.vetmic.2010.09.016
42. Lee D-H, Lee H-J, Lee Y-J, Kang H-M, Jeong O-M, Kim M-C, et al. DNA barcoding techniques for avian influenza virus surveillance in migratory bird habitats. J Wildl Dis. (2010) 46:649–54. doi: 10.7589/0090-3558-46.2.649
43. Jeong S, Lee D-H, Kwon J-H, Kim Y-J, Lee S-H, Cho AY, et al. Highly pathogenic avian influenza clade 2.3. 4.4 b subtype H5N8 virus isolated from mandarin duck in South Korea, 2020. Viruses. (2020) 2:1389. doi: 10.3390/v12121389
44. Ratnasingham S, Hebert PD. BOLD The Barcode of Life Data System (http://www.barcodinglifeorg). Mol Ecol Notes. (2007) 7:355–64. doi: 10.1111/j.1471-8286.2007.01678.x
45. Reck C, Menin Á, Canever MF, Miletti LC. Rapid detection of mycoplasma synoviae and avian reovirus in clinical samples of poultry using multiplex PCR. Avian Dis. (2013) 57:220–4. doi: 10.1637/10425-101712-Reg.1
46. Belkhir K. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations (2004). Available online at: https://kimura.univ-montp2.fr/genetix/
47. Kumar S, Stecher G, Li M, Knyaz C, Tamura K, MEGA X. molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:1547. doi: 10.1093/molbev/msy096
48. Strimmer K, Forslund K, Holland B, Moulton V. A novel exploratory method for visual recombination detection. Genome Biol. (2003) 4:R33. doi: 10.1186/gb-2003-4-5-r33
49. Docherty D, Converse K, Hansen W, Norman G. American woodcock (Scolopax minor) mortality associated with a reovirus. Avian Dis. (1994):899–904. doi: 10.2307/1592132
50. Sharafeldin TA, Mor SK, Sobhy NM, Xing Z, Reed KM, Goyal SM, et al. A newly emergent Turkey arthritis reovirus shows dominant enteric tropism and induces significantly elevated innate antiviral and t Helper-1 cytokine responses. PLoS ONE. (2015) 10:e0144085. doi: 10.1371/journal.pone.0144085
51. Takehara K, Chinen O, Jahangir A, Miyoshi Y, Ueno Y, Ueda S, et al. Ceramic powder made from chicken feces: anti-viral effects against avian influenza viruses. Avian Dis. (2009) 53:34–8. doi: 10.1637/8382-062008-Reg.1
52. Tamehiro CY, Alfieri AF, Médici C, Alfieri AA. Segmented double-stranded genomic RNA viruses in fecal samples from broiler chicken. Braz J Microbiol. (2003) 34:349–53. doi: 10.1590/S1517-83822003000400013
53. Egaña-Labrin S, Hauck R, Figueroa A, Stoute S, Shivaprasad H, Crispo M, et al. Genotypic characterization of emerging avian reovirus genetic variants in California. Sci Rep. (2019) 9:9351. doi: 10.1038/s41598-019-45494-4
54. Lin P-Y, Liu H-J, Chang C-D, Chang C-I, Hsu J-L, Liao M-H, et al. Avian reovirus S1133-induced DNA damage signaling and subsequent apoptosis in cultured cells and in chickens. Arch Virol. (2011) 156:1917. doi: 10.1007/s00705-011-1063-3
55. Zhang Z, Lin W, Li X, Cao H, Wang Y, Zheng SJ. Critical role of eukaryotic elongation factor 1 alpha 1 (EEF1A1) in avian reovirus sigma-C-induced apoptosis and inhibition of viral growth. Arch Virol. (2015) 160:1449–61. doi: 10.1007/s00705-015-2403-5
56. Shmulevitz M, Duncan R. A new class of fusion-associated small transmembrane (FAST) proteins encoded by the non-enveloped fusogenic reoviruses. EMBO J. (2000) 19:902–12. doi: 10.1093/emboj/19.5.902
57. Costas C, Martínez-Costas J, Bodelón G, Benavente J. The second open reading frame of the avian reovirus S1 gene encodes a transcription-dependent and CRM1-independent nucleocytoplasmic shuttling protein. J Virol. (2005) 79:2141–50. doi: 10.1128/JVI.79.4.2141-2150.2005
58. Ayalew LE, Gupta A, Fricke J, Ahmed KA, Popowich S, Lockerbie B, et al. Phenotypic, genotypic and antigenic characterization of emerging avian reoviruses isolated from clinical cases of arthritis in broilers in Saskatchewan, Canada. Sci Rep. (2017) 7:3565. doi: 10.1038/s41598-017-02743-8
59. Guo K, Dormitorio TV, Ou S-C, Giambrone JJ. Development of TaqMan real-time RT-PCR for detection of avian reoviruses. J Virol Methods. (2011) 177:75–9. doi: 10.1016/j.jviromet.2011.06.022
Keywords: avian reoviruses, wild bird, chicken, South Korea, σC-encoding gene, σNS-encoding genes
Citation: Kim S-W, Choi Y-R, Park J-Y, Wei B, Shang K, Zhang J-F, Jang H-K, Cha S-Y and Kang M (2022) Isolation and Genomic Characterization of Avian Reovirus From Wild Birds in South Korea. Front. Vet. Sci. 9:794934. doi: 10.3389/fvets.2022.794934
Received: 14 October 2021; Accepted: 06 January 2022;
Published: 28 January 2022.
Edited by:
Christina Leyson, United States Department of Agriculture, United StatesReviewed by:
Rodrigo A. Gallardo, University of California, Davis, United StatesLuiz Claudio Miletti, University of the State of Santa Catarina, Brazil
Copyright © 2022 Kim, Choi, Park, Wei, Shang, Zhang, Jang, Cha and Kang. 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: Min Kang, vet.minkang@gmail.com; Se-Yeoun Cha, seyeouncha@gmail.com
†These authors have contributed equally to this work