The sequence and phylogenetic analysis of avian reovirus genome segments M1, M2, and M3 encoding the minor core protein muA, the major outer capsid protein muB, and the nonstructural protein muNS

Whole Genomic Constellation of Avian Reovirus Strains Isolated from Broilers with Arthritis in North Carolina, USA

Avian reovirus (ARV) is an emerging pathogen which causes significant economic challenges to the chicken and turkey industry in the USA and globally, yet the molecular characterization of most ARV strains is restricted to a single particular gene, the sigma C gene. The genome of arthrogenic reovirus field isolates (R18-37308 and R18-38167), isolated from broiler chickens in North Carolina (NC), USA in 2018, was sequenced using long-read next-generation sequencing (NGS). The isolates were genotyped based on the amino acid sequence of sigma C (σC) followed by phylogenetic and amino acid analyses of the other 11 genomically encoded proteins for whole genomic constellation and genetic variation detection. The genomic length of the NC field strains was 23,494 bp, with 10 dsRNA segments ranging from 3959 bp (L1) to 1192 bp (S4), and the 5' and 3' untranslated regions (UTRs) of all the segments were found to be conserved. R18-37308 and R18-38167 were found to belong to genotype (G) VI based on the σC analysis and showed nucleotide and amino acid sequence identity ranging from 84.91-98.47% and 83.43-98.46%, respectively, with G VI strains. Phylogenetic analyses of individual genes of the NC strains did not define a single common ancestor among the available completely sequenced ARV strains. Nevertheless, most sequences supported the Chinese strain LY383 as a probable ancestor of these isolates. Moreover, amino acid analysis revealed multiple amino acid substitution events along the entirety of the genes, some of which were unique to each strain, which suggests significant divergence owing to the accumulation of point mutations. All genes from R18-37308 and R18-38167 were found to be clustered within genotypic clusters that included only ARVs of chicken origin, which negates the possibility of genetic pooling or host variation. Collectively, this study revealed sequence divergence between the NC field strains and reference ARV strains, including the currently used vaccine strains could help updating the vaccination regime through the inclusion of these highly divergent circulating indigenous field isolates.

First isolation and genomic characterization of avian reovirus from black swans (Cygnus atratus) in China

Identification and analysis of the avian reovirus from black swan. Isolation of the strain through the chorioallantoic membrane route of duck embryos, identified through transmission electron microscopy and RT-PCR based on the ARV S2 gene. The complete genome of the ARV strain was obtained using next-generation sequencing technology. The isolated strain of ARV was named CD200801 and was identified through transmission electron microscopy and RT-PCR based on the ARV S2 gene. Experimental infection with CD200801 resulted in the death of ducklings with serious spleen and liver focal necrosis. BLAST analysis of CD200801 sequences showed a 35.5 to 98.6% similarity to a novel duck reovirus that was isolated in recent years. Phylogenetic analysis revealed that CD200801 was closely related to ARV isolates YL, GX-Y7, and XT-18. We report the first avian reovirus infection in the black swan. This study provides important new insights into the evolutionary relationships among different ARV strains and highlights the need for continued surveillance and monitoring of these viruses in both domestic and wild bird flocks. These findings have significant implications for the development of effective strategies for disease prevention and control in the poultry industry.

Avian Reoviruses in Poultry Farms from Brazil

Avian reovirus (ARV) is highly disseminated in commercial Brazilian poultry farms, causing arthritis/tenosynovitis, runting-stunting syndrome, and malabsorption syndrome in different meat- and egg-type birds (breeders, broilers, grillers, and layers). In Brazil, ARV infection was first described in broilers in the 1970s but was not considered an important poultry health problem for decades. A more concerning outcome of field infections has been observed in recent years, including condemnations at slaughterhouses because of the unsightly appearance of chicken body parts, mainly the legs. Analyses of the performance of poultry flocks have further evidenced economic losses to farms. Genetic and antigenic characterization of ARV field strains from Brazil demonstrated a high diversity of lineages circulating in the entire country, including four of the five main phylogenetic groups previously described (I, II, III, and V). It is still unclear if all of them are associated with different diseases affecting flocks' performance in Brazilian poultry. ARV infections have been controlled in Brazilian poultry farms by immunization of breeders and young chicks with classical commercial live vaccine strains (S1133, 1733, 2408, and 2177) used elsewhere in the Western Hemisphere. However, genetic and antigenic variations of the field isolates have prevented adequate protection against associated diseases, so killed autogenous vaccines are being produced from isolates obtained on specific farms. In conclusion, ARV field variants are continuously challenging poultry farming in Brazil. Epidemiological surveillance combined with molecular biological analyses from the field samples, as well as the development of vaccine strains directed toward the ARV circulating variants, are necessary to control this economically important poultry pathogen.

Genetic characterization of a novel pheasant-origin orthoreovirus using Next-Generation Sequencing

