Characterization and pathogenicity of a novel variant infectious bursal disease virus in China

Infectious bursal disease (IBD) is a highly epidemic and immunosuppressive disease of 3- to 6-week-old chicks caused by infectious bursal disease virus (IBDV). Since 2017, there has been a notable increase in the isolation rates of novel variant IBDV strains in China, of which characteristic amino acid residues were different from those of early antigen variants. In this study, one IBDV strain was isolated from a farm with suspected IBD outbreak in Shandong Province, China, which was designated LY21/2. The strain LY21/2 could replicate in MC38 cells with previous culture adaption in SPF chick embryos. Phylogenetic analysis revealed that LY21/2 formed one branch with novel variant IBDVs and shared 96.8–98.6% nucleotide sequence identity with them. Moreover, LY21/2 serving as the major parent underwent the recombination event of a variant strain (19D69), while the minor parent was a very virulent strain (Harbin-1). SPF chicks inoculated with LY21/2 showed no gross clinic symptom, whereas bursal atrophy was exhibited and apoptosis was occurred in 55.21% of bursal cells. The results of histopathology and immunohistochemical staining showed that lymphocyte depletion and connective tissue hyperplasia and IBDV antigen-positive cells were observed in the bursa of LY21/2-infected chicks. Besides, DNA fragmentation was detected in the LY21/2-infected bursal tissue section by TUNEL assay. Collectivtely, these data presented analysis and evaluation of the genetic characteristics and pathogenicity of a novel variant IBDV strain. This study may help in the development of biosafety strategies for the prevention and control of IBDV in poultry.

Introduction

Infectious bursal disease (IBD) is caused by infectious bursal disease virus (IBDV), which is an immunosuppressive disease of young chicks. IBDV has two serotypes, including serotype I strains which consisted of classic strains, attenuated strains, variant strains, very virulent strains (vvIBDV), and serotype II strains. The genome of IBDV, a member of the family Birnaviridae and the genus Avibirnavirus (Azad et al., 1985), consists of two segments of double-strained RNA (Muller et al., 1979). The larger fragment A encodes a pVP2-VP3-VP4 polyprotein and a nonstructural protein VP5 in two partially overlapping reading frames, and the smaller fragment B encodes VP1 (Mahgoub et al., 2012). VP2 protein is the major structural protein related to virulence, cell tropism and antigenic variation (Brandt et al., 2001; Coulibaly et al., 2005; Letzel et al., 2007). The antigenic phenotype of IBDV is determined by the hypervariable region (HVR) of VP2 gene (Bayliss et al., 1990; Eterradossi et al., 1997, 1998), in which contains the amino acid (aa) residues 206–350.

In 1957, IBDV was first isolated in the United States (Muller et al., 2003) which was referred to as classic strain, and since then, the virus has spread worldwide and undergone a complex history of evolution. With the emergence of antigenic variant (Jackwood and Saif, 1987) and very virulent IBDVs (Chettle et al., 1989), new challenges on outbreak control and vaccine development have been come as follow. Since IBD first reported in 1979 in southern China (Cui et al., 2013), IBDV mainly of attenuated and very virulent strains infection became prevalent in several provinces and caused huge economic losses (He et al., 2014, 2019). Depending on the virulence, serotype I strains vary in clinical signs and mortality. Some of variant strains may not cause clinical signs or death but do cause bursal lesions. It might account for lack of monitoring and few studies reported about variant IBDVs. However, variant strains were reported and isolated increasingly after outbreaking since 2017 (Fan et al., 2019; Chen et al., 2022; Lian et al., 2022). As is known, IBDV causes immunosuppression mostly by attacking the bursa of Fabricius (BFs) particularly the developing B lymphocytes (Sharma et al., 2000). Moreover, IBDV affects the production of antibodies induced by B cells leading to increased susceptibility to other diseases and vaccination failure (Ranjbar et al., 2019; Fan et al., 2020). Therefore, in consideration of worldwide spread and immunosuppression induced by antigenic variation, it is of significance for preventing and controlling IBD to determine the dominant strains and the critical amino acid residues that may cause immune failure.

In this study, to explore the phenomenon of immunity failure of IBDV vaccine in chicks, we isolated and identified a novel variant IBDV strain (designated LY21/2) from an immunized chick flock with suspected IBD outbreak. The molecular characteristics and pathogenicity were analyzed and evaluated, respectively, by bioinformatics method and by conducting experiment on specific-pathogen-free (SPF) chicks. Our results can provide a reference for epidemic analysis of IBDV in China and studies on immune escape.

Materials and methods

Clinical samples

One clinical sample of bursa was collected from broiler in an immunized farm of Linyi, Shandong Province, China. The bursal tissue was homogenized and prepared to 10% (w/v) suspension in phosphate buffered saline (PBS, pH 7.2) with penicillin (100 U/ml) and streptomycin (0.1 mg/ml). The suspension was frozen and thawed three times followed by centrifugation at 10,000 × g for 10 min at 4°C, and then the supernatant was harvested for subsequent detection.

