key: cord-0711753-7cqu5oea
authors: Najimudeen, Shahnas M.; Hassan, Mohamed S. H.; Goldsmith, Dayna; Ojkic, Davor; Cork, Susan C.; Boulianne, Martine; Abdul-Careem, Mohamed Faizal
title: Molecular Characterization of 4/91 Infectious Bronchitis Virus Leading to Studies of Pathogenesis and Host Responses in Laying Hens
date: 2021-05-19
journal: Pathogens
DOI: 10.3390/pathogens10050624
sha: 72ef9d12288b24dd0d3fe3834c78ffd26b346976
doc_id: 711753
cord_uid: 7cqu5oea

Infectious bronchitis virus (IBV) initially establishes the infection in the respiratory tract and then spreads to other tissues depending on its virulence. During 2011–2018, the 4/91 IBV strain was isolated from poultry flocks affected by decreased egg production and quality in Eastern Canada. One of the Canadian 4/91 IBV isolates, IBV/Ck/Can/17-038913, was propagated in embryonated chicken eggs and molecularly characterized using whole genome sequencing. An in vivo study in laying hens was conducted to observe if IBV/Ck/Can/17-038913 isolate affects the egg production and quality. Hens were infected with IBV/Ck/Can/17-038913 isolate during the peak of egg lay, using a standard dose and routes maintaining uninfected controls. Oropharyngeal and cloacal swabs were collected at predetermined time points for the quantification of IBV genome loads. At 6 and 10 days post-infection, hens were euthanized to observe the lesions in various organs and collect blood and tissue samples for the quantification of antibody response and IBV genome loads, respectively. Egg production was not impacted during the first 10 days following infection. No gross lesions were observed in the tissues of the infected birds. The IBV genome was quantified in swabs, trachea, lung, proventriculus, cecal tonsils, kidney, and reproductive tissues. The serum antibody response against IBV was quantified in infected hens. In addition, histological changes, and recruitment of immune cells, such as macrophages and T cell subsets in kidney tissues, were measured. Overall, data show that IBV/Ck/Can/17-038913 isolate is not associated with egg production issues in laying hens infected at the peak of lay, while it demonstrates various tissue tropism, including kidney, where histopathological lesions and immune cell recruitments were evident.

Infectious bronchitis (IB) is a multi-system disease in chickens caused by infectious bronchitis virus (IBV), which is a Gammacoronavirus that belongs to the family Coronaviridae and order Nidovirales [1] . The disease was first recorded in the United States in the early 1930s [2] . Since then, outbreaks of IB in layer and broiler flocks have been reported globally with considerable economic losses due to poor growth performances in broilers, and egg production and quality defects in layers [3] . Although vaccine-mediated control measures are extensively employed, IB is still a major issue in the poultry industry in Canada and globally [4] [5] [6] [7] [8] .

The general time-reversible model with gamma distribution and 1000 bootstrap replicates was used for the construction of a maximum likelihood phylogenetic tree of the 87 complete S1 genes. IBV/Ck/Can/17-038913 isolate clustered within the GI-13 lineage along with IBVs isolated from different parts around the world that are known as 793B type (also known as 4/91 and CR88; Figure 1 ). The complete genome of IBV/Ck/Can/17-038913 isolate consists of 27,472 nucleotides, excluding the poly (A) tail. Twelve ORFs (5 -1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-3 ) were predicted (Table 1) . Sequence comparison analysis demonstrated that the genomic nucleotide identity between IBV/Ck/Can/17-038913 isolate and different IBV strains ranged from 86.1% (TW2575/98) to 92.0% (JMK) while it showed 87.2% identity with turkey coronavirus (TCoV) isolate (TCoV/IN-517/94). Among the different ORFs of the genome, the S gene was the most variable (70.9-96.7%), whereas it showed only 48.3% identity with TCoV (TCoV/IN-517/94). On the other hand, the 5b accessory protein was the most conserved (91.6-98.4%). IBV/Ck/Can/17-038913 isolate revealed higher sequence identities (93. .0%) of the 3b, E, M, 4b, 4c, 5a and 5b ORFs with the Canadian Massachusetts (Mass) type isolate (IBV_SES_15SK-02). The 3a and N had the highest nucleotide identity with the Conn46 1996 strain at 94.8% and 96.7%, respectively. ORF 1a had the highest nucleotide identity (93.5%) with TCoV/IN-517/94. The IBVs of 793B type, 4/91 vaccine, ck/CH/LSD/110857 and IBVUkr27-11, showed the highest identity (95.1-96.7%) with the S gene (Table 2) . Phylogenetic trees constructed based on different genomic regions showed different topologies and confirmed the results of the sequence comparison. The IBV/Ck/Can/17-038913 isolate did not cluster with 793B type strains outside of the S1 gene, whereas it clustered with the Mass type (1a and E genes) or Mass and non-Mass type strains (1b, M and N genes; Figure 2 ). 