A field isolate (Reo/SDWF /Pheasant/17608/20) of avian orthoreovirus (ARV), isolated from a flock of game-pheasants in Weifang, Shandong Province, was genetically characterized being a field variant or novel strain in our recent research studies in conducting whole genome sequencing by using Next-Generation Sequencing (NGS) technique on Illumina MiSeq platform. Among a total of 870,197 35-151-mer sequencing reads, 297,711 reads (34.21%) were identified as ARV sequences. The de novo assembly of the ARV reads resulted in generation of 10 ARV-related contigs with the average sequencing coverage from 1390× to 1977× according to 10 ARV genome segments. The complete genomes of this pheasant-origin ARV (Reo/SDWF /Pheasant/17608/20) were 23,495 bp in length and consist of 10 dsRNA segments ranged from 1192 bp (S4) to 3958 bp (L1) encoding 12 viral proteins. Sequence comparison between the SDWF17608 and classic ARV reference strains revealed that 58.1-100% nucleotide (nt) identities and 51.4-100% amino acid (aa) identities were in genome segment coding genes. The 10 RNA segments had conversed termini at 5' (5'-GCUUUU) and 3' (UCAUC-3') side, which were identical to the most published ARV strains. Phylogenetic analysis revealed that this pheasant ARV field variant was closely related with chicken ARV strains in 7 genome segment genes, but it possessed significant sequence divergence in M1, M3 and S2 segments. These findings suggested that this pheasant-origin field variant was a divergent ARV strain and was likely originated from reassortments between different chicken ARV strains.

Discovery of a novel recombinant avian orthoreovirus in China

In mid-2020, using next-generation sequencing (NGS) technology, we identified a recombinant cluster 2 avian orthoreovirus (ARV) variant named PHC-2020-0545, isolated from tendons of 33-day-old broilers with leg swelling in China. Complete genomic sequencing and analyses demonstrated that the isolate was genetically significantly distinct from known ARV strains in M1 and M3 genes and its σC coding gene had an extremely high variability, compared with the identified ARV strains grouped into other genotyping cluster. Further analysis showed that many base substitutions were silent and non-silent substitutions are most likely to occur in the first positions of codons. Multiple segmental recombination, intra-segmental recombination and accumulation of point mutations might contribute to the emergence of this isolate. The PHC-2020-0545 strain had a strong replication ability in 1-day-old broilers, and mainly affected the movement, digestion and metabolism of broilers. In addition, the infection route of the isolate is related to its pathogenicity to broilers. Therefore, combined with its unique genetic characteristics and potential origin, we determined that the PHC-2020-0545 field strain is a novel recombinant ARV strain, which has certain reference value for the preparation and evaluation of new vaccines.

Evidence of vertical transmission of novel duck orthoreovirus in ducks

Since 2017, duck spleen necrosis caused by new variant duck orthoreovirus (N-DRV) infection had been observed in several provinces in China. This disease retards the growth and development of ducks, thereby reducing feed return rate. N-DRV infection causes damage to duck spleen and other immune organs, leading to immunosuppression and susceptibility to other pathogens. In this study, we successfully constructed a breeding duck artificial infection model and found that N-DRV infection can cause pathologic changes, such as ovarian hemorrhage, follicle atrophy, and fallopian tube bleeding in breeding ducks, resulting in significantly reduced fertilization rate and egg hatching rate. Viral RNA was present in egg vitelline membrane, duck embryo, and duckling's spleen samples, as determined through quantitative polymerase chain reaction (qPCR). Autopsy revealed obvious pathologic changes in the spleen and other organs, although there were no obvious early clinical symptoms observed in ducklings. Sequence distance and phylogenetic analysis confirmed that N-DRV-SD19 re-isolated from the spleen samples of ducklings was consistent with the strain N-DRV-XT18 used for infecting breeding ducks. The findings in this study confirmed that N-DRV can be vertically transmitted through eggs, which provide an important reference for the disease prevention and control.

Genotypic characterization and molecular evolution of avian reovirus in poultry flocks from Brazil

Avian reovirus (ARV) is one of the main causes of infectious arthritis/tenosynovitis and malabsorption syndrome (MAS) in poultry. ARVs have been disseminated in Brazilian poultry flocks in the last years. This study aimed to genotype ARVs and to evaluate the molecular evolution of the more frequent ARV lineages detected in Brazilian poultry-producing farms. A total of 100 poultry flocks with clinical signs of tenosynovitis/MAS, from all Brazilian poultry-producing regions were positive for ARV by PCR. Seventeen bird tissues were submitted to cell culture and ARV RNA detection/genotyping by two PCRs. The phylogenetic classification was based on σC gene alignment using a dataset with other Brazilian and worldwide ARVs sequences. ARVs were specifically detected by both PCRs from the 17 cell cultures, and σC gene partial fragments were sequenced. All these sequences were aligned with a total of 451 ARV σC gene data available in GenBank. Phylogenetic analysis demonstrated five well-defined clusters that were classified into lineages I, II, III, IV, and V. Three lineages could be further divided into sub-lineages: I (I vaccine, Ia, Ib), II (IIa, IIb, IIc) and IV (IVa and IVb). Brazilian ARVs were from four lineages/sub-lineages: Ib (48.2%), IIb (22.2%), III (3.7%) and V (25.9%). The Bayesian analysis demonstrated that the most frequent sub-lineage Ib emerged in the world around 1968 and it was introduced into Brazil in 2010, with increasing spread soon after. In conclusion, four different ARV lineages are circulating in Brazilian poultry flocks, all associated with clinical diseases. RESEARCH HIGHLIGHTS One-hundred ARV-positive flocks were detected in all main poultry-producing regions from Brazil. A large dataset of 468 S1 sequences was constructed and divided ARVs into five lineages. Four lineages/sub-lineages (Ib, IIb, III and V) were detected in commercial poultry flocks from Brazil. Brazilian lineages shared a low identity with the commercial vaccine lineage (I vaccine). Sub-lineage Ib emerged around 1968 and was introduced into Brazil in 2010.

The genomic constellation of a novel duck reovirus strain associated with hemorrhagic necrotizing hepatitis and splenitis in Muscovy ducklings in Fujian, China