Animals, cells, and antibodies

SPF chick embryos were purchased from Lihua Agricultural Technology Co., Ltd. (Zhejiang, China) and were incubated at 37°C. Chicks were hatched from the embryos and housed in negative-pressure-filtered air isolators. All the animal experiments were approved by the Shanghai Veterinary Research Institute of the Chinese Academy of Agricultural Sciences and were performed in accordance with the animal ethics guidelines and approved protocols. The mouse colon cancer cell line MC38 was maintained in Dulbecco’s Modified Eagle’s medium (DMEM; BI, Israel) containing 10% foetal bovine serum (FBS; PAN, Germany) and cultured in an incubator at 37°C with 5% CO₂. Mouse monoclonal antibody against VP2 protein was prepared and stored in our lab. Alexa Fluor 488-conjugated and HRP-conjugated goat anti-mouse IgG antibodies were both stored in our lab (Thermo Fisher Scientific, United States).

Viral RNA extraction and RT-PCR assay

The total viral RNA was extracted from the bursal sample supernatant using the TRIzol^® Reagent (Invitrogen, United States) according to the manufacturer’s instructions. cDNA was reverse transcribed from the viral RNA with M-MLV (H-) Reverse Transcriptase (Vazyme, China) and random primers (Invitrogen, United States). The PCR was performed using a pair of primers, the forward primer (F: 5’-ATGACAAACCTGCAAGATCAAACCCAAC-3′) and the reverse primer (R: 5’-CCTTATGGCCCGGATTATGTCTTTGAAG-3′), which amplified one 1,356-bp fragment of the VP2 gene. A master mix was prepared by using 10 μl of 2 × Taq Master Mix (Vazyme, China), 10 μmol/l of forward primer and reverse primer, 1 μl of cDNA and adding ddH₂O to a total volume of 20 μl. The PCR was performed at 95°C for 5 min followed by 35 cycles at 95°C for 30 s, 56°C for 30 s and 72°C for 1 min and 15 s, and a final extension at 72°C for 10 min.

Sequencing and analysis

The RT-PCR product of VP2 was analyzed by 1% agarose gel electrophoresis before being purified with the FastPure^® Gel DNA Extraction Mini Kit (Vazyme, China) according to the manufacturer’s instructions. The purified product was cloned into pMD-19 T vector (Takara, Japan) and sequenced by Sangon Biotech Co., Ltd. (Shanghai, China). The complete sequence of VP2 gene was assembled using the DNASTAR software (version 5.0, United States). The alignment and phylogenetic analysis based on the nucleotide or amino acid sequences were performed using the Clustal W program of MEGA X software (version 6.4, New Zealand) and the confidence levels were assessed using 1,000 bootstrap replications. Forty-two sequences of very virulent strains, attenuated strains, classic strains, variant strains, and serotype II strains from GenBank were selected as the reference strains (Table 1). Potential recombination event and breakpoints were identified with RDP4 software (version 4.5, Martin et al., 2015). Further verification was performed with the SimPlot software (version 3.5, Lole et al., 1999).

TABLE 1

Table 1. Reference IBDV strains used in phylogenetic analysis and the sequence alignment.

Virus isolation

The supernatant above was inoculated into the chorioallantoic membranes (CAMs) of 9-day-old SPF chick embryos. Six serial embryo passages were performed. The main bodies of dead embryos were homogenized together with the CAMs in sterile PBS to make a suspension, and the debris was spun down by centrifugation at 10,000 × g for 10 min at 4°C. The supernatant was filtered by 0.22-μm membrane and then inoculated into MC38 cells. The cells were harvested by freezing and thawing three times if no cytopathic effect (CPE) shown during 5-day monitoring.

Indirect immunofluorescence assay

Approximately 1 × 10⁶ MC38 cells seeded on a 6-well plate (Nest Biotechnology, China) were infected with the isolate. The cells were fixed with 4% paraformaldehyde after washing three times with sterile PBS at 24 h post-infection. Cold methanol and 5% non-fat milk were, respectively, used to permeabilize for 15 min at −20°C and block the cells for 1 h at room temperature. Then, the cells were incubated with the mouse anti-VP2 monoclonal antibody and Alexa Fluor 488-conjugated goat anti-mouse IgG secondary antibody. After staining the cellular nuclei with DAPI (Beyotime, China), the images were obtained using an inverted fluorescent microscope.

Pathogenicity experiment

A total of fourteen 3-week-old SPF chicks were randomly divided into two groups. One group (n = 7) was inoculated with the isolated strain by intraperitoneal injection and the other (n = 7) was inoculated with sterile PBS as mock-infection. Food and water were provided ad randomly throughout the study. All chicks were observed for 10 days and the clinical symptoms were recorded daily. At 10 d post-infection, bursal tissues of two groups were collected and weighed. The BF:body weight index (BBIX) [BBIX = (BF:body weight ratio)/(BF:body weight ratio in negative group)] was calculated along with the standard deviation. Bursa with a BBIX less than 0.70 was considered as atrophy (Lucio and Hitchner, 1979).

