Prevalence, Genotype Diversity, and Distinct Pathogenicity of 205 Gammacoronavirus Infectious Bronchitis Virus Isolates in China during 2019–2023

Gammacoronavirus infectious bronchitis virus (IBV) causes a highly contagious disease in chickens and seriously endangers the poultry industry. The emergence and co-circulation of diverse IBV serotypes and genotypes with distinct pathogenicity worldwide pose a serious challenge to the development of effective intervention measures. In this study, we report the epidemic trends of IBV in China from 2019 to 2023 and a comparative analysis on the antigenic characteristics and pathogenicity of isolates among major prevalent lineages. Phylogenetic and recombination analyses based on the nucleotide sequences of the spike (S) 1 gene clustered a total of 205 isolates into twelve distinct lineages, with GI-19 as a predominant lineage (61.77 ± 4.56%) exhibiting an overall increasing trend over the past five years, and demonstrated that a majority of the variants were derived from gene recombination events. Further characterization of the growth and pathogenic properties of six representative isolates from different lineages classified four out of the six isolates as nephropathogenic types with mortality rates in one-day-old SPF chickens varying from 20–60%, one as a respiratory type with weak virulence, and one as a naturally occurring avirulent strain. Taken together, our findings illuminate the epidemic trends, prevalence, recombination, and pathogenicity of current IBV strains in China, providing key information for further strengthening the surveillance and pathogenicity studies of IBV.


Introduction
Coronaviruses are enveloped, single-stranded and positive-sensed RNA viruses in the order Nidovirales [1], and are classified into four genera, alpha-, beta-, gamma-, and deltacoronavirus.Infectious bronchitis virus (IBV), discovered in North Dakota, USA, in 1931, is the prototype coronavirus in the gammacoronavirus genus and mainly infects chickens with a wide range of tissue tropism, causing damage to the respiratory, genitourinary, and digestive systems [2].Chickens aged from one to four weeks are most susceptible, with a mortality rate as high as 20% in severe cases [3].Furthermore, IBV infection can also result in a deterioration of the performance in broilers and laying hens, leading to a decrease in the food conversion ratio as well as reductions in both egg production and quality.Some IBV strains are capable of causing intestinal, glandular, and muscular diseases.Currently, vaccination remains the most critical measure for the prevention and control of this disease, with live attenuated vaccines being the most widely utilized.Vaccine strains commonly samples with specific primers [15], and the positive PCR products were submitted for sequencing (Sangon Biotech, CO., LTD, Shanghai, China).

