key: cord-0850557-c1p2w375
authors: Fellahi, Siham; EL Harrak, Mehdi; Ducatez, Mariette; Loutfi, Chafiqa; Koraichi, Saad Ibn Souda; Kuhn, Jens H.; Khayi, Slimane; EL Houadfi, Mohammed; Ennaji, My Mustapha
title: Phylogenetic analysis of avian infectious bronchitis virus S1 glycoprotein regions reveals emergence of a new genotype in Moroccan broiler chicken flocks
date: 2015-08-04
journal: Virol J
DOI: 10.1186/s12985-015-0347-8
sha: d38f4068c4bf5d99a9f4fc20b1914a3332a42ead
doc_id: 850557
cord_uid: c1p2w375

BACKGROUND: Infectious bronchitis virus (IBV), a major pathogen of commercial poultry flocks, circulates in the form of several serotypes/genotypes. Only a few amino-acid changes in the S1 subunit of wild-type IBVS proteins may result in mutants unaffected by current vaccines. METHODS: Partial S1 gene sequences of 3 IBV isolates of the Moroccan Italy 02 genotype from vaccinated and unvaccinated broiler chicken flocks, located in southern and central regions of Morocco, were amplified by RT-PCR, sequenced, and aligned for phylogenetic and amino-acid similarity analyses. RESULTS: The three isolates were found genetically highly distant from known avian IBV based on partial sequences of their S1 genes: gammaCoV/chicken/Morocco/I01/2011(IBV/Morocco/01), gammaCoV/chicken/Morocco/I30/2010 (IBV/Morocco/30), and gammaCoV/chicken/Morocco/I38/2013 (IBV/Morocco/38), nucleotide sequence identities reached 89.5 % to 90.9 % among the three isolates. The deduced protein sequence identities ranged from 29.7 % (between IBV/Morocco/38 and Egypt SCU-14/2013-1) to 78.2 % (between IBV/Morocco/01 and Spain/05/866). Amino acid sequence comparison and phylogenetic analysis indicated the emergence of a new Moroccan genotype, clustering with regionally related isolates from Spain (Spain/05/866) and belonging to a new sub-genotype. CONCLUSION: Our sequencing results demonstrate a co-circulation of wild-type infectious bronchitis viruses in broiler chickens. These results justify permanent monitoring of circulating strains in order to rationally modify vaccination strategies to make them appropriate to the evolving field situation.

Infectious bronchitis (IB) is one of the most contagious diseases that affect poultry worldwide and is responsible for severe economic losses. Clinically, the disease causes respiratory distress, drop in egg production and quality in layers. Some strains of IBVare associated with nephritis [1] [2] [3] .

Infectious bronchitis virus (IBV) is a member of the species Avian coronavirus, genus Gammacoronavirus (Nidovirales: Coronaviridae: Coronavirinae) [4] . IBV has a linear, single-stranded RNA genome of positive polarity of approximately 27 kb in length and produces enveloped virions. IBV particles consist of three major structural proteins: spike (S) glycoprotein, membrane (M) protein, and nucleocapsid (N) protein. The spike of IBV is formed by post-translational cleavage of S into two separate polypeptide components, S1 and S2 [5, 6] . S1 mediates virion attachment to IBV host cells and is a major target of neutralizing antibodies in chickens. Genotype evolution of IBV is associated primarily with changes in the S1amino-acidsequence [7, 8] . Hence, the evolutionary characterization of IBV is mainly based on the analysis of the variable S1 gene or the expressed S1 protein [9, 10] .

Different IBV variants are distributed globally. Some of these variant are endemic only in particular regions, while others circulate worldwide [11] . More than 20 different IBV serotypes are differentiated worldwide that evolved from genomic insertions, deletions, substitutions, and/or RNA recombinations of the S1 gene [12] [13] [14] . This large diversity of serotypes is a major reason why commercial vaccines often fail or are only partially efficacious, and, therefore, new IB outbreaks continue to occur [1] .

The first isolation and characterization of IBV from poultry flocks in Morocco was reported by El-Houadfi et al. in 1986 [15] . Six isolates were obtained, of which isolates designated D, E, F, H, and M were found to be serologically related to the Massachusetts (Mass) serotype, whereas the sixth isolate, G, differed from the Mass serotype and other serotypes known at the time. Importantly, El-Houadfi et al. demonstrated that Mass-based IBV vaccines provided poor protection against infection with isolate G [15] . In 2004, Alarabi conducted a study to determine the relationship between IBV and nephropathogenic disease outbreaks observed in broiler flocks in Morocco between 1996 and 2000. Three different groups of IBV isolates were identified using RT-PCR coupled with restriction fragment length polymorphism (RFLP). Group I belonged to the Mass serotype, whereas groups II and III were distinct. Isolate 12/97 of group III, found to be closely related to isolate G, caused more severe kidney lesions and higher lethality in experimentally infected animals compared to isolate 7/97 of group II [16] . In 2005, El Bouqdaoui et al., while studying nephropathogenic IBV using RT-PCR and RFLP techniques, identified five genotypes, three of which differed from vaccine strains [17] .

