Abstract

Novel duck orthoreovirus can infect various hosts, mainly causing spleen necrosis and hemorrhagic necrotizing hepatitis, leading to immunosuppression and growth impairment. Infectious diseases characterized by swollen hock joints, movement disorders, and impaired growths have recently occurred in significant duck breeding areas in China. Infected ducks from Shandong Province were selected, and N-DRV-LY20 was isolated from the swollen joints. Next-generation sequencing (NGS) of N-DRV-LY20 and homology analysis showed that all ten fragments had the highest homology to N-DRV and the closest affinity to Duck/N-DRV-XT18/China/2018. After artificial inoculation with N-DRV-LY20 virus solution, Cherry Valley ducks showed the same symptoms of leg swelling and hock joint deformation, which were not explored in previous studies. A long history of evolution has led to significant variation in DRV and altered its pathogenicity. In this study, we confirmed the ability of N-DRV to cause viral arthritis in ducks and clarified this isolate’s genetic characteristics and level of variation, providing a basis for the control and variation detection of N-DRV.

1. Introduction

Avian orthoreovirus (ARV) was first successfully isolated in China in 1991 [1]. ARV belongs to the family Eutheroviridae, subfamily Spiny Eutherovirus, and genus Orthoreovirus [2], which has no envelope, and the viral particles are icosahedral symmetric sphere [3]. The virus induces cell fusion [4, 5], lacks hemagglutination activity [6], and leads to immunosuppression in poultry, secondary to other pathogenic infections [7]. A variety of avian species can be infected with ARV, such as chickens, turkeys, ducks, geese, pigeons, and wild birds [710], causing viral arthritis, tenosynovitis, growth retardation, and chronic respiratory diseases [11]. Based on sequence nucleotide homology and evolutionary analysis, waterfowl orthoreovirus (WRV) can be classified into two genotypes. Muscovy duck orthoreovirus (MDRV) and goose orthoreovirus (GRV) are classified as waterfowl orthoreovirus genotype 1 [12, 13]. Novel goose orthoreovirus (N-GRV) and novel duck orthoreovirus (N-DRV) are classified as waterfowl orthoreovirus genotype 2 [14].

Novel duck orthoreovirus is an infectious disease that mainly causes spleen necrosis in ducks, which is prevalent mainly in Hebei, Shandong, Jiangsu, Anhui, Zhejiang, Inner Mongolia, and Fujian [15, 16]. N-DRV has a wide range of natural hosts, can be transmitted horizontally or vertically, and has a high morbidity and mortality rate in ducklings. N-DRV has the same cultural characteristics as ARV that can be isolated and proliferated in various ways [17]. There are ten segments of double-stranded RNA in the genome of N-DRV, which are divided into three groups: large segments (L1, L2, and L3), medium segments (M1, M2, and M3), and small segments (S1, S2, S3, and S4) [18, 19]. The L fragment encodes λ proteins (λA, λB, and λC), the M fragment encodes μ proteins (μA, μB, and μNS), the S fragment encodes σ proteins (σC, σA, σB, and σNS), and the S1 gene additionally encodes proteins and [4, 2022]. Current studies of the N-DRV genome sequence have focused on the S-segment. The S1 gene encodes the σC protein, an adhesion protein located on the surface of the viral capsid [4], which attaches to host cells and is a significant site for binding cellular receptors, and also induces apoptosis [23]. In addition, the σC protein induces neutralizing antibodies against orthoreovirus [24]. The main functions of each protein encoded by each gene segment are shown in Table S1.