The recombination analysis identified two potential recombination events by the seven models employed by Recombination Detection Program (RDP) 4 software. The first recombinant region is located between 19,284 and 23,701 nucleotides that are found in the extreme 3 end of the replicase gene and the S gene. In this case, the minor parent belongs to the ck/CH/LSD/110857 strain (KP118885), and the major parent to the TCoV/IN-517/94 strain (GQ427175; Figure 3a ). The second recombinant region located between 24,125 and 25,927 nucleotides, containing E, 4 and 5 genes with a minor parent, was inferred as the IBV_SES_15SK-02 strain (MH539772), and the putative major parent was determined as the Conn46 1996 strain (FJ904716; Figure 3b ). 

Embryo passage resulted in IBV induced lesions such as stunting and curling.

A significant drop in egg production or changes in egg quality, such as egg shape, albumin thickness and eggshell structure, compared to the uninfected control group, were not observed during the study period (p = 0.8630; Supplementary Figure S1 ). No clinical signs were observed in the infected or the control group. At the necropsy examination performed at 6 and 10 dpi, no gross lesions were noted in any tissues of the infected or uninfected control hens.

Viral genome loads in oropharyngeal (OP) and cloacal (CL) swabs and tissue samples including trachea, lung, kidney, cecal tonsils, proventriculus, and reproductive tract (uterus, magnum and isthmus) in the infected group were quantified using SYBR Green qRT-PCR assay. A significant difference in the IBV genome load was observed between 3 and 10 dpi (p = 0.0452) in OP swabs ( Figure 4a ) and 3 and 6 dpi in CL swabs (p = 0.0472; Figure 4b ).

The viral loads in the tissues of different organs at 6 and 10 dpi are shown in Figure 4c ,d, respectively. At 6 dpi, the highest genome load was found in the trachea and isthmus, respectively, while the lowest amount was found in the lung. However, the highest viral load at 10 dpi was observed in the cecal tonsils, which is an organ of IBV persistence, followed by the lung. In contrast, the reproductive tract had the lowest amount of IBV genome load compared to other tissues. . IBV genome loads in swabs and tissues following IBV/Ck/Can/17-038913 isolate infection. The mean IBV genome loads in 100 ng of extracted RNA from OP swabs (a), CL swabs (b) at 3, 6, 10 dpi and from various tissues of infected hens at 6 dpi (c) and 10 dpi (d) are shown here. The IBV genome loads in swabs and tissues were compared using the Kruskal-Wallis test followed by Dunn's multiple comparison. The error bars represent the SD. Significance: * p < 0.05, ** p < 0.01, *** p < 0.001.

All infected hens developed anti-IBV antibody (Ab) at 10 dpi and two birds were positive at 6 dpi, whereas the hens from the control group remained negative at the two time points. A significant difference in anti-IBV Ab titer was observed between sera collected from the infected group at 6 and 10 dpi (p = 0.0043; Figure 5 ). 

A significant recruitment of CD8 + T cells, CD4 + T cells and macrophages at both 6 and 10 dpi was observed in the kidney tissues of the infected hens when compared to uninfected controls (Figure 6a -c). At both 6 and 10 dpi time points, the number of CD8 + cells exceeded that of CD4 + cells.