The Moroccan poultry industry has developed significantly under the framework "Plan Green Morocco" and meets the growing domestic consumer demand for poultry products, including meat and eggs The industry generates a turnover of around 30 billion dirhams/year, totaling mean investments of around 8.7 billion dirhams/year. The official data regarding the general poultry rearing systems in Morocco are: broiler production, 450 million/year (among 6,800 farms); layers: 20 million/year (250 farms); turkeys: 12 million/year (25 % imported d-old chicks); broiler breeders: 3.2 million/year (70 % imported from Spain and 30 % from other European countries); 48 hatcheries. The Interprofessional Federation of Poultry Sector (FISA) aims to develop exports of poultry products in particular to Northern and Western Africa. FISA is affiliated with Association of Moroccan Exporters (ASMEX) under the objective of promoting poultry product export in Africa, including chicks, hatching eggs, and compound feed for poultry. In 2012, 6.7 million hatching eggs and 1.9-million-d-old chicks were exported to Mauritania, Mali, Cameroon, Côte d'Ivoire, and Senegal, and 25,000 tons of compound feed were exported to Mauritania. Breeder chicks (182,000) were imported in January 2014 (an increase of 28 % compared to 2013). However, imports of layer and broiler breeder and turkey poults declined in recent years, mainly due the development of local production.

Three IBV vaccines are used in Morocco to protect the industry: 793B, Arkansas, and Mass. The Mass type has been used as of 1970 (H120, Ma5, and modified Mass strains), whereas 793B (4/91 and CR88 strains) and Arkansas vaccines were introduced in 2000 and 2013, respectively. New IBV serotypes and genotypes can emerge as a result of only very few changes in the amino acid sequence of the S protein. Vaccination programs largely rely on the use of the IBV Massachusetts strain, which is also the most commonly used IBV vaccine strain in Morocco. However, despite the use of this vaccine, the common presences of IB in vaccinated chickens continue to have a major negative economic impact [18] . In this study, IB was diagnosed between January 2010 and December 2013 in southern and central regions of Morocco. Suspected IBV infections were found in 47 commercial broiler chicken farms. Among these, three flocks showed severe clinical signs and mortality. The aim of the present study is to characterize a new emerging genotype of IBV from these outbreaks using molecular techniques, and to determine the relationship to reference IBV strains by nucleotide sequence analysis.

Infectious bronchitis was diagnosed between January 2010 and December 2013 in southern and central regions of Morocco. Suspicious IBV infections (prominent respiratory disease) were found in broiler chickens from 47 commercial farms. The flocks had been vaccinated against Newcastle disease and infectious bursal disease (Table 1) . Among these farms, three flocks had severe respiratory signs and experienced increased mortality. The first flock, located in southern Morocco (Marrakech), was 26 d old and was the only flock that had been vaccinated with the H120 live IB vaccine by spray on d 7 and 21. The second flock, located in central Morocco (Khamisset), was 24 d old, and the third flock, located also in central Morocco (Rabat) was 20 d old. Clinical presentation started with respiratory depression and distress, including nasal discharge, sneezing, coughing, and rales. Other signs included conjunctivitis and watery eyes. Within 10 d after the appearance of the disease, lethality increased from 10 % to 30 % of the total flocks. Post-mortem examination of dead birds revealed increased tracheal mucus, slight congestion, and presence of catarrhal exudates in the nasal turbinate and trachea.

We obtained a total of 47 Moroccan field isolates of IBV from the 47 affected farms and were able to distinguish them based on sequences of the hypervariable region (HVR) 3 of their S1 regions (nucleotides 700-1095). Of those isolates, three isolates, classified as the Moroccan Italy 02 genotype, underwent further sequencing to include all three hypervariable S1regions (nucleotides 1-1,095). These three IBV strains were isolated from trachea, lungs, and kidneys of broiler chickens. Between the first and fourth passage, virus propagation from the three suspect samples in specific pathogen-free eggs resulted in characteristic IBV lesions (e.g., curled and dwarfed embryos with head pressed over the feet and covered by a thickened amnion)in 90 %-100 % of the embryos. Results of both RT-PCR targeting the IBV N gene and RT-PCR using primers amplifying a partial sequence of the IBV S1 gene were positive, confirming the identity of the three isolates. The amplified partial S1 fragments (approximately 700 bp in length) of the three isolates were sequenced. The sequences were submitted to GenBank and assigned the accession numbers KJ701019, KJ701020 and KM594187 for isolates gammaCoV/chicken/ IBV/Morocco/I38/2013 (IBV/Morocco/38), gammaCoV/ chicken/IBV/Morocco/I30/2013 (IBV/Morocco/30), and gammaCoV/chicken/IBV/Morocco/I01/2013 (IBV/Morocco/ 01), respectively (Table 1) .