Viral arthritis in ducks has been emerging in the duck breeding areas of east-central China for the past two years, with the scope of infection expanding, characterized by swollen hock joints, movement disorders, decreased immune function, and impaired growth. Dissection revealed severe joint swelling and necrosis of the liver and spleen. To determine the pathogen causing the disease, we separately isolated viruses and bacteria from the swollen joints of diseased ducks; and finally, obtained the stably surviving virus strain. After whole genome sequencing and evolutionary genetic analysis, we identified the isolated virus as a novel orthoreovirus. The pathogenicity of this isolate was subsequently studied by artificially infecting Cherry Valley meat ducks. It was found that although the N-DRV-LY20 was genetically similar to the Duck/N-DRV-XT18/China/2018, the clinical symptoms caused were significantly different. This study improves our understanding of avian orthoreovirus’s evolutionary relationships and diversity. It provides a scientific basis for preventing and controlling new orthoreovirus infections in ducks and lays the foundation for further understanding the pathogenic mechanisms of orthoreovirus.

2. Materials and Methods

2.1. Ethics Statement

All the animal infection experiments were approved by the Animal Care and Use Committee of Shandong Agricultural University and conducted following the “Guidelines for Experimental Animals” of the Ministry of Science and Technology (Beijing, China).

2.2. Clinical Samples

Tendon samples of diseased ducks characterized by arthritis were used for virus isolation in Shandong Province, China. The existence of common avian pathogenic viruses such as avian influenza virus, duck virus hepatitis A types 1, duck virus hepatitis A types 3, Newcastle disease virus, avian orthoreovirus, and duck plague virus was excluded by the RT-PCR or PCR methods (Table S2) [25].

2.3. RNA Extraction and RT-PCR

According to the manufacturer’s instructions, total RNA from the tissue filtrate was extracted using the RNAprep Pure Tissue Kit (TIANGEN, Beijing, China). Extracted RNA samples were amplified using One-Step RT-PCR Kit (TaKaRa, Dalian, China) with specific primers (N-DRV-F 5′-ATG​GAT​CGC​AAC​GAG​GTG​ATA​C-3′, N-DRV-R 5′-CTA​GCC​CGT​GGC​GAC​GGT-3′) to a partial sequence of 966 bp in the N-DRV σC gene with the following reaction conditions: 42°C for 50 min, 94°C for 5 min, 94°C for 30 s, 54°C for 30 s, 72°C for 45 s, and 32 cycles. Beijing Genomics Institution (BGI) synthesized all primers. The PCR products were collected by electrophoresis in 1.0% agarose gel, and PCR-amplified fragments were sequenced and confirmed by Sanger sequencing (TsingKe Biotech, Beijing, China) for confirmation.

2.4. Viral Isolation

The collected tendon and swollen joint samples were processed by grinding and homogenizing after adding four times saline. The supernatant solution obtained after centrifugation was repeatedly frozen and thawed three times. We used LMH (Leghorn Male-chicken Hepatocellular-carcinoma, ATCC CRL-2117) cells inoculated with treated supernatant, Dulbecco’s Modified Eagle Medium (DMEM; Hyclone, Beijing, China) containing 2% fetal bovine serum (FBS; Hyclone, Beijing, China) was used as a cell growth maintenance solution, and the culture conditions were 37°C and 5% CO2. Virus cultures were collected and continuously re-cultured until stable CPE was observed. After sufficient dilution of the virus solution, the virus was purified using the plaque purification technology, and the final purified cell culture was frozen.

The titer of N-DRV-LY20 in LMH cells was determined according to the half-tissue culture infectious dose (TCID50). LMH cells were seeded in 96-well plates, inoculated with 105 cells per well, and incubated for 48 h at 37°C under 5% CO2. The viral solution was serially diluted 10-fold, and 100 μl was added to each well. Measure CPE 4 days after inoculation. Virus titers were determined according to the method of Reed and Muench [26].

2.5. Next-Generation Sequencing and Analysis