Flow cytometry

Bursal samples were cut into fragments in cold PBS and the resulting tissue suspensions were passed through 40-mm cell strainers. After washing with cold PBS two times, the cells were suspended in 1 × Binding Buffer to adjust the concentration to 1 × 10⁶ cells/mL. Then, 100 μl resulting cells were transferred to a 5-mL culture tube following by staining with 5 μl annexin-V/FITC (20 μg/ml; Solarbio, China) for 10 min and 5 μl PI (50 μg/ml, Solarbio, China) for 5 min at room temperature in darkness. Finally, 400 μl PBS was added to the stained cells and flow cytometry was preformed within 1 h.

Histopathology, immunohistochemical, and TUNEL staining

For analysis of light microscopy, the bursal tissue samples fixed with 4% paraformaldehyde were dehydrated in a graded series of ethanol, embedded in paraffin, and sectioned serially at a thickness of 4 μm. The obtained tissue sections were collected on slides and stained with hematoxylin and eosin (H&E). Furthermore, immunohistochemical (IHC) staining was subject to perform by using the mouse anti-VP2 monoclonal antibody. As for detection of DNA fragmentation, the sections were stained with the One Step TUNEL Apoptotic Assay Kit (Beyotime, China) according to the manufacturer’s instructions. Positive cells labeled red fluorescence were detected by inverted fluorescence microscope.

Statistical analysis

All statistical analyses were performed using GraphPad Prism (version 8.01, United States). Difference in this study was compared by unpaired t-test and statistical significance was set at p < 0.05. ** represents p < 0.01 and **** represents p < 0.0001.

Results

Isolation and identification of the novel variant IBDV isolate

The bursal sample from suspected IBD outbreak was tested positive for IBDV infection by amplifying a 1,356-bp specific fragment. After detection, the positive supernatant was inoculated into chick embryos for initial culture adaption. Embryos died within 72–120 h after inoculation with the isolate (designated LY21/2). Dead embryos showed lesions, such as incrassated CAMs and dwarfed carcass accompanied by congestion and haemorrhage (Figure 1A). The strain LY21/2 adapted to embryo culture was used to inoculate MC38 cells which attained shrinkage and agglomeration (Figure 1B). And the isolate was further identified by indirect immunofluorescence assay (IFA) using the mouse anti-VP2 monoclonal antibody. VP2 proteins (green) were present in the cytoplasm of LY21/2-infected cells, nevertheless, no positive signals were detected in negative cells (Figure 1C).

FIGURE 1

Figure 1. Isolation and identification of the variant IBDV strain LY21/2. (A) Isolation of the strain LY21/2 in 9-day-old chick embryo. Dwarfed carcass with congestion and haemorrhage was showed in dead embryo infected with LY21/2. (B) Isolation of the strain LY21/2 in cultured MC38 cells. The cells undergone shrinkage and aggregation at 24 h post infected with LY21/2. Scale bar, 100 μm. (C) Immunological identification of the strain LY21/2 VP2 proteins expressing in MC38 cells. MC38 cells with LY21/2 infection were fixed and incubated with the mouse anti-VP2 monoclonal antibody followed by an Alexa Fluor 488-conjugated goat anti-mouse IgG secondary antibody. MC38 cells expressing VP2 proteins (green) and nuclear staining (blue) with DAPI were merged. Scale bar, 200 μm.

Molecular characteristics of the VP2 gene

A phylogenetic tree based on the nucleotide sequences of the representative fragment of VP2 showed that the IBDV strains were distinctly divided into five major branches, including vvIBDV, attenuated strains, classic strains, variant strains, and serotype II strains. Phylogenetic analysis revealed that the isolate LY21/2 in this study belonged to the same branch as the variant strains (Figure 2A). Moreover, the variant strains can be further divided into two types, the early variant strains and the novel variant strains. The result of further phylogenetic analysis showed that LY21/2 was a novel variant strain (Figure 2A). VP2 gene of LY21/2 showed 93.8–98.6% nucleotide sequence identity with the corresponding regions from variant IBDV reference strains, and 96.8–98.6% identity with the novel variant reference strains. The HVR of VP2 was further analyzed by amino acid alignment, which showed that LY21/2 possessed the same characteristic amino acid residues as the reference novel variant strains, including 222 T, 249 K, 286I, and 318D (Figure 2B). However, two different amino acid residues (336A and 343S) were observed in the HVR of LY21/2 VP2 compared to other reference strains (336 T and 343P; Figure 2B).

FIGURE 2

Figure 2. Molecular Characteristics of the VP2 gene. (A) Phylogenetic analysis based on the nucleotide sequences of VP2 by MEGA X software v6.4. The phylogenetic tree represented the isolate LY21/2 formed one branch with novel variant IBDVs. The confidence levels were assessed using 1,000 bootstrap replications. The strain LY21/2 is indicated by red font. (B) The amino acid sequences in the HVR (206–350 aa) of LY21/2 VP2 were compared with variant strains, classic strains, attenuated strains and very virulent strains. The red boxes represent the characteristic amino acid residues (222 T, 249 K, 286I, and 318D). The blue boxes represent two different amino acid residues (336A and 343S) compared to other reference strains. LY21/2 is shown with a black bot.