The complete sequence of a reovirus, strain NP03 associated with necrotic focus formation in the liver and spleen of Muscovy ducklings in Fujian Province, China in 2009, was determined and compared with sequences of other waterfowl and chicken-origin avian reoviruses (ARVs). Sequencing of the complete genomes of strain NP03 showed that they consisted of 23,418 bp and were divided into 10 segments, ranging from 1191 bp (S4) to 3959 bp (L1) in length, and all segments contained conserved sequences in the 5' non-coding region (GCUUUU) and 3' non-coding region (UCAUC). Pairwise sequence comparisons demonstrated that NP03 strain showed the highest similarity with novel waterfowl origin reoviruses (WRVs). The genome analysis revealed that the S1 segment of novel WRV is a tricistronic gene, encoding the overlapping open reading frames (ORFs) for p10, p18, and σC, similar to the ARV S1 gene, but distinct from classical WRV S4 genome segment, which contained two overlapping ORFs encoding p10 and σC. Phylogenetic analyses of the nucleotide sequences of all 10 segments revealed that NP03 strain was clustered together with other novel WRVs and were distinct from classical WRVs and chicken-origin ARVs. The analyses also showed possible intra-segmental reassortment events in the segments encoding λA, λB, μB, μNS, σA, and σNS between novel and classical WRVs. Potential recombination events detection in segment L1 suggests that NP03 strain may be recombinants of novel WRVs. Based on our genetic analyses, multiple reassortment events, intra-segmental recombination, and accumulation of point mutations have possibly contributed to the emergence of this novel genotype of WRV, identified in China.

Partial Molecular Characterization and Pathogenicity Study of an Avian Reovirus Causing Tenosynovitis in Commercial Broilers

This study describes the molecular characterization of avian reoviruses (ARVs) isolated during an outbreak in commercial chickens between 2015 and 2016. In addition, a pathogenicity study of a selected ARV strain isolated from a field case of viral tenosynovitis in commercial broiler chickens was performed. On the basis of phylogenetic analysis of a 1088-bp fragment of the ARV S1 gene, the investigated sequences were differentiated into five distinct genotypic clusters (GCs), namely GC1, GC2, GC3, GC4, and GC6. Specific-pathogen-free (SPF) and commercial broiler chickens were challenged with the GC1 genetic type MK247011, at 14 days of age via the interdigital toe web. No significant effects in body weight gain and feed conversion were detected in both chicken types. The Δ interdigital web thickness was most severe at 4 days postchallenge (DPC) in both the SPF and broiler subgroups. The inflammation in SPF birds was slightly more severe compared with broilers. Neither mortality nor clinical signs occurred in the infected groups for the duration of the experiment, despite the presence of significant microscopic lesions in challenged birds. Microscopic changes of tenosynovitis became evident at 3 DPC, with the highest incidence and severity detected at 14 and 21 DPC, respectively. Seroconversion against ARV occurred 3 wk postchallenge, and the microscopic lesions detected in tendon and heart sections were highly compatible with those described in the field. Increased severity of tenosynovitis and epicarditis lesions were noted in the ARV-challenged groups compared with the control groups. Although SPF and broiler chickens showed comparable responses to the challenge with an ARV genetic variant, detected lesions were subclinical, denoting the limitations of our challenge approach. The age selected in this experiment possibly influenced the course of the infection. Data from this study highlight the genotypic diversity of isolates in California, and the outcome of the pathogenicity study can be used as a basis to improve protocols for pathogenicity studies to characterize ARV variants causing clinical disease in the field.

Isolation and characterization of a variant duck orthoreovirus causing spleen necrosis in Peking ducks, China

Since December 2017, an infectious disease has caused economic hardship for duck farms and breeding ducks in many regions of China. This disease characterized by spleen necrosis and swelling, is due to a variant strain of duck orthoreovirus (DRV) (Duck/N-DRV-XT18/China/2018), which we isolated from the spleen of diseased ducks. After isolating the virus, we used next-generation sequencing technology to determine the entire genomic of the virus. Our phylogenetic analysis of 10 genomic segments showed that the N-DRV-XT18 strain is closely related to orthoreovirus isolates derived from ducks and geese, with nucleotide sequence identities for 10 genomic fragments ranging between 49.8% and 99.3%. In contract, the nucleotide sequence of N-DRV-XT18 genomic fragments are only 38.6% to 78.8% similar to the chicken orthoreovirus isolate. Therefore, we determined that this pathogen, causing duck spleen necrosis, is a new variant of a duck orthoreovirus that is significantly different from any previously reported waterfowl-derived othoreovirus.

Genotypic Characterization of Emerging Avian Reovirus Genetic Variants in California

This study focuses on virus isolation of avian reoviruses from a tenosynovitis outbreak between September 2015 and June 2018, the molecular characterization of selected isolates based on partial S1 gene sequences, and the full genome characterization of seven isolates. A total of 265 reoviruses were detected and isolated, 83.3% from tendons and joints, 12.3% from the heart and 3.7% from intestines. Eighty five out of the 150 (56.6%) selected viruses for sequencing and characterization were successfully detected, amplified and sequenced. The characterized reoviruses grouped in six distinct genotypic clusters (GC1 to GC6). The most represented clusters were GC1 (51.8%) and GC6 (24.7%), followed by GC2 (12.9%) and GC4 (7.2%), and less frequent GC5 (2.4%) and GC3 (1.2%). A shift on cluster representation throughout time occurred. A reduction of GC1 and an increase of GC6 classified strains was noticed. The highest homologies to S1133 reovirus strain were detected in GC1 (~77%) while GC2 to GC6 homologies ranged between 58.5 and 54.1%. Over time these homologies have been maintained. Seven selected isolates were full genome sequenced. Results indicated that the L3, S1 and M2 genes, coding for proteins located in the virus capsid accounted for most of the variability of these viruses. The information generated in the present study helps the understanding of the epidemiology of reoviruses in California. In addition, provides insights on how other genes that are not commonly studied add variability to the reovirus genome.

