- 1Institute for Poultry Science and Health, Guangxi University, Nanning, China
- 2Guangxi Key Laboratory for Polysaccharide Materials and Modifications, School of Marine Sciences and Biotechnology, Guangxi Minzu University, Nanning, China
Novel variant infectious bursal disease virus (nvIBDV) is an emerging genotype (A2dB1b) that can cause severe and prolonged immunosuppression in young chickens. Despite current commercial vaccines being proven to lack complete protection against nvIBDV, it remains unclear whether the oil emulsion inactivated vaccines (OEVs) of the homologous and heterologous virus or booster immunization can provide effective protection. In this study, OEVs with two types of nvIBDV isolates QZ191002 (A-nv/B-nv) and YL160304 (A-nv/B-HLJ0504-like) were prepared and evaluated the protective effects of OEVs plus the booster immunizations with different current commercial vaccines against the challenge of nvIBDVs. The results from vaccination-challenge experiments showed that nvIBDV could break through the protection provided by only one immunization dose of the commercial vaccines, with the protection rates ranging from 40% to 60%. Interestingly, even with booster immunization with different commercial vaccines, the protection rates could only be increased to 60%–80%. As expected, only the OEVs of the homologous virus could provide 100% protection against the homologous nvIBDV, which could induce high-level specific antibodies, ameliorate target organ damage, and significantly reduce the viral load of the bursal in the challenged chickens. Notably, YL160304-OEV performed better than QZ191002-OEV, providing 100% protection not only against the challenge of homologous strain but also against that of heterologous QZ191002 strain. Antibody levels of the immunized chickens gradually increased after a short decline and reached the highest level on the age of 28 days. Similarly, the percentages of lymphocytes CD4+, CD8+ T, and B in peripheral blood lymphocytes (PBLs) were significantly increased on 21 d and 28 d. Notably, despite the nvIBDV, OEVs initially induced a delayed responses in the early stages but ultimately reach higher levels of CD4+ and CD8+ T lymphocytes. The results of study suggest that even booster immunization with different commercial vaccines cannot provide complete protection against nvIBDV, while the OEVs made by the nvIBDVs can provide full protection. Moreover, YL160304-OEV exhibits a broader protective spectrum against different nvIBDV strains, making it a potential candidate for the development of new vaccine.
Highlights
• This study finds that even the booster immunization with the current commercial vaccines still cannot provide complete protection against the nvIBDV according to the current vaccination procedure.
• The oil emulsion vaccine (OEV) made by two inactivated nvIBDV isolates (QZ191002 and YL160304) can induce high-level specific antibodies, ameliorate target organ damage, and significantly reduce the viral load of the challenged virus in the vaccinated chickens and provide complete protection against the nvIBDV.
• The enhanced protective spectrum of the OEV of YL160304 (A-nv/B-HLJ0504-like vvIBDV) compared to QZ191002 (A-nv/B-nv) highlights that the role of VP1 is not only involved in the virulence of the virus but also seems to have a certain impact on the antigenicity of the virus and indicated that YL160304 strain could be a potential candidate for the development of new vaccine against the nvIBDV.
1 Introduction
Infectious bursal disease (IBD), caused by infectious bursal disease virus (IBDV), is a worldwide epidemic immunosuppressive disease in young chickens that have lymphoid tissue, especially the bursa of Fabricius (BF) of the central immune organ as its primary target (1, 2). IBDV, a member of the family of Birnaviridae, is a non-enveloped double-stranded RNA virus with high resistance to harsh environmental conditions. The IBDV genome consists of two segments (A and B) (3). Segment A encodes VP5 protein and, in another and partially overlapping reading frame, a polyprotein that will yield outer capsid protein VP2, protease VP4, and scaffold protein VP3 upon co-translational cleavage by VP4 (4), whereas segment B encodes VP1, the viral RNA-dependent RNA polymerase (RdRp), which exists in the viral capsids both as a genome-linked and as a free protein (5). IBDV specifically infects and lyses B lymphocytes, resulting in lymphoid depletion of B cells and atrophy of the bursa, which culminates in severe and long-term immunosuppression (6). The immunosuppressive effect of the disease reduces the bird’s resistance to other infections and reduces its immune responsiveness to the vaccinations and nature infections.
Two serotypes of the virus (serotypes 1 and 2) have been described, and only serotype 1 virus is pathogenic to chickens (7). Serotype 1 virus can be further divided into several pathogenic types, including classical (cIBDV) (8), very virulent (vvIBDV) (9), and antigenic variant (avIBDV) strains (10). The infections of different pathogenic IBDVs in chickens resulted in different clinical symptoms and pathogenic changes. cIBDV infection is characterized by bursal inflammation, severe lymphocyte necrosis, and immune damage, resulting in a mortality rate of 20%–30% in chickens (11, 12). attIBDV, which does not cause obvious clinical symptoms, is usually a domesticated attenuated vaccine strain that can protect against homologous wild IBDV strains. The characteristic of vvIBDV infection is that it can break through the protection of maternal antibodies, with mortality rate as high as 30%–100% (13, 14). In contrast, avIBDV does not lead to high mortality, often manifesting as a subclinical infection. Especially, the Chinese novel variant virus (nvIBDV) does not cause clinical morbidity and death, but it can cause serious damage to the BF and have a certain impact on the weight gain of chickens (15–24). Notably, compared to vvIBDV, nvIBDV exhibits significant differences in both VP2 and VP1, with a total of nine amino acid (aa) mutations (N77D, D213N, Q221K, A222T, Q249K, G254N, T289I, G318D, and D323E) in the antigen-associated projection (P) domain of VP2, seven mutations in the N-terminal domain (V4I, A24V, I61V, V/L141I, T145N, N/G147D, and A163V), five mutations in the central polymerase domain (D240E, A287T, S511R, E515D, and S646G), and one mutation in the C-terminal domain (P687S) of VP1 (20). Most of these mutations on VP2 and VP1 are related to the cell tropism, antigenic variation, and virulence diversity among IBDV strains (25–33). This may be an important factor (significant differences in the protective antigen-associated genes) why current commercial vaccines (mainly targeting the epidemic vvIBDV strains) cannot provide complete protection against nvIBDV.
In addition to strict biosafety measures, massive immunization is the optimal strategy and effective method to prevent and control IBDV infection. The immune effect of vaccines is closely related to effectively stimulating the immune defense system of chickens (34, 35). At present, commercial IBD vaccines mainly include live-attenuated, immune complex, and whole-virus inactivated vaccines (36, 37). Live-attenuated IBD vaccines can replicate in the bursa and produce humoral and cellular immunity simultaneously (38). However, large-scale use of these vaccines in the poultry industry may cause reassortment and recombination risks with field strains, thereby accelerating the variation of IBDV (20, 39, 40). In addition, there is a risk of them returning to a virulent form and causing immunosuppression, resulting in diseases in normal or immune-impaired chickens (41). Furthermore, these vaccines contain intact pathogens, making it difficult to distinguish infected chickens from immunized chickens in disease management. Compared with live-attenuated IBD vaccines, inactivated vaccines are considered to be safer and more effective in stimulating and inducing IBDV-specific antibodies and T-cell immune responses (36, 42), but they cannot establish long-term immunity (often requiring repeated immunization) (36). Therefore, these vaccines rely on adjuvants that stimulate strong and long-term humoral and cellular immune responses, which are essential for the prevention of infectious viruses. Oil emulsion adjuvant, also known as water-in-oil emulsion, can retain the antigen at the injection site, leading to the continuous release of the antigen, and then inducing long-term efficacy in poultry (43), which is one of the most common and effective adjuvants used currently in the veterinary field.
Some research groups have demonstrated that current commercial vaccines cannot provide complete protection against the increasing challenge of nvIBDV infection (20, 22, 44, 45). However, these research studies were only based on one immunization dose, which is different from the twice immunization dose used in the current immunization practice. This raises the possibility that incomplete protection might be due to an insufficient immune response stimulated and induced by once immunization dose. Therefore, we adopted the twice immunization regimen, following current production practices, and explored whether booster immunization with different vaccines can enhance immune responses to better resist the challenges posed by nvIBDV. Concurrently, we assessed the immunization efficacy of OEVs prepared from homologous nvIBDV strains. In this study, we developed two safe and effective oil emulsion inactivated vaccines (OEVs) based on the recently circulating nvIBDV isolates (20) and evaluated the protective effects of the newly developed nvIBDV OEVs and the booster immunization with different commercial vaccines according to the current clinical immunization practice procedures for the first time.
2 Materials and methods
2.1 Viruses and the commercial vaccines
Two purified nvIBDV strains QZ191002 (genotype A2dB1b, A-nv/B-nv) and YL160304 (genotype A2dB3, A-nv/B-HLJ0504-like vvIBDV), whose segment A derived from the same origin while segment B derived from different origin, were isolated and identified in the natural outbreaks of IBD that involved previously immunized chickens in southern China (20). The commercial vaccines A (IBD VP2 subunit vaccine), B (IBD immunogenic complex vaccine), C (IBD intermediate-plus live vaccine), and B87 attenuated vaccine used in this study were purchased from the commercial vaccine market.
2.2 Virus propagation and titer determination
As mentioned earlier, two types of nvIBDV strains propagate and titrate on HD11 cells (chicken macrophages cell line) that has been proven to be suitable for nvIBDV proliferation (20, 46, 47). In short, the virus was taken out from the preservation solution in the −80°C refrigerator, and the viruses were inoculated onto monolayer cells of 80%–85% confluent and cultured at 37°C for 1–2 h for virus adsorption, respectively. Subsequently, the inoculum was discarded and replaced with Dulbecco’s Modified Eagle Medium (DMEM) containing 1% fetal bovine serum (Lonsera, Shanghai Shuangru Biology Science & Technology Co., Ltd.) incubated for 72 h. When the cytopathic effect (CPE) was over 80%, cell cultures were harvested and then freeze-thawed three times. The harvested suspension was used to titrate the virus concentration according to the method previously described (48).
2.3 Virus inactivation and the preparation of the OEVs
The viruses were continuously inactivated with a final concentration of 0.4% formaldehyde in a shaking incubator at 37°C for 24 h. After the inactivated virus was determined to be sterile by sterility test, the inactivation efficiency was checked by three passages of each inactivated virus in specific pathogen-free (SPF) chicken embryos via the chorioallantoic membrane (CAM) route. The OEVs of the QZ191002 and YL160304 strains were prepared as the previously described method (49). In brief, add 4 parts of the sterilized Tween 80 to 96 parts of the inactivated virus was used as the aqueous phase. The oil phase, consisted of 90 parts No. 10 mineral oil and 10 parts Span 80, was sterilized by autoclaving. The aqueous and oil phases are emulsified in a ratio of 40/60 (v/v) at high speed (6000–10,000 r/min) and stored at 4°C until use.
2.4 Chicken vaccination-challenge experiment
A total of 230 1-day-old chicks of three-yellow chicken were purchased in a local commercial hatchery. All these chicks without any prior immunization were randomly divided into 23 groups (n = 10/group) with unrestricted access to feed and water. All groups were immunized twice (on 1 d and 14 d, respectively) with a 2-week interval by intramuscularly inoculated with the volume of 500 μL of vaccine/phosphate-buffered saline (PBS). The chickens in each group were challenged on 28 d with the nvIBDV strains QZ191002 and YL160304 by 105 TCID50/0.2 mL per bird via the oral route according to the previously described (20, 50), respectively. The design and vaccination schedule of the experiments are shown in detail in Table 1.
After the challenge, the clinical signs and possible mortalities were recorded every day. At 7 days post-challenge (dpc), all chickens of each group were euthanized for necropsy. The bursa/body-weight index (BBIX) was calculated to measure the degree of the BF atrophy, and it was considered severe atrophy if BBIX<0.70 (51). Each bursa was divided into two parts for the following analysis: A part of the bursa samples was performed to analyze the viral load by quantitative RT-PCR (52), and the other part was fixed in 4% neutral buffered formalin and stained with HE for further histopathological examination. The results were quantified according to Skeeles’ scale for histopathological bursal lesion scores (HBLSs) (53).
2.5 Detection of IBDV-specific antibodies
Serum samples were collected from the wing veins of the immunized and unimmunized chickens at 1, 7, 14, 21, and 28 d, respectively, and the titers of IBDV-specific antibodies were determined by enzyme-linked immunosorbent assay (ELISA) (54). In brief, the purified VP2 protein was diluted with a coating buffer (1.59 g Na2CO3, 2.93 g NaHCO3, 1,000 mL ddH2O, pH 9.6) to a final concentration of 0.8 μg/mL. Each well of the ELISA plate was coated with 400 μL of the diluted VP2 protein and incubated at 37°C for 2 h. After washing with PBST buffer (PBS containing 0.05% Tween-20, pH 7.4) for three times to remove unbound antigens and impurities, the wells were blocked with 5% skim milk (BD Difco™, United States) at 37°C for 2 h and wash three times as above. Subsequently, 100 μL of chicken serum were added to each well at a 1:100 dilution and incubated at 37°C for 90 min. After washing, 100 μL Rabbit Anti-Chicken IgG Antibody (H + L)-HRP Conjugated (Biosynthesis Biotechnology Inc. Beijing, China), diluted 1:4000, was added and incubated for 30 min. Following another wash, the 100 μL volume of 3,3′,5,5′-tetramethylbenzidine (TMB) ELISA substrate solution (TransGen Biotech, Beijing, China) was added to visualize the reaction and incubated at 37°C for 10 min in the dark. Finally, the reaction was terminated by adding 50 μL of 2 M H2SO4 to each well, and absorbance values at 450 nm were measured using a Microplate Spectrophotometer.
2.6 Flow cytometry analysis of the lymphocytes of CD4+, CD8+ T, and B
Peripheral blood lymphocytes (PBLs) were isolated from the anticoagulation-treated (with sodium citrate) blood samples that were sampled from the wing veins of the immunized and unimmunized chickens at 1, 7, 14, 21, and 28 d, respectively, as the previously described method (55). In brief, the anticoagulation-treated blood samples were gently added above the liquid level of the peripheral blood lymphocyte isolation solution, followed by centrifugation at 800 × g for 20 min to collect the PBLs. PBLs were stained with Mouse Anti-Chicken CD4-APC, Mouse Anti-Chicken CD8a-FITC, and Mouse Anti-Chicken Bu-1-PE antibodies (SouthernBiotech, Birmingham, United States) according to the manufacturer’s instructions. Finally, the percentages of CD4+, CD8+ T lymphocytes, and B lymphocyte in all samples were evaluated by Invitrogen Attune™ NxT Flow Cytometer.
2.7 Statistical analysis
All data analyses were processed with GraphPad Prism 8.0 software and expressed as means ± standard deviation (SD). Statistical differences among the groups were evaluated using one-way analysis of variance (ANOVA). To compare the differences among treatments, Tukey’s multiple comparisons test was used. The results with p-values of <0.05 *, p < 0.01 **, p < 0.001 ***, and p < 0.0001 **** were considered statistically significant.
3 Results
3.1 The anti-IBDV antibody detection and B lymphocyte count revealed a sustained high level of humoral immune response in the immunized chickens
The antibody responses of chickens induced by different vaccines were detected by the anti-IBDV antibody and B lymphocyte count. Chickens vaccinated with nvIBDV OEVs and the commercial vaccines all produced different levels of antibodies (Figure 1A). The antibody titer of the unimmunized chickens decreased after hatching and reached levels of lower on 21 d and the lowest on 28 d. While those of the immunized chickens gradually increased over time after a short decline period and reached the highest level on the 28 d, indicating that all immunized chickens showed clear serological conversion (acquired immunity). In addition, the antibody levels of the groups that immunized twice with the nvIBDV OEVs were all the highest levels during the experiment, and the groups that with vaccine B (the immunogenic complex vaccine) as the prime and with nvIBDV OEVs as the booster immunizations were the second high. The number of B lymphocyte in all vaccinated groups continued to increase over time after immunization, significantly higher than those in the unimmunized groups (Figure 1B). It is worth noting that the anti-IBDV antibody level and the number of B lymphocyte of the nvIBDV OEV-immunized groups were the highest in all the immunized groups, and its stimulating effect on the humoral immunity seems to be stronger than the current commercial vaccines.
Figure 1. Anti-IBDV antibody detection and the percentages of B lymphocyte in PBLs of the experiment groups. (A) The anti-IBDV antibody levels at 1, 7, 14, 21, and 28 d, respectively, detected by ELISA. (B) The percentages of B lymphocyte in PBLs, which were isolated from the immunized chickens at 1, 7, 14, 21, and 28 d, respectively, analyzed by flow cytometry.
3.2 IBDV-specific T-cell immune responses in the immunized chickens
The T-cell immune response level of the immunized chickens was detected by the percentages of CD4+ and CD8+ T lymphocytes in the PBLs by flow cytometry (Figure 2). The results showed that the CD4+ and CD8+ T lymphocytes in all the immunized groups were increased in varying degrees from 7 to 28 d. On 21 d and 28 d, the values in each immunized group were significantly higher than those in the unimmunized group (p < 0.05). Most notably, the values of the nvIBDV OEV groups were increased relatively slower in the early stage (7–14 d) but reached a higher level later on 21 d and 28 d.
Figure 2. Percentages of CD4+ and CD8+ T lymphocytes in PBLs of the experiment groups. PBLs were isolated from immunized chickens at 0, 7, 14, 21, and 28 d, respectively, and analyzed by flow cytometry. (A) The percentages of CD4+ T lymphocytes. (B) The percentages of CD8+ T lymphocytes.
3.3 The booster immunization with different current commercial vaccines still cannot provide complete protection against the nvIBDV
Under the challenges of nvIBDV strain QZ191002 and its reassortant strain YL160304, no obvious clinical symptoms and death were observed in the chickens immunized with the commercial vaccines and the nvIBDV OEVs (Figure 3). Compared with the control group, severe atrophy of BF was observed in the unimmunized challenged group (BBIX<0.5). The booster immunization of nvIBDV OEVs based on the prime immunization with commercial vaccine B showed relatively mild bursal atrophy when compared to other groups that only immunized with the commercial vaccines (BBIX>0.7) (Figures 3A,D). Similarly, the same phenomenon can be observed in terms of HBLS (Figures 3B,E) and viral load (Figures 3C,F). All data indicated that the booster immunization of the commercial vaccines currently used still cannot provide complete protection against the nvIBDV. However, the booster immunization with the nvIBDV OEVs on the prime immunization with the commercial vaccine B seemed to have a better protection effect than others.
Figure 3. BBIX, HBLS, and viral loads of the groups with booster immunization with the different vaccines and challenged with the nvIBDV strains. (A–C) The BBIX, HBLS, and viral loads of the BF, respectively, in the booster immunization groups challenged with QZ191002. (D–F) The BBIX, HBLS, and viral loads, respectively, in the booster immunization groups challenged with YL160304.
3.4 The nvIBDV OEVs can provide complete protection against the homologous virus and YL160304-OEV performs better than QZ191002-OEV
In the groups vaccinated with nvIBDV OEVs, no significant morbidity or mortality was observed in chickens challenged with QZ191002 or YL160304, as in the control group (Figure 4). The unimmunized group challenged with YL160304 strain resulted in 10% (1/10) of the mortality rate, while the unimmunized group challenged with QZ191002 showed no obvious clinical symptoms or mortality. Further study on the values of BBIX, HBLS, and viral load showed that no significant bursal atrophy (BBIX>0.7) (Figures 4A,D) and histopathologic bursal lesions (HBLS<1) (Figures 4B,E) were observed in the nvIBDV OEV-immunized group. However, significant atrophy (BBIX<0.5) and severe histopathologic lesions in the BF (HBLS>3) can be observed in the unimmunized challenged group. The viral loads in the immunized groups were significantly lower than that in the unimmunized challenged group, which is consistent with the phenomenon observed by HBLS (Figures 4C,F). Interestingly, the data also revealed that the degree of damage to the BF is related to the replication efficiency of the virus (viral load) in the organ. It is worth noting that the YL160304-OEV (A-nv/B-HLJ0504-like) showed to be able to provide better protection against the challenge than those of the QZ191002-OEV (A-nv/B-nv), indicating that VP1 of IBDV is not only involved in the virulence of the virus but also has a certain impact on the antigenicity of the virus (Figures 4D–F).
Figure 4. BBIX, HBLS, and viral loads of the nvIBDV OEV-immunized groups challenged with the nvIBDV strains. (A–C)The BBIX, HBLS, and viral loads, respectively, in the nvIBDV OEV-immunized groups challenged with QZ191002 strain. (D–F) The BBIX, HBLS, and viral loads, respectively, of the nvIBDV OEV-immunized groups challenged with YL160304 strain.
3.5 Protective efficiency of the current commercial vaccines and the nvIBDV OEVs against the nvIBDVs challenge
The protective efficiency of nvIBDV OEVs and different commercial vaccines in booster immunization against the challenge of nvIBDVs were further clarified according to the BBIX, HBLS, and viral loads of the BF from the immunization-challenge experiments, as shown in Table 2. The results showed that nvIBDV could break through the protections provided by only one immunization dose of the current commercial vaccines, with the protection rate usually being 40%–60%. Interestingly, even with booster immunization with different current commercial vaccines, the protection rates could only be increased to 60%–80%, indicating complete protection could not be provided. As expected, only the homologous nvIBDV OEVs could provide 100% protections against the challenges of the homologous strains, which could induce high-level specific antibodies, ameliorate target organ damage, and significantly reduce viral load (Figure 3). Notably, YL160304-OEV performed better than QZ191002-OEV, providing 100% protection not only against the challenge of homologous strain but also against that of heterologous QZ191002 strain (Table 2).
Table 2. Protective efficiency of the current commercial vaccines and nvIBDV OEVs in the immunized groups challenged with the nvIBDV strains.
4 Discussion
Although a large-scale immunization program against IBDV infection has been carried out in the poultry industry, more and more different types of IBDV were isolated increasingly year by year from the vaccinated chickens (56). At present, many IBDVs with different path−/geno−/serotypes have been confirmed to be co-circulating in China, and the current commercial vaccines cannot provide sufficient protection, especially for the newly emerging nvIBDV (20–22, 45, 57, 58). It is reported that the epidemic of nvIBDV occurred in the important poultry-raising areas of China during 2016–2022, and the infection and prevalence of the nvIBDV strain were found in 17 provinces, including the major chicken production provinces like Shandong, Guangxi, and Hebei (17, 21, 57, 59–61). Due to the frequent commercial trade of live poultry and products, nvIBDV has spread to Southern China, Eastern China, and Northeast China, with Central China as the major transmission source, and it has even emerged and spread rapidly in immunized chickens in East Asia (Japan and South Korea) and Southeast Asia (Malaysia) (21), which seems to have a trend of further spreading to other regions of the world. To date, the raised calls for frequent outbreaks of nvIBDV among the vaccinated chickens in Asian countries/regions have raised the significance of thinking about vaccine effectiveness.
IBDV belongs to the double-stranded RNA viruses, and its conformation makes it more prone to variation. Previous studies have indicated that IBDV mainly evolves in diversity through gene mutation and reassortment to escape the host’s immune mechanisms and facilitate better survival of its own virus (20, 21, 27, 30, 62–64). Although China has made efforts to control the disease through regular immunization with vaccines against vvIBDV, the emergence of nvIBDV is raising new challenges. In this study, the twice immunization of current commercial vaccines (subunit vaccine, killed vaccine, or intermediate-plus live vaccine all mainly targeting the epidemic vvIBDV strains) has failed to provide complete protection against the nvIBDV (60%–80%), although they can provide better protection than previously reported once immunization (40%–60%) (44, 45). Leng et al. (65) introduced the “placeholder effect” as a novel criterion for assessing vaccine effectiveness, showing that the W2512 vaccine causes severe atrophy of the BF in SPF chickens and Yellow chickens, induces high levels of antibodies against IBDV, and protects chickens from infection with the novel variant strains via a placeholder effect. The results of this study indicate that developing an efficient vaccine against the increasing prevalence of nvIBDV is crucial for the prevention and control of the disease right now and in future.
Essentially, the difference in protection rates is probably related to the variation of VP2 between the field and vaccine strains. VP2 is considered as a crucial structural protein of IBDV and an important dominant antigen for vaccine development. However, many researchers have confirmed that the change of some specific amino acids in the VP2 region might lead to virus antigen drift, further reducing the ability of antibodies to capture the virus (57, 66–68). Compared to vvIBDV and attIBDV strains, changes in multiple antigenic sites have been found in nvIBDV (21, 57). This high mutation rate helps to enable it to escape the capture of the immune system induced by the current commercial vaccines and persist for long periods in the immunopotent chickens. As found in this study, twice immunizations (1, 14 d) of nvIBDV OEVs can provide good protection (90%–100%) against the epidemic strains of nvIBDV, which is superior to the protection (60%–80%) provided by the traditional commercial vaccines’ immunization (primary immunization) and nvIBDV (boost immunization) in the experiment. It suggests that to provide complete protection against nvIBDV strains, boost immunization with nvIBDV OEVs may be necessary. This poses new challenges to the IBD vaccination programs of regions or countries that are affected by the emergence of nvIBDV in current or in future. As in addition to the emerging nvIBDV, there are also co-circulations of vvIBDV (which are still dominant strains in many places) and/or cIBDV in the field. So, developing vaccines of multi-viral antigens that include multiple types of viruses is a possible solution.
More and more studies indicate that VP1 of IBDV is not only involved in the virulence of the virus (25, 26, 29, 32, 69, 70) but also seems to have a certain impact on the antigenicity of the virus. Chen et al. found that NN1005 (A-vv/B-new type) and JS7 (A-vv/B-c), two OEVs with the same VP2, exhibit different protection rates due to their different origins of VP1 (49). This has also been confirmed in this immunization-challenge experiment; the comparative study of two different genotypes of nvIBDV showed that the protective spectrum of YL160304-OEV (A-nv/B-HLJ0504-like vvIBDV) was superior to QZ191002-OEV (A-nv/B-nv), indicating that YL160304 strain could be a good candidate for the development of new vaccine against the nvIBDVs. Through such reassortants (VP1 originating from vvIBDV), the YL160304 strain might enhance its pathogenicity to chickens by increasing the virulence of the virus (20) and might also escape from the body’s immunity by altering the viral antigenicity, thereby increasing the chances of its survival and epidemic spread in the vaccinated chickens (49). Furthermore, the same results were also confirmed in the role of VP1 in the pathogenicity of the virus strain in this challenge experiment, which showed that the pathogenicity of YL160304 strain (10% mortality) with its VP1 originated from HLJ0504-like to three-yellow chicken was significantly higher than that of QZ191002 (0% mortality), indicating that VP1 reassortment contributes to enhance the virulence of nvIBDV. Therefore, in addition to be recognized and proven as the most important protective antigen of VP2, different origins of VP1 might also play an important role in the antigenicity of IBDV, although this still needs to be further verified by the construct of recombinant virus with different origins of VP2 and VP1 by the reverse genetics.
5 Conclusion
At present, many IBDVs with different path−/geno−/serotypes have been confirmed to be co-circulating in China and some other Asian countries, and the current commercial vaccines cannot provide sufficient protection, especially against the newly emerging nvIBDV strains. Our results suggest that even the booster immunization with the current commercial vaccines still cannot provide complete protection against the nvIBDVs, while the nvIBDV OEVs can provide complete protection. In addition, YL160304-OEV (A-nv/B-HLJ0504-like) can provide a broader protective spectrum against different nvIBDV strains. The findings of this study provide important insights in the current vaccination practice procedures of the nvIBDV and will assist in the development of programs for control strategies for these emerging viruses.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The animal study was approved by the Animal Experimental Ethical Committee of Guangxi University (No. Gxu-2021-113). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
WW: Conceptualization, Data curation, Methodology, Writing – original draft. YH: Data curation, Writing – review & editing. YZ: Data curation, Writing – review & editing. YQ: Data curation, Writing – review & editing. JS: Data curation, Writing – review & editing. JH: Formal analysis, Investigation, Writing – review & editing. TH: Formal analysis, Investigation, Writing – review & editing. TW: Formal analysis, Investigation, Writing – review & editing. MM: Formal analysis, Investigation, Writing – review & editing. XH: Writing – review & editing. PW: Conceptualization, Funding acquisition, Methodology, Project administration, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Natural Science Foundation of China (31660717, 32160824, and 31560706), the Guangxi Program for Modern Agricultural Industry Technical System Construction-Chicken Industry (nycytxgxcxtd-19-03), the Innovation Project of Guangxi Graduate Education (YCBZ2022034), and the Guangxi Special Funding on Science and Technology Research (AA17204057).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Berg, TP . Acute infectious bursal disease in poultry: a review. Avian Pathol. (2000) 29:175–94. doi: 10.1080/03079450050045431
2. Muller, H, Islam, MR, and Raue, R. Research on infectious bursal disease—the past, the present and the future. Vet Microbiol. (2003) 97:153–65. doi: 10.1016/j.vetmic.2003.08.005
3. Muller, H, Scholtissek, C, and Becht, H. The genome of infectious bursal disease virus consists of two segments of double-stranded RNA. J Virol. (1979) 31:584–9. doi: 10.1128/JVI.31.3.584-589.1979
4. Hudson, PJ, McKern, NM, Power, BE, and Azad, AA. Genomic structure of the large RNA segment of infectious bursal disease virus. Nucleic Acids Res. (1986) 14:5001–12. doi: 10.1093/nar/14.12.5001
5. Muller, H, and Nitschke, R. The two segments of the infectious bursal disease virus genome are circularized by a 90,000-Da protein. Virology. (1987) 159:174–7. doi: 10.1016/0042-6822(87)90363-1
6. Rodenberg, J, Sharma, JM, Belzer, SW, Nordgren, RM, and Naqi, S. Flow cytometric analysis of B cell and T cell subpopulations in specific-pathogen-free chickens infected with infectious bursal disease virus. Avian Dis. (1994) 38:16–21. doi: 10.2307/1591831
7. Jackwood, DJ, Saif, YM, and Hughes, JH. Characteristics and serologic studies of two serotypes of infectious bursal disease virus in turkeys. Avian Dis. (1982) 26:871–82. doi: 10.2307/1589875
8. Cosgrove, AS . An apparently new disease of chickens: avian nephrosis. Avian Dis. (1962) 6:385–9. doi: 10.2307/1587909
9. Chettle, N, Stuart, JC, and Wyeth, PJ. Outbreak of virulent infectious bursal disease in East Anglia. Vet Rec. (1989) 125:271–2. doi: 10.1136/vr.125.10.271
10. Jackwood, DH, and Saif, YM. Antigenic diversity of infectious bursal disease viruses. Avian Dis. (1987) 31:766–70. doi: 10.2307/1591028
11. Dulwich, KL, Asfor, A, Gray, A, Giotis, ES, Skinner, MA, and Broadbent, AJ. The stronger downregulation of in vitro and in vivo innate antiviral responses by a very virulent strain of infectious bursal disease virus (IBDV), compared to a classical strain, is mediated, in part, by the VP4 protein. Front Cell Infect Microbiol. (2020) 10:315. doi: 10.3389/fcimb.2020.00315
12. Kabell, S, Igyártó, B-Z, Magyar, A, Hajdú, Z, Biró, É, Bisgaard, M, et al. Impact of heterophil granulocyte depletion caused by 5-fluorouracil on infectious bursal disease virus infection in specific pathogen free chickens. Avian Pathol. (2006) 35:341–8. doi: 10.1080/03079450600821141
13. Tammiranta, N, Ek-Kommonen, C, Rossow, L, and Huovilainen, A. Circulation of very virulent avian infectious bursal disease virus in Finland. Avian Pathol. (2018) 47:520–5. doi: 10.1080/03079457.2018.1503642
14. Nooruzzaman, M, Hossain, I, Rahman, MM, Uddin, AJ, Mustari, A, Parvin, R, et al. Comparative pathogenicity of infectious bursal disease viruses of three different genotypes. Microb Pathogenesis. (2022) 169:105641. doi: 10.1016/j.micpath.2022.105641
15. Aliyu, HB, Hamisu, TM, Hair-Bejo, M, Omar, AR, and Aini, I. Comparative pathogenicity of Malaysian variant and very virulent infectious bursal disease viruses in chickens. Avian Pathol. (2022) 51:76–86. doi: 10.1080/03079457.2021.2006604
16. Fan, L, Wang, Y, Jiang, N, Chen, M, Gao, L, Li, K, et al. Novel variant infectious bursal disease virus suppresses Newcastle disease vaccination in broiler and layer chickens. Poultry Sci. (2020) 99:6542–8. doi: 10.1016/j.psj.2020.09.037
17. Lian, J, Wang, Z, Xu, Z, Pang, Y, Leng, M, Tang, S, et al. Pathogenicity and molecular characterization of infectious bursal disease virus in China. Poultry Sci. (2022) 101:101502. doi: 10.1016/j.psj.2021.101502
18. Myint, O, Suwanruengsri, M, Araki, K, Izzati, UZ, Pornthummawat, A, Nueangphuet, P, et al. Bursa atrophy at 28 days old caused by variant infectious bursal disease virus has a negative economic impact on broiler farms in Japan. Avian Pathol. (2021) 50:6–17. doi: 10.1080/03079457.2020.1822989
19. Thai, TN, Jang, I, Kim, HA, Kim, HS, Kwon, YK, and Kim, HR. Characterization of antigenic variant infectious bursal disease virus strains identified in South Korea. Avian Pathol. (2021) 50:174–81. doi: 10.1080/03079457.2020.1869698
20. Wang, W, Huang, Y, Zhang, Y, Qiao, Y, Deng, Q, Chen, R, et al. The emerging naturally reassortant strain of IBDV (genotype A2dB3) having segment a from Chinese novel variant strain and segment B from HLJ 0504-like very virulent strain showed enhanced pathogenicity to three-yellow chickens. Transbound Emerg Dis. (2022) 69:e566–79. doi: 10.1111/tbed.14336
21. Wang, W, He, X, Zhang, Y, Qiao, Y, Shi, J, Chen, R, et al. Analysis of the global origin, evolution and transmission dynamics of the emerging novel variant IBDV (A2dB1b): the accumulation of critical aa-residue mutations and commercial trade contributes to the emergence and transmission of novel variants. Transbound Emerg Dis. (2022) 69:e2832–51. doi: 10.1111/tbed.14634
22. Xu, A, Pei, Y, Zhang, K, Xue, J, Ruan, S, and Zhang, G. Phylogenetic analyses and pathogenicity of a variant infectious bursal disease virus strain isolated in China. Virus Res. (2020) 276:197833. doi: 10.1016/j.virusres.2019.197833
23. Legnardi, M, Poletto, F, Talaat, S, Selim, K, Moawad, MK, Franzo, G, et al. First detection and molecular characterization of novel variant infectious bursal disease virus (genotype A2dB1b) in Egypt. Viruses. (2023) 15:2388. doi: 10.3390/v15122388
24. Tomás, G, Techera, C, Marandino, A, Olivera, V, Williman, J, Panzera, Y, et al. Genomic characterization of infectious bursal disease virus in Argentina provides evidence of the recent transcontinental spread of Chinese genotype A2dB1b. Avian Pathol. (2024) 4:1–9. doi: 10.1080/03079457.2024.2355918
25. Escaffre, O, Le Nouen, C, Amelot, M, Ambroggio, X, Ogden, KM, Guionie, O, et al. Both genome segments contribute to the pathogenicity of very virulent infectious bursal disease virus. J Virol. (2013) 87:2767–80. doi: 10.1128/JVI.02360-12
26. Gao, L, Li, K, Qi, X, Gao, Y, Wang, Y, Gao, H, et al. N-terminal domain of the RNA polymerase of very virulent infectious bursal disease virus contributes to viral replication and virulence. Sci China Life Sci. (2018) 61:1127–9. doi: 10.1007/s11427-017-9297-8
27. He, X, Xiong, Z, Yang, L, Guan, D, Yang, X, and Wei, P. 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. (2014) 159:3279–92. doi: 10.1007/s00705-014-2195-z
28. Jackwood, DJ, Sreedevi, B, LeFever, LJ, and Sommer-Wagner, SE. Studies on naturally occurring infectious bursal disease viruses suggest that a single amino acid substitution at position 253 in VP2 increases pathogenicity. Virology. (2008) 377:110–6. doi: 10.1016/j.virol.2008.04.018
29. Nouen, CL, Toquin, D, Muller, H, Raue, R, Kean, KM, Langlois, P, et al. Different domains of the RNA polymerase of infectious bursal disease virus contribute to virulence. PLoS One. (2012) 7:e28064. doi: 10.1371/journal.pone.0028064
30. Qi, X, Gao, H, Gao, Y, Qin, L, Wang, Y, Gao, L, et al. 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. (2009) 84:225–33. doi: 10.1016/j.antiviral.2009.09.006
31. Schnitzler, D, Bernstein, F, Muller, H, and Becht, H. The genetic basis for the antigenicity of the VP2 protein of the infectious bursal disease virus. J Gen Virol. (1993) 74:1563–71. doi: 10.1099/0022-1317-74-8-1563
32. Wang, W, Huang, Y, Ji, Z, Chen, G, Zhang, Y, Qiao, Y, et al. The full region of N-terminal in polymerase of IBDV plays an important role in viral replication and pathogenicity: either partial region or single amino acid V4I substitution does not completely Lead to the virus attenuation to three-yellow chickens. Viruses. (2021) 13:107. doi: 10.3390/v13010107
33. Yamaguchi, T, Ogawa, M, Inoshima, Y, Miyoshi, M, Fukushi, H, and Hirai, K. Identification of sequence changes responsible for the attenuation of highly virulent infectious bursal disease virus. Virology. (1996) 223:219–23. doi: 10.1006/viro.1996.0470
34. Courtillon, C, Allée, C, Amelot, M, Keita, A, Bougeard, S, Härtle, S, et al. Blood B cell depletion reflects immunosuppression induced by live-attenuated infectious bursal disease vaccines. Front Vet Sci. (2022) 9:871549. doi: 10.3389/fvets.2022.871549
35. Trapp, J, and Rautenschlein, S. Infectious bursal disease virus' interferences with host immune cells: what do we know? Avian Pathol. (2022) 51:303–16. doi: 10.1080/03079457.2022.2080641
36. Muller, H, Mundt, E, Eterradossi, N, and Islam, MR. Current status of vaccines against infectious bursal disease. Avian Pathol. (2012) 41:133–9. doi: 10.1080/03079457.2012.661403
37. Ramon, G, Legnardi, M, Cecchinato, M, Cazaban, C, Tucciarone, CM, Fiorentini, L, et al. Efficacy of live attenuated, vector and immune complex infectious bursal disease virus (IBDV) vaccines in preventing field strain bursa colonization: a European multicentric study. Front Vet Sci. (2022) 9:978901. doi: 10.3389/fvets.2022.978901
38. Lauring, AS, Jones, JO, and Andino, R. Rationalizing the development of live attenuated virus vaccines. Nat Biotechnol. (2010) 28:573–9. doi: 10.1038/nbt.1635
39. Lu, Z, Zhang, L, Wang, N, Chen, Y, Gao, L, Wang, Y, et al. Naturally occurring reassortant infectious bursal disease virus in northern China. Virus Res. (2015) 203:92–5. doi: 10.1016/j.virusres.2015.04.003
40. Wei, Y, Li, J, Zheng, J, Xu, H, Li, L, and Yu, L. Genetic reassortment of infectious bursal disease virus in nature. Biochem Bioph Res Co. (2006) 350:277–87. doi: 10.1016/j.bbrc.2006.09.040
41. Mazariegos, LA, Lukert, PD, and Brown, J. Pathogenicity and immunosuppressive properties of infectious bursal disease "intermediate" strains. Avian Dis. (1990) 34:203–8. doi: 10.2307/1591353
42. Rautenschlein, S, Yeh, HY, and Sharma, JM. The role of T cells in protection by an inactivated infectious bursal disease virus vaccine. Vet Immunol Immunop. (2002) 89:159–67. doi: 10.1016/s0165-2427(02)00202-7
43. Herbert, WJ . The mode of action of mineral-oil emulsion adjuvants on antibody production in mice. Immunology. (1968) 14:301–18.
44. Fan, L, Wu, T, Wang, Y, Hussain, A, Jiang, N, Gao, L, et al. Novel variants of infectious bursal disease virus can severely damage the bursa of fabricius of immunized chickens. Vet Microbiol. (2020) 240:108507. doi: 10.1016/j.vetmic.2019.108507
45. Hou, B, Wang, CY, Luo, ZB, and Shao, GQ. Commercial vaccines used in China do not protect against a novel infectious bursal disease virus variant isolated in Fujian. Vet Rec. (2022) 191:e1840. doi: 10.1002/vetr.1840
46. Lam, KM . Alteration of chicken heterophil and macrophage functions by the infectious bursal disease virus. Microb Pathogenesis. (1998) 25:147–55. doi: 10.1006/mpat.1998.0224
47. Lee, CC, Wu, CC, and Lin, TL. Role of chicken melanoma differentiation-associated gene 5 in induction and activation of innate and adaptive immune responses to infectious bursal disease virus in cultured macrophages. Arch Virol. (2015) 160:3021–35. doi: 10.1007/s00705-015-2612-y
48. He, X, Chen, G, Yang, L, Xuan, J, Long, H, and Wei, P. Role of naturally occurring genome segment reassortment in the pathogenicity of IBDV field isolates in three-yellow chickens. Avian Pathol. (2016) 45:178–86. doi: 10.1080/03079457.2016.1139687
49. Chen, G, He, X, Yang, L, and Wei, P. Antigenicity characterization of four representative natural reassortment IBDVs isolated from commercial three-yellow chickens from southern China reveals different subtypes co-prevalent in the field. Vet Microbiol. (2018) 219:183–9. doi: 10.1016/j.vetmic.2018.04.024
50. Chen, R, Chen, J, Xiang, Y, Chen, Y, Shen, W, Wang, W, et al. Differential modulation of innate antiviral profiles in the intestinal Lamina Propria cells of chickens infected with infectious bursal disease viruses of different virulence. Viruses. (2022) 14:393. doi: 10.3390/v14020393
51. Lucio, B, and Hitchner, SB. Response of susceptible versus immune chicks to killed, live-modified, and wild infectious bursal disease virus vaccines. Avian Dis. (1979) 23:1037–50. doi: 10.2307/1589620
52. Wang, W, Shi, J, Zhang, Y, Qiao, Y, Zuo, W, Wang, Z, et al. Genetic and pathogenic characterizations of a naturally occurring reassortant and homologous recombinant strain of the classical infectious bursal disease virus re-emerging in chickens in southern China. Front Microbiol. (2023) 14:1293072. doi: 10.3389/fmicb.2023.1293072
53. Skeeles, JK, Lukert, PD, Fletcher, OJ, and Diseases, JDLJA. Immunization studies with a cell-culture-adapted infectious bursal disease virus. Avian Dis. (1979) 23:456–65. doi: 10.2307/1589576
54. Wang, W, Liu, T, Chen, G, Ji, Z, Zhao, Z, Huang, Y, et al. Prokaryotic expression of VP2, VP1 and VP2-VP1 multimer genes of infectious bursal disease virus and its preliminary clinical application. Animal Husband Vet Med. (2020) 52:73–80. (in Chinese)
55. Shi, MY, Li, M, Wang, WW, Deng, QM, Li, QH, Gao, YL, et al. The emergence of a vv + MDV can break through the protections provided by the current vaccines. Viruses. (2020) 12:1048. doi: 10.3390/v12091048
56. He, X, Wang, W, Chen, G, Jiao, P, Ji, Z, Yang, L, et al. 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. (2019) 239:108458. doi: 10.1016/j.vetmic.2019.108458
57. Fan, L, Wu, T, Hussain, A, Gao, Y, Zeng, X, Wang, Y, et al. Novel variant strains of infectious bursal disease virus isolated in China. Vet Microbiol. (2019) 230:212–20. doi: 10.1016/j.vetmic.2019.01.023
58. Zhang, W, Wang, X, Gao, Y, and Qi, X. The Over-40-years-epidemic of infectious bursal disease virus in China. Viruses. (2022) 14:2253. doi: 10.3390/v14102253
59. Chen, Z, Lian, J, Liang, Z, Leng, M, Lin, W, and Chen, F. Characterization and pathogenicity of infectious bursal disease virus in southern China. Poultry Sci. (2022) 101:102018. doi: 10.1016/j.psj.2022.102018
60. Huang, Y, Shu, G, Huang, C, Han, J, Li, J, Chen, H, et al. Characterization and pathogenicity of a novel variant infectious bursal disease virus in China. Front Microbiol. (2022) 13:1039259. doi: 10.3389/fmicb.2022.1039259
61. Jiang, N, Wang, Y, Zhang, W, Niu, X, Huang, M, Gao, Y, et al. Genotyping and molecular characterization of infectious bursal disease virus identified in important poultry-raising areas of China during 2019 and 2020. Front Vet Sci. (2021) 8:759861. doi: 10.3389/fvets.2021.759861
62. Legnardi, M, Franzo, G, Tucciarone, CM, Koutoulis, K, and Cecchinato, M. Infectious bursal disease virus in Western Europe: the rise of reassortant strains as the dominant field threat. Avian Pathol. (2023) 52:25–35. doi: 10.1080/03079457.2022.2130172
63. Legnardi, M, Poletto, F, Alam, S, Cherfane, A, Le-Tallec, B, Franzo, G, et al. Molecular epidemiology of infectious bursal disease virus in the near east and Persian gulf regions. Avian Pathol. (2024) 53:56–67. doi: 10.1080/03079457.2023.2270531
64. Wang, Q, Hu, H, Chen, G, Liu, H, Wang, S, Xia, D, et al. Identification and assessment of pathogenicity of a naturally reassorted infectious bursal disease virus from Henan. China Poultry Sci. (2019) 98:6433–44. doi: 10.3382/ps/pez498
65. Leng, M, Bian, X, Chen, Y, Liang, Z, Lian, J, Chen, M, et al. The attenuated live vaccine strain W2512 provides protection against novel variant infectious bursal disease virus. Arch Virol. (2023) 168:200. doi: 10.1007/s00705-023-05828-0
66. Jackwood, DJ, and Sommer-Wagner, SE. Amino acids contributing to antigenic drift in the infectious bursal disease Birnavirus (IBDV). Virology. (2011) 409:33–7. doi: 10.1016/j.virol.2010.09.030
67. Letzel, T, Coulibaly, F, Rey, FA, Delmas, B, Jagt, E, van Loon, AA, et al. Molecular and structural bases for the antigenicity of VP2 of infectious bursal disease virus. J Virol. (2007) 81:12827–35. doi: 10.1128/JVI.01501-07
68. Vakharia, VN, He, J, Ahamed, B, and Snyder, DB. Molecular basis of antigenic variation in infectious bursal disease virus. Virus Res. (1994) 31:265–73. doi: 10.1016/0168-1702(94)90009-4
69. Boot, HJ, Hoekman, AJ, and Gielkens, AL. The enhanced virulence of very virulent infectious bursal disease virus is partly determined by its B-segment. Arch Virol. (2005) 150:137–44. doi: 10.1007/s00705-004-0405-9
Keywords: novel variant infectious bursal disease virus (nvIBDV), oil emulsion inactivated vaccine (OEV), booster immunization, immune protection, antigenicity and pathogenicity
Citation: Wang W, Huang Y, Zhang Y, Qiao Y, Shi J, Huang J, Huang T, Wei T, Mo M, He X and Wei P (2024) The complete protections induced by the oil emulsion vaccines of the novel variant infectious bursal disease viruses against the homologous challenges indicating the important roles of both VP2 and VP1 in the antigenicity and pathogenicity of the virus. Front. Vet. Sci. 11:1466099. doi: 10.3389/fvets.2024.1466099
Edited by:
Chunhe Wan, Fujian Academy of Agricultural Sciences, ChinaReviewed by:
Dongdong Yin, Anhui Academy of Agricultural Sciences (CAAS), ChinaJianzhu Liu, Shandong Agricultural University, China
Copyright © 2024 Wang, Huang, Zhang, Qiao, Shi, Huang, Huang, Wei, Mo, He and Wei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiumiao He, hxm@gxmzu.edu.cn; Ping Wei, weiping@gxu.edu.cn
†These authors have contributed equally to this work