key: cord-269024-re0hyolh
authors: Bande, Faruku; Arshad, Siti Suri; Hair Bejo, Mohd; Kadkhodaei, Saeid; Omar, Abdul Rahman
title: Prediction and In Silico Identification of Novel B-Cells and T-Cells Epitopes in the S1-Spike Glycoprotein of M41 and CR88 (793/B) Infectious Bronchitis Virus Serotypes for Application in Peptide Vaccines
date: 2016-09-07
journal: Adv Bioinformatics
DOI: 10.1155/2016/5484972
sha: 
doc_id: 269024
cord_uid: re0hyolh

Bioinformatic analysis was used to predict antigenic B-cell and T-cell epitopes within the S1 glycoprotein of M41 and CR88 IBV strains. A conserved linear B-cell epitope peptide, YTSNETTDVTS(175–185), was identified in M41 IBV strains while three such epitopes types namely, VSNASPNSGGVD(279–290), HPKCNFRPENI(328–338), and NETNNAGSVSDCTAGT(54–69), were predicted in CR88 IBV strains. Analysis of MHCI binding peptides in M41 IBV strains revealed the presence of 15 antigenic peptides out of which 12 were highly conserved in 96–100% of the total M41 strains analysed. Interestingly three of these peptides, GGPITYKVM(208), WFNSLSVSI(356), and YLADAGLAI(472), relatively had high antigenicity index (>1.0). On the other hand, 11 MHCI binding epitope peptides were identified in CR88 IBV strains. Of these, five peptides were found to be highly conserved with a range between 90% and 97%. However, WFNSLSVSL(358), SYNISAASV(88), and YNISAASVA(89) peptides comparably showed high antigenicity scores (>1.0). Combination of antigenic B-cells and T-cells peptides that are conserved across many strains as approach to evoke humoral and CTL immune response will potentially lead to a broad-based vaccine that could reduce the challenges in using live attenuated vaccine technology in the control of IBV infection in poultry.

Infectious bronchitis virus (IBV) is a single stranded, enveloped RNA virus belonging to the family Coronaviridae, order Nidovirales [1] . The virus causes infectious bronchitis (IB), a contagious disease associated with huge economic loses in poultry industry worldwide [2] . Of major concern in the control of IB is the continued emergence of variant IBV strains that differ in terms of their tissue tropism, pathogenicity, and cross protection. Over the years, serological and molecular studies have been carried out extensively to determine the epidemiology of local IBV strains. Remarkably, both classical and variants IBV strains have been reported in different countries [3] . Among the widely identified IBV strains are IBV M41 (classical strain), originally recognized in USA [4, 5] , and CR88 (variant strain otherwise known as 793/B or 4/91) which was first reported in Europe [3, 6, 7] .

Over the years, control of IBV infection largely depends on vaccination using live attenuated and killed vaccines.

However, one of the challenges with live attenuated IB vaccines is that such vaccines are reported to encourage mutation and recombination, thus leading to the emergence of new variant strains. Live attenuated vaccines have also been linked with reversion to virulence, severe disease, and increased mortality rate [8] [9] [10] . On the other hand, killed vaccines induce humoral but not cell mediated immune (CMI) response and in most cases require adjuvants as well as repeated boosting especially in laying chickens and breeder flocks. These challenges therefore necessitate the needs for novel broad vaccines for the control of IB in poultry [11] . To achieve this, stimulation of both humoral (B-cell) and cell mediated immune (CMI) responses is considered very essential for any candidate vaccine [12] . Neutralizing antibodies are important in removing freely circulating IB virus, whereas cytotoxic T-lymphocytes (CTL) response is crucial for the control and clearance of virally infected cells. The latter is achieved through MHCI immune surveillance as well as antigen presentations which is the function of MHCI molecule and both have been associated with the epitope within the S1 glycoprotein [13, 14] . While much has been documented on MHC restricted allele in human and mouse models, little information is available on the biological functions of these molecules in poultry [15] . A large binding groove of BF2 * 2101 MHCI molecule identified in B21 chicken line is thought to confer conformational flexibility to the crucial Arg9 residue which allows remodeling of key peptide-binding sites and play a role in the resistance against poultry viral infections [16] . Chicken MHC B-F molecules have been structurally and functionally linked to mammalian MHC class I molecules and involved in antigen presentation to the CD8 + T lymphocytes, which is crucial in antiviral immune response [17] . Interestingly, the S1 glycoprotein of IBV (520 aa) contains different immune epitopes responsible for both antibodies and CTL-based immune responses, thus playing major protective role as viral antigenic determinant [18] .

Currently, the use of peptide based DNA vaccines represents a novel strategy for addressing challenges associated with the control of viral infections [19, 20] . This technology may employ the use of in silico analysis to predict novel B-cells and T-cells immune epitopes for further use in their in vivo applications. One of the innovations in using this technology is the ability to incorporate several epitope peptides directed against different viruses and/or multiple virus strains into one single delivery system with the view to induce broad and specific immune response in single administration [11] . To date, only few epitope based peptide vaccines have been developed and evaluated against IB [20, 21] . The objective of this study therefore is to identify novel B-cells and T-cell epitopes within the S1 glycoprotein of M41 and CR88 IBV strains. The antigenicity of the predicted peptides was also evaluated.