Partial S1 gene sequences of the three IBV isolates were sequenced. The nucleotide sequence and deduced amino acid sequences ( [19] when compared to non-Moroccan IBV strains ( Table 2 ). The alignment of the Moroccan IBV isolates amino-acid sequences with that of reference strain H120 identified the most variable regions in residues 62-77 and 117-131 (numbering in reference to the H120 strain, Fig. 1 ). Insertions in the partial S1 subunit of the three Moroccan strains were found at positions119, 120, 143, and 144 (Fig. 1) .

The partial S1 protein sequences of the Moroccan isolates were 69.4 % to 76. 8 (Table 2) .

Phylogenetic trees were constructed from the nucleotide and deduced amino-acid sequences of the partial S1 glycoprotein genes of Moroccan IBV isolates and non-Moroccan IBV referenced strains. The three Moroccan IBV isolates and Spanish isolates of the Italy 02 genotype formed a common branch (Fig. 2) , but they clearly clustered into 2 distinct sub-genotypes within the group. The sequence analysis of the partial S1 gene demonstrated that Moroccan isolates represent a unique subgenotype compared to other reference strains of various countries and are distantly related to African, European, and American strains.

IBV was detected in poultry in Morocco in the 1980s, and in more recent years, outbreaks of IB were reported repeatedly [15] [16] [17] . The present study suggests that a new genotype, represented by the three new IBV isolates described here, has been circulating in Morocco since at least 2010 (IBV/Morocco/30, accession number: KJ701020). Phylogenetic analyses comparing the partial S1 gene sequences of the three Moroccan field isolates with those of reference and variant strain sequences revealed that the Moroccan isolates were 69.4 % to 76.8 % homologous with Spanish isolates of the Italian 02 genotype [19] . This similarity can be explained by the geographical proximity and close transboundary of Morocco and Spain. Moreover, these two countries have an ancient history of commercial exchanges. In addition, Spain is the most frequent country that exports breeder chicks to Morocco (exports areunder the control of Office National de Sécurité Sanitaire des ProduitsAlimentaires, http://www.onssa.gov.ma/).

The deduced amino-acid sequences of Moroccan IBV isolate proteins demonstrated that the Moroccan isolates represent a unique genotype compared too ther reference strains of various countries. These Moroccan isolates are only distantly related to African, European, and American strains. The origin of this novel genotype remains to be identified. In addition to population diversity [20] [21] [22] , widespread use of various vaccines made from heterotypicIBVs in the field, recombination as a consequence of mixed infections, or a decrease of dominant serotypes as a result of vaccination may play important roles in the emergence of such novel genetic variants in Morocco.

Alignment of the partial S1 glycoprotein sequence of Moroccan IBV isolates with thatof published non-Moroccan IBV strain, H120, identified most sequence variations between residues 62-77 and 117-131 (Fig. 1) . Both of these sequence regions correspond to known HVR, mainly HVR1residues 50-69, and HVR2, residues 117-131, of the IBVS1 gene. HVR1 and HVR2 contain sequences that have been associated with specific IBV serotypes [7, 23] andserotype-specific neutralizing epitopes [8] . The nucleotide S1sequences of IBV serotypes commonly differ by 20 to 25 % [24] , but amino acid S1 sequences of some serotypes differ by as little as 2 % [25] . The HVRs are associated with two separate viral neutralizing and conformationally-dependent epitopes [7, 26] .

The insertions of four amino acids in the S1 subunit of the three Moroccan strains possibly havesome significance relatedto S1 functions (e.g., neutralization phenotypes, virus tropism) [7, 25, [27] [28] [29] . Further studies are needed to determine the role of theseinsertions in the S glycoprotein [30] , which might be related to reported changes in the clinical manifestations of IB.

In the phylogenetic tree shown in Fig. 2 , the Moroccan IBV isolates formed a distinct cluster withthe Spanish isolates of the Italy 02 genotype (supported by a high bootstrap value: 96). The Spanish IV genotype was previously established as a new genotype based on the mosaic nature of its S1 gene. Viral neutralization studies revealed Italy 02 and Spanish IV genotypes are subtypes within a common serotype [19] .

This molecular study has revealed that a new genotype can emerge as a result of only a few changes in the amino acid composition in the S1subunit of the virus Sprotein [25] . This is the first report of identification and genotyping of emergent IBV isolates in Morocco, which indicate the circulation of previously unknown IBV variants in both vaccinated and unvaccinated broiler chickens. Our data provide evidence of a new emerging genotype in Morocco. Three IBV field isolates were different from published referenced strains in nucleotide and amino acid sequences and clustered into a separate branch based onphylogenetic analysis. All Moroccan IBV variants presented specific nucleotide and amino acid sequence insertions in the S1 gene in comparison to the Massachusetts vaccine strain that might be associated with the occurrence of clinical disease in vaccinated flocks. Further investigation related to cross protection and pathogenicity of this new genotype detected in this country are needed.

Tissue samples from trachea, lungs, and kidneys were obtained from 10 birds from each flock, placed on ice, and transferred to the laboratory in sterile transport medium containing 5 % antibiotics (20,000 U/ml of penicillin, 10,000 μg/ml of streptomycin, and 5,000 μg/ml of kanamycin). These samples were used for virus isolation.

Virus isolation was performed by inoculation of 9-11-dold, embryonated, specific pathogen-free hen's eggs with 200 μl of 10 % tissue homogenates. Eggs were incubated at 37°C and candled daily. Embryos that died within 24 h after inoculation were discarded. On d 2 after inoculation, three eggs were chilled at 4°C for 2 h, and the allantoic fluids were harvested, centrifuged at 8,000 rpm (xxx x g) for 15 min, and stored at −80°C. The remaining eggs were further incubated for 5 more d and observed for typical embryonic lesions of IBV infection, such as dwarfing, stunting, and curling. Stored allantoic fluids were further passaged three times in specific pathogen-free eggs; the allantoic fluids were Fig. 2 Phylogenetic relationships of the obtained Moroccan isolates and selected reference strains based on partial S1 nucleotide sequences determined using MEGA 5.0 with the Clustal W method. Numbers along the branches refer to bootstrap values harvested and tested for the presence of IBV using rRT-PCR.

Viral RNA was extracted from infected allantoic fluids using the MagMax Express semi-automatic magnetic particle extractor (Life Technologies, Grand Island, NY) according to the manufacturer's instructions. Allantoic fluid samples and reagents were dispensed into X-well microplates prior to magnetic particle processing. After processing, each sample of extracted RNA recovered from the wells of the plates was subjected toRT-PCR amplification.

The oligonucleotides and the probe used for initial screening of IBV were described previously [31] . These oligonucleotides and the probe target the relatively conserved IBV N gene at nucleotide positions 741-1077 of the IBV Massachusetts H120 reference strain (GenBank accession number: AM260960). RNase free water and Newcastle disease virus strain (HB1) were included as negative controls. One-step RT-PCR was performed with 12.5 μl of 2× RT-PCR buffer mix, 0.5 μlof MgSO 4 (50 mmol/l), 0.5 μl of Rox reference dye (25 mmol/l, Life Technologies), 0.5 μl of Moloneymurine leukemia virus reverse transcriptase (200 U/μl), 0.5 μlof Taq DNA polymerase (Life Technologies), 0.5 μl of primers (0.2 μmol/ l), 0.25 μl of probe (0.1 μmol/l), and 5 μl of RNA template to make a final volume of 25 μl.The reaction was carried out using a StepOne Plus real-time PCR system (Smart Cycler, Cepheid, Sunnyvale, CA) at 50°C for 15 min, 95°C for 5 min, and followed by 40 cycles of 9°C for 15 s; 60°C for 45 s; 72°C for 30 s, and a final extension step of 74°C for 5 min. Amplifications were recorded and analyzed, and the threshold cycle (C t ) was determined with StepOne software (Smart Cycler).

Amplification of partialS1 gene sequences by RT-PCR was performed using the forward S1S15' mod (5'-TGA-AAA-CTG-AAC-AAA-AGA-3') and reverse S1CK2 (5'-CNG-TRT-TRT-AYT-GRC-A-3') primers [32] . RT-PCR was performed using the Applied Biosystems kit protocol (Life Technologies). The RT reaction was performed in a 20 μl reaction mixture containing 2 μl of buffer (10×), 2.5 μl of MgCl 2 (25 mmol/l), 2.5 μl of dNTP (10 mmol/l), 0.75 μl of each primer (10 μmol/l), 10.2 μl of sterilized water, 0.5 μl of RNAase inhibitor (20 U/μl), 0.3 μl of RT (50 U/μl), and 0.5 μl of Gold Taq polymerase (5 U/μl). The PCR products were analyzed on a 1 % agarose gel.

Partial RT-PCR products of the S1 gene (700 bp) containing HVR1, HVR2, and HVR3 were purified using the Gene Clean Kit (ExoSAP-IT, Affymetrix, Santa Clara, CA) according to the manufacturer's instructions. Purified PCR products were used as templates for sequencing using the BigDye® Terminator v1.1 Cycle Sequencing Kit (Life Technologies). The second purification step was performed using the Big Dye XTerminator Purification Kit (Life Technologies). Purified PCR products were sequenced from both directions using the same primers (S15' mod and CK2).

Assembly and analysis of sequence data were conducted using the BioEdit Software version 5.0.9 [33] . The opensource BLAST program (National Center for Biotechnology Information, Bethesda MD, http://blast.ncbi.nlm.nih.gov/ Blast.cgi) was used for sequence comparison. Nucleotide sequence and deduced amino acid sequences were aligned using ClustalW and MEGA software Version 6.0 [34] . Phylogenetic analyses and tree construction for partial S1 glycoprotein gene sequences were generated using the neighbor-joining method with 1,000 bootstrap replicates with MEGA. Author details

QX (AF193423)

Tunisia TN296/07 (FJ716133)

Isolate Q1 (AF286302)

Abbreviations IBV: Infectious bronchitis virus

Severe acute respiratory syndrome vaccine development: experiences of vaccination against avian infectious bronchitis coronavirus

Coronaviruses in poultry and other birds

Review of infectious bronchitis virus around the world

Completion of the sequence of the genome of the coronavirus avian infectious bronchitis virus

Coronavirus IBV: structural characterization of the spike protein

Coronavirus avian infectious bronchitis virus

Amino acids within hypervariable region 1 of avian coronavirus IBV (Massachusetts serotype) spike glycoprotein are associated with neutralization epitopes

Location of antigenic sites defined by neutralizing monoclonal antibodies on the S1 avian infectious bronchitis virus glycopolypeptide

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

Infectious bronchitis virus variants: a review of the history, current situation and control measures

Molecular characterization of avian infectious bronchitis virus strains isolated in Colombia during

Variant serotypes of infectious bronchitis virus isolated from commercial layer and broiler chickens

Evidence of genetic diversity generated by recombination among avian coronavirus IBV

The isolation and characterisation of six avian infectious bronchitis viruses isolated in Morocco

A field study of kidney disease among the broiler flocks in Morocco and its relationship to infectious bronchitis virus

Genetic grouping of nephropathogenic avian infectious bronchitis virus isolated in Morocco

Viral respiratory diseases (ILT, aMPV infections, IB): are they ever under control?

Molecular epidemiology and evolution of avian infectious bronchitis virus in Spain over a fourteen-year period

Antigenic and molecular characterization of isolates of the Italy 02 infectious bronchitis virus genotype

The quasispecies (extremely heterogeneous) nature of viral RNA genome populations: biological relevance-a review

The molecular biology of coronaviruses

Phylogeny of antigenic variants of avian coronavirus IBV

Identification of avian infectious bronchitis virus by direct automated cycle sequencing of the S-1 gene

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

Sequence evidence for RNA recombination in field isolates of avian coronavirus infectious bronchitis virus

Relationship between sequence variation in the S1 spike protein of infectious bronchitis virus and the extent of crossprotection in vivo

The S1 glycoprotein but not the N or M proteins of avian infectious bronchitis virus induces protection in vaccinated chickens

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

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

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

Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. World Organisation for Animal Health

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

MEGA6: Molecular Evolutionary Genetics Analysis version 6.0

We thank Laura Bollinger (Integrated Research Facility at Fort Detrick) for providing Technical Writing Services in critically editing this manuscript on behalf of Battelle Memorial Institute. The content of this publication does not necessarily reflect the views or policies of the US Department of Health 

The authors declare that they have no competing interests.Authors' contributions SF participated in the design of the study, performed the molecular genetic studies andthe phylogenetic sequence alignment,and drafted the manuscript. MD participated in the phylogenetic analysis and helped to draft the manuscript. CL participated in virus isolation. SIK participated in virus sequencing. JHK critically editedand corrected the manuscript. SK helped with sequence alignment. KK helped withvirus sampling and isolation. MEH participated in the design of the study, and helped to draft and improve the manuscript. MME conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.