Molecular characterization of two novel reoviruses isolated from Muscovy ducklings in Guangdong, China

Background

Novel Muscovy duck reovirus (N-MDRV), emerged in southeast China in 2002, which can infect a wide range of waterfowl and induces clinical signs and cytopathic effects that are distinct from those of classical MDRV, and continues to cause high morbidity and 5–50% mortality in ducklings. The present study aimed to investigate the characteristics of two novel reoviruses isolated from Muscovy ducklings in Guangdong, China.

Results

Two novel MDRV strains, designated as MDRV-SH12 and MDRV-DH13, were isolated from two diseased Muscovy ducklings in Guangdong province, China in June 2012 and September 2013, respectively. Sequencing of the complete genomes of these two viruses showed that they consisted of 23,418 bp and were divided into 10 segments, ranging from 1191 bp (S4) to 3959 bp (L1) in length, and all segments contained conserved sequences in the 5′ non-coding region (GCUUUU) and 3′ non-coding region (UCAUC). Pairwise sequence comparisons demonstrated that MDRV-SH12 and MDRV-DH13 showed the highest similarity with novel MDRVs. Phylogenetic analyses of the nucleotide sequences of all 10 segments revealed that MDRV-SH12 and MDRV-DH13 were clustered together with other novel waterfowl-origin reoviruses and were distinct from classical waterfowl-origin and chicken-origin reoviruses. The analyses also showed possible genetic re-assortment events in segment M2 between waterfowl-origin and chicken-origin reoviruses and the segments encoding λA, μA, μNS, σA, and σNS between classical and novel waterfowl-origin reoviruses. Potential recombination events detection in segment S2 suggests that MDRV-SH12 and MDRV-DH13 may be recombinants of classical and novel WRVs.

Conclusions

The results presented in this study, the full genomic data for two novel MDRV strains, will improve our understanding of the evolutionary relationships among the waterfowl-origin reoviruses circulating in China, and may aid in the development of more effective vaccines against various waterfowl-origin reoviruses.

Electronic supplementary material

The online version of this article (10.1186/s12917-019-1877-x) contains supplementary material, which is available to authorized users.

Complete genome sequence of a novel avian orthoreovirus isolated from gosling, China

Avian orthoreovirus (ARV) has been considered as a significant pathogen causing great infectious diseases to the avian, like broiler and waterfowl. The genome of this novel ARV(Reo/SDPY/Goose) was completely sequenced by next-generation sequencing. The complete genome was found to be 23517 bp in length with 10 segments. Although the Reo/SDPY/Goose was isolated from the gosling, it shares great similarity, no matter which segment within the genome, with those published as avian-origin reovirus. Genomic analysis revealed that this virus was distinct from published ARV strains and met criteria to become a novel ARV strain.

Dissection of mammalian orthoreovirus µ2 reveals a self-associative domain required for binding to microtubules but not to factory matrix protein µNS

Mammalian orthoreovirus protein μ2 is a component of the viral core particle. Its activities include RNA binding and hydrolysis of the γ-phosphate from NTPs and RNA 5´-termini, suggesting roles as a cofactor for the viral RNA-dependent RNA polymerase, λ3, first enzyme in 5´-capping of viral plus-strand RNAs, and/or prohibitory of RNA-5´-triphosphate-activated antiviral signaling. Within infected cells, μ2 also contributes to viral factories, cytoplasmic structures in which genome replication and particle assembly occur. By associating with both microtubules (MTs) and viral factory matrix protein μNS, μ2 can anchor the factories to MTs, the full effects of which remain unknown. In this study, a protease-hypersensitive region allowed μ2 to be dissected into two large fragments corresponding to residues 1-282 and 283-736. Fusions with enhanced green fluorescent protein revealed that these amino- and carboxyl-terminal regions of μ2 associate in cells with either MTs or μNS, respectively. More exhaustive deletion analysis defined μ2 residues 1-325 as the minimal contiguous region that associates with MTs in the absence of the self-associating tag. A region involved in μ2 self-association was mapped to residues 283-325, and self-association involving this region was essential for MT-association as well. Likewise, we mapped that μNS-binding site in μ2 relates to residues 290-453 which is independent of μ2 self-association. These findings suggest that μ2 monomers or oligomers can bind to MTs and μNS, but that self-association involving μ2 residues 283-325 is specifically relevant for MT-association during viral factories formation.

Detection and characterization of two co-infection variant strains of avian orthoreovirus (ARV) in young layer chickens using next-generation sequencing (NGS)