Although no gross lesions were observed during necropsy examination, histological changes were present in renal tissues of the infected hens at 10 dpi with no significant changes at 6 dpi ( Figure 7 ). Mild to moderate dilation of collecting ducts was observed with an increased luminal accumulation of urate spherules in all infected birds at 10 dpi ( Table 3 ). The dilated tubules were most prominent in the medullary area with occasional sloughing of epithelial cells and rare luminal hyper eosinophilic necrotic cellular debris. Mononuclear cell infiltration was observed in two birds of the infected group. Around the collecting ducts and ureteral branches, the subepithelial connective tissue was infiltrated by moderate numbers of lymphocytes and macrophages along with fewer plasma cells and leukocytes with large eosinophilic cytoplasmic vacuoles (suspected Mott cells). Small numbers of lymphocytes were observed extending into the lining epithelium and there was rare associated karyorrhectic debris. Mild collecting duct dilation was seen in one uninfected hen and no other notable histological changes were observed among uninfected kidney tissues. No significant histologic abnormalities within any of the examined kidney sections (infected and control) were seen at 6 dpi. Table 3 . Histopathological changes observed in the kidney of IBV/Ck/Can/17-038913 isolate infected and control hens (10 dpi).

Dilation of collecting duct 6/6 1/6 Increased urate spherules 6/6 0/6 Mononuclear cell infiltration 2/6 0/6 Sloughed epithelial cells in the lumen of collecting duct 4/6 0/6 Hyper eosinophilic necrotic cellular debris 2/6 0/6

The recovery of a significant percentage of 4/91 IBV isolates from Eastern Canadian poultry flocks, which had problems with egg production in the recent past [4] , led us to investigate the impact of Canadian 4/91 IBV infection in laying hens. The present study yielded several findings. First, based on S1 gene phylogenetic analysis, our IBV/Ck/Can/17-038913 isolate grouped within the GI-13 lineage along with 793B IBVs (also known as 4/91 and CR88). Second, the IBV/Ck/Can/17-038913 isolate leads to subclinical infection in laying hens with no impact on egg production. Third, the IBV genome load data suggest that IBV/Ck/Can/17-038913 isolate has a very broad tissue tropism. Fourth, when we focused our investigation on the kidney, we found that IBV/Ck/Can/17-038913 isolate induces histological changes in the kidney with significant immune cell recruitment.

In Canada, several IBV variants were detected in Canadian poultry operations over the last two decades [4, 37] . However, the published data that involved detailed molecular characterization of different Canadian IBV strains are scarce. Amongst the IBV strains isolated in Canada, there are only seven whole genome sequences available through GenBank [38, 39] . Hence, in the present study, the complete genome of the Canadian IBV isolate, IBV/Ck/Can/17-038913, was characterized to further understand the genetic backgrounds of IBV evolution in Canada. The results of comparative genomics showed high nucleotide identity between our isolate and TCoV/IN-517/94 in the replicase gene.

This finding was not surprising, whereas this TCoV was shown to be originated as a consequence of recombination that swapped the spike gene of a Mass IBV strain [40] . Additionally, another Mass type isolate, IBV-SES-15SK-02, shared high nucleotide sequence identities with our isolate in genes 4 and 5 and most of gene 3. The role of Mass type strains in the evolution of our isolate was further elucidated by identifying IBV-SES-15SK-02 isolate as a minor parent for a recombinant region in the 3 end of the genome. Over the last two decades, a marked proportion of the IBVs detected in the Canadian poultry operations has shown to be of Mass type [4, 37] . Despite the high similarity of different genomic fragments with Mass type strains, another recombination event occurred during the origin and evolution of the IBV/Ck/Can/17-038913 isolate and caused it to cluster with the 793/B group. Within the S gene, natural recombination events associated with shifts in the serotype and/or genotype of the circulating IBVs were described in the past [39, 41] .

The 4/91 IBV strain was introduced into Eastern Canada in mid-2011, and around 50% of IBV isolates that was detected from the flocks affected with respiratory and egg production problems during 2012-2013 were identified as 4/91 type [4] . According to Animal Health Laboratory data, the recovery of 4/91 IBVs has declined towards 2018 in Canada, although they are still a significant problem in other countries [7, 27, 36, 42, 43] .

Although 4/91 IBV has shown to be a highly pathogenic virus strain causing severe clinical signs and mortality [42, 44] , some studies have reported only mild respiratory signs [21, 45] . In the present study, we observed neither clinical signs nor gross lesions in any of the tissues throughout the experimental period. It is possible that this variant is less virulent in laying hens. However, the results could be different in field conditions as the severity of the disease is also influenced by other factors such as coinfections, breed, and immunity of the host [18, 21, 24, 25] , which are controlled in our experimental design.