Phylogenetic, Recombination, Deletion and Insertion Analysis of S1 Gene
The S1 gene of all isolates was sequenced and the complete S1 sequences used in this study are listed in Supplementary File S1.The nucleotide sequences of the S1 gene from all isolates and nine reference strains were used to construct a phylogenetic tree using the neighbor-joining method in MEGA version 7.0, and bootstrap values were determined based on 1000 replicates of the original data [16].The phylogenetic tree was further enhanced using the online software iTOL v6 (https://itol.embl.de/itol.cgi).Homology analysis of S1 nucleotide (S1-N) and amino acid (S1-A) sequences between representative strains of different lineages and vaccine strain H120 was carried out as described [17].Briefly, amino acid sequences of analyzed strains were aligned by the ClustalW method in Mega software 7.0, saved as Fasta files, and were subsequently analyzed by the ClustalW method in Megalign software 7.10.As cleavage of the S protein by furin facilitates viral entry, promotes syncytium formation, and enhances viral infectivity [18], the presence of a furin restriction site in IBV lineages was analyzed by multiple alignment and comparative analysis using the online software Weblogo 3 (http://weblogo.berkeley.edu/logo.cgi).

Growth Characteristics of Isolates in 9-Day-Old SPF Chicken Embryos
The virus solutions were ten-fold diluted in sterilized PBS, ranging from 10 −1 to 10 −8 in five replicates.The mortality rate and the presence of lesions (dwarfism and bleeding) in chicken embryos were observed and recorded, followed by calculation of EID 50 using the Reed-Muench method [19].
Nine-day-old SPF chicken embryos were inoculated with virus at a dose of ~10 3 EID 50 in duplicate.Following inoculation, embryos were cultured at 37 • C, and allantoic fluids (200 µL each) were collected at 0, 12, 24, 36, 48, and 60 h post-inoculation, respectively.Viral RNA copy numbers were determined by real-time qPCR (RT-qPCR) and calculated at each time point, and standard curves were constructed.

Pathogenicity Evaluation of Isolates
A total of 35 one-day-old chickens were randomly divided into seven groups, A, B, C, D, E, F, and G, with the control group designated as Group A and other groups infected and labeled as CZ_GI-7 (Group B), GZ_GI-13 (Group C), KP_GI-19 (Group D), KM-1_GI-22 (Group E), QY_GI-28 (Group F), and YC_GVI-1 (Group G).Five Chickens in each group were infected with ~10 5.5 EID 50 each of one isolate by the nasal-ocular route, respectively.The infected chicks were observed daily for clinical symptoms, such as tracheal rales, wheezing, nasal discharge, or death, for 7 days.The surviving chicks in each group were euthanized at 7 days post-challenge, and necropsies were performed immediately postmortem.
Trachea, lung, and kidney tissues were collected and fixed in 10% neutral buffered formalin for further histopathological and immunohistochemical analysis, respectively.The immunohistochemical staining was conducted following previously reported protocols [20], using a monoclonal antibody against IBV N as the primary antibody.Meanwhile, total RNAs were extracted from the trachea, lung, and kidney tissues harvested at 7 days post-infection, and the viral RNA copy numbers were quantified by RT-qPCR.

Statistical Analysis
Data were analyzed using the GraphPad software 9 package, and two-way analysis of variance (ANOVA) was used to analyze significant differences between the indicated samples and the respective control samples.Significance levels are presented by the p-value (ns, non-significant; * p < 0.05; ** p < 0.01; *** p < 0.001; and **** p < 0.0001).
Analysis of the prevalent patterns of these isolates in the past five years demonstrated that the prevalence of GI-1 remained relatively stable, prevalence trends of GI-19 and GI-13 isolates were consistently upward, and GI-7 was in a declining trend (Figure 2B).The prevalence of the GVI-1 variant slightly declined in the last two years, but was over 8.97% in the previous year, ranking second only to the predominant GI-19 strain (Figure 2B).
The furin motif sequences in the S protein from isolates with different lineages were then analyzed by Weblogo 3, and the following relatively conserved patterns were observed: RR(F/S)RR for GI-7, (R/H)R(S/R)RR for GI-13, ( for GI-22, and RRFRR for both GI-1 and GI-28 (Figure 2C).Due to the availability of either a limited number of isolates (fewer than three isolates) or insufficient isolates with the full-length S sequences (fewer than three), the furin motif sequences for the GI-7-like, GI-22-like, GVI-1, and GVI-1-like strains could not be subjected to statistical analysis and were excluded.These results demonstrate that, although there is no significant conservation in the furin motif sequences of different genotypic isolates, they are comparatively more conserved in the GI-1 and GI-28 lineages.Afte conducting a comparative analysis of the S1 sequences between the 205 isolates and referenc strains, an evolutionary tree was constructed by the neighbor-joining method using the MEGA software, and the resulting evolutionary tree was further beautified utilizing the iTOL onlin software.
Figure 1.Geographical distribution and evolution analysis of IBV isolates.(A) Regional distribution of 205 IBV isolates from 2019 to 2023 in China.(B) Evolution analysis of IBV isolates.After conducting a comparative analysis of the S1 sequences between the 205 isolates and reference strains, an evolutionary tree was constructed by the neighbor-joining method using the MEGA7 software, and the resulting evolutionary tree was further beautified utilizing the iTOL online software.
The furin motif sequences in the S protein from isolates with different lineages were then analyzed by Weblogo 3, and the following relatively conserved patterns were observed: RR(F/S)RR for GI-7, (R/H)R(S/R)RR for GI-13, (H/R)R(R/F)(R/K)R for GI-13-like,

Analysis of Recombination, Deletion, and Insertion on S1 Genes of Different Isolates
Recombination analysis was conducted on the S1 gene sequences of all 214 isolated and reference strains to identify putative parental sequences with statistical significance set at p-values <0.05 using RDP 4.0.Table 1 lists the breakpoint positions and specific p-values for the analysis of each recombination event in five variant-lineages.The results showed that 80% (4/5) variant-lineages were involved in potential recombination events (Table 1).Determination of the potential recombination events and breakpoints in these four variant-lineages by SimPlot analysis showed that the occurrence rates of gene recombination events were 83.33% (10/12), with lineages GI-19, GI-22, GI-28, GI-13-like, and GI-19-like identified as major parental strains, and GI-7, GI-13, GI-19, GI-22, GI-19-like, and GI-22-like as minor parental strains, respectively, in ten events (Figure 3A-E,G,J-L).It suggests that gene recombination may serve as a primary factor contributing to the emergence of these variants.
Deletion and insertion events of five variant-lineages were determined by Megalign analysis.Compared to the reference strains GI-7, GI-13, GI-19, GI-22, and GVI-1, variants belonging to the GI-7-like, GI-13-like, and GI-19-like lineages exhibited a significantly higher frequency (ranging from 67 to 400) of deletion and insertion events in the S1 gene (Figure 4).Among variants in the GI-7-like lineage, two distinct variants displayed a characteristic three-base deletion (AGG) at positions 285-287, and variants in the GI-22-like lineage featured a distinctive five-base deletion (GCCATT) at positions 70-75 (Figure 4).No specific mutations or insertions were detected in the S1 gene of variants in the GI-13-like, GI-19-like, and GVI-1-like lineages (Figure 4).These results suggest that deletion and insertion events are also significant factors contributing to the emergence of these variants.Deletion and insertion events of five variant-lineages were determined by Megalign analysis.Compared to the reference strains GI-7, GI-13, GI-19, GI-22, and GVI-1, variants  (Figure 4).Among variants in the GI-7-like lineage, two distinct variants displayed a characteristic three-base deletion (AGG) at positions 285-287, and variants in the GI-22like lineage featured a distinctive five-base deletion (GCCATT) at positions 70-75 (Figure 4).No specific mutations or insertions were detected in the S1 gene of variants in the GI-13-like, GI-19-like, and GVI-1-like lineages (Figure 4).These results suggest that deletion and insertion events are also significant factors contributing to the emergence of these variants.
The autopsy results of the dead chickens in the CZ_GI-7, KP_GI-19, KM-1_GI-22, and QY_GI-28 groups showed no visible pathological changes in the trachea and lungs, but with urate deposition in the kidneys of all infected birds.Among them, severe urate deposition and typical nephritis of swollen and pale kidneys were observed in chickens infected with CZ_GI-7, KP_GI-19, and KM-1_GI-22 (Figure 5D).As all chickens in the YC_GVI-1-and GZ_GI-13 groups showed no obvious clinical symptoms and death within 7 days after challenge, three chickens each in these two groups were randomly euthanized and necropsied at 7 days post-challenge, showing no obvious macroscopic lesions in lungs and kidneys, but mild tracheal hemorrhages in chickens infected with GI-13 (Figure 5D).

Histopathological Characterization of Chickens Infected with the Six Representative Isolates
Further evaluation of tissue damage by routine HE sections was conducted to investigate the pathological features of these isolates and to confirm the macroscopic lesions observed.Compared to the control group, chickens in the CZ_GI-7, KP_GI-19, and KM1_GI-22 groups displayed extensive necrosis and exfoliation of mucosal epithelial cells in the trachea (black arrows), along with moderate lamina propria edema, loose connective tissue arrangement, and a small amount of lymphocyte infiltration (blue arrows) (Figure 6).Chickens in the GZ_GI-13, QY_GVI-28, and YC_GVI-1 groups showed irregularly arranged tracheal cilia, slightly thickened tracheal mucosa (green arrows), and the presence of cell debris in the tracheal cavity (purple arrows) (Figure 6).Some necrosis of mucosal epithelial cells was also found in chickens in the GZ_GI-13 and QY_GVI-28 groups (red arrows) (Figure 6).

Determination of Viral Loads in Different Organs in Chickens Infected with the Six Representative Isolates
Viral loads in the trachea, lungs, and kidneys from three randomly selected chickens in each group were assessed at 7 days post-challenge (Figure 7A).The average viral loads in all three organs were first determined by RT-qPCR.With the exception of YC_GV-1, significant higher viral loads (p < 0.0001) were detected in chickens infected with the other five isolates than in the control group.The average viral loads in all three organs as well as in individual organs from chickens infected with YC_GVI-1 were at similar levels as in the control group, with no statistical difference between the infected and control groups, Chickens in the CZ_GI-7, GZ_GI-13, KP_GI-19, and KM1_GI-22 groups exhibited irregular arrangement of glomeruli and renal tubules in the kidneys, with a significant number of renal tubular epithelial cells that were necrotic and exfoliated (red arrows), and, to a lesser extent, localized lymphocyte infiltration (blue arrows) (Figure 6).Renal tubular expansion was also specifically detected in chickens infected with CZ_GI-7 (green arrows) (Figure 6).In contrast, chickens in the QY_GVI-28 and YC_GVI-1 groups showed minimal renal tubular epithelium degeneration (black arrows), along with focal lymphocyte infiltration within the interstitium only in the QY_GVI-28 group (blue arrows) (Figure 6).

Determination of Viral Loads in Different Organs in Chickens Infected with the Six Representative Isolates
Viral loads in the trachea, lungs, and kidneys from three randomly selected chickens in each group were assessed at 7 days post-challenge (Figure 7A).The average viral loads in all three organs were first determined by RT-qPCR.With the exception of YC_GV-1, significant higher viral loads (p < 0.0001) were detected in chickens infected with the other five isolates than in the control group.The average viral loads in all three organs as well as in individual organs from chickens infected with YC_GVI-1 were at similar levels as in the control group, with no statistical difference between the infected and control groups, indicating that this isolate is a naturally occurring avirulent strain, and it was excluded in the following analysis.indicating that this isolate is a naturally occurring avirulent strain, and it was excluded in the following analysis.The viral loads in the kidneys of the infected chickens were determined and calculated (Figure 7B).Very similar levels of average viral loads in the CZ_GI-7, KM1_GI-22, and QY_GI-28 groups as those in chickens infected with the predominant lineage representative strain KP_GI-19 were detected, suggesting that these three strains were nephropathogenic.The average viral loads in the GZ_GI-13 group were significantly lower than those in the KP_GI-19 group (**** p < 0.0001).KM1_GI-22 displayed a significantly higher viral load in the lungs (* p < 0.05), but no significant difference was observed in the trachea from chickens infected with the five isolates.
Comparison of the viral loads in the three organs infected with individual isolates was then conducted (Figure 7C), showing that viral loads in the kidneys of the QY_GI-28 group were significantly higher than those in both trachea (* p < 0.05) and lungs (**** p < 0.0001).Viral loads in both kidneys (*** p < 0.001) and trachea (** p < 0.01) of the KP_GI-19 The viral loads in the kidneys of the infected chickens were determined and calculated (Figure 7B).Very similar levels of average viral loads in the CZ_GI-7, KM1_GI-22, and QY_GI-28 groups as those in chickens infected with the predominant lineage representative strain KP_GI-19 were detected, suggesting that these three strains were nephropathogenic.The average viral loads in the GZ_GI-13 group were significantly lower than those in the KP_GI-19 group (**** p < 0.0001).KM1_GI-22 displayed a significantly higher viral load in the lungs (* p < 0.05), but no significant difference was observed in the trachea from chickens infected with the five isolates.
Comparison of the viral loads in the three organs infected with individual isolates was then conducted (Figure 7C), showing that viral loads in the kidneys of the QY_GI-28 group were significantly higher than those in both trachea (* p < 0.05) and lungs (**** p < 0.0001).Viral loads in both kidneys (*** p < 0.001) and trachea (** p < 0.01) of the KP_GI-19 group were significantly higher than in the lungs.No significant differences were observed in chickens infected with GZ_GI-13, KM1_GI-22, and QY_GVI-1, respectively, in all three organs.These findings further support the strong renal tropism of CZ_GI-7, QY_GI-28, and KP_GI-19.
Immunohistochemistry results (Figure 8) demonstrated that isolates CZ_GI-7, KM1_GI-22, and QY_GI-28 had a similar significant renal tropism as KP_GI-19, and also exhibited a broad tissue tropism to the trachea and lung tissues.GZ_GI-13 displayed a limited tropism to the trachea and lung tissues, with no associated virions detected in the kidneys.Careful observation also revealed that the main cell types primarily infected by these isolates were the tracheal mucosal epithelial cells, the bronchial and parabronchial epithelial cells in the lungs, and the renal tubular epithelial cells.These findings further support the classification of CZ_GI-7, KM1_GI-22, QY_GI-28, and KP_GI-19 as nephropathogenic strains, GZ_GI-13 a respiratory type, and YC_GVI-1 an avirulent strain.
Viruses 2024, 16, x FOR PEER REVIEW 15 of 19 limited tropism to the trachea and lung tissues, with no associated virions detected in the kidneys.Careful observation also revealed that the main cell types primarily infected by these isolates were the tracheal mucosal epithelial cells, the bronchial and parabronchial epithelial cells in the lungs, and the renal tubular epithelial cells.These findings further support the classification of CZ_GI-7, KM1_GI-22, QY_GI-28, and KP_GI-19 as nephropathogenic strains, GZ_GI-13 a respiratory type, and YC_GVI-1 an avirulent strain.and kidney autopsies of dead chickens or chickens randomly selected were conducted, and lung and kidney sections were immunohistochemically stained with a monoclonal antibody against IBV N protein.

Discussion
The co-existence of multiple genotypes of IBV strains with varying antigenicity and pathogenicity in many countries challenges the current surveillance and prevention strategies.The current classification scheme classifies IBV into seven genotypes (GI-GVII) and a total of 35 distinct lineages [8].QX-type (GI-19) strains are the most prevalent in six Asian countries, including Japan, Korea, India, Indonesia, Thailand, and China, and 793B-type (GI-13) strains are the second most prevalent lineage in these countries except Korea [21].However, during 2013-2015, 89.8% isolates were the QX-type (GI-19), TW I-type (GI-7), and 4/91-type (GI-13) [21], with the QX-type (GI-19) as the most common isolates (46.7%) in China [22].In 2016-2018, an epidemiological survey conducted in Henan, Hubei, and Hunan provinces showed that the prevalence of GI-19 (QX) was 48.2% in 56 isolates [23].The GVI-1 (TC07-2) lineage was first detected in China in 2007 [24], and has significantly increased in recent years.This lineage is currently prevalent in more than five countries, including China, Colombia, Japan, Korea, and Vietnam [25].In this study, we report the isolation and identification of seven predominant lineages and five variant-lineages circulating in China over the past five years.The GI-19 (QX or LX4) lineage has emerged as the most prevalent strain, accounting for 61.77 ± 4.56%, and the GVI-1 lineage as the second, accounting for 6.88 ± 1.84%, supporting that the GI-19 (QX-type) lineage has become the dominant strain in China.These findings are largely consistent with previous reports.It is worth noting that GVI-1 genotype (TC07-2) strains have been isolated and identified, in this study, in Sichuan in Southwest China, Liaoning in Northeast China, Jiangsu in East China, and Shandong in North China, apart from isolates in South China.This reveals a nationwide epidemic trend of GVI-1 in China.As the existing vaccines may provide poor immune protection against GVI-1 [26], these findings highlight the pressing needs to update IBV vaccine strains and develop novel broad-spectrum vaccines.
The coronavirus replicase gene contains an ExoN domain in non-structural protein (nsp) 14, playing a crucial role in proofreading and repair activities.Despite this unique feature, the estimated substitution rate of IBV genome at 10−4-10−5 substitutions/site/year remains remarkably high [27], ensuring the introduction of mutations and the emergence of multiple variants.A more important force that may drive the continuous evolution of IBV and emergence of variants would be through recombination [28].In recent years, isolation and identification of variants have been increasingly reported [16,29,30].In this study, potential recombination events showed an incidence rate of 66.67% (12/18).A novel recombinant strain of IBV was reported to emerge from three attenuated live vaccine strains (H120, 4/91, and LDT3-A) [31], but no such recombination events were found from any of the 205 isolates in this study.Variants of GI-7-like, GI-13-like, G-19-like, and GI-22-like are most likely generated via genetic recombination, deletion, and insertion.On the other hand, the GI-1 and GI-28 lineages exhibited remarkable stability, with no associated mutations detected across all isolates.Isolates in the two lineages also share highly conserved furin motif sequences (RRFRR).The sustained stability of the GI-1 lineage under the pressure of widespread use of vaccines based on the same GI-1 lineage (H120 and H52 strains) highlights the continuous usefulness of these vaccines for preventing and controlling outbreaks caused by IBV strains of the GI-1 lineage.Meanwhile, the relative stability of GI-28 and the emergence of multiple other lineages and variants would call for the development of novel bi-and multivalent vaccines.
Furin, a proprotein convertase, is located in the trans-Golgi network and activated by acidic PH.Furin can cleave precursor proteins with specific motifs to produce mature proteins with biological activity [32].Cleavage of the S protein by furin is thought to play an important role in promoting productive IBV infection in cultured cells [33].In this study, the furin motif sequences of each genotype were analyzed systematically, reveal-ing some relatively conserved patterns: RR(F/S)RR for GI-7, (R/H)R(S/R)RR for GI-13, (H/R)R(R/F)(R/K)R for GI-13-like, (H/R)R(R/H/F)RR for GI-19, (H/R)R(R/H)(R/K)R for GI-19-like, (H/R)R(R/F/L)(K/R)R for GI-22, and RRFRR for both GI-1 and GI-28.The furin motif, R-X-X-R↓X (X: any amino acid; ↓: cleavage site), displays a significant level of conservation.Interestingly, we find that the furin cleavage site in the S protein of three out of four nephropathogenic isolates (CZ_GI-7, KM1_GI-22, and QY_GI-28) is RRFRR, exhibiting renal tropism with urate deposition.The same motif in GZ_GI-13, an isolate with strong tracheal tropism, is RRSRR.As this motif contains physical properties critical for furin cleavage and S protein-mediated fusion efficiency [34], this sequence difference may play a role in determining the tissue tropism and pathogenicity of IBV.In addition, the relative abundance of furin proteases in different tissues and cells may also determine the IBV tissue and cell tropism.Cellular furin abundance has been strongly associated with the susceptibility of cells to IBV infection [33,35].Subtle variations in furin expression across tissues and cell types may account for differential susceptibility to IBV infection, despite its ubiquitous expression [36,37].Further studies are warranted to clarify these issues.
The virulence and mortality rates of prevalent IBV isolates in different lineages were reported to vary markedly.Some QX-like strains recently isolated in China exhibited a mortality rate ranging from 10% to 50% among infected chickens [38,39].As observed in this study, the up to 60% fatality rate in chickens infected with the KP_GI-19 strain is significantly higher than the approximate 20% fatality rate in chickens infected by other prevalent strains of different lineages (CZ_GI-7, KM1_GI-22, and QY_GI-28).Chickens infected with the two attenuated strains (GZ_GI-13 and YC_GVI-1) showed no mortality.However, it was recently reported that the pathogenicity of an isolate of the GVI-1 lineage (SX/2204) was significantly increased, with a mortality rate up to 60%, compared to previously isolated GVI-1 strains [25].It appears that the pathogenicity of isolates of the GI-19 and GVI-1 lineages exhibit an upward trend, as evidenced by this and several other recent studies.
Although the clinical and pathological outcomes of IBV infection are heavily influenced by the infecting strains, it is noteworthy that multiple physiological systems can be concurrently infected.The target organs susceptible to IBV infection encompass mainly the trachea, lung, kidney, oviduct, and gut [3,40].In this study, immunohistochemistry results were consistent with their viral loads in the three organs, and we showed that four out of six representative isolates (CZ_GI-7, KP_GI-19, KM1_GI-22, and QY_GI-28) could infect the tracheal, lung, and renal tissues of SPF chickens, with typical 'macular nephropathy' pathological changes in the kidney autopsy.However, the other two isolates (GZ_GI-13 and YC_GVI-1) exhibited no to mild clinical symptoms and pathological alterations.Our results support that GI-19 is more pathogenic than other genotypes, but variation in virulence was also noted in the same lineage.
In summary, this study has demonstrated that the GI-19 lineage is currently the predominant IBV strain, with high pathogenicity and an escalating mortality rate.The prevalence of this lineage and the second most prevalent strains of the GVI-1 lineage, along with their related variants, are steadily increasing, and the primary driving force for the emergence of variants is by genetic recombination.These observations urge the development of novel bi-and multivalent vaccines targeting the prevalent strains of the GI-19 and GVI-1 lineages.In addition, the possibility of developing the naturally occurring avirulent YC_GVI-1 strain as a live attenuated vaccine warrants further investigation to assess its immunogenicity and protection efficacy.

Figure 1 .
Figure 1.Geographical distribution and evolution analysis of IBV isolates.(A) Regional distributio of 205 IBV isolates from 2019 to 2023 in China.(B) Evolution analysis of IBV isolates.Afte conducting a comparative analysis of the S1 sequences between the 205 isolates and referenc strains, an evolutionary tree was constructed by the neighbor-joining method using the MEGA software, and the resulting evolutionary tree was further beautified utilizing the iTOL onlin software.

Figure 2 .
Figure 2. Analysis of proportion and prevalence of different lineages of IBV isolates, as well as the furin motif sequences.(A) Percentages of different lineages of 205 isolates from 2019 to 2023 in China.(B) The annual proportion of different lineages in 2019-2023.(C) Analysis of furin motif sequences in diverse lineages of isolates.The furin motif sequences were aligned using Mega sequencing, saved as Fasta files, and subsequently subjected to analysis using WebLogo 3 online software.The amino acid proportion at a certain position with sequence variation is represented by the relative height of the single amino acid letter.

Figure 2 .
Figure 2. Analysis of proportion and prevalence of different lineages of IBV isolates, as well as the furin motif sequences.(A) Percentages of different lineages of 205 isolates from 2019 to 2023 in China.(B) The annual proportion of different lineages in 2019-2023.(C) Analysis of furin motif sequences in diverse lineages of isolates.The furin motif sequences were aligned using Mega sequencing, saved as Fasta files, and subsequently subjected to analysis using WebLogo 3 online software.The amino acid proportion at a certain position with sequence variation is represented by the relative height of the single amino acid letter.

Viruses 2024 , 19 Figure 5 .
Figure 5. Homology, growth kinetics, and pathogenicity analyses of six representative isolates.(A) Two software programs, MEGA 7.0 and Megalign, were successively used to analyze the homology of S1 sequences between the vaccine strain H120 and six representative isolates with different lineages.S1-N and S1-A represent homology of nucleic acid and amino acid sequences, respectively.(B) Growth curves of six representative isolates in chicken embryos.Nine-day-old chicken embryos were infected at an EID50 ~10 3, respectively.The copy numbers of IBV in allantoic fluids were determined and calculated by RT-qPCR.(C) Surviving curves of six representative isolates in 1-dayold SPF chickens.Chickens were infected with ~10 5.5 EID50 each of the six representative isolates by the nasal-ocular route, respectively.The number of deaths occurring within 7 days after the challenge was recorded, and a survival curve was generated using the Graphpad Prism 9 program.(D) Examination of pathological lesions in the trachea, lung, and kidney autopsies.At 7 days postchallenge, the trachea, lung, and kidney autopsies of dead chickens or chickens randomly selected

Figure 5 .
Figure 5. Homology, growth kinetics, and pathogenicity analyses of six representative isolates.(A) Two software programs, MEGA 7.0 and Megalign, were successively used to analyze the homology of S1 sequences between the vaccine strain H120 and six representative isolates with different lineages.S1-N and S1-A represent homology of nucleic acid and amino acid sequences, respectively.(B) Growth curves of six representative isolates in chicken embryos.Nine-day-old chicken embryos were infected at an EID 50 ~10 3, respectively.The copy numbers of IBV in allantoic fluids were determined and calculated by RT-qPCR.(C) Surviving curves of six representative isolates in 1-day-old SPF chickens.Chickens were infected with ~10 5.5 EID 50 each of the six representative isolates by the nasal-ocular route, respectively.The number of deaths occurring within 7 days after the challenge was recorded, and a survival curve was generated using the Graphpad Prism 9 program.(D) Examination of

Figure 6 .
Figure 6.Histopathological examination of the trachea, lungs, and kidneys in chickens infected with each of the six representative isolates.At 7 days post-challenge, the trachea, lung, and kidney autopsies of dead chickens or chickens randomly selected from each experimental group were conducted and examined by microscopy.

Figure 6 .
Figure 6.Histopathological examination of the trachea, lungs, and kidneys in chickens infected with each of the six representative isolates.At 7 days post-challenge, the trachea, lung, and kidney autopsies of dead chickens or chickens randomly selected from each experimental group were conducted and examined by microscopy.

Figure 7 .
Figure 7. Viral loads in the trachea, lung, and kidney tissues from chickens infected with six representative isolates.At 7 days post-challenge, three chickens from each experimental group were randomly selected, the trachea, lung, and kidney tissues were collected, and the viral RNA copy numbers were determined and calculated by RT-qPCR.Significance levels are presented by the pvalue (ns, non-significant; * p < 0.05; ** p < 0.01; *** p < 0.001; and **** p < 0.0001).(A) Statistical analysis of the total viral loads in the three tissues within each challenging group and the control group.(B) Statistical analysis of viral load variations in different tissues from chickens challenged with six representative isolates.(C) Statistical analysis of viral load variations in the same organs from chickens challenged with the other five isolates and the predominant lineage strain KP_GI-19, respectively.

Figure 7 .
Figure 7. Viral loads in the trachea, lung, and kidney tissues from chickens infected with six representative isolates.At 7 days post-challenge, three chickens from each experimental group were randomly selected, the trachea, lung, and kidney tissues were collected, and the viral RNA copy numbers were determined and calculated by RT-qPCR.Significance levels are presented by the p-value (ns, non-significant; * p < 0.05; ** p < 0.01; *** p < 0.001; and **** p < 0.0001).(A) Statistical analysis of the total viral loads in the three tissues within each challenging group and the control group.(B) Statistical analysis of viral load variations in different tissues from chickens challenged with six representative isolates.(C) Statistical analysis of viral load variations in the same organs from chickens challenged with the other five isolates and the predominant lineage strain KP_GI-19, respectively.

Figure 8 .
Figure 8. Immunohistochemical examination of IBV replication in the trachea, lungs, and kidneys from chickens infected with six representative isolates.At 7 days post-challenge, the trachea, lung, and kidney autopsies of dead chickens or chickens randomly selected were conducted, and lung and kidney sections were immunohistochemically stained with a monoclonal antibody against IBV N protein.

Figure 8 .
Figure 8. Immunohistochemical examination of IBV replication in the trachea, lungs, and kidneys from chickens infected with six representative isolates.At 7 days post-challenge, the trachea, lung,

Table 2 .
EID 50 of the six representative isolates.