Using next-generation sequencing (NGS) for full genomic characterization studies of the newly emerging avian orthoreovirus (ARV) field strains isolated in Pennsylvania poultry, we identified two co-infection ARV variant strains from one ARV isolate obtained from ARV-affected young layer chickens. The de novo assembly of the ARV reads generated 19 contigs of two different ARV variant strains according to 10 genome segments of each ARV strain. The two variants had the same M2 segment. The complete genomes of each of the two variant strains were 23,493 bp in length, and 10 dsRNA segments ranged from 1192 bp (S4) to 3958 bp (L1), encoding 12 viral proteins. Sequence comparison of nucleotide (nt) and amino acid (aa) sequences of all 10 genome segments revealed 58.1–100% and 51.4–100% aa identity between the two variant strains, and 54.3–89.4% and 49.5–98.1% aa identity between the two variants and classic vaccine strains. Phylogenetic analysis revealed a moderate to significant nt sequence divergence between the two variant and ARV reference strains. These findings have demonstrated the first naturally occurring co-infection of two ARV variants in commercial young layer chickens, providing scientific evidence that multiple ARV strains can be simultaneously present in one host species of chickens.

Avian reovirus S1133-induced apoptosis is associated with Bip/GRP79-mediated Bim translocation to the endoplasmic reticulum

In this study the mechanism of avian reovirus (ARV) S1133-induced pathogenesis was investigated, with a focus on the contribution of ER stress to apoptosis. Our results showed that upregulation of the ER stress response protein, as well as caspase-3 activation, occurred in ARV S1133-infected cultured cells and in SPF White Leghorn chicks organs. Upon infection, Bim was translocated specifically to the ER, but not mitochondria, in the middle to late infectious stages. In addition, ARV S1133 induced JNK phosphorylation and promoted JNK-Bim complex formation, which correlated with the Bim translocation and apoptosis induction that was observed at the same time point. Knockdown of BiP/GRP78 by siRNA and inhibition of BiP/GRP78 using EGCG both abolished the formation of the JNK-Bim complex, caspase-3 activation, and subsequent apoptosis induction by ARV S1133 efficiently. These results suggest that BiP/GRP78 played critical roles and works upstream of JNK-Bim in response to the ARV S1133-mediated apoptosis process.

Critical role of eukaryotic elongation factor 1 alpha 1 (EEF1A1) in avian reovirus sigma-C-induced apoptosis and inhibition of viral growth

Avian reovirus (ARV) causes viral arthritis, chronic respiratory diseases, retarded growth and malabsorption syndrome. It is well established that the ARV sigma-C protein induces apoptosis in host cells. However, the underlying molecular mechanism of this induction is still unclear. We report here the identification of eukaryotic elongation factor 1 alpha 1 (EEF1A1) as the interacting partner of σC. We found that σC-induced apoptosis in DF-1 cells could be completely abolished by knockdown of EEF1A1 by siRNA. Furthermore, knockdown of EEF1A1 markedly reduced ARV-induced apoptosis associated with decreased caspase-9 and -3 activation and cytochrome C release, leading to increased ARV growth in host cells. Thus, EEF1A1 plays a critical role in σC-induced apoptosis and inhibition of viral growth.

Genomic characterization of a turkey reovirus field strain by Next-Generation Sequencing

The genome of a turkey arthritis reovirus (TARV) field strain (Reo/PA/Turkey/22342/13), isolated from a turkey flock in Pennsylvania (PA) in 2013, has been sequenced using Next-Generation Sequencing (NGS) on the Illumina MiSeq platform. The genome of the PA TARV field strain was 23,496bp in length with 10 dsRNA segments encoding 12 viral proteins. The lengths of the genomic segments ranged from 1192bp (S4) to 3959bp (L1). The 5' and 3' conserved terminal sequences of the PA TARV field strain were similar to the two Minnesota (MN) TARVs (MN9 and MN10) published recently and avian orthoreovirus (ARV) reference strains. Phylogenetic analysis of the nucleotide sequences of all 10 genome segments revealed that there was a low to significant nucleotide sequence divergence between the PA TARV field strain and reference TARV and ARV strains. Analysis of the PA TARV sequence indicates that this PA TARV field strain is a unique strain and is different from the TARV MN9 or MN10 in M2 segment genes and ARV S1133 vaccine strain.

Phylogenetic analysis, genomic diversity and classification of M class gene segments of turkey reoviruses

From 2011 to 2014, 13 turkey arthritis reoviruses (TARVs) were isolated from cases of swollen hock joints in 2-18-week-old turkeys. In addition, two isolates from similar cases of turkey arthritis were received from another laboratory. Eight turkey enteric reoviruses (TERVs) isolated from fecal samples of turkeys were also used for comparison. The aims of this study were to characterize turkey reovirus (TRV) based on complete M class genome segments and to determine genetic diversity within TARVs in comparison to TERVs and chicken reoviruses (CRVs). Nucleotide (nt) cut off values of 84%, 83% and 85% for the M1, M2 and M3 gene segments were proposed and used for genotype classification, generating 5, 7, and 3 genotypes, respectively. Using these nt cut off values, we propose M class genotype constellations (GCs) for avian reoviruses. Of the seven GCs, GC1 and GC3 were shared between the TARVs and TERVs, indicating possible reassortment between turkey and chicken reoviruses. The TARVs and TERVs were divided into three GCs, and GC2 was unique to TARVs and TERVs. The proposed new GC approach should be useful in identifying reassortant viruses, which may ultimately be used in the design of a universal vaccine against both chicken and turkey reoviruses.

Characterization of a novel orthoreovirus isolated from fruit bat, China

Background

In recent years novel human respiratory disease agents have been described for Southeast Asia and Australia. The causative pathogens were classified as pteropine orthoreoviruses with a strong phylogenetic relationship to orthoreoviruses of bat origin.

Results