Recombination event of the novel variant IBDV isolate

Potential recombination events between different IBDV strains were detected using the RDP4 based on VP2 nucleotide sequences. The result revealed that the novel variant IBDV strain LY21/2 as the major parent was involved in a recombination event of the variant 19D69 strain (MT550878, South Korea), and the vvIBDV Harbin-1 strain (EF517528, China) was the minor parent (Figure 3A). SimPlot was used to further confirm the result and to identify the most likely recombination breakpoints. As Figure 3B shown, the most likely recombination breakpoints were at nt 637 and nt 1,269 in VP2 sequence, which were correspond to amino acid residues 213 aa and 423 aa.

FIGURE 3

Figure 3. Recombination event of the variant IBDV strain LY21/2. (A) Comparisons of the similarity of 19D69 VP2 gene scales with LY21/2 (cyan) and Harbin-1 (purple) were calculated by RDP4 v4.5. The strain LY21/2 was the major parent and the strain Harbin-1 was the minor parent. (B) The recombination breakpoints were determined by using Simplot v3.5 and located at nt 637 and nt 1,269, respectively. The major parent (LY21/2) is shown in green, and the minor parent (Harbin-1) is shown in blue. The recombination breakpoints are both exhibited as vertical lines.

Pathogenicity of the novel variant IBDV isolate

To evaluate the pathogenicity of the novel variant IBDV strain LY21/2, 3-week-old chicks were inoculated with LY21/2 or sterile PBS by enterocoelia. Chicks inoculated with LY21/2 showed no gross clinical symptoms or mortality compared to the mock-infected chicks. The body weight of LY21/2-infected chicks was negatively influenced (p < 0.01; Figure 4B). Bursal atrophy was observed in LY21/2-infected chicks (Figure 4A), and as expected, the BBIX of infection group (0.219 ± 0.027) was lower than 0.7 (Figure 4C). After separating the bursal tissues into individual cells, the apoptosis rates of two groups were detected by flow cytometry (Figure 4D). Compared with the mock-infection group, the mean apoptosis rate was 49.46% higher in LY21/2-infected group (p < 0.0001; Figure 4E). Bursal tissues from the LY21/2-infected group and from the mock-infected group were subjected to pathological sections followed by H&E staining, IHC staining and TUNEL staining. According to the result of H&E staining, the numbers of lymphocytes were decreased and lymphoid nodule atrophy as well as connective tissue hyperplasia were observed in LY21/2-infected bursal tissue (Figure 4F). No obvious pathological changes showed in mock-infected bursal tissue (Figure 4F). IBDV antigen-positive cells were observed in the bursal tissue of LY21/2-infected chicks, while no antigen-positive cells showed in the bursa of mock-infected chicks (Figure 4G). To detect whether DNA fragmentation occurred in apoptotic bursal cells, TUNEL staining was used to perform in bursal tissue sections. Red TUNEL-positive cells were observed in LY21/2-infected bursal tissue, nevertheless, no red signal was detected in mock-infected group (Figure 4H).

FIGURE 4

Figure 4. Pathogenicity of the novel variant IBDV strain LY21/2. 3-week-old SPF chicks were inoculated with the IBDV strain LY21/2 and euthanatized at 10 days post-infection. Bursal atrophy was observed in infected chick (A).The mean body weight (B) and the BBIX values (BF:body weight index) (C) of LY21/2-infected chicks were both negatively influenced. The apoptosis rates of bursal cells of two groups were detected, and the mean rate of LY21/2-infected cells were 49.46% higher than that of mock-infected cells (D,E). **p < 0.01, ****p < 0.0001; an unpaired t-test was used to compare differences between the LY21/2-infected and mock-infected group. (F) Histological observation of bursal tissues from the mock- and LY21/2-infected chicks. The bursal section of the mock-infected group showed no apparent histological lesions. However, lymphocyte depletion and connective tissue hyperplasia were observed in the bursa of LY21/2-infected chicks. Arrows indicate atrophy of the follicle and lymphocyte depletion in the bursa. Empty triangle indicates a proliferation of connective tissue. (G) IHC staining of bursal tissues from the non-infected and infected chicks. No staining was observed in the bursal tissue from non-infected chicks. The bursal sample of LY21/2-infected chicks showed IBDV antigen-positive cells using the anti-VP2 monoclonal antibody. (H) TUNEL staining of bursal sections. Infection of LY21/2 induced DNA fragmentation in apoptotic bursal cells (red), whereas, positive signal was hardly detected in normal bursal cells. Scale bar, 100 μm.

Discussion