Infected cell cultures of the N-DRV-LY20 strain were extracted from viral RNA using the TruSeq Stranded Total RNA Sample Prep Kit (Illumina, Shanghai, China) to construct sequenced RNA libraries. The library fragments were purified using the Ampre XP (Beckman-Coulter) and quantified using the Bioanalyzer 2100 system [7, 27]. The prepared cDNA libraries were loaded onto the Illumina MiSeq platform to assemble the viral genome, reads of nontarget genes were removed, and the raw reads of NGS data were analyzed. Assemble all pure reads using QIAGEN software and align sequences. Virus (ORF) prediction, amino acid (aa) translation, and nucleotide (nt) sequence alignment were performed using Editseq software in DNASTAR software and Megalign software [28, 29]. The BLAST program on the NCBI website was used to search for published related isolates and further analyze homology relationships. Information on the viral strains used in the analysis is presented in Table S3. Genetic evolution trees were constructed using the Clustal W program in MEGA 11 software.

2.6. Animals and Samples Collection

One hundred and twenty healthy Cherry Valley ducklings (1 day old) were purchased from Liuhe Jingwei Agricultural and Livestock Company in Tai’an, Shandong Province. The ducklings were tested negative for N-DRV and other pathogens in blood and cloacal swabs using qPCR and negative for antibodies to N-DRV in serum using ELISA before the experiment. They were randomly divided into four groups of thirty ducks each. Three experimental groups (A, B, and C) were set up to be inoculated with the virus by footpad, intramuscular, and oral administration. Each duckling was inoculated with 0.2 ml of N-DRV-LY20 virus (TCID50: 10−5.8/0.1 ml), and Group D was inoculated with 0.2 ml of saline as the control group. Later, all ducklings were continuously observed daily and recorded clinical symptoms, weight changes, and mortality following infection.

Three ducks in each group were euthanized by CO2 on days 1, 3, 5, 7, 10, 14, 18, and 22 following infection, and blood and cloacal swabs were collected. Various tissue specimens were collected from Fabricius, the thymus, and the tendons’ liver, spleen, lung, kidney, and bursa. One part of the tissue sample set was also fixed in 10% neutral buffered formalin to make paraffin sections. Sections were examined microscopically with hematoxylin and eosin staining.

2.7. Genomic RNA Extraction and Detection of Viral Load

The extraction of total RNA from all samples was performed using the MiniBEST Universal RNA Extraction Kit (TIANGEN, Beijing, China) according to the instructions, and the concentrations were measured using a DeNovix DS-11 spectrophotometer (Nanodrop, USA) before being stored at −80°C. The viral load was detected using the N-DRV TaqMan probe fluorescent quantitative RT-PCR method established previously in our laboratory [25, 30]. The primers and the TaqMan probe (Table S2) used to detect the viral load were designed based on the S1 gene of the N-DRV strain.

One-Step PrimeScript™ RT-PCR kit (TaKaRa, Dalian, China) was used to perform qPCR reactions. The reaction system contained 2 × One-Step RT-PCR Buffer III 10 μl, 0.4 μl each of TaKaRa Ex Taq HS, PrimeScript RT Enzyme Mix II, PCR Forward primer and PCR Reverse primer, total RNA 80 ng. Finally, RNase-free H2O was added to make a total volume of 25 μl. Roche LightCycler 96 Real-Time PCR System (Roche, Switzerland) system was used with a reaction program of 42°C for 5 min, 95°C for 10 s, 95°C for 5 s, and 60°C for 20 s for 40 cycles, and fluorescence signals were collected.

2.8. Antibody Level Detection

We used an indirect ELISA method to detect N-DRV antibody levels. The blood was first incubated at 37°C for 1 h and centrifuged at 3000 r/min for 10 min to collect the serum. The 96-well plate was precoated with 2000-fold dilution N-DRV-σC protein in 0.05 M carbonate buffer , then washed three times with PBST (0.01 M PBS, , 0.05% tween 20) for 3 min each time the next day. Add 200 μl blocking solution (5% skim milk powder) and leave at 37°C for 1 h. Dilute serum 1 : 20 times with 5% skim milk powder, add 100 μl per well and leave at 37°C for 1 h. Dilute HRP-labeled antiduck secondary antibody (KPL) (Solarbio, China) 1 : 500 times with 5% skim milk powder, then add 100 μl per well and leave at 37°C for 1 h. Add 100 μl of TMB color development solution under light-proof conditions and let it stand for 15 min at 37°C before adding 50 μl of termination solution (2 M sulfuric acid). The final reading was taken at an OD value of 450 using a microplate reader.