In this report, we isolated a novel Melaka-like reovirus (named “Cangyuan virus”) from intestinal content samples of one fruit bat residing in China’s Yunnan province. Phylogenetic analysis of the whole Cangyuan virus genome sequences of segments L, M and S demonstrated the genetic diversity of the Cangyuan virus. In contrast to the L and M segments, the phylogenetic trees for the S segments of Cangyuan virus demonstrated a greater degree of heterogeneity.

Conclusions

Phylogenetic analysis indicated that the Cangyuan virus was a novel orthoreovirus and substantially different from currently known members of Pteropine orthoreovirus (PRV) species group.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0293-4) contains supplementary material, which is available to authorized users.

Sequencing and phylogenetic analysis of an avian reovirus genome

Avian reovirus infection causes considerable economic loss to the commercial poultry industry. Live-attenuated vaccine strain S1133 (v-S1133, derived from parent strain S1133) is considered the safest and most effective vaccine and is currently used worldwide. To identify the genes responsible for its attenuation, DNA sequences of open reading frames (ORF) of S1133 and its parent strains S1133, 1733, 526, and C78 along with three field isolates (GuangxiR1, GuangxiR2, and GX110058) and one isolate (GX110116) from a vaccinated chicken were performed. The sequence data were compared with available sequences in nucleotide sequence databases of American (AVS-B, 138, 176) and Chinese (C-98 and T-98) origin. Sequence analysis identified that several v-S1133 specific nucleotide substitutions existed in the ORFs of λA, λB, λC, μA, μB, μNS, σA, σB, and σNS genes. The v-S1133 strain could be differentiated from the field-isolated strains based on single nucleotide polymorphisms. Phylogenetic analysis revealed that v-S1133 shared the highest sequence homologies with S1133 and reovirus isolates from China, grouped together in one cluster. Chinese isolates were clearly more distinct from the American reovirus AVS-B strain, which is associated with runting-stunting syndrome in broilers.

Identification of a linear B-Cell epitope on avian reovirus protein sigmaC

SigmaC (σC) protein, which mediates virus attachment to target cells, is the most variable proteins of avian reovirus (ARV). It is responsible for inducing protective antibody immune responses in animals. To understand the antigenic determinants of σC protein, a set of partially overlapping and consecutive peptides spanning σC were expressed and then screened with the monoclonal antibody (mAb) 2B5 directed against σC. The mAb 2B5 recognized peptides with the σC motif (45)ELLHRSISDISTTV(58). Further identification of the displayed B-cell epitope was conducted with a set of truncated peptides expressed as GST fusion proteins. The Western blot and ELISA results indicated that (45)ELLHRSISDI(54) was the minimal determinant of the linear B-cell epitope. Using sequences analysis, we found that this epitope was not a common motif shared among the other members of the ARV and DRV groups. Furthermore, cross reactivity analysis showed that the associated coding motif of other ARV and DRV groups was not recognized by 2B5. These data suggested that (45)ELLHRSISDI(54) was a type-specific linear B-cell epitope of avian reovirus. The results in this study may have potential applications in the development of diagnostic techniques and epitope-based marker vaccines against ARV, which is prevalent in China.

Complete sequence of a reovirus associated with necrotic focus formation in the liver and spleen of Muscovy ducklings

The complete sequence of a reovirus, strain 815-12 associated with necrotic focus formation in the liver and spleen of Muscovy ducklings in China, was determined and compared with sequences of other duck-, goose-, and chicken-origin reoviruses. The 815-12 genome comprised 22,969 bp with 10 dsRNA segments ranging from 1125 bp (S4) to 3958 bp (L1), all of which (except S4) were almost identical in length to the cognate segments of other waterfowl and chicken isolates. Detailed analyses revealed that 815-12 and other waterfowl isolates contained the conserved 3'-terminal pentanucleotide sequence (UCAUC-3') of the orthoreoviruses and 5'-terminal hexanucleotide sequence (5'-GCUUUU) of avian orthoreoviruses (ARVs), and conserved functional motifs previously identified in ARV proteins. Several notable differences, including organization of the polycistronic genome segments and genomic coding assignments of the S segments, existed between viruses represented by 815-12 and the waterfowl reoviruses emerging in China in recent years; the latter was somewhat similar to chicken isolates. Pairwise sequence comparisons demonstrated extensive sequence diversity among the various waterfowl isolates and between waterfowl and chicken isolates. Phylogenetic analyses identified two genetic groups for waterfowl reoviruses, and potential genetic reassortment of segment M2 between waterfowl and chicken reoviruses and segments encoding for λA, λB, μA, μNS and σA between waterfowl reoviruses. Taken together, it was suggested that common designation ARV-Wa should be used to represent ARV isolates from different waterfowl species and that the two ARV-Wa genotypes should be considered as two separate groups distinct from chicken isolates within the species Avian orthoreovirus.

Identification of a new reovirus causing substantial losses in broiler production in France, despite routine vaccination of breeders

Numerous cases of tenosynovitis appeared in France causing high morbidity in free-range and standard broilers. The main clinical findings were lameness, stunting and non-uniform bodyweights. Although the natural mortality was low, the economic losses due to birds that had to be removed from the flock prematurely, downgrading of carcases and lower average weights at slaughter were substantial. Postmortem examinations, bacteriological, virological and serological examination confirmed the aetiology of avian orthoreovirus (ARV)-induced tenosynovitis. The isolated ARVs were analysed serologically and genetically. Sequencing of σC RT-PCR products and phylogenetic analysis revealed a new type of ARV. The virus was not neutralised in serum neutralisation test using monovalent sera from vaccinated chickens. Together with the flock data, epidemiology of these recent reovirus outbreaks in France was reconstructed. It is concluded that these reovirus isolates differ serologically and genetically from the well described reovirus isolates used in commercial vaccines which were not capable of preventing the disease. The outbreaks resulted in substantial losses in broilers from vaccinated breeders.