After classic IBDV first reported in 1957 (Muller et al., 2003), one antigenic variant IBDV emerged in the 1980s, which led a failure of the immune protection against the classic IBDV (Jackwood and Saif, 1987). In eastern China, it has been reported that subclinical infection of IBD caused by variant strains in some poultry farms since 2017 resulted in considerable economic losses (Coulibaly et al., 2005; Letzel et al., 2007; Fan et al., 2019). Because of the resulting severe bursal damage and immunosuppression, the disease induced by variant IBDV has been regard as an economically significant disease worldwide. It is noticeable that quite a few variant strain cases were detected among very virulent strains, therefore, seasonable monitoring and accurate control are of importance.

IBDV strains can replicate in different systems from avian cell lines such as DF-1 (Chen et al., 2020) or mammalian cell lines such as Vero (Rasool and Hussain, 2006). Nevertheless, the majority of IBDV strains from clinical isolates do not replicate in those systems without adaptation (Soubies et al., 2018; Wang et al., 2020). In addition, there has been an exploration of the possibility of IBDV strains replicating in primary bursal cells (Soubies et al., 2018) or tumor cell lines utilizing the apoptotic characteristics (Delgui et al., 2009). The result of culture in MC38 cells, a kind of tumor cell line from mouse colon, showed that the isolate LY21/2 could achieve virus replication in this cell line for the first time and provided a new reference for in vitro systems to replicate IBDV.

As is well-known to all, VP2 is usually used for the identification and analysis of IBDV because of its antigenic determinant (Becht and Muller, 1991; Brandt et al., 2001). In this study, the IBDV strain LY21/2 formed one branch with variant IBDVs according to the phylogenetic analysis based on the VP2 gene nucleotide sequences. We found that LY21/2 shared 93.8–98.6% nucleotide sequence identity with the variant reference stains. The antigenic phenotype of IBDV is determined by the HVR (206–350 aa) of VP2 gene. It was showed by further analysis that the HVR of VP2 gene of LY21/2 contained identical characteristic residues compared to variant IBDVs, including 222 T, 249 K, 286I, and 318D. Besides, the same distinct amino acid sites (221 K, 252I and 299S) were also observed in LY21/2 as the novel variant strains compared to the early variant strains. These characteristic residues were, respectively, located in four loop structures in the HVR which were designated P_BC (204–236 aa), PDE (240–265 aa), P_FG (270–293 aa) and P_HI (305–337 aa; Coulibaly et al., 2005). As reported, the residue 222 T (in P_BC) is closely related to the antigenicity of virus (Jackwood and Sommer-Wagner, 2011) and replication efficiency (Qi et al., 2016). Besides, residues 249 K (in P_DE), 286I (in P_FG) and 318D (in P_HI) are involved in the reactivity of different monoclonal antibodies or contribute to the antigenic drift of IBDV (Vakharia et al., 1994; Qi et al., 2013; Fan et al., 2022). Nevertheless, previous studies have confirmed that exchange of a single amino acid at the HVR was sufficient for altering the neutralizing properties, replication efficiency or virulence (Qi et al., 2009, 2013, 2016). In this study, we found that one different amino acid residue, T336A, was observed in the HVR of LY21/2 VP2 compared to other reference strains. In spite of that seemed a rare or unique mutation, further experiments are needed to be confirmed.

Homologous recombination contributes to change on antigenic variation and viral pathogenicity. In some previous studies, recombination among classic, very virulent, attenuated and variant IBDV strains in VP2 resulted in the exchange of some antigenic epitopes and change on pathogenicity (Wu et al., 2020; Feng et al., 2022; Guzman et al., 2022). The result of recombination analysis showed that LY21/2 as the major parent was involved in the recombination event of the novel variant 19D69 strain. The recombination breakpoints were at nt 637 and nt 1,269, respectively, corresponding to the residues 213 aa and 423 aa. Notably, the recombination fragment covered the HVR of VP2 hinting that there might be differences in viral antigenicity or pathogenicity after recombination event happened even if they were both classified as the novel variant strains. Moreover, it is still possible of LY21/2 to be a recombination parent of other strains during the process of nature evolution.

Different virulence of IBDV may lead to different pathological changes, such as enlarged and hemorrhagic bursa caused by very virulent strains, while the variant strains cause bursal atrophy. To evaluate the pathogenicity of the novel variant LY21/2 strain, 3-week-old SPF chicks were infected in our study. No gross clinical symptoms or death were present caused by LY21/2 during 10-day observation. However, we weighed all the chicks and found that the mean body weight of IBDV-infected group were negatively influenced. Experiment in which bursectomized chicks survived IBDV infection proves that bursa is the target organ for serotpye I strains (Kaufer and Weiss, 1980). The degree of damage to bursa differs in diverse virulence which has been used as an index of virulence for subclinical cases induced by the variant strains. In this present study, bursal atrophy was observed by necropsy in the LY21/2-infected chicks and correspondingly the BBIX decreased to 0.219 ± 0.027, when it is generally acknowledged that bursal atrophy occurs when the BBIX is less than 0.7. The result of histopathology further verified that after infected with LY21/2, bursal cells depletion was observed in the bursal tissue. Although the novel variant IBDVs do not cause death in chicks, it is likely to weaken the central immune function due to severe atrophy occurred in the bursa which may result in immunosuppression even vaccine failure and secondary or complicating diseases (Tomas et al., 2019; Aliyu et al., 2022). Many factors are involved in tissue or organ atrophy, such as apoptosis. It is reported that IBDV targets the precursors of antibody-producing B lymphocytes in bursa which leads to lymphoid depletion of B cells and apoptosis is responsible for that depletion (Cubas-Gaona et al., 2018). The application of double-staining technique in our study showed that apoptosis in cells of bursa was induced by LY21/2 infection, and correspondingly the apoptosis rate was significantly increased compared to mock-infected group. Besides, DNA fragmentation was detected in bursal tissue of LY21/2-infected chicks which is one of the main features of apoptosis. It might further account for infection by novel variant IBDV strains resulting in bursal atrophy and B lymphocytes depletion in young chicks.

