The prediction flow adopted is shown in Fig. 1.
Download sequences
Twelve F protein sequences and twelve HN protein sequences from NDV genotype XI from Madagascar were downloaded in FASTA format from the NCBI database (http://www.ncbi.nlm.nih.gov/protein). Details of these sequences are given in Table 1.
Table 1
List of sequences used in this study
Strains | Identification Number | Collection Date | References |
F proteins | HN proteins |
MG-1992 | ADQ64395 | ADQ64396 | 1992 | [19] |
MG-725 | ADQ64389 | ADQ64390 | 2008 |
MG-MEOLA | ADQ64398 | ADQ64399 | 2008 |
MG-39 | ADQ64401 | ADQ64402 | 2008 |
MGMNJ | AGC23353 | AGW43230 | 2009 | [20] |
MGF003C | AGC23347 | AGW43224 | 2010 |
MGF082T | AGC23349 | AGW43226 | 2010 |
MGF015C | AGC23348 | AGW43225 | 2011 |
MGF120T | AGC23350 | AGW43227 | 2011 |
MGF166 | AGC23351 | AGW43228 | 2011 |
MGF192C | AGC23352 | AGW43229 | 2011 |
MGS1130T | AGC23354 | AGW43231 | 2011 |
Prediction of continuous B epitopes
Continuous B-cell epitopes were predicted using the Bepipred2.0 Linear Epitope Prediction server (http://tools.iedb.org/bcell/) with a default threshold of 0.35. This server combines two different methods, including propensity scales for amino acids and the hidden Markov model [7]. In the results, only continuous B epitopes containing at least 6 amino acids are retained.
CTL epitope prediction
The NetCTLv1.2 server (https://www.cbs.dtu.dk/services/NetCTL/) was used to predict the cytotoxic T lymphocyte (CTL) epitopes present in the target proteins using the default threshold value of 0.75. The prediction is made using HLA alleles, as the immunocomputer servers do not possess chicken MHC alleles. In addition, three epitope characteristics including TAP transport efficiency, C-terminal cleavage and MHC-I binding peptides were combined in the prediction [8]. CTL epitopes with a combined score ≤ 1.0 percentile rank were selected.
HTL epitope prediction
The NetMHCIIpan4.0 server (http://www.cbs.dtu.dk/services/NetMHCIIpan/) was used to predict T helper lymphocyte (HTL) epitopes in F and HN proteins. This server predicted peptide binding to all HLA-II alleles (HLA-DR, HLA-DQ, HLA-DP), using artificial neural networks [9]. Peptide length was set at 15-mer and other parameters were used by default.
Analysis of antigenicity, allergenicity and membrane topology of epitopes
The antigenicity of predicted epitopes was analyzed with the VaxiJen v2.0 server (http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html) using a threshold of 0.4 [10]. In addition, the allergenicity and membrane topology of epitopes found to be antigenic are respectively analyzed in the AllerTOP v2.0 server (https://www.ddg-pharmfac.net/AllerTOP/) [11] and TMHMM v.20 (https://services.healthtech.dtu.dk/services/TMHMM-2.0/) [12].
Results and interpretation
Identification of continuous B epitopes
BepiPred 2.0 proposed 150 continuous B epitopes, including 72 for F proteins and 78 for HN proteins. However, only 18 of these 150 predicted continuous B epitopes (5 for F proteins and 13 for HN proteins) were found to be antigenic, non-allergenic and localized in the exo-membrane region (Table 2). Of these, the epitopes E1 and E15 are respectively the most antigenic in F and HN proteins. The continuous B epitopes identified in HN proteins are more numerous than those identified in F proteins.
Table 2
Continuous B epitopes in NDV F and HN proteins
Proteins | Epitope | Position | Peptide | Score VaxiJen | Allergenicity | Membrane topology |
F protein | E1 | 29–36 | TSSLDGRP | 1,3312 | NA | EM |
E2 | 222–232 | GPQITSPALTQ | 0,4337 | NA | EM |
E3 | 222–233 | GPQITSPALTQL | 0,4637 | NA | EM |
E4 | 250–262 | LTKLGVGNNQLSS | 0,7468 | NA | EM |
E5 | 349–357 | RVVTFPMSP | 0,6366 | NA | EM |
HN protein | E6 | 52–79 | TPGDLAGIPTEISRAEERITSTLSSNRD | 0,4949 | NA | EM |
E7 | 52–79 | TPGDLAGIPTEITRAEERITSALSSNR* | 0,5048 | NA | EM |
E8 | 52–80 | TPGDLAGIPTEISRAEERITSTLSSNRDV | 0,4563 | NA | EM |
E9 | 52–80 | TPGDLAGIPTEISRAEERITSTLSSNRDI* | 0,4749 | NA | EM |
E10 | 108–170 | TSLSYQINGAANNSGCGAPVHDPDYIGGIGKELIVDDTSEVTSFYPSAFQEHLNFIPAPTTGS | 0,4673 | NA | EM |
E11 | 108–170 | TSLSYQINGAANTSGCGAPVHDPDYIGGIGKELIIDDTSEVTSFYPSAFQEHLNFIPAPTTGS* | 0,4377 | NA | EM |
E12 | 108–171 | TSLSYQINGAANNSGCGAPVHDPDYIGGIGKELIVDDTSEVTSFYPSAFQEHLNFIPAPTTGSG | 0,4415 | NA | EM |
E13 | 109–171 | SLSYQINGAANNSGCGAPVHDPDYIGGIGKELIVDDTSEVTSFYPSAFQEHLNFIPAPTTGSG | 0,4226 | NA | EM |
E14 | 161–170 | NFIPAPTTGS | 0,7259 | NA | EM |
E15 | 280–292 | QYHEKDLDVITLF | 0,7614 | NA | EM |
E16 | 321–364 | KPNSPSDIAQEGRYVIYKRYNDTCPDEQDYQIRMAKSSYKPGRF | 0,6749 | NA | EM |
E17 | 321–364 | KPNSLSDIAQEGRYVIYKRYNDTCPDEQDYQIRMAKSSYKPERF* | 0,7293 | NA | EM |
E18 | 380–385 | TSLGED | 0,6376 | NA | EM |
(*) = Mutant epitopes, NA = Non-allergenic, EM = Exo-membrane
Identification of CTL epitopes
NetCTLv1.2 predicted 364 CTL epitopes, including 161 in F proteins and 203 in HN proteins. Of these, only 74 CTL epitopes (35 in F proteins and 39 in HN proteins) are antigenic, non-allergenic, non-toxic and located in the exo-membrane region (Table 3). Of these, the E epitopes 44 and E71 are the most antigenic in F and HN proteins respectively. Thus, CTL epitopes identified in HN proteins are more numerous than those identified in F proteins.
Table 3
CTL epitopes in NDV F and HN proteins
Proteins | Epitope | Position | Peptide | Score VaxiJen | Allergenicity | Membrane topology |
F protein | E19 | 12–20 | FPMPTVLIV | 0,4136 | NA | EM |
E20 | 19–27 | IVLALGCVR | 1,2102 | NA | EM |
E21 | 20–28 | VLALGCVRL | 1,3961 | NA | EM |
E22 | 35–43 | RPLAAAGIV | 0,4797 | NA | EM |
E23 | 117–125 | FVGAVIGSV | 0,4429 | NA | EM |
E24 | 119–127 | GAVIGSVAL | 0,7104 | NA | EM |
E25 | 119–127 | GAIIGSVAL* | 0,5387 | NA | EM |
E26 | 121–129 | VIGSVALGV | 1,1096 | NA | EM |
E27 | 121–129 | IIGSVALGV* | 0,9459 | NA | EM |
E28 | 205–213 | QVGVELSLY | 0,8613 | NA | EM |
E29 | 212–220 | LYLTELTTV | 0,8958 | NA | EM |
E30 | 227–235 | SPALTQLTI | 0,9003 | NA | EM |
E31 | 240–248 | NLAGSNMDY | 0,8480 | NA | EM |
E32 | 245–253 | NMDYLLTKL | 0,7282 | NA | EM |
E33 | 300–308 | ATYLETLSV | 0,4129 | NA | EM |
E34 | 329–337 | VIEELDTSY | 0,5085 | NA | EM |
E35 | 379–387 | LTTPYLTLK | 0,7754 | NA | EM |
E36 | 383–391 | YLTLKGSVI | 0,6714 | NA | EM |
E37 | 427–435 | LSLDGITLR | 0,8224 | NA | EM |
E38 | 428–436 | SLDGITLRL | 0,5295 | NA | EM |
E39 | 432–440 | ITLRLSGEF | 0,7007 | NA | EM |
E40 | 504–512 | YIVLTVVSL | 1,2955 | NA | EM |
E41 | 505–513 | IVLTVVSLV | 1,1583 | NA | EM |
E42 | 506–514 | VLTVVSLVF | 0,9971 | NA | EM |
E43 | 509–517 | VVSLVFGML | 1,0571 | NA | EM |
E44 | 511–519 | SLVFGMLSL | 1,6161 | NA | EM |
E45 | 512–520 | LVFGMLSLV | 1,4689 | NA | EM |
E46 | 513–521 | VFGMLSLVL* | 0,9617 | NA | EM |
E47 | 513–521 | VFGVLSLVL | 0,8891 | NA | EM |
E48 | 516–524 | MLSLVLICY* | 1,2434 | NA | EM |
E49 | 516–524 | VLSLVLICY | 1,0715 | NA | EM |
E50 | 517–525 | LSLVLICYL | 1,2555 | NA | EM |
E51 | 518–518 | SLVLICYLM | 0,8804 | NA | EM |
E52 | 519–527 | LVLICYLMY | 0,7853 | NA | EM |
E53 | 536–544 | LLWLGNNTL | 0,4064 | NA | EM |
HN protein | E54 | 22–30 | RLAFRIAIL | 1,2391 | NA | EM |
E55 | 22–30 | RLAFRISIL* | 1,3115 | NA | EM |
E56 | 23–31 | LAFRIAILL | 1,1265 | NA | EM |
E57 | 23–31 | LAFRISILL* | 1,1814 | NA | EM |
E58 | 25–33 | FRIAILLLI | 0,6820 | NA | EM |
E59 | 25–33 | FRISILLLI* | 0,7761 | NA | EM |
E60 | 29–37 | ILLLITIIL | 0,4397 | NA | EM |
E61 | 29–37 | ILLLITITL* | 1,0588 | NA | EM |
E62 | 30–38 | LLLITIILA | 0,4607 | NA | EM |
E63 | 30–38 | LLLITITLA* | 1,1094 | NA | EM |
E64 | 31–39 | LLITIILAL | 0,4530 | NA | EM |
E65 | 31–39 | LLITITLAL* | 0,9878 | NA | EM |
E66 | 36–44 | ILALSTATL | 0,5396 | NA | EM |
E67 | 36–44 | TLALSTAAL* | 0,4753 | NA | EM |
E68 | 95–103 | ALLSTESVI | 0,5043 | NA | EM |
E69 | 132–140 | YIGGIGKEL | 0,7564 | NA | EM |
E70 | 148–156 | VTSFYPSAF | 0,4192 | NA | EM |
E71 | 177–185 | SFDISATHY | 1,7398 | NA | EM |
E72 | 203–211 | HQYLAIGVL | 1,4544 | NA | EM |
E73 | 203–211 | HQYLALGVL* | 1,4158 | NA | EM |
E74 | 205–213 | YLALGVLRT | 1,1717 | NA | EM |
E75 | 273–281 | GRLGFDGQY | 1,9611 | NA | EM |
E76 | 274–282 | RLGFDGQYH | 1,3727 | NA | EM |
E77 | 281–289 | YHEKDLDVI | 1,3509 | NA | EM |
E78 | 299–307 | YPGVGGGSF | 1,0615 | NA | EM |
E79 | 312–320 | VWFPVYGGL | 0,4965 | NA | EM |
E80 | 321–329 | KPNSPSDIA | 0,5389 | NA | EM |
E81 | 321–329 | KPSSPSDIA* | 0,7347 | NA | EM |
E82 | 390–398 | VSPNTVTLM | 0,6229 | NA | EM |
E83 | 397–405 | LMGAEGRVL | 0,4415 | NA | EM |
E84 | 405–412 | LTVGTSHFL | 0,6614 | NA | EM |
E85 | 406–413 | TVGTSHFLY | 0,5265 | NA | EM |
E86 | 418–426 | SSYFSPALL | 0,7036 | NA | EM |
E87 | 419–427 | SYFSPALLY | 0,6322 | NA | EM |
E88 | 468–476 | GVYTDPYPL | 0,6599 | NA | EM |
E89 | 469–477 | VYTDPYPLV | 0,6099 | NA | EM |
E90 | 470–478 | YTDPYPLVF | 1,0589 | NA | EM |
E91 | 551–559 | TLFGEFRIV | 0,4567 | NA | EM |
E92 | 557–565 | RIVPLLVEI | 1,1852 | NA | EM |
(*) = Mutant epitopes, NA = Non-allergenic, EM = Exo-membrane
Identification of HTL epitopes
NetMHCIIpan 4.0 predicted 157 HTL epitopes, including 72 in F proteins and 85 in HN proteins. However, only 28 HTL epitopes (12 in F proteins and 16 in HN proteins) are antigenic, non-allergenic, non-toxic and localized in exo-membrane regions (Table 4). Of these, the E epitopes 93 and E111 are the most antigenic in F and HN proteins respectively. Thus, the HTL epitopes identified in HN proteins are more numerous than those identified in F proteins.
Tableau 4 : Les épitopes HTL dans les protéines F et HN de VMN
Proteins | Epitope | Position | Peptide | Score VaxiJen | Allergenicity | Membrane topology |
F protein | E93 | 8–22 | WIPISPMPTILIALA | 0,8658 | NA | EM |
E94 | 8–22 | WIPIFPMPTVLIVLA* | 0,4612 | NA | EM |
E95 | 34–48 | GRPLAAAGIVVTGDK | 0,4569 | NA | EM |
E96 | 35–49 | RPLAAAGIVVTGDKA | 0,5427 | NA | EM |
E97 | 118–132 | VGAIIGSVALGVATA | 0,8164 | NA | EM |
E98 | 126–139 | ALGVATAAQITAAAA | 0,8498 | NA | EM |
E99 | 235–249 | IQALYNLAGSNMDYL | 0,4999 | NA | EM |
E100 | 318–332 | LVPKVVTQVGSVIEE | 0,4500 | NA | EM |
E101 | 323–337 | VTQVGSVIEELDTSH | 0,6443 | NA | EM |
E102 | 452–466 | DSQVIVTGNLDISTE | 0,5541 | NA | EM |
E103 | 452–466 | DSQVIVTGNLDISAE* | 0,7456 | NA | EM |
E104 | 532–546 | QQKTLLWLGNNTLDQ | 0,4601 | NA | EM |
HN protein | E105 | 43–57 | ALIYSMEASTPGDLA | 0,4307 | NA | EM |
E106 | 53–67 | PGDLAGIPTEITRAE | 0,4359 | NA | EM |
E107 | 129–143 | DPDYIGGIGKELIVD | 0,5995 | NA | EM |
E108 | 129–143 | DPDYIGGIGKELIID* | 0,6788 | NA | EM |
E109 | 130–144 | PDYIGGIGKELIIDD | 0,4315 | NA | EM |
E110 | 157–171 | QEHLNFIPAPTTGSG | 0,9029 | NA | EM |
E111 | 202–216 | SHQYLAIGVLRTSAT | 0,9224 | NA | EM |
E112 | 202–216 | SHQYLALGVLRTSA* | 0,8896 | NA | EM |
E113 | 218–232 | RVLFSTLRSINLDDN | 0,6104 | NA | EM |
E114 | 218–232 | RVFFSTLRSINLDDN* | 0,4580 | NA | EM |
E115 | 221–235 | FSTLRSINLDDNQNR | 0,6869 | NA | EM |
E116 | 317–331 | YGGLKPNSPSDIAQE | 0,5425 | NA | EM |
E117 | 379–393 | STSLSEDPVLTVSPN | 0,6043 | NA | EM |
E118 | 475–489 | PLVFHRNHTLRGVFG | 0,5237 | NA | EM |
E119 | 499–513 | LNPVSAVFDNIARSR | 0,5124 | NA | EM |
E120 | 555–569 | EFRIVPLLVEILKDG | 0,7710 | NA | EM |
(*) = Mutant epitopes, NA = Non-allergenic, EM = Exo-membrane