Nucleotide sequences containing complete S1-gene of homologous to M41 and CR88 IBV strains were retrieved from NCBI database and translated into deduced amino acid sequences for further analysis. Consensus amino acid sequences from the two datasets (M41 and CR88) were used to predict B-and Tcells epitope peptides within S1 glycoprotein. The criteria used for the selection of M41 and CR88 were simply based on epidemiological relevance of the two strains in terms of their wide geographic distribution as classical and variant IBV strains, respectively. All sequences were assembled with Geneious5 software R7 version [22] . The names and accession numbers of IBV sequences used in this study are presented in Supplementary File 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/5484972.

Epitopes associated with Bcells were predicted using the BepiPred epitope prediction server version 1.0 [23] . Similarly, antigenicity index of the predicted epitopes was analysed using VaxiJen v2.0 online antigen prediction tool http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html [22] . Prediction criteria were set at 75% classifier specificity and 20 epitope lengths with an overlap. The location of the predicted epitope within transmembrane regions was evaluated using TMHMM [24] . Only epitope present at the surface of the membrane was selected and further analysed for antigenicity. For this search, target organism was narrowed to viruses and only scores above 0.4 thresholds were considered as a good antigenic epitope. Further, conservancy analysis was carried out using IEDB tools available http://tools.immuneepitope.org/tools/ conservancy/iedb input [25] .

The T-cell epitope prediction was carried out using ProPred-1 software which covers 47 MHC class-I alleles. Both proteasome and immuneproteasome were set at 5% threshold by MHC class I and only peptides with proteosomal cleavage site at the C terminus were considered [26] .

Immune epitopes were analysed for conservancy and variability using IEDB conservation analysis tool: http://tools.immuneepitope.org/tools/ conservancy/iedb input [25] .

A total of 86 M41 and 72 CR88 IBV isolates were considered in this study. However, to build the consensus sequence, one isolate from each data set was removed due to some unidentified amino acid at one or two positions. The M41 IBV isolates were found to have an average nucleotides sequence length of 1614 (537 aa) while CR88 isolates had an average length of 1653 (553 aa).

Analysis of M41 associated linear B-cell epitopes revealed the presence of six unique B-cells peptides within the S1 glycoprotein. Of these, only YTSNETTDVTS 175-185 peptide was found to be antigenic (Table 1) . Analysis of transmembrane potentials within the M41 IBV S1 gene showed that the predicted linear B-cell epitopes interact with the outside surface membrane especially with maximum threshold occurring at amino acid positions 280 to 290 (Figure 1) .

Interestingly, conservancy analysis showed that YTSNETTDVTS 175-185 peptide is common in 82.35% of the M41 strains. On the other hand, a total of six antigenic epitope peptides were predicted in CR88 IBV strains; however, three of these peptides, NETNNAGSVSDCTAGT 54-69 , VSNASPNSGGVD 279-290 , and HPKCNFRPENI 328-338 were demonstrated to be antigenic (Table 2) . Further, conservancy analysis showed that the three CR88 associated antigenic B-cell peptides had 39.44%, 84.5%, and 43.66% conservancy rates, respectively.

As with the M41 IBV strains, most of these epitopes are shown to interact with the surface at varying threshold (yellow) as depicted in Figure 2 . 

In M41 IBV strains, a total of 21 MHCI peptides were predicted at different amino acid positions; however, 15 of these peptides were found to be antigenic at various conservancy rates (98.82-100%). Three epitope peptides, GGPITYKVM 208 , WFNSLSVSI 356 , and YLADAGLAI 472 relatively have high antigenicity index (>1.0) as compared to other peptides (Table 3 ).

In the case of CR88 IBV strains, 18 MHCI epitope binding peptides were predicted. Of these, 11 peptides were found to be antigenic with most of the epitopes having conservancy rate ranging from 52.11% to 94.37% except peptide occurring at position S69 which is conserved in only 8.45%. Interestingly, WFNSLSVSL 358 , SYNISAASV 88 , and YNISAASVA 89 demonstrated high antigenicity index compared to other predicted peptides (Table 4 ). 

The use of bioinformatic and/or in silico analyses to understand infectious diseases as well as to predict novel vaccine candidate has been recently extended to poultry [27, 28] . Most of the epitopes responsible for virus neutralization have been mapped to be located within the first and third quarters of the linear S1 glycoprotein [13, 16] . The present study identified novel B-cells and T-cells epitopes presence in the S1 glycoprotein of M41 and CR88 IBV strains. Predicted antigenic B-cell epitopes were found to be highly conserved and demonstrated strong transmembrane potentials which probably predicts regions of virus-cell interaction. Linear Bcells epitopes located within the S1 region have been reported to play a role in virus neutralization. Using a phage display library, Zou et al. [29] identified two neutralizing linear Bcells epitopes within the S1 glycoprotein.

