key: cord-0748822-gkm4g6co
authors: Rohaim, Mohammed A.; El Naggar, Rania F.; Hamoud, Mohamed M.; Bazid, Abdel-Hamid I.; Gamal, Abdulrhman M.; Laban, Samah E.; Abdel-Sabour, Mohammed A.; Nasr, Shimaa A. E.; Zaki, Manal M.; Shabbir, Muhammad Z.; Zahran, Osama K.; Munir, Muhammad
title: Emergence and genetic analysis of variant pathogenic 4/91 (serotype 793/B) infectious bronchitis virus in Egypt during 2019
date: 2019-08-01
journal: Virus Genes
DOI: 10.1007/s11262-019-01693-9
sha: 31047a34b6e7cfe9ebd619ed6f518b5e983fd650
doc_id: 748822
cord_uid: gkm4g6co

Infectious bronchitis virus (IBV) affects both vaccinated and unvaccinated flocks worldwide, with a significant impact on the poultry industry. The aim of the present study is to characterize an emerging variant pathogenic IBV originating from field outbreaks in vaccinated Egyptian layer flock. Samples were collected from disease-suspected flock with a history of administration of live and inactivated IBV vaccines (Ma5 type). Virus propagation in embryonated chicken eggs (ECEs), after three successive passages, revealed typical IBV lesions such as curling and dwarfism. The reported isolate was identified by a real-time reverse transcriptase PCR assay targeting nucleocapsid (N) gene and, further characterized by full-length spike (S1) gene sequencing. Phylogenetic analysis revealed clustering of the isolated virus within 4/91 genotype of GI-13 lineage. Deduced amino acid sequences identity revealed 75–76% and 88–90% similarity with the currently used classic (H120, Ma5, and M41) and variant vaccine strains (4/91 and CR88) in Egypt, respectively. Recombination analysis gave an evidence for distinct patterns of origin for the studied isolate providing another example of intra-genotypic recombination among IBVs and the first example of recombination within the GI-13 lineage in the Egyptian field. The studied isolate (IBV/CK/EG/Fadllah-10/2019) emerged as a result of recombination between the variant group (Egy/var I genotype, GI-23 lineage) as a major parent and the CR88 variant vaccine strain (4/91 genotype, GI-13 lineage) as minor parent. Our data suggest that both mutation and recombination may be contributing to the emergence of IBV variants which ascertain the importance of disease monitoring in vaccinated flocks as well as re-appropriation for the current vaccine strategies. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s11262-019-01693-9) contains supplementary material, which is available to authorized users.

. IBV belongs to genus Gammacoronavirus within the Coronaviridae family [3] , has positive-sense single-stranded RNA genome (27.6 kb) which encodes four major structural proteins; spike glycoprotein (S), membrane glycoprotein (M), envelope protein (E), and phosphorylated nucleocapsid protein (N) [4] . The S1 subunit of spike (S) glycoprotein gene contains epitopes for virus neutralization, cell attachment, and serotype specificity [5] [6] [7] .

Genetic classification and evolutionary analysis of IBV variants are mainly based on S1 gene sequencing [8, 9] . The latest classification for IBV recognized six main genotypes (GI-GVI), 32 viral lineages , and a number of interlineage recombinants worldwide [10] . IBV was first reported in Egypt in 1954 [11] . By then, several reports have confirmed the disease circulation with continuous genetic diversity and evolution [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] . Despite mass vaccination strategies in Egypt utilizing Mass-type (H120, M41, and Ma5) and 793B-type(CR88 and 4/91) of vaccine, variant IBVs still devastate the poultry industry periodically. During the years 2010 to 2018, majority of the circulating IBV strains in the Egyptian poultry industry were clustered within GI-23 lineage that were further sub-categorized into Egy/variant I and Egy/variant II of IS/1494 and IS/885 origin, respectively [10, 23, 24] . Indeed, the IBVs pose a significant economic impact on poultry industry and, therefore, disease consequences are devastating the national economy of Egypt. The aim of the present study was to (1) characterize newly emerged variant IBV or escape mutant from a breeder layer flock and (2) determine the genetic divergence between this circulating field variant and the currently used classic and variant vaccine strains in Egypt.