Since 4/91 IBV was associated with kidney tissues [21, 27, 46] , we focused our studies on the kidney. Although gross lesions were not observed in IBV/Ck/Can/17-038913 isolate infected kidneys, histopathological changes were seen. At 10 dpi, infected kidney tissues were comprised of tubular dilation with luminal urate accumulation, desquamation, and mononuclear cell infiltration with rare epithelial necrosis. These findings were similar to the histological lesions caused by nephropathogenic IBV strains, such as QX-like, Aust T and TW-like IBV [18, 27, 47] . Moreover, nephropathogenic lesions and mortality induced by nephropathogenic IBV strains are more prominent in young chickens than adults [26, 29] . Therefore, our 4/91 IBV isolate may be a nephropathogenic strain and this requires to be confirmed by infecting young chickens.

In agreement with the previous literature [21, 27, 45] , the Canadian 4/91 IBV isolate, IBV/Ck/Can/17-038913, had a broad tissue tropism, including cecal tonsils, trachea, lung, kidney, proventriculus, and reproductive tract. At 10 dpi, the highest IBV genome load was observed in cecal tonsils, which represent an organ of IBV persistence [16, 43, 44, 48, 49] , followed by the trachea. In agreement with a study by Raj and Jones [21] , our data show high IBV genome loads in kidney tissues at 6 dpi followed by 10 dpi.

The egg production issues caused by IBV can vary from only mild discoloration of the eggshell to severe production failure with the deterioration of egg quality [31, 50] . IBV infection in laying hens can result in a 20-90% drop in egg production [51] . However, a significant drop in egg production or deterioration of the internal or external quality of the eggs were not observed in our study. A reduction in the length of the oviduct and regressed ovaries, which are features of damage caused by IBV [24, 52] , were also not observed in all 12 IBV/Ck/Can/17-038913 isolate infected hens.

In summary, based on S1 gene phylogenetic analysis and comparative genomics, the IBV/Ck/Can/17-038913 isolate was assigned to a GI-13 lineage and showed varying similarities with different reference strains across its different ORFs. In laying hens, the Canadian 4/91 IBV isolate, IBV/Ck/Can/17-038913, had a very broad tissue tropism including viral replication in the reproductive tract, while it did not cause egg production or egg quality problems with the dose of infection used in the current study. However, the IBV/Ck/Can/17-038913 isolate led to histological changes in the kidney with significant recruitment of T cell subsets and macrophages. Further studies are necessary to determine the impact of IBV/Ck/Can/17-038913 isolate infection in young chickens on renal and reproductive tissues.

Specific pathogen-free (SPF) eggs for virus propagation and titration and 25-weeks-old SPF layer chickens for experimental infection were purchased from the Canadian Food Inspection Agency (CFIA), Ottawa. Ethical approval for the proposed work was obtained from the Health Science Animal Care Committee (HSACC) of the University of Calgary (Protocol number: AC19-0011).

The Canadian 4/91 IBV isolate, designated as IBV/Ck/Can/17-038913, used in the present study was recovered from a tracheal swab collected from a 2-week-old layer flock presented with respiratory manifestations in Ontario, Canada. The archived swab was obtained from Animal Health Laboratory, University of Guelph, Guelph, ON, Canada. Virus isolation in 9-11-days-old SPF embryonated chicken eggs was carried out as previously described [39] . The virus was titrated using 9-days-old SPF eggs and the titer was determined to be 1 × 10 6 EID 50 / mL.

After three passages in embryonated chicken eggs, RNA was extracted from allantoic fluid using Trizol LS ® reagent (Ambion, Invitrogen Canada Inc., Burlington, ON, Canada) according to the manufacturer's protocol. The genome was transcribed to complimentary deoxyribonucleic acid (cDNA) using random primers (High Capacity Reverse Transcription Kit, Applied Biosystems™, Invitrogen Canada Inc., Burlington, ON, Canada) according to the manufacturer's instructions. Whole genome sequencing on the cDNA library was performed at the Faculty of Veterinary Medicine, University of Montreal, Montreal QC, Canada, using a Miseq platform (Illumina corp, San Diego, CA, USA). The genome was deposited in GenBank under the accession number MT665806.