Recent studies have shown that chicken MHCI genes are categorized into MHCI associated genes (B-F) and MHCII (B-L) associated B-G genes, probably found only in poultry. Remarkably, this study has identified 21 and 18 MHCI binding peptides located within the S1 glycoproteins of M41 and CR88 IBV strains, respectively. Chickens with MHC homozygous B12 and B19 are reported to be more susceptible to infection with IBV-Grey strains as compared with B2/B2 and B5/B5 haplotype which often resist infection [30] . Recently, Tan et al. [21] used BF2 * 4, BF2 * 12, BF2 * 15, and BF2 * 19 chicken MHCI haplotypes to predict 21 CTL-peptide candidates in Massachusetts and Australian T strains. The study revealed that the constructed poly-CTL-epitope DNA vaccine was capable of inducing protection in 90% of the vaccinated chicken following challenge with SH1208 IBV strain. In a different study, Tian et al., 2008 [31] have identified seven T-and B-cells epitopes within the IBV S1 and S2 and NP proteins chimeric DNA vaccine derived from these epitope was found to be associated with 80% protection rate. Interestingly, the regions through which B-cells and T-cells epitope were predicted in this study correspond to the IBV receptor binding domain (RBD) which play a role in sufficient binding as well as entry of IBV viruses [32] . Similarly, epitopes peptides found within RBD have been reported in different study to induce neutralizing antibody response, for example, against human coronavirus and mouse hepatitis virus [33, 34] .

This study predicted novel antigenic B-cells and T-cell epitopes within the S1 glycoprotein of M41 and CR88 IBV strains. The characteristic antigenicity index as well as epitope conservancy rates demonstrates potentials of the identified epitope peptides as polyvalent synthetic or DNA-based peptide vaccine for application in the control of IBV infection. The use of such vaccines will likely reduce the challenges associated with live attenuated vaccines and allow broad coverage of the target IBV strains.

Nidovirales: a new order comprising Coronaviridae and Arteriviridae

Coronavirus avian infectious bronchitis virus

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

Genetic and antigenic diversity in avian infectious bronchitis virus isolates of the 1940s

Detection of Massachusetts and Arkansas serotypes of infectious bronchitis virus in broilers

Characterisation of an infectious bronchitis virus isolated from vaccinated broiler breeder flocks

Proceedings of the Angers 1eres Journee De La Recherche Avicole, Centre De Congres d' Angers

Reversion to virulence of chicken-passaged infectious bronchitis vaccine virus

Attenuated live vaccine usage affects accurate measures of virus diversity and mutation rates in avian coronavirus infectious bronchitis virus

Pathogenesis and diagnostic approaches of avian infectious bronchitis

Progress and challenges toward the development of vaccines against avian infectious bronchitis

Priming with a DNA vaccine and boosting with an inactivated vaccine enhance the immune response against infectious bronchitis virus

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

Cytotoxic T lymphocytes are critical in the control of infectious bronchitis virus in poultry

Alleledependent association of chicken MHC class I molecules with the invariant chain

Structures of an MHC class I molecule from B21 chickens illustrate promiscuous peptide binding

Peptide motifs of the single dominantly expressed class I molecule explain the striking MHC-determined response to Rous sarcoma virus in chickens

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

Prediction and identification of T cell epitopes in the H5N1 influenza virus nucleoprotein in chicken

The protective immune response against infectious bronchitis virus induced by multiepitope based peptide vaccines

Prediction and identification of novel IBV S1 protein derived CTL epitopes in chicken

Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data

Predicting linear B-cell epitopes using string kernels

Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes

Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines

In silico prediction of epitopes in virulence proteins of mycobacterium tuberculosis H37Rv for diagnostic and subunit vaccine design

Bioinformatics and evolutionary insight on the spike glycoprotein gene of QXlike and Massachusetts strains of infectious bronchitis virus

Prediction of CTL epitope, in silico modeling and functional analysis of cytolethal distending toxin (CDT) protein of Campylobacter jejuni

Two novel neutralizing antigenic epitopes of the s1 subunit protein of a QX-like avian infectious bronchitis virus strain Sczy3 as revealed using a phage display peptide library

Association of the chicken MHC B haplotypes with resistance to avian coronavirus

The immunoreactivity of a chimeric multi-epitope DNA vaccine against IBV in chickens

Mapping of the receptorbinding domain and amino acids critical for attachment in the spike protein of avian coronavirus infectious bronchitis virus

Localization of neutralizing epitopes and the receptor-binding site within the amino-terminal 330 amino acids of the murine coronavirus spike protein

Identification of a receptor-binding domain of the spike glycoprotein of human coronavirus HCoV-229E

The authors would like to thank the Universiti Putra Malaysia and Ministry of Science, Technology and Innovation (MOSTI), Project no. 02-01-04-SF1070 for funding the research.

The authors declared no conflict of interests regarding the present publication.