2.9. Statistical Analysis

Data on body weight, viral load, and blood antibody levels were obtained and statistically analyzed using GraphPad Prism software 9.0 (GraphPad Software Inc.) and were expressed as means ± SD. The data from all experimental and control groups were compared in the same graph. was considered a statistically significant difference.

3. Results

3.1. N-DRV Isolation and Identification

Third-generation cell cultures showed obvious cytopathic effects (CPE), with clear syncytial features at 24 h, and 90% of cells showed cell shedding at about 72 h (Figure 1). ARV causes rounding and shedding of LMH cells, forming a characteristic giant or bloom-like subrounded syncytial CPE [31]. Using RT-PCR amplification (with F/R primers based on the S1 gene design), we confirmed that the cell culture tested N-DRV positive. We named the isolate Duck/N-DRV-LY20/China/2021. Meanwhile, RT-PCR or PCR was used to detect other common pathogens, and the results showed that all other viruses tested negative. According to the method of Reed and Muench, the virus titer was determined as TCID50 = 10−5.8/0.1 ml.

3.2. NGS Analysis

After viral genome assembly, we obtained 10 contigs that covered the full length of ARV. The entire genome of the N-DRV-LY20 strain was 23,419 bp in length and contained 10 genomic fragments ranging in size from 1191 bp (S4) to 3959 bp (L1). Ten segments of the whole genome sequence of the N-DRV-LY20 strain were uploaded to GenBank with the accession numbers ON040907 (L1) to ON040916 (S4). We determined ORFs by predicting each segment’s ORF start and stop codons. The S1 gene contained three ORFs for the (98aa), (163aa), and σC (322aa) proteins, and the remaining gene segments contained only one ORF, each encoding one protein (Table 1). Sequence analysis of the 5′ and 3′ untranslated regions (UTR) showed that the nucleotide sequences of the fragments of the N-DRV-LY20 isolate remained conserved. They were consistent with the published sequences for ARV, DRV, and GRV: the identical GCUUUUU sequence at the 5′ UTR and the UCAUC sequence at the 3′ UTR. No other pathogens were found during the analysis of the samples.

3.3. Phylogenetic Analysis

Nucleotide sequences of the N-DRV-LY20 strain were analyzed for nucleotide (nt) and amino acid (aa) homology with avian-derived, goose-derived, and duck-derived orthoreoviruses published on GenBank (Table S4). The results showed that strain LY20 had the highest homology with N-DRV among all fragments and was most closely related to strain N-DRV-XT18. The S1 gene had the lowest homology with other sequences, a relatively small genetic distance among strains, and a close genetic evolutionary relationship with the ARVbranch because the S1 gene encodes the σC protein of ARV and DRV, whereas the σC protein of MDRV is encoded by the S4 gene [4]. The homology comparison showed that the S1 gene of the N-DRV-LY20 had higher homology with the N-DRV-XT18 (nt: 98.3%; aa: 97.7%) and lower homology with the ARV/S1133 (nt: 45.2%; aa: 11.3%) and the MDRV/815-12 (nt: 33.4%; aa: 14.5%). Similarly, the σC gene of N-DRV-LY20 had the highest homology with the N-DRV-XT18 (nt: 98.9%; aa: 99.1%) and lower homology with the MDRV/815-12 (nt: 49.6%; aa: 42.2%) and the ARV/S1133 (nt: 39.6%; aa: 28.3%).

The results of the genetic evolution analysis showed that the avian orthoreovirus formed three evolutionary branches, namely, the ARV branch, MDRV branch, and N-DRV branch, in the genetic evolution analysis of 10 segments (Figure 2). In the evolutionary genetic analysis tree of L1, L2, L3, M1, M3, S2, and S4 genes, N-DRV was more closely related to MDRV, and ARV was in a separate branch. In the evolutionary genetic analysis tree of the M2, S1, and S3 genes, N-DRV was more closely related to ARV, and MDRV was in a separate evolutionary branch.