The S1 gene of the current isolate, in addition to the S1 gene sequences of 86 IBV reference strains of different genotypes (Supplementary Table S1 ) [53] , was selected and used for phylogenetic analysis. Clustal Omega [54] and RAxML [55] plug-ins in the Geneious ® software v10.2.6 (https://www.geneious.com/) were used for multiple sequence alignments and phylogenetic tree construction, respectively. The coding sequences and ORFs prediction was carried out in https://www.ncbi.nlm.nih.gov/orffinder/. Both the complete genome and nucleotide sequences of different ORFs of IBV/Ck/Can/17-038913 isolate were used for BLAST searching of the National Center for Biotechnology Information (NCBI) database and were compared to 23 reference IBV sequences of North America (USA and Canada), Asia (China, Korea and Taiwan), Africa (Nigeria and Egypt), Europe (Sweden, Poland, Ukraine, Belgium and Italy) and Australia. Phylogenetic analyses of the nucleotide sequences of 1a, 1b, E, M and N genes were also performed with the selected 23 reference IBV sequences.

The complete genome sequence of the Canadian 4/91 IBV isolate, IBV/Ck/Can/17-038913 and 23 reference IBV sequences mentioned in the previous section were aligned using multiple alignments with fast Fourier transformation (MAFFT) [56] and used for recombination analysis conducted in the RDP4 (RDP v.4.101) [57] . Specific models (RDP, Genecov, Bootscan, Maxchi, Chimaera, Siscan, and 3Seq) embedded in RDP4 were used to identify putative recombination events in the genome of our isolate using default settings for most models. Changes from default settings were the window size for RDP (60) and Bootscan and SiScan (500) to reduce masking of recombination signals. Only the recombination events identified by at least five of these models were considered.

Twenty-four 25-weeks-old SPF hens (at the peak of egg lay) were purchased from CFIA, Ottawa, and allowed to acclimatize for 19 days. Hens were placed in Animal Resource Center (ARC), Foothill campus, University of Calgary. Feed and water were provided ad libitum. The laying hens were kept under a daily light cycle of 16-h light:

At the end of the acclimatization period, 12 hens were infected with IBV/Ck/Can/17-038913 isolate at a dose of 1× 10 6 EID 50 per chicken under isoflurane anesthesia. We used the standard dose of IBV used in adult chickens as was performed previously [58, 59] . A total volume of 1 mL was inoculated via intraocular (100 µL), intranasal (100 µL) and intratracheal (800 µL) routes. The infected group was placed in high containment poultry isolators in Prion Virology facility (PVF). The remaining 12 hens were mock-infected with phosphate-buffered saline (PBS) and acted as uninfected controls. Egg production, egg quality (egg weight, shape index, albumin thickness) and clinical signs were monitored daily. Six hens from the infected and control groups were euthanized at 6 and10 days post-infection (dpi) to observe the gross lesions in various organs and to collect samples.

OP and CL swabs were collected at 3, 6 and 10 dpi into tubes containing viral transport medium (Puritan ® UniTranz-RT ® Media Transport Systems, VWR, Edmonton, AB, Canada). Tissue samples of trachea, lung, cecal tonsils, kidney, oviduct (magnum, isthmus, uterus), and proventriculus were collected individually in RNA Save (Biological Industries, FroggaBio, Toronto ON, Canada), and preserved at −20 • C, for virus genome load quantification by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). In addition, kidney samples were snap frozen in Optimum Cutting Temperature (OCT) medium (Tissue-Tek ® , Sakura Finetek USA Inc., Torrance, CA, USA) for immunofluorescence studies and fixed in 10% neutral buffered formalin (VWR International, West Chester, PA, USA) for examining histological changes. Blood samples were collected into plain tubes on the day of euthanasia (6 and 10 dpi) to monitor the anti-IBV Ab response. The blood samples were kept at room temperature for 30 min, centrifuged at 2000× g for 10 min at 4 • C, the serum was separated, and aliquots were stored at −20 • C.

The tubes containing OP and CL swabs were vortexed and 250 µL of the medium was used to extract the RNA using Trizol LS ® reagent (Ambion, Invitrogen Canada Inc., Burlington, ON, Canada) according to manufacturer's guidelines. Total RNA was also extracted from the trachea, lung, kidney, reproductive tract (uterus, isthmus and magnum), cecal tonsils, and proventriculus using Trizol reagent (Ambion, Invitrogen Canada Inc., Burlington, ON, Canada) following the manufacturer's protocol. The amount of the extracted RNA was quantified using Nanodrop 1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA) at 260 nm wavelength. The quality of the RNA was determined based on the 260/280 absorbance ratio. Approximately 1000 ng from swab samples and 2000 ng of tissue samples were used to synthesize cDNA using RT random primers (High capacity cDNA reverse transcriptase kit, Invitrogen Life Technologies, Carlsbad, CA, USA) according to manufacturer's guidelines.