Evolutionary relationship of the L- and M-class genome segments of bat-borne fusogenic orthoreoviruses in Malaysia and Australia

We previously described three new Malaysian orthoreoviruses designated Pulau virus, Melaka virus and Kampar virus. Melaka and Kampar viruses were shown to cause respiratory disease in humans. These viruses, together with Nelson Bay virus, isolated from Australian bats, are tentatively classified as different strains within the species Pteropine orthoreovirus (PRV), formerly known as Nelson Bay orthoreovirus, based on the small (S) genome segments. Here we report the sequences of the large (L) and medium (M) segments, thus completing the whole-genome characterization of the four PRVs. All L and M segments were highly conserved in size and sequence. Conserved functional motifs previously identified in other orthoreovirus gene products were also found in the deduced proteins encoded by the cognate segments of these viruses. Detailed sequence analysis identified two genetic lineages divided into the Australian and Malaysian PRVs, and potential genetic reassortment among the M and S segments of the three Malaysian viruses.

The genomic constellation of a novel avian orthoreovirus strain associated with runting-stunting syndrome in broilers

Avian orthoreoviruses (ARVs) are responsible for considerable economic losses in broiler chickens; yet, the genetic characterization of most ARV strains is limited to a few genes, and the full coding region has been determined for only S1133 and 138, two ARV strains associated with tenosynovitis. Recently, in parts of the United States, ARVs with novel neutralization antigen type were isolated from chickens afflicted with runting-stunting syndrome. One such strain, AVS-B, was selected for full genome sequencing and phylogenetic analysis. The complete genome was 23,494 bp in size and included 12 open reading frames. The lengths of the coding regions, as well as those of the 5' and 3' ends, were fairly well conserved between AVS-B and other reference strains. In pairwise comparisons to the S1133 and 138 strains, the AVS-B strain shared a wide range of sequence identities along each genome segment, i.e., a range of 54-55% for the σC coding region of S1 genome segment and 91-93% for the S2 genome segment. Phylogenetic analyses of individual genes of AVS-B did not identify any single common ancestor among more completely characterized ARV strains for which sequence data are available. One exception to this lack of identity was strain 138, which shared 90-93% nt identity with AVS-B along seven of ten genome segments; only M2, M3, and S1 segments of these strains shared lower sequence identities. Collectively, our analyses indicated that multiple reassortment events and strong divergence caused by the accumulation of point mutations could have led to the observed assortment and genetic heterogeneity of the AVS-B genome.

Characterization of interleukin-1beta mRNA expression in chicken macrophages in response to avian reovirus

Inhibitors of viral disassembly or RNA and protein synthesis, viral disassembly intermediates (infectious subviral particles, ISVP), binary ethylenimine-inactivated virions, and viral particles lacking genomic double-stranded (ds) RNA (empty particles) were used to assess the expression of interleukin-1beta (IL-1beta) mRNA in chicken (chIL-1beta) macrophages in response to avian reovirus. The results demonstrate that two distinct expression patterns of chIL-1beta mRNA mediated by different steps in viral replication were found. Viral disassembly was required for the induction of a rapid, transient expression pattern of chIL-1beta mRNA that was rapidly induced at 30 min, with maximal levels reached by 2 h, and fell to a low level within 6 h post-inoculation, while viral RNA synthesis rather than protein translation, which was subsequent to membrane penetration, was required to induce a stable, sustained expression pattern of chIL-1beta mRNA that occurred at and after 6 h post-inoculation. In addition, the induction of chIL-1beta mRNA expression by the empty particles and ISVP was extremely weak, compared with the active dsRNA(+) virions or binary ethylenimine-inactivated virions, suggesting that the presence of dsRNA, even if transcriptionally inactive, may be an important factor in this response.

Formation of the factory matrix is an important, though not a sufficient function of nonstructural protein mu NS during reovirus infection

Genome replication of mammalian orthoreovirus (MRV) occurs in cytoplasmic inclusion bodies called viral factories. Nonstructural protein microNS, encoded by genome segment M3, is a major constituent of these structures. When expressed without other viral proteins, microNS forms cytoplasmic inclusions morphologically similar to factories, suggesting a role for microNS as the factory framework or matrix. In addition, most other MRV proteins, including all five core proteins (lambda1, lambda2, lambda3, micro2, and sigma2) and nonstructural protein sigmaNS, can associate with microNS in these structures. In the current study, small interfering RNA targeting M3 was transfected in association with MRV infection and shown to cause a substantial reduction in microNS expression as well as, among other effects, a reduction in infectious yields by as much as 4 log(10) values. By also transfecting in vitro-transcribed M3 plus-strand RNA containing silent mutations that render it resistant to the small interfering RNA, we were able to complement microNS expression and to rescue infectious yields by ~100-fold. We next used microNS mutants specifically defective at forming factory-matrix structures to show that this function of microNS is important for MRV growth; point mutations in a C-proximal, putative zinc-hook motif as well as small deletions at the extreme C terminus of microNS prevented rescue of viral growth while causing microNS to be diffusely distributed in cells. We furthermore confirmed that an N-terminally truncated form of microNS, designed to represent microNSC and still able to form factory-matrix structures, is unable to rescue MRV growth, localizing one or more other important functions to an N-terminal region of microNS known to be involved in both micro2 and sigmaNS association. Thus, factory-matrix formation is an important, though not a sufficient function of microNS during MRV infection; microNS is multifunctional in the course of viral growth.