3.4. Pathogenicity of N-DRV-LY20 in Ducklings

Thirty days of clinical sign monitoring were used to study the pathogenicity of the N-DRV-LY20 isolate in Cherry Valley ducklings. After ducklings were infected with N-DRV, the clinical signs caused by different inoculation methods were similar, with the most severe signs in the footpad injection group. The morbidity rate in the infected group was 100%, and mortality rates differed, with 30% in the footpad injection group, 20% in the intramuscular injection group, and 17% in the oral group. All deaths of ducklings were concentrated within 7 dpi (Figure 3(a)). The weight monitoring results indicated that artificially infected ducklings’ weight gain was suppressed (Figure 3(b)). The symptoms were mainly mental lethargy, decreased feeding and drinking, loose white stools, slow growth, and a sparse coat of the affected ducks. In addition, diseased ducks showed severe movement disorders with obvious deformation of joints or paralysis (Figures 3(c)3(f)).

It was observed that the symptoms caused by different inoculation methods were consistent with the naturally infected ducks. Liver necrosis existed in all groups of ducklings (Figure 4(a)). Hemorrhagic diffuse necrotic spots of variable size with black necrotic foci visible at the margins were seen on the liver. The spleen was severely lesioned (Figure 4(b)), swollen at 1 dpi, appeared black necrotic spots at 2 dpi, and had lost its original texture and color, turning gray-black at 3 dpi dissection of the joints showed subcutaneous hemorrhage, suppuration of the joint cavity, and tibial deformation (Figure 4(c)). Other tissues and organs were damaged to different degrees, mainly showing intestinal hemorrhage and thinning of the intestinal wall; bleeding and black necrotic spots in the lungs; swollen and hemorrhage in the thymus; and a swollen hemorrhagic kidney.

3.5. Pathological Changes of N-DRV-LY20 in Ducklings

Different virus inoculation methods caused no significant difference in pathological changes; no abnormal changes were seen in the control group. The spleen was extensively necrotic and subsequently replaced by fibrous tissue. The demarcation between the cortex and medulla was blurred, with massive erythrocyte filling, increased macrophages, and inflammatory cell infiltration (Figure 5(a)). It showed steatosis, cellular hypertrophy or atrophy, disorganized hepatic cords, hepatic sinus hemorrhage, neutrophil infiltration, lymphocytosis, and foci of coagulative necrosis in the liver (Figure 5(b)). The bursa of Fasciola showed interstitial edema with reduced lymphocytes, and scattered inflammatory cells in the necrotic area (Figure 5(c)). Tendon fibers were widened, broken, and lost their original fibrous structure. Inflammatory cells proliferated in the tendon sheath (Figure 5(d)). Widening and degeneration of interstitial myocardial fibers with inflammatory cell infiltration. The renal tubules and the glomeruli are necrotic and atrophic, with neutrophil infiltration, epithelial swelling of the renal tubules, narrowing of the tubular lumen, and interstitial vasodilatation or congestion. The thymus was hemorrhagic with erythrocyte sludge. The alveolar tissue is destroyed or atrophied, and lymphocytes predominantly infiltrate the inflammatory cells.

3.6. Viral Replication and Shedding Pattern

The viral replication pattern in various tissues and the shedding pattern of Cherry Valley ducks infected with the N-DRV-LY20 can be derived using TaqMan probe qPCR at different times and organs. Virus levels and their changes were more consistent with different inoculation methods but not similar between different organs (Figures 6(a)6(g)). Viral content in the spleen was significantly higher than in other organs. Positive data could be monitored as early as 1 dpi in all organs and continued to rise for several days before reaching a peak, followed by a gradual decrease, and then maintained at a relatively stable level. The spleen, lung, and liver maintained high viral levels at 3–7 dpi. Differently, the bursa of Fabricius, the thymus, and kidney had elevated levels at 3–5 dpi, while a small peak occurred at 22 dpi in the bursa of Fabricius and the thymus. When comparing the three infection groups, it was found that the virus replication rate of the organs in the oral group rose slower than in the other two groups.

