key: cord-0998475-b5fhgxud authors: Mousavi, Fatemeh Sadat; Ghalyanchilangeroudi, Arash; Hosseini, Hossein; Nayeri Fasaei, Bahar; Ghafouri, Seyed Ali; Abdollahi, Hamed; Fallah-Mehrabadi, Mohammad Hosein; Sadri, Naser title: Complete genome analysis of Iranian IS-1494 like avian infectious bronchitis virus date: 2018-06-09 journal: VirusDisease DOI: 10.1007/s13337-018-0462-4 sha: de49d38fbe22e6078e967a7b15ab27a977609250 doc_id: 998475 cord_uid: b5fhgxud The nephropathogenic infectious bronchitis virus (IBV) strain IS-1494 like (variant-2; GI-23) was first isolated in the Middle East (1998). Despite intensive vaccinations, IS-1494 like IBVs are still circulating in Iran (the dominant genotype) and spread to other countries. Here, the full-length genome of this Iranian IS-1494 like IBV was (Mahed) determined to understand its evolutionary relationships. The genome consists of 27,652 nucleotides, with mutations in most of the structural genes. Thirteen open reading frames (ORFs) were predicted in the Mahed isolate (5′ UTR-1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3′ UTR). ORFs 4b, 4c, and 6b, which has rarely been reported, were present in the Mahed genome. According to phylogenetic analysis of the full-length genome, 1a, S2, M, E, N protein, Mahed isolate clustered with the QX type strain. Based on the partial 1b, S1, Mahed clustered with the Q1 strain. The full-length genome of Mahed isolate shared the highest sequence homology with Gray and JMK (90.06–90.07%) and was least related to the Vic-s (86.21%). These data show that evolutionary variation because of recombination in IBV plays a major role in the adaptation and origin of IBV leading to new genetic and types of the virus strain. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s13337-018-0462-4) contains supplementary material, which is available to authorized users. ) and five structural proteins on the remaining: spike (S), envelope protein (E), membrane protein (M), and nucleocapsid protein (N), besides the accessory proteins 3a, 3b, 5a, and 5b, with an untranslated region (UTR) at both the 5 0 and 3 0 ends [4, 8] . Two-thirds of the 5 0 region in the IBV genome encodes the 1ab or 1b polyproteins, which are proteolytically cleaved by two virus-encoded replicase proteins (the papain-like and 3Clike proteinases) into 15 nonstructural proteins (nsp2-nsp16). The 5 0 and 3 0 untranslated regions (UTRs) usually harbor important structural elements that are involved in replication and translation [22] . Mutations and recombination in the IBV genome have resulted in viruses that have different tissue tropism, increased virulence, and an increased ability to persist in the chicken host [13] . Currently, more than 60 IBV serotypes have been reported because of mutation in its genome [10] . The isolation of IS/ 1494/06, one of the Israeli variant-2 (GI-23) isolates, and information on its S1 gene sequence in GenBank (Accession number: EU780077) has first been reported by Meier & Maher from Israel [12] . Strains of this lineage were first identified in Israel in 1998. The IS-1494 like IBVs is the main genotype in the Middle East, and African countries such as Iran, Jordan, Egypt, Turkey and recently reported from Europe [9] . The IS-1494 like IBV is known to be nephropathogenic, and it also affects the respiratory system [15] . Till now, seven genotypes of IBV (Massachusetts (Mass), 793/B, IS-1494, IS-720, QX, IR-1, and IR-2) are reported in Iran [16] . Contrary to the application of various IB vaccines (Mass & 793/B types) in Iran, regular outbreaks have been reported and a recent study showed IS-1494 like is the dominant (* 70%) IBV genotype in Iran [6] . This study aimed to determine the complete genome analysis of an Iranian IS-1494 like IBV isolate (Mahed;). The Iranian Is-1494 like IBV lineage (Mahed) was identified in Iran in December 2015 and occurred in a commercial Ross broiler flock vaccinated with IBV Mass type at day old one. Virus isolation has been done according to OIE manual. For genotyping, positive samples were submitted to a nested PCR reaction amplifying of S1 gene (Data has not been shown). To extract viral genomes from allantoic fluid (5th SF Egg passage) RNAeasy mini kit (Qiagen) was used according to manufacturer's instructions. For cDNA synthesis, random hexamer primer (SinaClon, Iran) was used and did according previous jobs [19] . Twenty-five primer pairs were synthesized based on the published LX4 strain genome (AY338732) [23] . PCR was performed with 3 ul cDNA as template in a total of 25 ul reaction volume containing 0.5 ul of each primer and 12.5 ul 29 PCR Mix (Sinaclon, Iran). To determine the complete genomic sequence, the extreme 5 0 and 3 0 termini of the isolates were identified by amplification of cDNA ends, using a 3 0 /5 0 RACE kit (Takara Bio Inc., Shiga, Japan). PCR products were inserted into a pTG19-T vector (Sinaclon, Iran) and transformed into a Top10 (Sinaclon, Iran) competent cell according to the TA cloning kit (Sinaclon, Iran) manufacturer's instructions.The positive clone have been selected based on blue white screening teat and clocny PCR. The nucleotide (nt) sequences of the positive clones were determined using T7 primers on a commercial service (Bioneer Co., Korea). Chromatograms were evaluated with ChromasPro (ChromasPro Version 1.5). Multiple alignments of different genes were achieved using Clustal W using the MEGA 7 software. The phylogenetic tree was constructed by using the same software with the neighbor-joining method (bootstrap replicates: 1000) [20] . The Mahed strain has submitted to gene bank under Accession number: MG233398. The complete genome sequence of the Mahed isolate is 27,652 nucleotides (nt) in length, excluding the poly (A) tail, and includes 10 open reading frames (ORFs). Thirteen ORFs (5 0 -1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3 0 ) were predicted in the Mahed genome (Supplementary material 1). ORFs 4b, 4c, and 6b were predicted in the Mahed genome and most of the GenBank IBV genomes. The length of nucleotide acid, amino acid and ORFs (Start and End) has been shown in (Supplementary material 1). Mahed Strain had the highest nucleotide sequence identity (90.07%) to Gray strain and the lowest nucleotide sequence identity (86.21%) to Vic-S strain. The nucleotide identity among strain Mahed and the H120 and 4/91 vaccine strain 89.41 and 88.39%, respectively ( Table 1 ). The phylogenetic tree showed the genetic relationships of the fulllength sequences of strain Mahed with other strains representing the majority of IBV genotypes showed that Mahed strain was closely related to QX (Fig. 1) . Based on the phylogenetic analysis of the replicase genes, Mahed clustered with Q1 and B1648 (1a) and QX (1b). Among the 16 NSPs of polyprotein 1ab, NSP15 (85.8-96.85%) and NSP16 (84.07-96.32%) were most variable, whereas the other NSPs (Especially NSP 4, 7, 13) were generally more conserved. Phylogenetic trees were constructed from the S1, and S2 genes of Mahed shows that the virus related to Q1 and QX, respectively (Supplementary material 2). We observed 81% nucleotide sequence identity (The highest) for the S1 gene between our isolate with Q1 strain. Based on the S2 analysis, the Mahed has high homology with QX (93.62%) ( Table 1 ). The amino acid (aa) sequence of cleavage site on the spike was RRTRR. The S2 fusion aa site in S2 of mahed is PSGRS that unique to QX and Vics strain. According to a phylogenetic study of the full-length genome, 1a, M, E, N protein, Mahed strain clustered with the QX type strain. Based on the partial 1b, Mahed clustered with the Q1 strain (Supplementary material 2). There were 7 probable accessory proteins in the Mahed strain, such as 3a (57 aa), 3b (64 aa), 4b (94 aa), 4c (56 aa) 5a (65 aa), 5b (82 aa) and 6b (74 aa). Over the past 15 years, the nephropathogenic IBV strains have been emerging as most predominant IBV strains in the poultry industry, especially in Asian and Middle East countries [17] . Hosseini et al. [7] reported IS-1494 like the virus with the frequency of 17.2% during their molecular surveillance of IBV genotypes involved in outbreaks between ''2010-2014''. In the study done by Najafi et al. (2015) , IS-1494 like IBV with the total prevalence of 34% in Iranian chicken flocks. Regarding the pathogenicity of Iranian IS-1494 like a virus, Najdfi et al. [14] demonstrated variant-2 has high tropism to the respiratory tract, digestive system, and renal tissue that is due to its epitheliotropic nature. The Mahed genome organization was: (5 0 -1a-1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b-3 0 ). 4b, 4c, and 6b were additional ORFs present in the Mahed genome. However, 4b, 4c, and 6b were present in the most of the IBV genomes, they have sometimes been reported. ORFs 4b, 4c, and 6b have also been reported in a turkey coronavirus. The detailed reason for the infrequent reports Complete genome analysis of Iranian IS-1494 like avian infectious bronchitis virus 391 of these ORFs (4b, 4c, and 6b) in most of the IBV genomes is not known. It could be that 4b, 4c, and 6b ORFs are present in most of the IBV genomes, but that the success of their identification depends on the algorithms of ORF prediction software that was used [17] . Lately, Bentley et al. [1] have confirmed and indicated 4b as a 5th secondary protein in IBV, as well 3a, 3b, 5a and 5b. The 4b homolog of Middle East respiratory syndrome coronavirus (MERS-CoV) and the 6b homolog of SARS coronavirus were found to be an interferon antagonist and apoptosis inducer [21] but for IBV, additional research is necessary to show the production of 4b, 4c, and 6b proteins and to recognize their actions in pathogenesis [17] . The S1 gene of Mahed IBV isolate with published sequences reference strains in GenBank, the isolate shared 95.12% nucleotide identity with attenuated IBVAR2-06, 95.22% homology with both IS/1494/06, Eg/CLEVB-2/ IBV/012 and 95.42% nucleotide similarity with TR8 (Turkish IBV) strains. Mahed shared amino acid sequence similarities of 92.16% with Attenuated IBVAR2-06 and 92.36% with both IS/1494/06 and Eg/CLEVB/2/IBV_012. Three hypervariable regions (HVR) of S1 were determined. The most variation was seen in HVR2. The S glycoprotein of IBV is translated as a precursor protein (S0), and then cleaved into two subunits (S1 and S2). The studies reveal that the cleavage recognition site sequence, which consists of five basic amino acid residues, does not appear to correlate with increased cleavability, host cell range and increased virulence as it does with envelope glycoproteins in orthomyxoviruses and paramyxoviruses, but correlates with IBVs in different geographic regions. In this study, the S1 cleavage site of strain Mahed was R-R-T-R-R (Arg-Arg-Thr-Arg-Arg). The proteolytic cleavages of the spike protein are essential for virus infection of host cells. The cleavage site between S1 and S2 contained the furin- Mahed strain is marked with a black circle recognition motif RRTRRS and was shared by all of the compared viruses [3] .The 3c gene encodes the envelope protein that has a critical role in assembly and budding of the virus [2] . Mahed IBV based on 3c gene is located in QX cluster in the phylogenetic tree based on E gene and has 90.07% homology to QX-like IB virus.The Matrix proteins are diverged considerably across genera but moderately well conserved within each CoV genus and are the most abundant component of CoV virions [11] . Variant-2 and Mass strains shared 99.69% homologies based on M gene sequences which contrast with the approximately 20% difference of the S1 gene. Our results agree with published data by Ren et al. [18] who characterized the M gene of IBV Chinese strain HH06. The M gene and protein share 83.9-97.9 and 83.6-96.5% homologous identities compared with IBV reference strains derived from different regions or countries, respectively. In our study, Mahed IBV had high homology with 4/91 (91.00%) and JMK (91.16%) IBVs. The effects of the single amino acid insertions in the ER motif of E and ectodomain of M proteins are unknown. However, the different amino acid sequences of the M protein ectodomains, identified when examining recurrent viruses, encourages further study of a role for IBV antigenic variation [5] . Finally, based on Nucleoprotein gene, Mahed is cluster with QX and 4/91. The close phylogenetic relationship of Mahed IBV (Iranian IBV) and QX are consistent with the hypothesis that their ancestors shared hosts during their evolutionary history and thus with the hypothesis of recombination between these two genotypes. Phylogenic analysis, it was shown that Mahed was a chimeric strain. The data obtained from this study suggests that this variant of IBV may have arisen through genetic modifications, which conferred changes in pathogenicity and antigenicity. In conclusion, our data suggested that both mutation and recombination were involving the emergence of IBV variants in the field. So, periodic monitoring is really necessary to track the emergence of new variants and is required to develop an effective vaccination plan to control IB. Identification of a noncanonically transcribed subgenomic mRNA of infectious bronchitis virus and other gammacoronaviruses The production of recombinant infectious DI-particles of a murine coronavirus in the absence of helper virus Severe acute respiratory syndrome coronavirus (SARS-CoV) infection inhibition using spike protein heptad repeat-derived peptides Complete genome sequence of a Brazil-type Avian coronavirus detected in a chicken Monoclonal antibodies to the matrix (El) glycoprotein of mouse hepatitis virus protect mice from encephalitis Genotyping of Avian infectious bronchitis viruses in Iran (2015-2017) reveals domination of IS-1494 like virus Epidemiology of avian infectious bronchitis virus genotypes in Iran Molecular evolution and emergence of avian gammacoronaviruses First characterization of a Middle-East GI-23 lineage (Var2-like) of infectious bronchitis virus in Europe Complete genome sequence analysis of a predominant infectious bronchitis virus (IBV) strain in China Identification of a novel nephropathogenic infectious bronchitis virus in Israel Genotypic and phenotypic characterization of the California 99 (Cal99) variant of infectious bronchitis virus Pathogenicity characteristics of an Iranian variant-2 (IS-1494) like infectious bronchitis virus in experimentally infected SPF chickens Pathogenicity characteristics of an Iranian variant-2 (IS-1494) like infectious bronchitis virus in experimentally infected SPF chickens Molecular characterization of infectious bronchitis viruses isolated from broiler chicken farms in Iran Genetic characterization of the Belgian nephropathogenic infectious bronchitis virus (NIBV) reference strain B1648 Characterization and membrane genebased phylogenetic analysis of avian infectious bronchitis virus Chinese strain HH06 Genotyping of infectious bronchitis viruses from broiler farms in Iraq during 2014-2015 MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony Complete genome analysis of Iranian IS-1494 like avian infectious bronchitis virus 393 methods A SARS-CoV protein, ORF-6, induces caspase-3 mediated, ER stress and JNK-dependent apoptosis Serotype shift of a 793/B genotype infectious bronchitis coronavirus by natural recombination Analysis of a QX-like avian infectious bronchitis virus genome identified recombination in the region containing the ORF 5a, ORF 5b, and nucleocapsid protein gene sequences