Avian reovirus core protein μA expressed in Escherichia coli possesses both NTPase and RTPase activities

Analysis of the amino acid sequence of core protein μA of avian reovirus has indicated that it may share similar functions to protein μ2 of mammalian reovirus. Since μ2 displayed both nucleotide triphosphatase (NTPase) and RNA triphosphatase (RTPase) activities, the purified recombinant μA ( μA) was designed and used to test these activities. μA was thus expressed in bacteria with a 4.5 kDa fusion peptide and six His tags at its N terminus. Results indicated that  μA possessed NTPase activity that enabled the protein to hydrolyse theβ–γphosphoanhydride bond of all four NTPs, since NDPs were the only radiolabelled products observed. The substrate preference was ATP>CTP>GTP>UTP, based on the estimatedkcatvalues. Alanine substitutions for lysines 408 and 412 (K408A/K412A) in a putative nucleotide-binding site of  μA abolished NTPase activity, further suggesting that NTPase activity is attributable to protein  μA. The activity of  μA is dependent on the divalent cations Mg2+or Mn2+, but not Ca2+or Zn2+. Optimal NTPase activity of  μA was achieved between pH 5.5 and 6.0. In addition,  μA enzymic activity increased with temperature up to 40 °C and was almost totally inhibited at temperatures higher than 55 °C. Tests of phosphate release from RNA substrates with  μA or K408A/K412A  μA indicated that  μA, but not K408A/K412A  μA, displayed RTPase activity. The results suggested that both NTPase and RTPase activities of  μA might be carried out at the same active site, and that protein μA could play important roles during viral RNA synthesis.

Genetic variation of the lambdaA and lambdaC protein encoding genes of avian reoviruses

Sequence and phylogenetic analysis of lambdaA and lambdaC protein encoding genes of 12 avian reoviruses is described. The sequence of lambdaA possesses a variable region (residues 19-51) located within a conserved hydrophilic region (residues 1-110) and a C(2)H(2) zinc-binding motif (residues 182-202). lambdaC shows the two conserved K residues at positions 169 and 188 indicative of guanylyltransferase activity, an ATP/GTP-binding site motif A (residues 379-386), and a conserved S-adenosyl-l-methionine-binding motif (residues 822-830). Pairwise sequence comparisons show that the mean sequence identities of lambdaA encoding genes and lambdaA proteins are 92% and 98%, respectively, and those of lambdaC encoding genes and lambdaC proteins are 91% and 95%, respectively. Phylogenetic analysis of lambdaA and lambdaC encoding genes reveals that both encoding genes have diverged into three distinct lineages, respectively, and that there is no correlation between lineages and viral serotypes or pathotypes. Also, reassortment of gene segments L1 and L3 has been observed between viruses.

Characterization of M-class genome segments of muscovy duck reovirus S14

This report documents the first sequence analysis of the entire M1, M2, and M3 genome segments of the muscovy duck reovirus (DRV) S14. The complete sequence of each of the three M gene segments was determined. The M1 genome segment was 2283 nucleotides in length and was predicted to encode muA protein of 732 residues. The Escherichia coli expressed M1 transcripts generated a 108kDa protein, as expected for muA. A cleavage product of muA, muA1, could be detected by Western blotting with duck anti-reovirus and mouse anti-muA polyclonal serum. muA was distributed diffusely in the cytoplasma and nucleus of transfected Vero cells, which provides evidence that muA might be functional related to the mammalian reovirus (MRV) mu2. The M2 gene was 2155 nucleotides in length and was predicted to encode muB major outer capsid protein of 676 amino acids. The M3 genome segment was 1996 nucleotides in length and was predicted to encode a muNS protein of 635 amino acids. It was unexpectedly found that 5'-termini of the M1 and M2 genes ended with 5'-ACUUUU and 5'-UCUUUU, respectively, instead of 5'-GCUUUU, which is present on most mRNAs of other avian reoviruses (ARV). The UCAUC 3'-terminal sequences of the S14 M1, M2, and M3 genome segments are shared by DRV, ARV, and MRV. Alignment of the DRV muA-, muB-, and muNS-encoding genes with ARV revealed 72.9-73.9%, 67.1-69.6%, and 69.4-70.8% nucleotide identity, respectively. The amino acid sequence homology between DRV and ARV ranged from 85.3 to 86.2% (muA), 75.0 to 76.5% (muB), and 78.4 to 79.8% (muNS). Phylogenetic analyses of the M1, M2, M3, and S-class [Kuntz-Simon, G., Le Gall-Recule, G., de Boisseson, C., Jestin, V., 2002. Muscovy duck reovirus sigmaC protein is a typically encoded by the smallest genome segment. J. Gen. Virol. 83, 1189-1200; Zhang, Y., Liu, M., Hu, Q.L., Ouyang, S.D., Tong, G.Z., 2006a. Characterization of the sigmaC-encoding gene from muscovy duck reovirus. Virus Genes 36, 169-174; Zhang, Y., Liu, M., Ouyan, S.D., Hu, Q.L., Guo, D.C., Han, Z., 2006b. Detection and identification of avian, duck, and goose reoviruses by RT-PCR: goose and duck reoviruses aggregated the same specified genogroup in Orthoreovirus Genus II. Arch. Virol. 151, 1525-1538] genome segments suggests that DRV and ARV share a recent common ancestor and that the two lineages have subsequently undergone host dependent evolution.