It is rapid for the virus to enter the bloodstream and excrete from the body, with similar trends in virus levels and changes between various infection groups. High amounts of the virus could be detected in blood and cloacal swab samples at 1 dpi (Figures 6(h)6(i)). Virus levels in the blood increased rapidly over 3 dpi, and then slowly decreased until a rebound occurred at 14 dpi. There was also a tendency for the virus number in the cloaca to rising then fall. The peak in groups A and B occurred at 5 dpi, while group C occurred at 3 dpi, and the virus content then gradually decreased but still showed two slight rises at 14 dpi and 22 dpi.

3.7. Antibody Response in Affected Ducks

The ELISA method measured the antibody levels in the serum of infected ducks to understand the protective effect produced against the virus. All infected groups differed from the control group regarding antibody levels and response effects. A comparison of the different experimental groups showed that the antibody levels changed fastest in the footpad injection group, while they were produced least in the oral group, but still higher than the control group. The changes in antibody levels in the footpad injection group were similar to those in the control group, with the highest antibody levels at 10 dpi (Figure 6(j)). For the intramuscular group, antibody levels tended to increase after inoculation, reaching a maximum of 18 dpi, where the changes were more moderate (Figure 6(k)). While for the oral group, there is a distinct peak at 10 dpi and two less distinct peaks at 1 dpi and 18 dpi (Figure 6(l)).

4. Discussion

Avian orthoreovirus (ARV) was first isolated in 1954 and reported in China in 1985, which frequently causes viral arthritis and tenosynovitis in various avian species [32]. Several waterfowl orthoreovirus have been isolated successively in China, showing significant differences in host and pathogenicity, separately causing white spot disease, hemorrhagic necrotizing hepatitis, and spleen necrosis [13, 14]. The novel duck orthoreovirus was first reported to cause duck spleen necrosis in 2016 [33]. Over the past decade, the disease has become a severe disease that seriously affects the duck breeding industry, evolving with increasing lethality and different clinical features [25, 34]. Diseased ducks rapidly contaminate feed and drinking water through the digestive and respiratory tracts. Tolerance of the virus to the environment allows it to survive for a long time, resulting in a wider spread of the pathogen. In the past two years, we have found a frequent occurrence of viral arthritis in major duck areas in China, mainly characterized by swollen joints, unstable standing, and impaired growth. We finally confirmed that it was caused by novel duck orthoreovirus through pathogen isolation and identification, which differed from the previously reported symptoms. The connection between the novel duck orthoreovirus and viral arthritis was not discussed previously.

After successfully isolating Duck/N-DRV-LY20/China/2021, we sequenced the entire genome using NGS and compared it with other orthoreovirus species at the nucleotide and amino acid levels. Data analysis showed that the 10 genomic fragments of N-DRV-LY20 were the same length as the homologous fragments of the previous N-DRV isolates, with highly conserved sequences at the 5′ and 3′ ends, and showed the closest affinity to N-DRV-XT18. The genetic relationship of the L segment is stable. Still, the significant amino acid variation suggests that most of the nucleotides of the L gene should be synonymous mutations, excluding the variable region. The σC protein has a significant role in viral attachment and becomes the most variable part of all ARV proteins under frequent selection pressure [23]. However, the variation was not significant in N-DRV-LY20 compared with N-DRV-XT18 between σC gene fragments (nt: 98.9%; aa: 99.1%) but more prominent between M1 fragments (nt: 96.1%; aa: 90.7%), suggesting that a new highly variable region may exist in M1 gene fragments.