The qRT-PCR assay was carried out using a CFX 96-c1000 Thermocycler (Bio-Rad Laboratories, Mississauga, ON, Canada) to quantify the IBV genome loads in cDNA samples. The assay was performed using Fast SYBR ® Green Master Mix (Invitrogen, Burlington, ON, Canada) maintaining the final reaction volume at 20 µL. Each reaction volume consisted of 10 µL of SYBR Green master mix, 100 ng of respective samples, 0.5 µL of forward Fw-5 GACGGAGGACCTGATGGTAA-3 and reverse Re-5 CCCTTCTTCTGCTGATCCTG-3 specific primers targeting the conserved IBV N gene [60] . Thermal cycling conditions were initial denaturation at 95 • C for 20 seconds (s); followed by 40 cycles of final denaturation at 95 • C for 3 s, and annealing at 60 • C for 30 s, a melting curve analysis was conducted between 65 • C to 95 • C with an increment of 0.5 • C at every 5 s.

The kidney samples fixed in 10% neutral buffered formalin were sent to the Histopathology Diagnostic Services Unit at the University of Calgary, Faculty of Veterinary Medicine for further processing, sectioning, and hematoxylin-eosin (H&E) staining.

Anti-IBV Abs in serum samples (1:500 dilution) were quantified using IDEXX ELISA kit (IDEXX Laboratories, Westbrook, ME, USA) according to manufacturer's guidelines. This ELISA kit does not differentiate different Ab isotypes. The absorbance values were taken at 650 nm using the SpectraMax M2 microplate reader (Molecular Devices, Sunnyvale, NS, Canada) and the titers were calculated according to the manufacturer's instruction. Titers greater than 396 were considered positive for anti-IBV Abs.

OCT preserved frozen tissues were cut into sections of 5 µm thickness, adhered to positively charged slides (VWR, Missisauga, ON, Canada), and preserved at −20 • C until stained for CD4 + T cells, CD8 + T cells and macrophages. On the day of immunostaining, the sections were air-dried for 20 min (min) at room temperature and fixed with cold acetone for 5 min. The tissues were blocked with 5% goat serum diluted in Trizma buffered saline (TBS) buffer (Trizma base: 2.42 g; NaCl: 8 g in 1 L of distilled water; pH 7.6) and the mouse monoclonal Ab specific for chicken macrophages/monocytes (KUL01, Southern Biotech, Birmingham, AL, USA) or CD8α (CT-8, Southern Biotech, Birmingham, Alabama, USA) were used in a 1:200 dilution in 5% goat serum. The secondary Ab, goat anti-mouse IgG (H + L) conjugated with Dylight ® 550 (red fluorescence; Bethyl Laboratories Inc., Montgomery, TX, USA) was used in 1:500 dilution in 5% goat serum. The sections were incubated with primary Ab for 30 min and with secondary Ab for 60 min.

Immunostaining of CD4 + T cells consists of an additional 15 min of blocking step with avidin and biotin (Vector Laboratories, Inc., Burlingame, CA, USA) before blocking with 5% goat serum in TBS for 30 min. Each step was followed by 3 min of washing twice with Trisbuffered saline containing 0.1% Tween 20 (TBST) and once with PBS. Following blocking, the sections were incubated with 1:200 dilution of mouse monoclonal Ab specific to chicken CD4 + T cells (CT-4, Southern Biotech, Birmingham, Alabama, USA) in blocking buffer and incubated for 30 min. Next, the tissues were incubated for 30 min with 1:250 dilution of biotinylated goat anti-mouse IgG (H + L; Southern Biotech, Birmingham, AL, USA) and 15:1000 dilution of DyLight ® 488 (green fluorescence) streptavidin (Vector Laboratories Inc., Burlingame, CA, USA), respectively. Each incubation step was followed by 5 min of washing with TBST buffer twice and once with PBS. Finally, the slides were mounted with 4 ,6-diamidino-2-phenylindole dihydrochloride (DAPI; DAPI, Vector Laboratories Inc., Burlingame, CA, USA; Blue fluorescence) and cover slipped.