Despite a history of vaccination using live attenuated and inactivated Ma5 vaccine (Table 1) , clinical infectious bronchitis disease was suspected in a 85-week-old layer flock (White Bovans breeds) exhibiting respiratory signs along with high mortality and morbidity. The affected flock was raised in cages at Ismailia governorate, Egypt. Respiratory tract tissues (trachea) collected on 15th January 2019, showed congestion with caseous exudates, cloudiness, and turbidity of air sacs. Furthermore, degeneration of the ovary and swollen oviducts associated with egg peritonitis was observed. For virus isolation and propagation, a volume of 0.2 ml of tracheal homogenate was inoculated into the allantoic cavity of specific pathogen-free embryonated chicken eggs (SPF ECEs) (9-to 11-day old) [27] , kindly obtained from Kom Oshim farm, Fayoum, Egypt. Inoculated eggs were monitored daily by candling for embryonic mortality. Two days after inoculation, the allantoic fluid was harvested and further passaged into another two successive passages.

Viral RNA was extracted from the allantoic fluid of 3rd passage using a QIAamp viral RNA mini kit (Qiagen, Valencia, CA) according to the manufacturer's instructions. Detection of IBV was conducted using a real-time reverse transcriptase PCR based on the highly conserved N gene, as described previously [28] . Amplification of full-length spike (S1) gene was conducted as described previously [25] using the Mx3005Pe system (Stratagene, Santa Clara, CA). PCR products were analyzed by agarose gel electrophoresis (1%) and then purified using a QIAquick Gel Extraction Kit (Qiagen) following the manufacturer's instructions. Sequencing was carried out using BigDye Terminator v3.1 sequencing kit (Applied Biosystems) and an automated sequencer (ABI, 3130, Applied Biosystems, Foster City, CA). The quality of obtained sequence (S1 gene) was checked, assembled, edited using BioEdit software version 7.0.4.1 [29] , and submitted to GenBank using BankIt tool of the GenBank (http://www. ncbi.nlm.nih.gov/WebSu b/?tool=genba nk), and obtained the accession number KP729419.

Deduced amino acid sequences alignment for the S1 subunit of the spike protein was performed by comparing with commonly utilized commercial vaccine strains in Egypt (Ma5, Ma41, H120, CR88, and 4/91). Phylogenetic analysis was conducted to explain the phylogenetic relationships with high-level of clustering pattern dependent on the full-length S1 gene between the reported isolate and recently described strains from the Middle East including Egypt and other parts of the world. For maximum-likelihood analysis of the phylogenetic relationship, Kimura 2-parameter model of nucleotide substitution was used with gamma rate distribution, which was the best-fitting substitution model overall for the three clades, as determined by the substitution model-fitting function in MEGA 7 [30] . Trees were finally visualized and annotated using FigTree v1.4.2 software (http://tree. bio.ed.ac.uk/softw are/figtr ee/). For recognition of recombination occasions, sequence alignments were conducted, in barrel with different available viruses resembling genotypes GI-1, GI-13, GI-19, and GI-23. Recombination Detection Program 4 (RDP4) software suite [31] were used to identify recombination events in the full-length S1 gene of IBV/CK/EG/Fadllah-10/2019 isolate through detection of breaking points using specific algorithms implemented in RDP4; RDP, Genecov, Bootscan, Maxchi, Chimaera, Siscan, and 3Seq with the highest acceptable P value adjusted to 0.05.

Here, we present the isolation and integrative genetic analysis that maps the emergence of 4/91 (793B) genotype in Egypt during 2019. Virus propagation in embryonated chicken eggs (ECEs), after three successive passages, revealed typical IBV lesions such as curling and dwarfism. After identification of IBV using a N gene-targeted real-time RT-PCR assay (Ct 12.22), the isolated IBV strain (IBV/CK/EG/Fadllah-10/2019) was genotyped based on sequencing and subsequent sequence analysis of the fulllength S1 gene. The obtained nucleotide (1617bp) and corresponding amino acid sequences were aligned and compared with reference and vaccine strains (H120, Ma5, M41, 4/91, and CR88) including previously described IBV variants in Egypt and neighboring countries.

Compared to currently used vaccines in Egypt, the characterized virus shared different levels of nucleotide and amino acid sequence identities which were between 75 and 90% (Tables S1 and S2). The S1 protein contained hypervariable regions (HVRs) associated with serotype specificity and virus-neutralizing epitopes and located within the amino acid residues 38-67, 91-141, and 274-387 [32, 33] . A significant clustering of substitutions were detected in IBV/CK/ EG/Fadllah-10/2019 in the three HVRs as compared to 4/91 and other vaccine strains. The prominent amino acid substitutions within the three HVRs in comparison with 4/91 vaccine were as follows at HVR1; D60F, V62A, S63G, D64Q, T69S, F70I, Y71H, E72W, Y74K, I76F, A79S. While at HVR2, the mutations were F115Y, S117N, Q118G, N131D, I135R, R139M, S141Y, F143I, and at HVR3 was only single-mutation (A281P) ( Table 2) . It is well known that even a small change in the amino acid sequence of the spike protein can result in generation of novel genotypes and/or serotypes that differ antigenically from the existing classic and variant vaccine strains [24, 25, 34] . In addition, it has been shown that changes as little as 5% in the S1 gene are able to alter the protective ability of a vaccine [33] . Moreover, residues N38, H43, P63, and T69 have been described to be critical for binding of the IBV spike protein to the chicken respiratory tract [35] . In the analyzed sequence of the isolate under study, residue substitutions associated with virus tropism were identified at positions N38P and G63Q in the receptorbinding site of the variant strains.

Previously Massachusetts-type strains were used for IBV vaccination programs in Egypt. Recently, the use of other strain types (e.g., 793/B serotype) has been officially authorized. However, the antigenic relatedness between GI-1, GI-16, and GI-23 lineages has not yet been analyzed, even though it may be useful to understand IBV epidemiology and to design specific control programs. Phylogenomic analysis based on S1 gene of the IBVs collected from Egypt and rest of the world indicated that Egyptian virus reported here was clustered within 4/91 genotype of GI-13 lineage (Fig. 1) . Unique mechanism of RNA synthesis involving polymerase jumping and discontinuous transcription might be responsible for recombination in coronaviruses [36] . Studies reported that occurrence of recombination among different IBV field strains results in the emergence Table 2 Sequence alignment of HVRs amino acid sequences of the reported isolate in this study compared to other reference strains/genotypes of new variants [37] [38] [39] . Many recombination events have been reported in different IBV strains, not only between field (wild type) and vaccine viruses but also among field viruses either within the same genotype (intra-genotypic) or between different genotypes (inter-genotypic) [40, 41] , giving rise to new IBV genotypes [42] . Different available viruses resembling genotype GI-1, GI-12, GI-13, and the two subgroups of GI-23 were used in the analysis. The employed recombination detection methods revealed that IBV/CK/EG/Fadllah-10/2019 strain has undergone intergenotypic recombination and different breakpoints within the S1 gene which might be a precursor for further GI-13 lineage evolution in Egypt in the near future. So, future studies are required to determine the pathobiological and clinical features of this virus in a chicken model. Also, the results showed that the reported isolate IBV/CK/EG/ Fadllah-10/2019 is a recombinant virus which probably emerged from at least two different genotypes, including the Egy/Variant I genotype (GI-23 lineage) as a major parent and the CR88 vaccine strain (4/91 genotype of GI-13 lineage) as minor parents (Fig. 2) .

To recapitulate, the current study revealed the isolation and molecular identification of pathogenic 4/91 IBV collected from Egyptian vaccinated layer flock during 2019. Amino acids alignment revealed distinct mutations within HVRs compared to the commercially used vaccines. Likewise, phylogeny revealed that clustering of IBV/CK/EG/ Fadllah-10/2019 strain within GI-13 lineage. Further analysis based on recombination detection analysis revealed the isolated strain undergoes distinct spots of recombination. Thus, the identified variant pathogenic 4/91 IBV may express as a new genotype that requires re-evaluation of vaccination strategies employed in Egypt. In conclusion, owing to high mutation rate and subsequent residue substations, there is an ongoing evolution among the circulating IBV variants in Egypt that results in frequent vaccination failure. A continuous disease monitoring and surveillance is required not only to elucidate sequence characteristics of prevailing strains but also to revise or appropriate vaccine strains and strategies accordingly. 

GI-13

FN182243.1_IBV_strain_NGA/B401/2006 FN182272

JX173488.1_IBV_isolate_Eg

EU780077.2_IBV_isolate_IS_1494_06

Egy/Variant 2 GI-23

GU938413.1_IBV_strain_CK/CH/Chongqing/0908 AY606323

AF151958.1_IBV_strain_K43 AF151955.1_IBV_strain_C AF151960

KC577397.1_IBV_isolate_66GD-98VI HQ018890.1_IBV_strain_CK/CH/JX

AJ311317.1_IBV_strain_Beaudette_CK DQ487085.1_IBV_isolate_Egypt_F_03 AY561711.1_IBV_strain_M41 GU393335.1_IBV_serotype_H120 AY561713.1_IBV_strain_Ma5 FJ904719.1_IBV_strain_Conn46_1991

AY296745.1_IBV_strain_JP8443 AY296744.1_IBV_strain_JP8127 KC577388.1_IBV_isolate_48SD-96VI KC577391

Coronavirus avian infectious bronchitis virus

Avian infectious bronchitis. Manual of diagnostic tests and vaccines for terrestrial animals

Family Coronaviridae

Coronavirus proteins: biogenesis of avian infectious bronchitis virus virion proteins

Monoclonal antibodies to the S1 spike and membrane proteins of avian infectious bronchitis coronavirus strain Massachusetts M41

Coronavirus IBV: removal of spike glycopolypeptide S1 by urea abolishes infectivity and haemagglutination but not attachment to cells

Location of the amino acid differences in the S1 spike glycoprotein subunit of closely related serotypes of infectious bronchitis virus

Typing of field isolates of infectious bronchitis virus based on the sequence of the hypervariable region in the S1 gene

Evolutionary implications of genetic variations in the S1 gene of infectious bronchitis virus

S1 gene-based phylogeny of infectious bronchitis virus: an attempt to harmonize virus classification

Incidence and treatment of some infectious viral respiratory diseases of poultry in Egypt

Present status of infectious bronchitis in Egypt

Isolation and identification of Egypt/Beni-Seuf/01 a novel genotype of infectious bronchitis virus

S1 gene sequence analysis of a nephropathogenic strain of avian infectious bronchitis virus in Egypt

Isolation and characterization of nephropathogenic strain of infectious bronchitis virus in Egypt

Emergence of a novel genotype of avian infectious bronchitis virus in Egypt

Molecular characterization of infectious bronchitis viruses isolated from broiler and layer chicken farms in Egypt during 2012

Sequence analysis of infectious bronchitis virus IS/1494 like strain isolated from broiler chicken co-infected with Newcastle disease virus in Egypt during 2012

Molecular diversity between field isolates and vaccinal strains of avian infectious bronchitis virus in Egypt

Isolation and characterization of new variant strains of infectious bronchitis virus in Northern Egypt

Characterization of infectious bronchitis viruses circulating in Egyptian chickens during 2012 and 2013

Molecular studies on infectious bronchitis virus isolated from broiler chickens in Damietta Governorate

Genotyping and pathotyping of diversified strains of infectious bronchitis viruses circulating in Egypt

The sequence of the full spike S1 glycoprotein of infectious bronchitis virus circulating in Egypt reveals evidence of intra-genotypic recombination

Isolation and molecular characterization of novel infectious bronchitis virus variants from vaccinated broiler flocks in Egypt

Genetic diversity and phylodynamics of Avian coronaviruses in Egyptian wild birds

Manual diagnostic tests vaccines terrestrial animals (mammals, birds bees)

Development of a real-time TaqMant RT-PCR assay for the detection of infectious bronchitis virus in chickens, and comparison of RT-PCR and virus isolation

BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT

MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets

RDP4: detection and analysis of recombination patterns in virus genomes

Identification of amino acids involved in a serotype and neutralization specific epitope within the S1 subunit of avian infectious bronchitis virus

Relationship between serotypes and genotypes based on the hypervariable region of the S1 gene of infectious bronchitis virus

Complete genome sequences of two avian infectious bronchitis viruses isolated in Egypt: evidence for genetic drift and genetic recombination in the circulating viruses

Detection of infectious bronchitis virus

Recombination in large rna viruses: coronaviruses

Complete genomic sequence analysis of infectious bronchitis virus ark dpi strain and its evolution by recombination

Characterization and analysis of the full-length genome of a strain of the european qx-like genotype of infectious bronchitis virus

Genomic and single nucleotide polymorphism analysis of infectious bronchitis coronavirus

Complete genome sequence and recombination analysis of infectious bronchitis virus attenuated vaccine strain H120

Altered pathogenicity of a tl/CH/ LDT3/03 genotype infectious bronchitis coronavirus due to natural recombination in the 5'-17 kb region of the genome

Emergence of a group 3 coronavirus through recombination

We would like to thank all co-workers and col-