To better understand the pathogenicity of N-DRV-LY20, we successfully conducted animal regression experiments in Cherry Valley ducks, and the virus was inoculated by footpad, intramuscular, and oral injection routes. All groups consistent with the natural disease state showed severe swelling of the leg joints and marked necrosis of the liver and spleen. Ducks inoculated via footpad were sometimes more susceptible to N-DRV than intramuscularly, especially when the duck’s flippers were ruptured and infected by the contaminated surroundings [35]. The footpad injection group caused the most severe symptoms, showing significant mobility impairment, leg swelling, joint abscesses and ulcers at 10 dpi, and paralysis at 18 dpi. The N-DRV-LY20 could also cause splenic necrosis disease rapidly like the N-DRV-XT18, and severe liver necrotic spots similar to the NP03 strain [15, 34, 36]. Furthermore, diseased ducks exhibit other immune organ damage, such as a swollen bursa of Fabricius and thymic hemorrhage. Wang recently reported that the N-DRV-XT18 strain has extensive tissue orientation and may cause severe damage to immune organs [37, 38]. In addition, N-DRV can cause decreased immune function and is more likely to cause secondary infections with other pathogens [38].

There was 100% morbidity in the infected group, but mortality and weight change were significantly influenced by the route of vaccination and age at infection. The highest mortality was in the footpad injection group (30%), followed by the intramuscular injection group (20%), and the oral group (17%), which all appeared within 7 dpi different DRV strains isolated in China in recent years differed widely in the mortality rates, with rates of 10%–15% in Pekin ducks earlier discovered, rising to 40% in 2011 for the TH11 strain and up to 60% in 2013 for the variant HN5d strain [39]. The pattern of pathogenicity was consistent, but the mortality was higher compared to the N-DRV-XT18. Body weight changes indicate that ducklings infected with the N-DRV virus by any route showed weight gain suppression and poor growth, especially in the footpad-inoculated group.

The virus was detected in all tissues after inoculating the viral solution, with higher viral loads in the spleen and bursa of Fabricius, suggesting that the spleen and bursa may be our organs of choice for clinical testing of N-DRV. N-DRV destroys the spleen and bursa of Fasciola, reduces immune function, accelerates the invasion of other pathogens, and causes secondary infections. More minor ages of infection cause more severe disease reactions, and when over two weeks, it tends to be more easily tolerated and then becomes a stiff duck. Infected duck excrement contaminates feed and drinking water, which can easily cause mass transmission in the duck flock. When the disease is found in production, disinfection measures should be done, and disinfection efforts should be increased to prevent the spread of pathogens.

Overall, the N-DRV-LY20 isolate was able to cause duck viral arthritis with altered pathogenicity in Cherry Valley ducks compared to previously reported strains. Long-term evolution has caused significant variation in DRV. Detailed exploration of the genome and pathogenicity of the N-DRV-LY20 isolate provides a basis for further investigation of the pathogenesis of the newly identified N-DRV variant. To control the continued occurrence and prevalence of N-DRV, persistent monitoring of the genes and pathogenicity of the virus is needed, as well as efficient vaccines for prevention and treatment. Further studies on the mechanisms of enhanced virulence and genome reorganization of the duck orthoreovirus are also needed.

Data Availability

All data in this study are available on request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

YT, YD, and BL conceived and designed the experiments. BL, XJ, WJ, QW, and YC performed the experiments. BL, XJ, SZ, and QW analyzed the data. BL, WJ, and QW contributed to reagents/materials/analysis tools. BL wrote the paper. Bing Li and Xiaoning Jiang contributed equally to this work.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31872500 and 32172845); Natural Science Foundation of Shandong Province (ZR2020MC183); China Agriculture Research System of MOF and MARA (CARS-42-19); and The Higher Education Support Program of Youth Innovation and Technology of Shandong Province, China (2019KJF022).

Supplementary Materials

Supplementary Table S1. Main functions of proteins encoded by orthoreovirus genes. Supplementary Table S2. Primers for PCR in this study. Supplementary Table S3. GenBank accession numbers of avian orthoreoviruses for sequence alignment. Supplementary Table S4. Homology analysis between LY20 strain and other strain (%). (Supplementary Materials)