In summary, we isolated and investigated the molecular characteristics and pathogenicity of a novel variant IBDV strain LY21/2. Our results revealed that LY21/2 strain could replicate in MC38 cells, shared a high level of nucleotide sequence identity and several characteristic amino acid residues in the HVR of VP2 gene with the novel variant IBDVs. Besides, LY21/2 as the major parent was involved in a recombination event of other one variant IBDV stain. Interestingly, LY21/2 had no obvious clinical signs in 3-week-old SPF chicks, while it caused bursal atrophy and lesions, as well as increased apoptosis rate. Findings of our study indicate that it is urgent to take action on further molecular analysis and pathogenicity tests for developing a new vaccine program for IBDV.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found at: https://www.ncbi.nlm.nih.gov/nuccore/OQ566936.

Ethics statement

The animal study was reviewed and approved by Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences.

Author contributions

YH and ZC conceived and designed the experiments. YH, CH, HC, and JH performed the experiments. YH, JL, HC, and ZC analyzed the data. HC and GS contributed to reagents/materials/analysis tools. YH, GS, and ZC wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by National Key R&D Program of China (2022YFD1801000).

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.

The reviewer KL declared a past collaboration with the authors to the handling editor at the time of review.

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.

References

Aliyu, H. B., Hamisu, T. M., Hair Bejo, M., Omar, A. R., and Ideris, A. (2022). Comparative pathogenicity of Malaysian variant and very virulent infectious bursal disease viruses in chickens. Avian Pathol. 51, 76–86. doi: 10.1080/03079457.2021.2006604

PubMed Abstract | CrossRef Full Text | Google Scholar

Azad, A. A., Barrett, S. A., and Fahey, K. J. (1985). The characterization and molecular cloning of the double-stranded RNA genome of an Australian strain of infectious bursal disease virus. Virology 143, 35–44. doi: 10.1016/0042-6822(85)90094-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayliss, C. D., Spies, U., Shaw, K., Peters, R. W., Papageorgiou, A., Muller, H., et al. (1990). A comparison of the sequences of segment a of four infectious bursal disease virus strains and identification of a variable region in VP2. J. Gen. Virol. 71Placeholder Text:1303–1312. doi: 10.1099/0022-1317-71-6-1303

PubMed Abstract | CrossRef Full Text | Google Scholar

Becht, H., and Muller, H. (1991). Infectious bursal disease--B cell dependent immunodeficiency syndrome in chickens. Behring Inst. Mitt. 89, 217–225.

Google Scholar

Brandt, M., Yao, K., Liu, M., Heckert, R. A., and Vakharia, V. N. (2001). Molecular determinants of virulence, cell tropism, and pathogenic phenotype of infectious bursal disease virus. J. Virol. 75, 11974–11982. doi: 10.1128/JVI.75.24.11974-11982.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Lian, J., Liang, Z., Leng, M., Lin, W., and Chen, F. (2022). Characterization and pathogenicity of infectious bursal disease virus in southern China. Poult. Sci. 101:102018. doi: 10.1016/j.psj.2022.102018

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Qin, Y., Qian, K., Shao, H., Ye, J., and Qin, A. (2020). HSC70 is required for infectious bursal disease virus (IBDV) infection in DF-1 cells. Virol. J. 17:65. doi: 10.1186/s12985-020-01333-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chettle, N., Stuart, J. C., and Wyeth, P. J. (1989). Outbreak of virulent infectious bursal disease in East Anglia. Vet. Rec. 125, 271–272. doi: 10.1136/vr.125.10.271

PubMed Abstract | CrossRef Full Text | Google Scholar