The number of IBV genome copies in 100 ng of RNA was calculated based on a standard curve obtained by serial dilution of plasmids. For quantification of immune cells, 5 areas (x20 magnification) of each tissue with the highest signals were captured using an epifluorescent microscope (Olympus IX51, Center Valley, PA, USA) along with applicable nuclear stained (DAPI) areas. The density of the cells in each area was quantified using Image-J software (National Institute of Health, Bethesda, MD, USA).

The difference in egg production between the infected and control groups was analyzed using a Mann-Whitney U test. The Mann-Whitney U test was also used to compare the anti-IBV Ab responses in serum of the infected group at two time points, 6 and 10 dpi. IBV genome load in OP, CL swabs and tissues were compared using the Kruskal-Wallis test followed by Dunn's multiple comparison. Differences of CD4 + T cells, CD8 + T cells, and macrophages at 6 and 10 dpi were detected using the Mann-Whitney U test. All the statistical analysis was performed in GraphPad Prism Software (GraphPad Prism 9 Software, San Diego, CA, USA).

The following are available online at https://www.mdpi.com/article/10 .3390/pathogens10050624/s1, Figure S1 : Egg production following infection with IBV/Ck/Can/17-038913 isolate, Table S1 : List of reference IBV sequences used in the S1 sequence analysis. 

Data Availability Statement: The datasets used and/or analyzed within the frame of the study can be provided by the corresponding author upon reasonable request.

Coronaviruses in Avian Species-Review with Focus on Epidemiology and Diagnosis in Wild Birds

An Apparently New Respiratory Disease of Baby Chicks

Infectious Bronchitis Virus Variants: Molecular Analysis and Pathogenicity Investigation

Genotyping of Infectious Bronchitis Viruses Identified in Canada between

Global Distributions and Strain Diversity of Avian Infectious Bronchitis Virus: A Review

Prevalence of Avian Influenza, Newcastle Disease, and Infectious Bronchitis Viruses in Broiler Flocks Infected with Multifactorial Respiratory Diseases in Iran

Epidemiology and Characterization of Avian Infectious Bronchitis Virus Strains Circulating in Southern China during the Period from

Current Situation, Genetic Relationship and Control Measures of Infectious Bronchitis Virus Variants Circulating in African Regions

The Molecular Biology of Coronaviruses

The Replicase Gene of Avian Coronavirus Infectious Bronchitis Virus Is a Determinant of Pathogenicity

Coronavirus Avian Infectious Bronchitis Virus

Gene 5 of the Avian Coronavirus Infectious Bronchitis Virus is not Essential for Replication

A Contemporary View of Coronavirus Transcription

Binding of Avian Coronavirus Spike Proteins to Host Factors Reflects Virus Tropism and Pathogenicity

Ultrastructural Pathology of Trachea in Chicken Experimentally Infected with Infectious Bronchitis Virus-MH-5365/95

Early Pathogenesis in Chicks of Infection with an Enterotropic Strain of Infectious Bronchitis Virus

Productive Replication of Nephropathogenic Infectious Bronchitis Virus in Peripheral Blood Monocytic Cells, a Strategy for Viral Dissemination and Kidney Infection in Chickens

The Pathogenesis of Nephritis in Chickens Induced by Infectious Bronchitis Virus

Pathogenesis of Oviduct Lesions in Immature Chickens Following Exposure to Infectious Bronchitis Virus at One Day Old

A 'New' Strain of Infectious Bronchitis Virus Infecting Domestic Fowl in Great Britain

Immunopathogenesis of Infection in SPF Chicks and Commercial Broiler Chickens of a Variant Infectious Bronchitis Virus of Economic Importance

Infectious Bronchitis

Nephropathogenicity of Infectious Bronchitis Virus

Pathogenicity of Four Serotypes of Avian Infectious Bronchitis Virus for the Oviduct of Young Chickens of Various Ages

Comparison of the Susceptibility to Avian Infectious Bronchitis Virus Infection of Two Inbred Lines of White Leghorn Chickens

Comparison of the Nephropathogenicity of Four Strains of Infectious Bronchitis Virus