Coulibaly, F., Chevalier, C., Gutsche, I., Pous, J., Navaza, J., Bressanelli, S., et al. (2005). The birnavirus crystal structure reveals structural relationships among icosahedral viruses. Cells 120, 761–772. doi: 10.1016/j.cell.2005.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Cubas-Gaona, L. L., Diaz-Beneitez, E., Ciscar, M., Rodriguez, J. F., and Rodriguez, D. (2018). Exacerbated apoptosis of cells infected with infectious Bursal disease virus upon exposure to interferon alpha. J. Virol. 92:e00364-18. doi: 10.1128/JVI.00364-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, P., Ma, S. J., Zhang, Y. G., Li, X. S., Gao, X. Y., Cui, B. A., et al. (2013). Genomic sequence analysis of a new reassortant infectious bursal disease virus from commercial broiler flocks in Central China. Arch. Virol. 158, 1973–1978. doi: 10.1007/s00705-013-1682-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Delgui, L., Gonzalez, D., and Rodriguez, J. F. (2009). Infectious bursal disease virus persistently infects bursal B-lymphoid DT40 cells. J. Gen. Virol. 90, 1148–1152. doi: 10.1099/vir.0.008870-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Eterradossi, N., Arnauld, C., Toquin, D., and Rivallan, G. (1998). Critical amino acid changes in VP2 variable domain are associated with typical and atypical antigenicity in very virulent infectious bursal disease viruses. Arch. Virol. 143, 1627–1636. doi: 10.1007/s007050050404

PubMed Abstract | CrossRef Full Text | Google Scholar

Eterradossi, N., Toquin, D., Rivallan, G., and Guittet, M. (1997). Modified activity of a VP2-located neutralizing epitope on various vaccine, pathogenic and hypervirulent strains of infectious bursal disease virus. Arch. Virol. 142, 255–270. doi: 10.1007/s007050050075

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, L., Wang, Y., Jiang, N., Chen, M., Gao, L., Li, K., et al. (2020). Novel variant infectious bursal disease virus suppresses Newcastle disease vaccination in broiler and layer chickens. Poult. Sci. 99, 6542–6548. doi: 10.1016/j.psj.2020.09.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, L., Wang, Y., Jiang, N., Gao, Y., Niu, X., Zhang, W., et al. (2022). Residues 318 and 323 in capsid protein are involved in immune circumvention of the atypical epizootic infection of infectious bursal disease virus. Front. Microbiol. 13:909252. doi: 10.3389/fmicb.2022.909252

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, L., Wu, T., Hussain, A., Gao, Y., Zeng, X., Wang, Y., et al. (2019). Novel variant strains of infectious bursal disease virus isolated in China. Vet. Microbiol. 230, 212–220. doi: 10.1016/j.vetmic.2019.01.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, X., Zhu, N., Cui, Y., Hou, L., Zhou, J., Qiu, Y., et al. (2022). Characterization and pathogenicity of a naturally reassortant and recombinant infectious bursal disease virus in China. Transbound. Emerg. Dis. 69, e746–e758. doi: 10.1111/tbed.14347

PubMed Abstract | CrossRef Full Text | Google Scholar

Guzman, M., Cadiz, L., Guerrero-Moncayo, A., Caceres, F., Vidal, S., Lapierre, L., et al. (2022). Molecular characterization of infectious Bursal disease virus isolated in Chile reveals several mutations in VP2 coding region and a reassortment in its genome. Vet. Res. Commun. 46, 1281–1289. doi: 10.1007/s11259-022-09956-x

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Wang, W., Chen, G., Jiao, P., Ji, Z., Yang, L., et al. (2019). Serological study reveal different antigenic IBDV strains prevalent in southern China during the years 2000-2017 and also the antigenic differences between the field strains and the commonly used vaccine strains. Vet. Microbiol. 239:108458. doi: 10.1016/j.vetmic.2019.108458

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Xiong, Z., Yang, L., Guan, D., Yang, X., and Wei, P. (2014). Molecular epidemiology studies on partial sequences of both genome segments reveal that reassortant infectious bursal disease viruses were dominantly prevalent in southern China during 2000-2012. Arch. Virol. 159, 3279–3292. doi: 10.1007/s00705-014-2195-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackwood, D. H., and Saif, Y. M. (1987). Antigenic diversity of infectious bursal disease viruses. Avian Dis. 31, 766–770. doi: 10.2307/1591028

CrossRef Full Text | Google Scholar

Jackwood, D. J., and Sommer-Wagner, S. E. (2011). Amino acids contributing to antigenic drift in the infectious bursal disease Birnavirus (IBDV). Virology 409, 33–37. doi: 10.1016/j.virol.2010.09.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaufer, I., and Weiss, E. (1980). Significance of bursa of Fabricius as target organ in infectious bursal disease of chickens. Infect. Immun. 27, 364–367. doi: 10.1128/iai.27.2.364-367.1980

PubMed Abstract | CrossRef Full Text | Google Scholar

Letzel, T., Coulibaly, F., Rey, F. A., Delmas, B., Jagt, E., van Loon, A. A., et al. (2007). Molecular and structural bases for the antigenicity of VP2 of infectious bursal disease virus. J. Virol. 81, 12827–12835. doi: 10.1128/JVI.01501-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Lian, J., Wang, Z., Xu, Z., Pang, Y., Leng, M., Tang, S., et al. (2022). Pathogenicity and molecular characterization of infectious bursal disease virus in China. Poult. Sci. 101:101502. doi: 10.1016/j.psj.2021.101502

PubMed Abstract | CrossRef Full Text | Google Scholar

Lole, K. S., Bollinger, R. C., Paranjape, R. S., Gadkari, D., Kulkarni, S. S., Novak, N. G., et al. (1999). Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol. 73, 152–160. doi: 10.1128/JVI.73.1.152-160.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucio, B., and Hitchner, S. B. (1979). Response of susceptible versus immune chicks to killed, live-modified, and wild infectious bursal disease virus vaccines. Avian Dis. 23, 1037–1050. doi: 10.2307/1589620

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahgoub, H. A., Bailey, M., and Kaiser, P. (2012). An overview of infectious bursal disease. Arch. Virol. 157, 2047–2057. doi: 10.1007/s00705-012-1377-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, D. P., Murrell, B., Golden, M., Khoosal, A., and Muhire, B. (2015). RDP4: detection and analysis of recombination patterns in virus genomes. Virus Evol 1:vev003. doi: 10.1093/ve/vev003

PubMed Abstract | CrossRef Full Text | Google Scholar

Muller, H., Islam, M. R., and Raue, R. (2003). Research on infectious bursal disease--the past, the present and the future. Vet. Microbiol. 97, 153–165. doi: 10.1016/j.vetmic.2003.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Muller, H., Scholtissek, C., and Becht, H. (1979). The genome of infectious bursal disease virus consists of two segments of double-stranded RNA. J. Virol. 31, 584–589. doi: 10.1128/jvi.31.3.584-589.1979

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, X., Gao, H., Gao, Y., Qin, L., Wang, Y., Gao, L., et al. (2009). Naturally occurring mutations at residues 253 and 284 in VP2 contribute to the cell tropism and virulence of very virulent infectious bursal disease virus. Antivir. Res. 84, 225–233. doi: 10.1016/j.antiviral.2009.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, X., Gao, X., Lu, Z., Zhang, L., Wang, Y., Gao, L., et al. (2016). A single mutation in the PBC loop of VP2 is involved in the in vitro replication of infectious bursal disease virus. Sci. China Life Sci. 59, 717–723. doi: 10.1007/s11427-016-5054-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, X., Zhang, L., Chen, Y., Gao, L., Wu, G., Qin, L., et al. (2013). Mutations of residues 249 and 256 in VP2 are involved in the replication and virulence of infectious Bursal disease virus. PLoS One 8:e70982. doi: 10.1371/journal.pone.0070982

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranjbar, V. R., Mohammadi, A., and Dadras, H. (2019). Infectious bursal disease virus suppresses H9N2 avian influenza viral shedding in broiler chickens. Br. Poult. Sci. 60, 493–498. doi: 10.1080/00071668.2019.1621991

PubMed Abstract | CrossRef Full Text | Google Scholar

Rasool, M. H., and Hussain, I. (2006). Preparation and evaluation of Vero-cell infectious bursal disease vaccine in Pakistan. Vaccine 24, 2810–2814. doi: 10.1016/j.vaccine.2006.01.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, J. M., Kim, I. J., Rautenschlein, S., and Yeh, H. Y. (2000). Infectious bursal disease virus of chickens: pathogenesis and immunosuppression. Dev. Comp. Immunol. 24, 223–235. doi: 10.1016/S0145-305X(99)00074-9

CrossRef Full Text | Google Scholar

Soubies, S. M., Courtillon, C., Abed, M., Amelot, M., Keita, A., Broadbent, A., et al. (2018). Propagation and titration of infectious bursal disease virus, including non-cell-culture-adapted strains, using ex vivo-stimulated chicken bursal cells. Avian Pathol. 47, 179–188. doi: 10.1080/03079457.2017.1393044

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomas, G., Marandino, A., Courtillon, C., Amelot, M., Keita, A., Pikula, A., et al. (2019). Antigenicity, pathogenicity and immunosuppressive effect caused by a south American isolate of infectious bursal disease virus belonging to the "distinct" genetic lineage. Avian Pathol. 48, 245–254. doi: 10.1080/03079457.2019.1572867

PubMed Abstract | CrossRef Full Text | Google Scholar

Vakharia, V. N., He, J., Ahamed, B., and Snyder, D. B. (1994). Molecular basis of antigenic variation in infectious bursal disease virus. Virus Res. 31, 265–273. doi: 10.1016/0168-1702(94)90009-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Fan, L., Jiang, N., Gao, L., Li, K., Gao, Y., et al. (2020). Naturally occurring cell-adapted classic strain of infectious bursal disease virus. Vet. Microbiol. 243:108620. doi: 10.1016/j.vetmic.2020.108620

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Wang, Y., Li, H., Fan, L., Jiang, N., Gao, L., et al. (2020). Naturally occurring homologous recombination between novel variant infectious bursal disease virus and intermediate vaccine strain. Vet. Microbiol. 245:108700. doi: 10.1016/j.vetmic.2020.108700

PubMed Abstract | CrossRef Full Text | Google Scholar