Comparative Histopathology and Immunohistochemistry of QX-like, Massachusetts and 793/B Serotypes of Infectious Bronchitis Virus Infection in Chickens

Comparative Histopathology of Two Serotypes of Infectious Bronchitis Virus (T and n1/88) in Laying Hens and Cockerels

Pathogenesis of H13 Nephropathogenic Infectious Bronchitis Virus

Differences in the Tissue Tropism to Chicken Oviduct Epithelial Cells between Avian Coronavirus IBV Strains QX and B1648 are not Related to the Sialic Acid Binding Properties of Their Spike Proteins

Newly Isolated Serotypes of Infectious Bronchitis Virus: Their Role in Disease

Effects of Infectious Bronchitis in Baby Chicks

Persistence of Virus in the Tissues and Development of the Oviduct in the Fowl Following Infection at Day Old with Infectious Bronchitis Virus

Co-circulation of Four Types of Infectious Bronchitis Virus (793/B, 624/I, B1648 and Massachusetts)

Europe: History, Current Situation and Control Measures for Infectious Bronchitis

The First Evidence of a New Genotype of Nephropathogenic Infectious Bronchitis Virus Circulating in Vaccinated and Unvaccinated Broiler Flocks in Algeria

Using Phylogenetic Analysis to Examine the Changing Strains of Infectious Bronchitis Virus Infections in Ontario Over Time

Comparative Features of Infections of Two Massachusetts (Mass) Infectious Bronchitis Virus (IBV) Variants Isolated from Western Canadian Layer Flocks

DMV/1639) Infectious Bronchitis Virus (IBV) Variants Isolated in Eastern Canada Show Evidence of Recombination

Emergence of a Group 3 Coronavirus through Recombination

Serotype Shift of a 793/B Genotype Infectious Bronchitis Coronavirus by Natural Recombination

Emergence of Avian Infectious Bronchitis in a Non-vaccinating Country

Co-circulation of Three Clusters of 793/B-like Avian Infectious Bronchitis Virus Genotypes in Iranian Chicken Flocks

Comparison of the Pathogenicity of QX-like, M41 and 793/B Infectious Bronchitis Strains from Different Pathological Conditions

Pathogenesis and Tissue Distribution of Avian Infectious Bronchitis Virus Isolate IRFIBV32 (793/B serotype) in Experimentally Infected Broiler Chickens

Experimental Histopathologic Study of the Lesions Induced by Serotype 793/B (4/91) Infectious Bronchitis Virus

Pathogenicity of a TW-Like Strain of Infectious Bronchitis Virus and Evaluation of the Protection Induced against It by a QX-Like Strain

Isolation of Avian Infectious Bronchitis Virus from Experimentally Infected Chickens

A Long-term Study of the Pathogenesis of Infection of Fowls with Three Strains of Avian Infectious Bronchitis Virus

Egg Quality and Production Following Infectious Bronchitis Virus Exposure at One Day Old

The Influence of Infectious Bronchitis Virus on Egg Production, Fertility, Hatchability and Mortality Rate in Chickens (Author's Transl.). Nord. Vet

Effects of Infectious Bronchitis on the Reproductive Tracts, Egg Production, and Egg Quality of Laying Chickens

Monne, I. S1 Gene-based Phylogeny of Infectious Bronchitis Virus: An Attempt to Harmonize Virus Classification

Scalable Generation of High-quality Protein Multiple Sequence Alignments Using Clustal Omega

RAxML Version 8: A Tool for Phylogenetic Analysis and Post-analysis of Large Phylogenies

Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability

RDP4: Detection and Analysis of Recombination Patterns in Virus Genomes

Shell-Less Egg Syndrome (SES) Widespread in Western Canadian Layer Operations Is Linked to a Massachusetts (Mass) Type Infectious Bronchitis Virus (IBV) Isolate. Viruses

Comparative Nephropathogenicity of Different Strains of Infectious Bronchitis Virus in Chickens

Induction of Innate Immune Response Following Infectious Bronchitis Corona Virus Infection in the Respiratory Tract of Chickens

We would like to thank the University of Guelph for kindly providing us the IBV-infected samples. We acknowledge the help of the staff at the Prion-Virology Animal Facility and Animal Resource Center at Foothill campus, University of Calgary for the animal management.

The authors declare no conflict of interest.